Hybrid MM/CG Webserver: Automatic Set Up of Molecular Mechanics/Coarse-Grained Simulations for Human G Protein-Coupled Receptor/Ligand Complexes

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.

Introduction

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% hGPCRs¹ (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 residues² (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

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).

TABLE 1

Table 1. Available online resources for MD simulations of GPCRs.

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 webserver³ 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

Input

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.

Setup Modes

The Hybrid MM/CG webserver offers two system setup modes:

Automatic Preparation

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 webserver³ (Sandal et al., 2013). All simulation parameters are set according to default values that are defined in the documentation⁴.

Interactive Preparation

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 webserver⁴.

Workflow

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

Figure 2. Webserver workflow. Interactive mode with manual intervention is presented on the left and fast automatic mode on the right.

FIGURE 3

Figure 3. Sample 3D visualizations created by the webserver using as example the hA_2AR/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.

Output

Results can be viewed and downloaded for 2 weeks (renewal possible) by bookmarking the link or alternatively by using the corresponding ID⁵. 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).

Server Architecture

Front-End

The web interface was developed with the DJANGO 3.0 Web framework⁶ and designed with the Bootstrap 4 front-end open-source toolkit⁷. Interactive 3D visualizations were implemented with the NGL library for molecular visualization (Rose and Hildebrand, 2015) and the interactive graphs were created using the Highcharts SVG-based JavaScript charting library⁸.

Back-End

For historical reasons, different programming and scripting languages are used for the workflow such as Bash, Python, Perl, Tcl, AWK, and JavaScript. Furthermore, different software packages are used in the preparation: a patched GROMACS 4.5.1 version (Berendsen et al., 1995; Hess et al., 2008), VMD 1.9.3 (Humphrey et al., 1996), PPM (Lomize et al., 2012), lovoalign (Martínez et al., 2007), ACPYPE Rev: 10101 (Sousa da Silva and Vranken, 2012), Antechamber (Wang et al., 2001) from AmberTools16 and OpenBabel 2.3.2 (O’Boyle et al., 2011). All parts were combined using the DJANGO 3.0 Web framework⁶.

Results

As an application case, we considered the human adenosine 2A receptor (hA_2AR) 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 hA_2AR/CFF complex used for the Example can be downloaded from the same page. The demonstration can be run by simply starting the preparation workflow⁹ without uploading any PDB file. Precomputed results for this example case can be found on the server¹⁰ and selected parts of these results can be seen in Figure 3. Figures 3A–C show the whole hA_2A 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).

Conclusion

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 hA_2AR/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 download¹¹) 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.

Author Contributions

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.

Funding

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.

Acknowledgments

We thank Emiliano Ippoliti for the setup of the server and Jonas Goßen for programming advice.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2020.576689/full#supplementary-material

Abbreviations

AA, all-atom; CFF, caffeine; CG, coarse-grained; FDA, United States Food and Drug Administration; hA_2AR, human adenosine 2A receptor; hGPCRs, human G protein-coupled receptors; MD, molecular dynamics; MM, molecular mechanics; MM/CG, molecular mechanics/coarse-grained.

Footnotes

1. ^ https://gpcrdb.org/structure/statistics as of June 2020.
2. ^ the residues whose mutation is known experimentally to affect binding
3. ^ gomodo.grs.kfa-juelich.de
4. ^ mmcg.grs.kfa-juelich.de/documentation
5. ^ mmcg.grs.kfa-juelich.de/results/+ID
6. ^ djangoproject.com
7. ^ getbootstrap.com
8. ^ highcharts.com
9. ^ https://mmcg.grs.kfa-juelich.de/
10. ^ mmcg.grs.kfa-juelich.de/example_results/
11. ^ mmcg.grs.kfa-juelich.de/download/install

References

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The Protein Data Bank. Nucl. Acids Res. 28, 235–242. doi: 10.1093/nar/28.1.235

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Case, D. A., Belfon, K., Ben-Shalom, I. Y., Brozell, S. R., Cerutti, D. S., Cheatham, T.E. III., et al. (2020). Amber 2020. San Francisco, CL: University of California.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Gō, N., and Abe, H. (1981). Noninteracting local-structure model of folding and unfolding transition in globular proteins. I. Formulation. Biopolymers 20, 991–1011. doi: 10.1002/bip.1981.360200511

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Hauser, A. S., Chavali, S., Masuho, I., Jahn, L. J., Martemyanov, K. A., Gloriam, D. E., et al. (2018). Pharmacogenomics of GPCR drug targets. Cell 172, 41–54.e19. doi: 10.1016/j.cell.2017.11.033

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Humphrey, W., Dalke, A., and Schulten, K. (1996). VMD: Visual molecular dynamics. J. Mol. Graph 14, 33–38. doi: 10.1016/0263-7855(96)00018-5

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Jo, S., Kim, T., Iyer, V. G., and Im, W. (2008). CHARMM-GUI: a web-based graphical user interface for CHARMM. J. Comput. Chem. 29, 1859–1865. doi: 10.1002/jcc.20945

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, G. R., and Seok, C. (2016). Galaxy7TM: flexible GPCR–ligand docking by structure refinement. Nucleic Acids Res. 44, W502–W506. doi: 10.1093/nar/gkw360

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez, L., Andreani, R., and Martínez, J. M. (2007). Convergent algorithms for protein structural alignment. BMC Bioinformatics 8:306. doi: 10.1186/1471-2105-8-306

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Boyle, N. M., Banck, M., James, C. A., Morley, C., Vandermeersch, T., and Hutchison, G. R. (2011). Open babel: an open chemical toolbox. J. Cheminform. 3:33. doi: 10.1186/1758-2946-3-33

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Pons, J.-L., and Labesse, G. (2009). @ TOME-2: a new pipeline for comparative modeling of protein–ligand complexes. Nucleic Acids Res. 37(Suppl._2), W485–W491. doi: 10.1093/nar/gkp368

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Rose, A. S., and Hildebrand, P. W. (2015). NGL Viewer: a web application for molecular visualization. Nucleic Acids Res. 43, W576–W579. doi: 10.1093/nar/gkv402

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandal, M., Duy, T. P., Cona, M., Zung, H., Carloni, P., Musiani, F., et al. (2013). GOMoDo: a GPCRs online modeling and docking webserver. PLoS One 8:e74092. doi: 10.1371/journal.pone.0074092

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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].

CrossRef Full Text | PubMed Abstract | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Sousa da Silva, A. W., and Vranken, W. F. (2012). ACPYPE - AnteChamber PYthon Parser interfacE. BMC Res. Notes 5:367. doi: 10.1186/1756-0500-5-367

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wang, W., Kollman, P. A., and Case, D. A. (2001). Antechamber: an accessory software package for molecular mechanical calculations. J. Am. Chem. Soc. 222:U403.

Google Scholar

Wang, J., Wolf, R. M., Caldwell, J. W., Kollman, P. A., and Case, D. A. (2004). Development and testing of a general amber force field. J. Comput. Chem. 25, 1157–1174. doi: 10.1002/jcc.20035

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar