1
|
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: 0] [Impact Index Per Article: 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
|
2
|
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
|
3
|
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
|
4
|
Abrol R, Serrano E, Santiago LJ. Development of enhanced conformational sampling methods to probe the activation landscape of GPCRs. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2021; 128:325-359. [PMID: 35034722 DOI: 10.1016/bs.apcsb.2021.11.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
G protein-coupled receptors (GPCRs) make up the largest superfamily of integral membrane proteins and play critical signal transduction roles in many physiological processes. Developments in molecular biology, biophysical, biochemical, pharmacological, and computational techniques aimed at these important therapeutic targets are beginning to provide unprecedented details on the structural as well as functional basis of their pleiotropic signaling mediated by G proteins, β arrestins, and other transducers. This pleiotropy presents a pharmacological challenge as the same ligand-receptor interaction can cause a therapeutic effect as well as an undesirable on-target side-effect through different downstream pathways. GPCRs don't function as simple binary on-off switches but as finely tuned shape-shifting machines described by conformational ensembles, where unique subsets of conformations may be responsible for specific signaling cascades. X-ray crystallography and more recently cryo-electron microscopy are providing snapshots of some of these functionally-important receptor conformations bound to ligands and/or transducers, which are being utilized by computational methods to describe the dynamic conformational energy landscape of GPCRs. In this chapter, we review the progress in computational conformational sampling methods based on molecular dynamics and discrete sampling approaches that have been successful in complementing biophysical and biochemical studies on these receptors in terms of their activation mechanisms, allosteric effects, actions of biased ligands, and effects of pathological mutations. Some of the sampled simulation time scales are beginning to approach receptor activation time scales. The list of conformational sampling methods and example uses discussed is not exhaustive but includes representative examples that have pushed the limits of classical molecular dynamics and discrete sampling methods to describe the activation energy landscape of GPCRs.
Collapse
Affiliation(s)
- Ravinder Abrol
- Department of Chemistry and Biochemistry, California State University, Northridge, CA, United States.
| | - Erik Serrano
- Department of Chemistry and Biochemistry, California State University, Northridge, CA, United States
| | - Luis Jaimes Santiago
- Department of Chemistry and Biochemistry, California State University, Northridge, CA, United States
| |
Collapse
|
5
|
Wang J, Miao Y. Recent advances in computational studies of GPCR-G protein interactions. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2019; 116:397-419. [PMID: 31036298 PMCID: PMC6986689 DOI: 10.1016/bs.apcsb.2018.11.011] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Protein-protein interactions are key in cellular signaling. G protein-coupled receptors (GPCRs), the largest superfamily of human membrane proteins, are able to transduce extracellular signals (e.g., hormones and neurotransmitters) to intracellular proteins, in particular the G proteins. Since GPCRs serve as primary targets of ~1/3 of currently marketed drugs, it is important to understand mechanisms of GPCR signaling in order to design selective and potent drug molecules. This chapter focuses on recent advances in computational studies of the GPCR-G protein interactions using bioinformatics, protein-protein docking and molecular dynamics simulation approaches.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States.
| |
Collapse
|
6
|
Saleh N, Kleinau G, Heyder N, Clark T, Hildebrand PW, Scheerer P. Binding, Thermodynamics, and Selectivity of a Non-peptide Antagonist to the Melanocortin-4 Receptor. Front Pharmacol 2018; 9:560. [PMID: 29910730 PMCID: PMC5992272 DOI: 10.3389/fphar.2018.00560] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2018] [Accepted: 05/10/2018] [Indexed: 11/22/2022] Open
Abstract
The melanocortin-4 receptor (MC4R) is a potential drug target for treatment of obesity, anxiety, depression, and sexual dysfunction. Crystal structures for MC4R are not yet available, which has hindered successful structure-based drug design. Using microsecond-scale molecular-dynamics simulations, we have investigated selective binding of the non-peptide antagonist MCL0129 to a homology model of human MC4R (hMC4R). This approach revealed that, at the end of a multi-step binding process, MCL0129 spontaneously adopts a binding mode in which it blocks the agonistic-binding site. This binding mode was confirmed in subsequent metadynamics simulations, which gave an affinity for human hMC4R that matches the experimentally determined value. Extending our simulations of MCL0129 binding to hMC1R and hMC3R, we find that receptor subtype selectivity for hMC4R depends on few amino acids located in various structural elements of the receptor. These insights may support rational drug design targeting the melanocortin systems.
Collapse
Affiliation(s)
- Noureldin Saleh
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Berlin, Germany.,Computational Modelling and Dynamics of Molecular Complexes, Berlin, Germany
| | - Gunnar Kleinau
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Berlin, Germany.,Group Protein X-ray Crystallography and Signal Transduction, Berlin, Germany
| | - Nicolas Heyder
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Berlin, Germany.,Group Protein X-ray Crystallography and Signal Transduction, Berlin, Germany
| | - Timothy Clark
- Computer-Chemie-Centrum, Department of Chemistry and Pharmacy, Friedrich-Alexander University Erlangen-Nürnberg, Erlangen, Germany
| | - Peter W Hildebrand
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Berlin, Germany.,Computational Modelling and Dynamics of Molecular Complexes, Berlin, Germany.,Institute of Medical Physics and Biophysics, Leipzig University, Leipzig, Germany
| | - Patrick Scheerer
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Berlin, Germany.,Group Protein X-ray Crystallography and Signal Transduction, Berlin, Germany
| |
Collapse
|
7
|
Mechanism of the G-protein mimetic nanobody binding to a muscarinic G-protein-coupled receptor. Proc Natl Acad Sci U S A 2018; 115:3036-3041. [PMID: 29507218 DOI: 10.1073/pnas.1800756115] [Citation(s) in RCA: 84] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Protein-protein binding is key in cellular signaling processes. Molecular dynamics (MD) simulations of protein-protein binding, however, are challenging due to limited timescales. In particular, binding of the medically important G-protein-coupled receptors (GPCRs) with intracellular signaling proteins has not been simulated with MD to date. Here, we report a successful simulation of the binding of a G-protein mimetic nanobody to the M2 muscarinic GPCR using the robust Gaussian accelerated MD (GaMD) method. Through long-timescale GaMD simulations over 4,500 ns, the nanobody was observed to bind the receptor intracellular G-protein-coupling site, with a minimum rmsd of 2.48 Å in the nanobody core domain compared with the X-ray structure. Binding of the nanobody allosterically closed the orthosteric ligand-binding pocket, being consistent with the recent experimental finding. In the absence of nanobody binding, the receptor orthosteric pocket sampled open and fully open conformations. The GaMD simulations revealed two low-energy intermediate states during nanobody binding to the M2 receptor. The flexible receptor intracellular loops contribute remarkable electrostatic, polar, and hydrophobic residue interactions in recognition and binding of the nanobody. These simulations provided important insights into the mechanism of GPCR-nanobody binding and demonstrated the applicability of GaMD in modeling dynamic protein-protein interactions.
Collapse
|
8
|
Accelerated structure-based design of chemically diverse allosteric modulators of a muscarinic G protein-coupled receptor. Proc Natl Acad Sci U S A 2016; 113:E5675-84. [PMID: 27601651 DOI: 10.1073/pnas.1612353113] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Design of ligands that provide receptor selectivity has emerged as a new paradigm for drug discovery of G protein-coupled receptors, and may, for certain families of receptors, only be achieved via identification of chemically diverse allosteric modulators. Here, the extracellular vestibule of the M2 muscarinic acetylcholine receptor (mAChR) is targeted for structure-based design of allosteric modulators. Accelerated molecular dynamics (aMD) simulations were performed to construct structural ensembles that account for the receptor flexibility. Compounds obtained from the National Cancer Institute (NCI) were docked to the receptor ensembles. Retrospective docking of known ligands showed that combining aMD simulations with Glide induced fit docking (IFD) provided much-improved enrichment factors, compared with the Glide virtual screening workflow. Glide IFD was thus applied in receptor ensemble docking, and 38 top-ranked NCI compounds were selected for experimental testing. In [(3)H]N-methylscopolamine radioligand dissociation assays, approximately half of the 38 lead compounds altered the radioligand dissociation rate, a hallmark of allosteric behavior. In further competition binding experiments, we identified 12 compounds with affinity of ≤30 μM. With final functional experiments on six selected compounds, we confirmed four of them as new negative allosteric modulators (NAMs) and one as positive allosteric modulator of agonist-mediated response at the M2 mAChR. Two of the NAMs showed subtype selectivity without significant effect at the M1 and M3 mAChRs. This study demonstrates an unprecedented successful structure-based approach to identify chemically diverse and selective GPCR allosteric modulators with outstanding potential for further structure-activity relationship studies.
Collapse
|
9
|
Multiscale design of coarse-grained elastic network-based potentials for the μ opioid receptor. J Mol Model 2016; 22:227. [DOI: 10.1007/s00894-016-3092-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2016] [Accepted: 08/08/2016] [Indexed: 01/10/2023]
|
10
|
G-protein coupled receptors: advances in simulation and drug discovery. Curr Opin Struct Biol 2016; 41:83-89. [PMID: 27344006 DOI: 10.1016/j.sbi.2016.06.008] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2016] [Accepted: 06/07/2016] [Indexed: 11/21/2022]
Abstract
G-protein coupled receptors (GPCRs), the largest family of human membrane proteins, mediate cellular signaling and represent primary targets of about one third of currently marketed drugs. GPCRs undergo highly dynamic structural transitions during signal transduction, from binding of extracellular ligands to coupling with intracellular effector proteins. Molecular dynamics (MD) simulations have been utilized to investigate GPCR signaling mechanisms (such as pathways of ligand binding and receptor activation/deactivation) and to design novel small-molecule drug candidates. Future research directions point towards modeling cooperative binding of multiple orthosteric and allosteric ligands to GPCRs, GPCR oligomerization and interactions of GPCRs with different intracellular signaling proteins. Through methodological and supercomputing advances, MD simulations will continue to provide important insights into GPCR signaling mechanisms and further facilitate structure-based drug design.
Collapse
|
11
|
Abstract
Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.
Collapse
|
12
|
Chen Y, Tang H, Seibel W, Papoian R, Li X, Lambert NA, Palczewski K. A High-Throughput Drug Screening Strategy for Detecting Rhodopsin P23H Mutant Rescue and Degradation. Invest Ophthalmol Vis Sci 2015; 56:2553-67. [PMID: 25783607 DOI: 10.1167/iovs.14-16298] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
PURPOSE Inherent instability of the P23H mutant opsin accounts for approximately 10% of autosomal dominant retinitis pigmentosa cases. Our purpose was to develop an overall set of reliable screening strategies to assess if either stabilization or enhanced degradation of mutant rhodopsin could rescue rod photoreceptors expressing this mutant protein. These strategies promise to reveal active compounds and clarify molecular mechanisms of biologically important processes, such as inhibition of target degradation or enhanced target folding. METHODS Cell-based bioluminescence reporter assays were developed and validated for high-throughput screening (HTS) of compounds that promote either stabilization or degradation of P23H mutant opsin. Such assays were further complemented by immunoblotting and image-based analyses. RESULTS Two stabilization assays of P23H mutant opsin were developed and validated, one based on β-galactosidase complementarity and a second assay involving bioluminescence resonance energy transfer (BRET) technology. Moreover, two additional assays evaluating mutant protein degradation also were employed, one based on the disappearance of luminescence and another employing the ALPHA immunoassay. Imaging of cells revealed the cellular localization of mutant rhodopsin, whereas immunoblots identified changes in the aggregation and glycosylation of P23H mutant opsin. CONCLUSIONS Our findings indicate that these initial HTS and following assays can identify active therapeutic compounds, even for difficult targets such as mutant rhodopsin. The assays are readily scalable and their function has been proven with model compounds. High-throughput screening, supported by automated imaging and classic immunoassays, can further characterize multiple steps and pathways in the biosynthesis and degradation of this essential visual system protein.
Collapse
Affiliation(s)
- Yuanyuan Chen
- Department of Pharmacology Case Western Reserve University, Cleveland, Ohio, United States
| | - Hong Tang
- Drug Discovery Center, College of Medicine, University of Cincinnati, Cincinnati, Ohio, United States
| | - William Seibel
- Drug Discovery Center, College of Medicine, University of Cincinnati, Cincinnati, Ohio, United States
| | - Ruben Papoian
- Drug Discovery Center, College of Medicine, University of Cincinnati, Cincinnati, Ohio, United States
| | - Xiaoyu Li
- Department of Pharmacology Case Western Reserve University, Cleveland, Ohio, United States
| | - Nevin A Lambert
- Department of Pharmacology and Toxicology, Georgia Regents University, Augusta, Georgia, United States
| | - Krzysztof Palczewski
- Department of Pharmacology Case Western Reserve University, Cleveland, Ohio, United States
| |
Collapse
|
13
|
Meng XY, Mezei M, Cui M. Computational approaches for modeling GPCR dimerization. Curr Pharm Biotechnol 2015; 15:996-1006. [PMID: 25307013 DOI: 10.2174/1389201015666141013102515] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2014] [Revised: 07/09/2014] [Accepted: 08/18/2014] [Indexed: 11/22/2022]
Abstract
Growing experimental evidences suggest that dimerization and oligomerization are important for G Protein- Coupled Receptors (GPCRs) function. The detailed structural information of dimeric/oligomeric GPCRs would be very important to understand their function. Although it is encouraging that recently several experimental GPCR structures in oligomeric forms have appeared, experimental determination of GPCR structures in oligomeric forms is still a big challenge, especially in mimicking the membrane environment. Therefore, development of computational approaches to predict dimerization of GPCRs will be highly valuable. In this review, we summarize computational approaches that have been developed and used for modeling of GPCR dimerization. In addition, we introduce a novel two-dimensional Brownian Dynamics based protein docking approach, which we have recently adapted, for GPCR dimer prediction.
Collapse
Affiliation(s)
| | | | - Meng Cui
- Institute of Quantitative Biology and Medicine, Jiangsu Provincial Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China.
| |
Collapse
|
14
|
Miao Y, Nichols SE, McCammon JA. Free energy landscape of G-protein coupled receptors, explored by accelerated molecular dynamics. Phys Chem Chem Phys 2014; 16:6398-406. [PMID: 24445284 PMCID: PMC3960983 DOI: 10.1039/c3cp53962h] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2013] [Accepted: 01/14/2014] [Indexed: 11/21/2022]
Abstract
G-protein coupled receptors (GPCRs) mediate cellular responses to various hormones and neurotransmitters and are important targets for treating a wide spectrum of diseases. They are known to adopt multiple conformational states (e.g., inactive, intermediate and active) during their modulation of various cell signaling pathways. Here, the free energy landscape of GPCRs is explored using accelerated molecular dynamics (aMD) simulations as demonstrated on the M2 muscarinic receptor, a key GPCR that regulates human heart rate and contractile forces of cardiomyocytes. Free energy profiles of important structural motifs that undergo conformational transitions upon GPCR activation and allosteric signaling are analyzed in detail, including the Arg(3.50)-Glu(6.30) ionic lock, the Trp(6.48) toggle switch and the hydrogen interactions between Tyr(5.58)-Tyr(7.53).
Collapse
Affiliation(s)
- Yinglong Miao
- Howard Hughes Medical Institute , University of California at San Diego , La Jolla , CA 92093 , USA .
| | - Sara E. Nichols
- Department of Chemistry and Biochemistry , University of California at San Diego , La Jolla , CA 92093 , USA .
- Department of Pharmacology , University of California at San Diego , La Jolla , CA 92093 , USA
| | - J. Andrew McCammon
- Howard Hughes Medical Institute , University of California at San Diego , La Jolla , CA 92093 , USA .
- Department of Chemistry and Biochemistry , University of California at San Diego , La Jolla , CA 92093 , USA .
- Department of Pharmacology , University of California at San Diego , La Jolla , CA 92093 , USA
| |
Collapse
|
15
|
Senćanski M, Došen-Mićović L. In Silico Study of the Structurally Similar ORL1 Receptor Agonist and Antagonist Pairs Reveal Possible Mechanism of Receptor Activation. Protein J 2014; 33:231-42. [DOI: 10.1007/s10930-014-9555-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
16
|
Campomanes P, Neri M, Horta BAC, Röhrig UF, Vanni S, Tavernelli I, Rothlisberger U. Origin of the Spectral Shifts among the Early Intermediates of the Rhodopsin Photocycle. J Am Chem Soc 2014; 136:3842-51. [DOI: 10.1021/ja411303v] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Pablo Campomanes
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Marilisa Neri
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Bruno A. C. Horta
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Ute F. Röhrig
- Molecular Modeling
Group, Swiss Institute of
Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Stefano Vanni
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Ivano Tavernelli
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|
17
|
Arnatt CK, Zhang Y. G Protein-Coupled Estrogen Receptor (GPER) Agonist Dual Binding Mode Analyses toward Understanding of its Activation Mechanism: A Comparative Homology Modeling Approach. Mol Inform 2013; 32:647-658. [PMID: 26229572 PMCID: PMC4517684 DOI: 10.1002/minf.201200136] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
G protein-coupled estrogen receptor (GPER) has been shown to be important in several disease states such as estrogen sensitive cancers. While several selective ligands have been identified for the receptor, little is known about how they interact with GPER and how their structures influence their activity. Specifically, within one series of ligands, whose structure varied only at one position, the replacement of a hydrogen atom with an acetyl group changed a potent antagonist into a potent agonist. In this study, two GPER homology models were constructed based on the x-ray crystal structures of both the active and inactive β2-adrenergic receptors (β2AR) in an effort to characterize the differences of binding modes between agonists and antagonists to the receptor, and to understand their activity in relation to their structures. The knowledge attained in this study is expected to provide valuable information on GPER ligands structure activity relationship to benefit future rational design of potent agonists and antagonists of the receptor for potential therapeutic applications.
Collapse
Affiliation(s)
- Christopher K Arnatt
- Department of Medicinal Chemistry, Virginia Commonwealth University, 800 E. Leigh Street, Richmond, VA 23289, USA, ; Tel: 804-828-7625
| | - Yan Zhang
- Department of Medicinal Chemistry, Virginia Commonwealth University, 800 E. Leigh Street, Richmond, VA 23289, USA, ; Tel: 804-828-7625
| |
Collapse
|
18
|
Cascella M, Bärfuss S, Stocker A. Cis-retinoids and the chemistry of vision. Arch Biochem Biophys 2013; 539:187-95. [PMID: 23791723 DOI: 10.1016/j.abb.2013.06.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2013] [Revised: 06/04/2013] [Accepted: 06/10/2013] [Indexed: 11/17/2022]
Abstract
We discuss here principal biochemical transformations of retinoid molecules in the visual cycle. We focus our analysis on the accumulating evidence of alternate pathways and functional redundancies in the cycle. The efficiency of the visual cycle depends, on one hand, on fast regeneration of the photo-bleached chromophores. On the other hand, it is crucial that the cyclic process should be highly selective to avoid accumulation of byproducts. The state-of-the-art knowledge indicates that single enzymatically active components of the cycle are not strictly selective and may require chaperones to enhance their rates. It appears that protein-protein interactions significantly improve the biological stability of the visual cycle. In particular, synthesis of thermodynamically less stable 11-cis-retinoid conformers is favored by physical interactions of the isomerases present in the retina with cellular retinaldehyde binding protein.
Collapse
Affiliation(s)
- Michele Cascella
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland.
| | | | | |
Collapse
|
19
|
Parrill AL. Computational Design and Experimental Characterization of GPCR Segment Models. Methods Enzymol 2013; 522:81-95. [DOI: 10.1016/b978-0-12-407865-9.00005-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|