ORIGINAL RESEARCH article
Hybrid MM/CG Webserver: Automatic Set Up of Molecular Mechanics/Coarse-Grained Simulations for Human G Protein-Coupled Receptor/Ligand Complexes
- 1Computational Biomedicine, Institute for Advanced Simulation IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, Germany
- 2JARA-Institute: Molecular Neuroscience and Neuroimaging, Institute for Neuroscience and Medicine INM-11/JARA-BRAIN Institute JBI-2, Forschungszentrum Jülich GmbH, Jülich, Germany
- 3Department of Physics, RWTH Aachen University, Aachen, Germany
- 4Department of Biotechnology, University of Verona, Verona, Italy
- 5Medical Faculty, Cécile and Oskar Vogt Institute for Brain Research, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
Hybrid Molecular Mechanics/Coarse-Grained (MM/CG) simulations help predict ligand poses in human G protein-coupled receptors (hGPCRs), the most important protein superfamily for pharmacological applications. This approach allows the description of the ligand, the binding cavity, and the surrounding water molecules at atomistic resolution, while coarse-graining the rest of the receptor. Here, we present the Hybrid MM/CG Webserver (mmcg.grs.kfa-juelich.de) that automatizes and speeds up the MM/CG simulation setup of hGPCR/ligand complexes. Initial structures for such complexes can be easily and efficiently generated with other webservers. The Hybrid MM/CG server also allows for equilibration of the systems, either fully automatically or interactively. The results are visualized online (using both interactive 3D visualizations and analysis plots), helping the user identify possible issues and modify the setup parameters accordingly. Furthermore, the prepared system can be downloaded and the simulation continued locally.
Human G protein-coupled receptors (hGPCRs) are the largest drug target superfamily (Hauser et al., 2017). One third of FDA-approved drugs target ∼14% hGPCRs (Hauser et al., 2018) and this protein class has a further, untapped pharmacological potential. Unfortunately, rational drug design is hampered by the lack of experimental structures for more than 90% hGPCRs1 (Munk et al., 2019; Qu et al., 2020). Structural insights into ligand poses for these cases can be obtained by computational modeling. Reliable predictions can be made by docking approaches on homology models, based on templates sharing overall sequence identity >35–40% and/or high conservation of binding site residues (Beuming and Sherman, 2012; Kufareva et al., 2014). For lower resolution models, however, the uncertainty in the structure (and particularly in the orientation of side chains) decreases the accuracy of the docking predictions and, thus, follow-up molecular dynamics (MD) simulations are recommended (Kufareva et al., 2014; Cavasotto and Palomba, 2015; Esguerra et al., 2016; Heifetz et al., 2016; Fierro et al., 2017; Lupala et al., 2018; Rodríguez-Espigares et al., 2020). In an effort at addressing this issue, we have developed a Hybrid Molecular Mechanics/Coarse-Grained (MM/CG) simulation approach (Neri et al., 2005, 2008; Leguèbe et al., 2012; Marchiori et al., 2013; Sandal et al., 2015; Capaldi et al., 2018; Alfonso-Prieto et al., 2019; Fierro et al., 2019). The receptor/ligand interactions are described in atomistic detail, including explicit water molecules in the binding site (MM region), while the rest of the receptor is coarse-grained (CG region) (Schneider et al., 2018). The all-atom force fields used (Schneider et al., 2020) for the MM part of the protein and water are the Amber14SB (Maier et al., 2015) and TIP3P (Jorgensen et al., 1983), respectively, whereas the ligand can be described using either GAFF or GAFF2 (Wang et al., 2004; Case et al., 2020). The CG region is described by a Gō-like (Gō and Abe, 1981) potential. A region at the interface between the MM and CG parts couples the two levels of resolution (Figure 1). The membrane is described implicitly by introducing five potential walls (Leguèbe et al., 2012; Schneider et al., 2018). Two planar walls coincide with the height of the head groups of the membrane lipids, two hemispheric walls cap the extracellular and intracellular ends of the protein and prevent water evaporation, and the last wall follows the initial shape of the interface between protein and membrane, mimicking the effect of the lipid acyl tails (Figure 1). The scheme can be used also for GPCRs other than those from Homo sapiens. This approach turned out to be able to reproduce the ligand poses for four different hGPCRs (Schneider et al., 2020). These include the adenosine 2A receptor in complex with caffeine (i.e., the example case of the webserver), the human bitter receptor 16 in complex with phenyl-β-D-glucopyranoside, the β2-adrenergic receptor with adrenaline, and the dopamine D3 receptor with eticlopride. The structures of these hGPCR/ligand complexes were either experimentally determined, taken from all-atom MD trajectories or predicted based on templates with decreasing resolution, up to extremely low sequence identity (<15%) (Schneider et al., 2020). Retrospective validation against available X-ray structures and mutagenesis data confirmed that the MM/CG approach can predict correct ligand poses and identify experimentally determined binding residues2 (Schneider et al., 2020), regardless of the model resolution. In addition, the MM/CG simulations can provide insights into the flexibility of receptor–ligand interactions and hydration of the binding cavity, at a lower computational cost than all-atom molecular dynamics simulations.
Figure 1. Hybrid MM/CG scheme. The three regions (MM, I, and CG) are framed with different background colors and the five potential walls (upper and lower hemispheres, upper and lower membrane planes, and membrane surface) are indicated with black lines.
The increasing number of applications and requests to use the MM/CG approach has prompted us to develop the Hybrid MM/CG webserver, a publicly accessible web interface aimed at preparing and running short MM/CG simulations of hGPCR/ligand complexes (available since December 2019). To the best of our knowledge, this is the only webserver dedicated to this task. It complements other excellent online resources (Table 1) aimed at preparing all-atom or coarse-grained molecular dynamics (MD) simulations of GPCRs and GPCR/ligand complexes (or in general membrane proteins).
The Hybrid MM/CG webserver requires only the coordinates of the receptor/ligand complex (as PDB file). These may come from experimental structures, simulation snapshots, or computational models generated with other webservers, such as the ones listed in the “Input” section. Furthermore, we have linked our GOMoDo webserver3 for modeling and docking of GPCRs to the Hybrid MM/CG webserver, so that the complexes generated with GOMoDo can be directly transferred. The user is first guided through a set of user-friendly forms for preparing the setup files. The procedure requires only a few minutes. Then, the user runs the initial MM/CG simulation steps (up to 10 ns) directly on the server. The resulting files can be downloaded and the MM/CG simulation can be continued using local resources.
Materials and Methods
The input PDB file can be obtained from experimental databases, such as the Protein Data Bank (Berman et al., 2000; Rose et al., 2016), GPCRdb (Pándy-Szekeres et al., 2018) or GPCRmd (Rodríguez-Espigares et al., 2020), as well as computational services like GOMoDo (Sandal et al., 2013), GPCR-ModSim (Esguerra et al., 2016), GPCR-SSFE (Worth et al., 2017), GPCRM (Miszta et al., 2018), Galaxy7TM (Lee and Seok, 2016), GPCRautomodel (Launay et al., 2012), @TOME (Pons and Labesse, 2009) and others (reviewed in Busato and Giorgetti, 2016). The GOMoDo webserver (Sandal et al., 2013), which can be used for homology modeling of hGPCRs and subsequent docking of ligands, is linked directly with the Hybrid MM/CG webserver. Several structures of the same hGPCR/ligand complex obtained under different membrane compositions and/or in different activation states can be funneled into the webserver to indirectly (and very approximately) account for the influence of explicit lipids and of large conformational changes of the receptor, respectively.
The Hybrid MM/CG webserver offers two system setup modes:
Here, the user just needs to upload a PDB file of the hGPCR in complex with its ligand or transfer it directly from the GOMoDo webserver3 (Sandal et al., 2013). All simulation parameters are set according to default values that are defined in the documentation4.
Using the interactive preparation method, the same steps as in the automatic preparation are carried out to set up the system. The advantage is that several parameters, such as the position of the interface between MM and CG regions and the size of the hemispheric potential walls (see Figure 1), can be adjusted. Furthermore, intermediate results can be visualized, such as the input structure, aligned orientation in membrane, solvation (water drop), level of coarse-graining, and position of the wall potentials. More details about the individual preparation steps and parameters are explained in the Documentation section of the webserver4.
The steps carried out by the webserver for MM/CG system preparation and short simulation (Figure 2) are the following:
• File upload. The preparation starts with a PDB file of the hGPCR/ligand complex. This file can be obtained from one of the databases and webservers mentioned in the “Input” section and uploaded by the user into the Hybrid MM/CG webserver. Alternatively, it can be transferred directly from GOMoDo (Sandal et al., 2013).
• Check Input. The input PDB file is checked for missing residues and the numbering is corrected if possible (i.e., no residues are missing). The ligand name is determined by comparing the list of residue names in the input PDB file and in the Amber14SB force field. Known amino acids and capping groups are discarded and the remaining residue name is considered to be the ligand.
• Alignment. The orientation and position of the receptor with respect to the hydrophobic core of the lipid bilayer are determined using the PPM tool (Lomize et al., 2012). The initial PDB file of the hGPCR/ligand complex is then superimposed to the PPM-aligned structure of the receptor using lovoalign (Martínez et al., 2007).
• Build topology: Protein. The program pdb2gmx (Berendsen et al., 1995; Hess et al., 2008) is used to build the receptor topology using the Amber14SB force field (Maier et al., 2015). Protonation states of receptor titratable residues are determined automatically by pdb2gmx. In case the user wants to bypass the automatic assignment, different protonation states can be enforced with minimal manual editing of the input PDB file by changing the corresponding residue name (e.g., from “HIS” to “HIP” for a doubly protonated histidine).
• Build topology: Ligand. Protonation of the ligand at pH = 7 is determined using OpenBabel (O’Boyle et al., 2011). The ligand topology is built with ACPYPE (Sousa da Silva and Vranken, 2012) and Antechamber (Wang et al., 2001), using AM1-BCC charges (Jakalian et al., 2004) and the general Amber force field. GAFF (Wang et al., 2004) is used in the automatic preparation, but the user can choose between GAFF and GAFF2 (Case et al., 2020) in the interactive preparation mode.
• Solvation. The simulation box is created and the system is solvated using the TIP3P (Jorgensen et al., 1983) force field. Water molecules below the upper membrane plane are deleted. Hemisphere sizes are defined (default radii: 50 Å) and water outside the upper hemisphere is deleted, so that only a water “drop” solvating the extracellular part of the receptor is kept (see Figure 1).
• Coarse-graining. The regions of different resolutions are defined by two cutoff values, one defining the boundary between the MM and interface regions (default is centered between the two membrane planes) and the other between the interface and CG regions (6 Å below the first cutoff). Coarse-graining is performed on residues below the second cutoff level using a Gō-like (Gō and Abe, 1981) model.
• Minimization. A simple minimization of the system is done using the steepest descent method until a maximum force of 1,000 kJ mol–1 nm–1 is reached.
• Add wall potentials. Wall potentials are added to the system according to the aforementioned heights and radii. Wall grid points that are too close to the ligand are deleted according to the chosen cutoff (default 7 Å).
• Visualization. The user is redirected to the “Results” section where preparation results can be inspected using interactive 3D visualizations and graphs.
• Equilibration. An optional short equilibration (2 ns) and short MD simulation (2 ns) can be started from the “Results” section as well.
• MD. The simulation can be extended up to 10 ns, in 2 ns increments, in order to offer the possibility to check the intermediate results.
• Visualization. The minimization, equilibration, and MD are analyzed automatically in terms of potential energy and temperature after completion of the individual step (see Supporting Information). Interactive 3D visualizations are available for the results of the preparation, equilibration, and MD (see Figure 3).
Figure 2. Webserver workflow. Interactive mode with manual intervention is presented on the left and fast automatic mode on the right.
Figure 3. Sample 3D visualizations created by the webserver using as example the hA2AR/CFF complex. (A–C) Full view of the hGPCR/ligand complex and (D–F) close-up view of the ligand and its interactions with protein and water molecules. Images were generated with the screenshot tool in the NGL viewer (Jakalian et al., 2004) used in the “Results” section of the webserver. Atoms are colored using the default NGL color code and ligand interactions are determined using the default NGL definition.
Results can be viewed and downloaded for 2 weeks (renewal possible) by bookmarking the link or alternatively by using the corresponding ID5. The full output of the preparation can be downloaded as a compressed archive file including the input, output, and log files of all preparation and simulation steps. The downloaded files can be used to continue the MM/CG simulations locally. Experienced users have the possibility to download the prepared system and tune the partial charges, as well as other parameters, before running the simulation on their local computer. The Download Code section of the server provides an installation script, which contains the links to the source code for both GROMACS 4.5.1 (from the GROMACS website) and the Hybrid MM/CG patch (from our server website).
As an application case, we considered the human adenosine 2A receptor (hA2AR) in complex with its antagonist caffeine (CFF). CFF was proposed to have a protective and therapeutic effect against Parkinson’s disease (Chen et al., 2001; Sonsalla et al., 2012; Petzer and Petzer, 2015; Nazario et al., 2017). The input PDB file of the hA2AR/CFF complex used for the Example can be downloaded from the same page. The demonstration can be run by simply starting the preparation workflow9 without uploading any PDB file. Precomputed results for this example case can be found on the server10 and selected parts of these results can be seen in Figure 3. Figures 3A–C show the whole hA2A receptor in complex with CFF, with and without water solvating the binding site and the extracellular loops, as well as the walls mimicking the membrane and preventing water evaporation. Figures 3D–F show the ligand and its direct and water-mediated interactions with residues in the binding cavity, as determined by the NGL viewer (Jakalian et al., 2004). The system can be inspected in 3D after the preparation, short equilibration, and production runs. In addition, plots showing the time evolution of the potential energy and temperature during the different simulation steps are generated (see Supplementary Figures 1–3).
The hybrid MM/CG protocol has been successfully used to predict ligand poses in a variety of hGPCRs (Leguèbe et al., 2012; Marchiori et al., 2013; Sandal et al., 2015; Fierro et al., 2017; Capaldi et al., 2018; Fierro et al., 2019). Moreover, the Amber-based MM/CG poses improved significantly relative to those obtained by simple docking, especially for low resolution starting models (Schneider et al., 2020). However, the setup of such simulations has been time-consuming and system-dependent. The webserver presented here automatically prepares all files needed to start the equilibration, in a short time: the hA2AR/CFF example requires only 6 min. An additional advantage of the Hybrid MM/CG webserver over most of the online services for simulation of GPCRs (Table 1) is the automatic ligand parameterization. Although other webservers can also include ligands (see Table 1), their parameterization usually has to be done externally. The results of the preparation steps can be checked on-the-fly (e.g., ligand structure and full solvation of extracellular loops). No additional software (except the patch for the GROMACS code, which is available for download11) is needed to extend the simulations locally. This makes the method useful for people acquainted with MD simulations and not necessarily familiar with every single MM/CG preparation step.
The server could prospectively be used for setting up MM/CG simulations of different ligands in complex with the same hGPCR. Capitalizing on the limited computational cost of the MM/CG approach, this could pave the way for low-throughput virtual screening efforts. Furthermore, comparative studies of hGPCRs in multiple states, e.g., active/inactive, can be accomplished by setting up several simulations starting from different initial structures. Modeling of allosteric effects is currently not possible, because of the use of the Gō-like model for the part of the receptor pointing toward the intracellular side. We also plan to integrate the Amber-based MM/CG approach used by the webserver with the recently developed open boundary MM/CG for grand canonical simulations (Tarenzi et al., 2017, 2019), aiming at calculating ligand binding free energies. Further planned developments include the extension to other MD codes, the implementation of the CHARMM36m force field (Huang et al., 2017) for the MM/I regions, as well as the adaptation of the scheme to other membrane protein families, such as ion channels.
Data Availability Statement
The Hybrid MM/CG Webserver, available at https://mmcg.grs.kfa-juelich.de/, is free and open to all. There is no login requirement.
JS, MA-P, PC, and AG contributed to workflow design and documentation. JS developed the server. RR designed the interface on the GOMoDo webserver. JS, RR, MA-P, and AG tested the webserver. All authors participated in writing the article.
This work was supported by the “Ernesto Illy Foundation” (Trieste, Italy) and by the EBRAINS research infrastructure, funded from the European Union’s Horizon 2020 Framework Programme for Research and Innovation under the Specific Grant Agreement No. 945539 (Human Brain Project SGA3). Funding for the open access charge was provided by the Central Library of Forschungszentrum Jülich.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Emiliano Ippoliti for the setup of the server and Jonas Goßen for programming advice.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2020.576689/full#supplementary-material
AA, all-atom; CFF, caffeine; CG, coarse-grained; FDA, United States Food and Drug Administration; hA2AR, human adenosine 2A receptor; hGPCRs, human G protein-coupled receptors; MD, molecular dynamics; MM, molecular mechanics; MM/CG, molecular mechanics/coarse-grained.
- ^ https://gpcrdb.org/structure/statistics as of June 2020.
- ^ the residues whose mutation is known experimentally to affect binding
- ^ gomodo.grs.kfa-juelich.de
- ^ mmcg.grs.kfa-juelich.de/documentation
- ^ mmcg.grs.kfa-juelich.de/results/+ID
- ^ djangoproject.com
- ^ getbootstrap.com
- ^ highcharts.com
- ^ https://mmcg.grs.kfa-juelich.de/
- ^ mmcg.grs.kfa-juelich.de/example_results/
- ^ mmcg.grs.kfa-juelich.de/download/install
Alfonso-Prieto, M., Navarini, L., and Carloni, P. (2019). Understanding ligand binding to G-protein coupled receptors using multiscale simulations. Front. Mol. Biosci. 6:29. doi: 10.3389/fmolb.2019.00029
Berendsen, H. J. C., van der Spoel, D., and van Drunen, R. (1995). GROMACS: a message-passing parallel molecular dynamics implementation. Comput. Phys. Commun. 91, 43–56. doi: 10.1016/0010-4655(95)00042-e
Beuming, T., and Sherman, W. (2012). Current assessment of docking into GPCR crystal structures and homology models: successes, challenges, and guidelines. J. Chem. Inf. Model 52, 3263–3277. doi: 10.1021/ci300411b
Busato, M., and Giorgetti, A. (2016). Structural modeling of G-protein coupled receptors: An overview on automatic web-servers. Int. J. Biochem. Cell Biol. 77(Pt. B), 264–274. doi: 10.1016/j.biocel.2016.04.004
Capaldi, S., Suku, E., Antolini, M., Di Giacobbe, M., Giorgetti, A., and Buffelli, M. (2018). Allosteric sodium binding cavity in GPR3: a novel player in modulation of Aβ production. Sci. Rep. 8:11102. doi: 10.1038/s41598-018-29475-7
Cavasotto, C. N., and Palomba, D. (2015). Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models. Chem Commun. 51, 13576–13594. doi: 10.1039/c5cc05050b
Chen, J.-F., Xu, K., Petzer, J. P., Staal, R., Xu, Y.-H., Beilstein, M., et al. (2001). Neuroprotection by caffeine and A2A adenosine receptor inactivation in a model of Parkinson’s disease. J. Neurosci. 21, RC143–RC143. doi: 10.1523/JNEUROSCI.21-10-j0001.2001
Damre, M., Marchetto, A., and Giorgetti, A. (2019). MERMAID: dedicated web server to prepare and run coarse-grained membrane protein dynamics. Nucleic Acids Res. 47, W456–W461. doi: 10.1093/nar/gkz416
Esguerra, M., Siretskiy, A., Bello, X., Sallander, J., and Gutiérrez-de-Terán, H. (2016). GPCR-ModSim: a comprehensive web based solution for modeling G-protein coupled receptors. Nucleic Acids Res. 44, W455–W462. doi: 10.1093/nar/gkw403
Fierro, F., Giorgetti, A., Carloni, P., Meyerhof, W., and Alfonso-Prieto, M. (2019). Dual binding mode of “bitter sugars” to their human bitter taste receptor target. Sci. Rep. 9:8437. doi: 10.1038/s41598-019-44805-z
Fierro, F., Suku, E., Alfonso-Prieto, M., Giorgetti, A., Cichon, S., and Carloni, P. (2017). Agonist binding to chemosensory receptors: a systematic bioinformatics analysis. Front. Mol. Biosci. 4:63. doi: 10.3389/fmolb.2017.00063
Hauser, A. S., Attwood, M. M., Rask-Andersen, M., Schiöth, H. B., and Gloriam, D. E. (2017). Trends in GPCR drug discovery: new agents, targets and indications. Nat. Rev. Drug Discov. 16, 829–842. doi: 10.1038/nrd.2017.178
Heifetz, A., James, T., Morao, I., Bodkin, M. J., and Biggin, P. C. (2016). Guiding lead optimization with GPCR structure modeling and molecular dynamics. Curr. Opin. Pharmacol. 30, 14–21. doi: 10.1016/j.coph.2016.06.004
Hess, B., Kutzner, C., van der Spoel, D., and Lindahl, E. (2008). GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 4, 435–447. doi: 10.1021/ct700301q
Huang, J., Rauscher, S., Nawrocki, G., Ran, T., Feig, M., de Groot, B. L., et al. (2017). CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods 14, 71–73. doi: 10.1038/nmeth.4067
Jakalian, A., Jack, D. B., and Bayly, C. I. (2004). Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation. J. Comput. Chem. 23, 1623–1641. doi: 10.1002/jcc.10128
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., and Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926–935. doi: 10.1063/1.445869
Kufareva, I., Katritch, V., Stevens, R. C., and Abagyan, R. (2014). Advances in GPCR modeling evaluated by the GPCR Dock 2013 assessment: meeting new challenges. Structure 22, 1120–1139. doi: 10.1016/j.str.2014.06.012
Launay, G., Téletchéa, S., Wade, F., Pajot-Augy, E., Gibrat, J.-F., and Sanz, G. (2012). Automatic modeling of mammalian olfactory receptors and docking of odorants. Protein Eng. Des. Sel. 25, 377–386. doi: 10.1093/protein/gzs037
Lee, J., Cheng, X., Swails, J. M., Yeom, M. S., Eastman, P. K., Lemkul, J. A., et al. (2016). CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field. J. Chem. Theory Comput. 12, 405–413. doi: 10.1021/acs.jctc.5b00935
Leguèbe, M., Nguyen, C., Capece, L., Hoang, Z., Giorgetti, A., and Carloni, P. (2012). Hybrid molecular mechanics/coarse-grained simulations for structural prediction of G-protein coupled receptor/ligand complexes. PLoS One 7:e47332. doi: 10.1371/journal.pone.0047332
Lomize, M. A., Pogozheva, I. D., Joo, H., Mosberg, H. I., and Lomize, A. L. (2012). OPM database and PPM web server: resources for positioning of proteins in membranes. Nucleic Acids Res. 40, D370–D376. doi: 10.1093/nar/gkr703
Lupala, C. S., Rasaeifar, B., Gomez-Gutierrez, P., and Perez, J. J. (2018). Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling. J. Biomol. Struct. Dyn. 36, 2436–2448. doi: 10.1080/07391102.2017.1357503
Maier, J. A., Martinez, C., Kasavajhala, K., Wickstrom, L., Hauser, K. E., and Simmerling, C. (2015). ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 11, 3696–3713. doi: 10.1021/acs.jctc.5b00255
Marchiori, A., Capece, L., Giorgetti, A., Gasparini, P., Behrens, M., Carloni, P., et al. (2013). Coarse-grained/molecular mechanics of the TAS2R38 bitter taste receptor: experimentally-validated detailed structural prediction of agonist binding. PLoS One 8:e64675. doi: 10.1371/journal.pone.0064675
Miszta, P., Pasznik, P., Jakowiecki, J., Sztyler, A., Latek, D., and Filipek, S. (2018). GPCRM: a homology modeling web service with triple membrane-fitted quality assessment of GPCR models. Nucleic Acids Res. 46, W387–W395. doi: 10.1093/nar/gky429
Munk, C., Mutt, E., Isberg, V., Nikolajsen, L. F., Bibbe, J. M., Flock, T., et al. (2019). An online resource for GPCR structure determination and analysis. Nat. Methods 16, 151–162. doi: 10.1038/s41592-018-0302-x
Nazario, L. R., da Silva, R. S., and Bonan, C. D. (2017). Targeting adenosine signaling in parkinson’s disease: from pharmacological to non-pharmacological approaches. Front. Neurosci. 11:658. doi: 10.3389/fnins.2017.00658
Neri, M., Anselmi, C., Cascella, M., Maritan, A., and Carloni, P. (2005). Coarse-grained model of proteins incorporating atomistic detail of the active site. Phys. Rev. Lett. 95:218102. doi: 10.1103/PhysRevLett.95.218102
Neri, M., Baaden, M., Carnevale, V., Anselmi, C., Maritan, A., and Carloni, P. (2008). Microseconds dynamics simulations of the outer-membrane protease T. Biophys. J. 94, 71–78. doi: 10.1529/biophysj.107.116301
Newport, T. D., Sansom, M. S. P., and Stansfeld, P. J. (2019). The MemProtMD database: a resource for membrane-embedded protein structures and their lipid interactions. Nucleic Acids Res. 47, D390–D397. doi: 10.1093/nar/gky1047
Pándy-Szekeres, G., Munk, C., Tsonkov, T. M., Mordalski, S., Harpsøe, K., Hauser, A. S., et al. (2018). GPCRdb in 2018: adding GPCR structure models and ligands. Nucleic Acids Res. 46, D440–D446. doi: 10.1093/nar/gkx1109
Petzer, J. P., and Petzer, A. (2015). Caffeine as a lead compound for the design of therapeutic agents for the treatment of Parkinson’s disease. Curr. Med. Chem. 22, 975–988. doi: 10.2174/0929867322666141215160015
Qu, X., Wang, D., and Wu, B. (2020). “Progress in GPCR structure determination,” in GPCRs: Structure, Function, and Drug Discovery, ed. S.-H. Beata Jastrzebska Paul (Cambridge, MA: Academic Press). doi: 10.1016/C2018-0-00025-0
Rodríguez-Espigares, I., Torrens-Fontanals, M., Tiemann, J. K., Aranda-García, D., Ramírez-Anguita, J. M., Stepniewski, T. M., et al. (2020). GPCRmd uncovers the dynamics of the 3D-GPCRome. Nat. Methods 17, 777–787. doi: 10.1038/s41592-020-0884-y
Rose, P. W., Prlić, A., Altunkaya, A., Bi, C., Bradley, A. R., Christie, C. H., et al. (2016). The RCSB protein data bank: integrative view of protein, gene and 3D structural information. Nucleic Acids Res. 45:gkw1000. doi: 10.1093/nar/gkw1000
Sandal, M., Behrens, M., Brockhoff, A., Musiani, F., Giorgetti, A., Carloni, P., et al. (2015). Evidence for a transient additional ligand binding site in the TAS2R46 bitter taste receptor. J. Chem. Theory Comput. 11, 4439–4449. doi: 10.1021/acs.jctc.5b00472
Schneider, J., Korshunova, K., Musiani, F., Alfonso-Prieto, M., Giorgetti, A., and Carloni, P. (2018). Predicting ligand binding poses for low-resolution membrane protein models: perspectives from multiscale simulations. Biochem. Biophys. Res. Commun. 498, 366–374. doi: 10.1016/j.bbrc.2018.01.160
Schneider, J., Korshunova, K., Si Chaib, Z., Giorgetti, A., Alfonso-Prieto, M., and Carloni, P. (2020). Ligand pose predictions for human G protein-coupled receptors: insights from the amber-based hybrid molecular mechanics/coarse-grained approach. J. Chem. Inf. Model. doi: 10.1021/acs.jcim.0c00661 [Epub ahead of print].
Schott-Verdugo, S., and Gohlke, H. (2019). PACKMOL-Memgen: a simple-to-use, generalized workflow for membrane-protein-lipid-bilayer system building. J. Chem. Inf. Model. 59, 2522–2528. doi: 10.1021/acs.jcim.9b00269
Sonsalla, P. K., Wong, L.-Y., Harris, S. L., Richardson, J. R., Khobahy, I., Li, W., et al. (2012). Delayed caffeine treatment prevents nigral dopamine neuron loss in a progressive rat model of Parkinson’s disease. Exp. Neurol. 234, 482–487. doi: 10.1016/j.expneurol.2012.01.022
Tarenzi, T., Calandrini, V., Potestio, R., and Carloni, P. (2019). Open-boundary molecular mechanics/coarse-grained framework for simulations of low-resolution G-protein-coupled receptor–ligand complexes. J. Chem. Theory Comput. 15, 2101–2109. doi: 10.1021/acs.jctc.9b00040
Tarenzi, T., Calandrini, V., Potestio, R., Giorgetti, A., and Carloni, P. (2017). Open boundary simulations of proteins and their hydration shells by Hamiltonian adaptive resolution scheme. J. Chem. Theory Comput. 13, 5647–5657. doi: 10.1021/acs.jctc.7b00508
Worth, C. L., Kreuchwig, F., Tiemann, J. K. S., Kreuchwig, A., Ritschel, M., Kleinau, G., et al. (2017). GPCR-SSFE 2.0-a fragment-based molecular modeling web tool for Class A G-protein coupled receptors. Nucleic Acids Res. 45, W408–W415. doi: 10.1093/nar/gkx399
Keywords: MM/CG, molecular mechanics, coarse-grained, hybrid methods, webserver, G protein-coupled receptor, ligand, molecular dynamics simulation
Citation: Schneider J, Ribeiro R, Alfonso-Prieto M, Carloni P and Giorgetti A (2020) Hybrid MM/CG Webserver: Automatic Set Up of Molecular Mechanics/Coarse-Grained Simulations for Human G Protein-Coupled Receptor/Ligand Complexes. Front. Mol. Biosci. 7:576689. doi: 10.3389/fmolb.2020.576689
Received: 26 June 2020; Accepted: 13 August 2020;
Published: 04 September 2020.
Edited by:Sergio Decherchi, Italian Institute of Technology (IIT), Italy
Reviewed by:Yong Wang, University of Copenhagen, Denmark
Giovanni Grazioso, University of Milan, Italy
Copyright © 2020 Schneider, Ribeiro, Alfonso-Prieto, Carloni and Giorgetti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.