1
|
Zhang M, Chen T, Lu X, Lan X, Chen Z, Lu S. G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery. Signal Transduct Target Ther 2024; 9:88. [PMID: 38594257 PMCID: PMC11004190 DOI: 10.1038/s41392-024-01803-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 02/19/2024] [Accepted: 03/13/2024] [Indexed: 04/11/2024] Open
Abstract
G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.
Collapse
Affiliation(s)
- Mingyang Zhang
- Key Laboratory of Protection, Development and Utilization of Medicinal Resources in Liupanshan Area, Ministry of Education, Peptide & Protein Drug Research Center, School of Pharmacy, Ningxia Medical University, Yinchuan, 750004, China
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Ting Chen
- Department of Cardiology, Changzheng Hospital, Affiliated to Naval Medical University, Shanghai, 200003, China
| | - Xun Lu
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Xiaobing Lan
- Key Laboratory of Protection, Development and Utilization of Medicinal Resources in Liupanshan Area, Ministry of Education, Peptide & Protein Drug Research Center, School of Pharmacy, Ningxia Medical University, Yinchuan, 750004, China
| | - Ziqiang Chen
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, 200433, China.
| | - Shaoyong Lu
- Key Laboratory of Protection, Development and Utilization of Medicinal Resources in Liupanshan Area, Ministry of Education, Peptide & Protein Drug Research Center, School of Pharmacy, Ningxia Medical University, Yinchuan, 750004, China.
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China.
| |
Collapse
|
2
|
Jardón-Valadez E, Ulloa-Aguirre A. Tracking conformational transitions of the gonadotropin hormone receptors in a bilayer of (SDPC) poly-unsaturated lipids from all-atom molecular dynamics simulations. PLoS Comput Biol 2024; 20:e1011415. [PMID: 38206994 PMCID: PMC10807830 DOI: 10.1371/journal.pcbi.1011415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 01/24/2024] [Accepted: 12/15/2023] [Indexed: 01/13/2024] Open
Abstract
Glycoprotein hormone receptors [thyrotropin (TSHR), luteinizing hormone/chorionic gonadotropin (LHCGR), and follicle stimulating hormone (FSHR) receptors] are rhodopsin-like G protein-coupled receptors. These receptors display common structural features including a prominent extracellular domain with leucine-rich repeats (LRR) stabilized by β-sheets and a long and flexible loop known as the hinge region (HR), and a transmembrane (TM) domain with seven α-helices interconnected by intra- and extracellular loops. Binding of the ligand to the LRR resembles a hand coupling transversally to the α- and β-subunits of the hormone, with the thumb being the HR. The structure of the FSH-FSHR complex suggests an activation mechanism in which Y335 at the HR binds into a pocket between the α- and β-chains of the hormone, leading to an adjustment of the extracellular loops. In this study, we performed molecular dynamics (MD) simulations to identify the conformational changes of the FSHR and LHCGR. We set up a FSHR structure as predicted by AlphaFold (AF-P23945); for the LHCGR structure we took the cryo-electron microscopy structure for the active state (PDB:7FII) as initial coordinates. Specifically, the flexibility of the HR domain and the correlated motions of the LRR and TM domain were analyzed. From the conformational changes of the LRR, TM domain, and HR we explored the conformational landscape by means of MD trajectories in all-atom approximation, including a membrane of polyunsaturated phospholipids. The distances and procedures here defined may be useful to propose reaction coordinates to describe diverse processes, such as the active-to-inactive transition, and to identify intermediaries suited for allosteric regulation and biased binding to cellular transducers in a selective activation strategy.
Collapse
Affiliation(s)
- Eduardo Jardón-Valadez
- Departamento de Recursos de la Tierra, Unidad Lerma, Universidad Autónoma Metropolitana, Lerma de Villada, Estado de México, Mexico
| | - Alfredo Ulloa-Aguirre
- Instituto Nacional de Ciencias Medicas y Nutrición “Salvador Zubiran”. Mexico City, Mexico
- Red de Apoyo a la Investigación, Universidad Nacional Autónoma de México. Mexico City, Mexico
| |
Collapse
|
3
|
Do HN, Wang J, Miao Y. Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors. JACS AU 2023; 3:3165-3180. [PMID: 38034960 PMCID: PMC10685416 DOI: 10.1021/jacsau.3c00503] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/26/2023] [Revised: 10/10/2023] [Accepted: 10/18/2023] [Indexed: 12/02/2023]
Abstract
G-protein-coupled receptors (GPCRs) make up the largest superfamily of human membrane proteins and represent primary targets of ∼1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon the binding of positive and negative allosteric modulators (PAMs and NAMs). The mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy prOfiling Workflow (GLOW). GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence and absence of the modulator. DL and free energy calculations revealed significantly reduced dynamic fluctuations and conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G-protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes. Therefore, GPCR allostery exhibits a dynamic "conformational selection" mechanism. In the absence of available modulator-bound structures as for most current GPCRs, it is critical to use a structural ensemble of representative GPCR conformations rather than a single structure for compound docking ("ensemble docking"), which will potentially improve structure-based design of novel allosteric drugs of GPCRs.
Collapse
Affiliation(s)
| | - Jinan Wang
- Computational Biology Program
and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | | |
Collapse
|
4
|
Yang PC, Rose A, DeMarco KR, Dawson JRD, Han Y, Jeng MT, Harvey RD, Santana LF, Ripplinger CM, Vorobyov I, Lewis TJ, Clancy CE. A multiscale predictive digital twin for neurocardiac modulation. J Physiol 2023; 601:3789-3812. [PMID: 37528537 PMCID: PMC10528740 DOI: 10.1113/jp284391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 07/11/2023] [Indexed: 08/03/2023] Open
Abstract
Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the β-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.
Collapse
Affiliation(s)
- Pei-Chi Yang
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | - Adam Rose
- Department of Mathematics, University of California Davis, Davis, CA
| | - Kevin R. DeMarco
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | - John R. D. Dawson
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | - Yanxiao Han
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | - Mao-Tsuen Jeng
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | | | - L. Fernando Santana
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | | | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
| | - Timothy J. Lewis
- Department of Mathematics, University of California Davis, Davis, CA
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA
- Center for Precision Medicine and Data Science, University of California Davis, Sacramento, CA
| |
Collapse
|
5
|
Han Y, Dawson JRD, DeMarco KR, Rouen KC, Bekker S, Yarov-Yarovoy V, Clancy CE, Xiang YK, Vorobyov I. Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein. Proc Natl Acad Sci U S A 2023; 120:e2215916120. [PMID: 36853938 PMCID: PMC10013855 DOI: 10.1073/pnas.2215916120] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2022] [Accepted: 01/27/2023] [Indexed: 03/01/2023] Open
Abstract
G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-β2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and β2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.
Collapse
Affiliation(s)
- Yanxiao Han
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
| | - John R. D. Dawson
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
- Biophysics Graduate Group, University of California, Davis, CA95616
| | - Kevin R. DeMarco
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
| | - Kyle C. Rouen
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
- Biophysics Graduate Group, University of California, Davis, CA95616
| | - Slava Bekker
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
- Department of Science and Engineering, American River College, Sacramento, CA95841
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
- Department of Anesthesiology and Pain Medicine, University of California, Davis, CA95616
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
- Department of Pharmacology, University of California, Davis, CA95616
| | - Yang K. Xiang
- Department of Pharmacology, University of California, Davis, CA95616
- VA Northern California Health Care System, Mather, CA95655
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California, Davis, CA95616
- Department of Pharmacology, University of California, Davis, CA95616
| |
Collapse
|
6
|
Do H, Wang J, Miao Y. Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors. RESEARCH SQUARE 2023:rs.3.rs-2543463. [PMID: 36865316 PMCID: PMC9980202 DOI: 10.21203/rs.3.rs-2543463/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
Abstract
G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~ 1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.
Collapse
|
7
|
Do HN, Wang J, Miao Y. Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.15.524128. [PMID: 36711515 PMCID: PMC9882226 DOI: 10.1101/2023.01.15.524128] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.
Collapse
Affiliation(s)
- Hung N. Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| |
Collapse
|
8
|
Elgeti M, Hubbell WL. DEER Analysis of GPCR Conformational Heterogeneity. Biomolecules 2021; 11:778. [PMID: 34067265 PMCID: PMC8224605 DOI: 10.3390/biom11060778] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2021] [Revised: 05/18/2021] [Accepted: 05/19/2021] [Indexed: 02/06/2023] Open
Abstract
G protein-coupled receptors (GPCRs) represent a large class of transmembrane helical proteins which are involved in numerous physiological signaling pathways and therefore represent crucial pharmacological targets. GPCR function and the action of therapeutic molecules are defined by only a few parameters, including receptor basal activity, ligand affinity, intrinsic efficacy and signal bias. These parameters are encoded in characteristic receptor conformations existing in equilibrium and their populations, which are thus of paramount interest for the understanding of receptor (mal-)functions and rational design of improved therapeutics. To this end, the combination of site-directed spin labeling and EPR spectroscopy, in particular double electron-electron resonance (DEER), is exceedingly valuable as it has access to sub-Angstrom spatial resolution and provides a detailed picture of the number and populations of conformations in equilibrium. This review gives an overview of existing DEER studies on GPCRs with a focus on the delineation of structure/function frameworks, highlighting recent developments in data analysis and visualization. We introduce "conformational efficacy" as a parameter to describe ligand-specific shifts in the conformational equilibrium, taking into account the loose coupling between receptor segments observed for different GPCRs using DEER.
Collapse
Affiliation(s)
- Matthias Elgeti
- Jules Stein Eye Institute and Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| | - Wayne L. Hubbell
- Jules Stein Eye Institute and Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| |
Collapse
|
9
|
Fleetwood O, Carlsson J, Delemotte L. Identification of ligand-specific G protein-coupled receptor states and prediction of downstream efficacy via data-driven modeling. eLife 2021; 10:60715. [PMID: 33506760 PMCID: PMC7886328 DOI: 10.7554/elife.60715] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Accepted: 01/27/2021] [Indexed: 12/22/2022] Open
Abstract
Ligand binding stabilizes different G protein-coupled receptor states via a complex allosteric process that is not completely understood. Here, we have derived free energy landscapes describing activation of the β2 adrenergic receptor bound to ligands with different efficacy profiles using enhanced sampling molecular dynamics simulations. These reveal shifts toward active-like states at the Gprotein-binding site for receptors bound to partial and full agonists, and that the ligands modulate the conformational ensemble of the receptor by tuning protein microswitches. We indeed find an excellent correlation between the conformation of the microswitches close to the ligand binding site and in the transmembrane region and experimentally reported cyclic adenosine monophosphate signaling responses. Dimensionality reduction further reveals the similarity between the unique conformational states induced by different ligands, and examining the output of classifiers highlights two distant hotspots governing agonism on transmembrane helices 5 and 7.
Collapse
Affiliation(s)
- Oliver Fleetwood
- Science for Life Laboratory, Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Lucie Delemotte
- Science for Life Laboratory, Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
| |
Collapse
|
10
|
|
11
|
Renault P, Giraldo J. Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors. Int J Mol Sci 2020; 22:ijms22010187. [PMID: 33375427 PMCID: PMC7795036 DOI: 10.3390/ijms22010187] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2020] [Revised: 12/22/2020] [Accepted: 12/24/2020] [Indexed: 01/14/2023] Open
Abstract
G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.
Collapse
Affiliation(s)
- Pedro Renault
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain;
- Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 08193 Bellaterra, Spain
| | - Jesús Giraldo
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain;
- Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 08193 Bellaterra, Spain
- Correspondence:
| |
Collapse
|
12
|
Salmaso V, Jacobson KA. In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors. Biomolecules 2020; 10:E812. [PMID: 32466404 PMCID: PMC7356333 DOI: 10.3390/biom10060812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 05/21/2020] [Accepted: 05/22/2020] [Indexed: 12/14/2022] Open
Abstract
Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.
Collapse
Affiliation(s)
| | - Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA;
| |
Collapse
|
13
|
Ribeiro JML, Filizola M. Insights From Molecular Dynamics Simulations of a Number of G-Protein Coupled Receptor Targets for the Treatment of Pain and Opioid Use Disorders. Front Mol Neurosci 2019; 12:207. [PMID: 31507375 PMCID: PMC6716474 DOI: 10.3389/fnmol.2019.00207] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2019] [Accepted: 08/07/2019] [Indexed: 01/20/2023] Open
Abstract
Effective treatments for pain management remain elusive due to the dangerous side-effects of current gold-standard opioid analgesics, including the respiratory depression that has led to skyrocketing death rates from opioid overdoses over the past decade. In an attempt to address the horrific opioid crisis worldwide, the National Institute on Drug Abuse has recently proposed boosting research on specific pharmacological mechanisms mediated by a number of G protein-coupled receptors (GPCRs). This research is expected to expedite the discovery of medications for opioid overdose and opioid use disorders, leading toward a safer and more effective treatment of pain. Here, we review mechanistic insights from recent all-atom molecular dynamics simulations of a specific subset of GPCRs for which high-resolution experimental structures are available, including opioid, cannabinoid, orexin, metabotropic glutamate, and dopamine receptor subtypes.
Collapse
Affiliation(s)
- João Marcelo Lamim Ribeiro
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States
| | - Marta Filizola
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States
| |
Collapse
|
14
|
Editorial overview: Theory and simulation: demystifying GPCRs – structure, function and drug design. Curr Opin Struct Biol 2019; 55:vi-viii. [DOI: 10.1016/j.sbi.2019.07.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|