1
|
Wong AR, Hung A, Yang AWH, Gill H, Lenon GB. Poria cocos compounds targeting neuropeptide Y1 receptor (Y1R) for weight management: A computational ligand- and structure-based study with molecular dynamics simulations identified beta-amyrin acetate as a putative Y1R inhibitor. PLoS One 2023; 18:e0277873. [PMID: 37390097 PMCID: PMC10313034 DOI: 10.1371/journal.pone.0277873] [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: 02/14/2022] [Accepted: 10/09/2022] [Indexed: 07/02/2023] Open
Abstract
Poria cocos (PC) is a medicinal herb frequently used in weight-loss clinical trials, however the mechanisms by which its compounds target orexigenic receptors including the neuropeptide Y1 receptor (Y1R) remain largely unknown. This study aimed to screen PC compounds for favourable pharmacokinetics profiles and examine their molecular mechanisms targeting Y1R. Forty-three PC compounds were systematically sought from pharmacological databases and docked with Y1R (PDB: 5ZBQ). By comparing the relative binding affinities, pharmacokinetics and toxicity profiles, we hypothesised that compounds designated PC1 3,4-Dihydroxybenzoic acid, PC8 Vanillic acid, PC40 1-(alpha-L-Ribofuranosyl)uracil, could be potential antagonists as they contact major residues Asn283 and Asp287, similar to various potent Y1R antagonists. In addition, PC21 Poricoic acid B, PC22 Poricoic acid G and PC43 16alpha,25-Dihydroxy-24-methylene-3,4-secolanosta-4(28),7,9(11)-triene-3,21-dioic acid, contacting Asn299, Asp104 and Asp200 proximal to the extracellular surface could also interfere with agonist binding by stabilising the extracellular loop (ECL) 2 of Y1R in a closed position. Owing to their selective interaction with Phe302, an important residue in binding of selective Y1R antagonists, PC12 beta-Amyrin acetate, PC26 3-Epidehydrotumulosic acid and PC27 Cerevisterol were proposed as putative antagonists. Following the consensus approach, PC12 beta-Amyrin acetate, PC26 3-Epidehydrotumulosic acid and PC27 Cerevisterol were identified as candidate compounds due to their high affinities (-12.2, -11.0 and -10.8 kcal, respectively), high drug-likeness and low toxicity profiles. Trajectory analyses and energy contributions of PC12-Y1R complex further confirmed their structural stability and favourable binding free energies, highlighting the feasibility and possible development of PC12 beta-Amyrin acetate as a future Y1R inhibitor.
Collapse
Affiliation(s)
- Ann Rann Wong
- School of Health and Biomedical Sciences, RMIT University, Bundoora, Victoria, Australia
| | - Andrew Hung
- School of Science, RMIT University, Melbourne, Victoria, Australia
| | - Angela Wei Hong Yang
- School of Health and Biomedical Sciences, RMIT University, Bundoora, Victoria, Australia
| | - Harsharn Gill
- School of Science, RMIT University, Melbourne, Victoria, Australia
| | - George Binh Lenon
- School of Health and Biomedical Sciences, RMIT University, Bundoora, Victoria, Australia
| |
Collapse
|
2
|
Palma J, Pierdominici-Sottile G. On the Uses of PCA to Characterise Molecular Dynamics Simulations of Biological Macromolecules: Basics and Tips for an Effective Use. Chemphyschem 2023; 24:e202200491. [PMID: 36285677 DOI: 10.1002/cphc.202200491] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2022] [Revised: 08/24/2022] [Indexed: 01/20/2023]
Abstract
Principal Component Analysis (PCA) is a procedure widely used to examine data collected from molecular dynamics simulations of biological macromolecules. It allows for greatly reducing the dimensionality of their configurational space, facilitating further qualitative and quantitative analysis. Its simplicity and relatively low computational cost explain its extended use. However, a judicious implementation of PCA requires the knowledge of its theoretical grounds as well as its weaknesses and capabilities. In this article, we review these issues and discuss several strategies developed over the last years to mitigate the main PCA flaws and enhance the reproducibility of its results.
Collapse
Affiliation(s)
- Juliana Palma
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes.,Consejo Nacional de Investigaciones Científicas y Técnicas
| | - Gustavo Pierdominici-Sottile
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes.,Consejo Nacional de Investigaciones Científicas y Técnicas
| |
Collapse
|
3
|
Turan Yücel N, Evren AE, Kandemir Ü, Can ÖD. Antidepressant-like effect of tofisopam in mice: A behavioural, molecular docking and MD simulation study. J Psychopharmacol 2022; 36:819-835. [PMID: 35638175 DOI: 10.1177/02698811221095528] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
BACKGROUND Depression is a disease that affects millions of people worldwide, and the discovery and development of effective and safe antidepressant drugs is one of the important topics of psychopharmacology. OBJECTIVES In this study, it was aimed to investigate the antidepressant-like activity potential of tofisopam, an anxiolytic drug with 2,3-benzodiazepine structure, and to elucidate the pharmacological mechanisms mediating this effect. METHODS The antidepressant-like activity of tofisopam was investigated using tail suspension and modified forced swimming tests. Possible interactions of tofisopam with µ- and δ-opioid receptor subtypes were clarified by pharmacological antagonism, molecular docking and molecular dynamics simulation studies. RESULTS Tofisopam (50 and 100 mg/kg) significantly shortened the immobility time of mice in both the tail suspension and the modified forced swimming tests. The drug, at the same doses, prolonged the duration of swimming and climbing behaviours measured in modified forced swimming tests. A dosage of 25 mg/kg was ineffective. Mechanistic studies showed that the pretreatment with p-chlorophenylalanine methyl ester (serotonin synthesis inhibitor; 4 consecutive days, 100 mg/kg), α-methyl-para-tyrosine methyl ester (catecholamine synthesis inhibitor; 100 mg/kg), naloxonazine (selective µ-opioid receptor blocker, 7 mg/kg) and naltrindole (a selective δ-opioid receptor blocker, 0.99 mg/kg) abolished the anti-immobility effect induced by the 50 mg/kg dose of tofisopam in the tail suspension tests. Our in silico studies supported the behavioural findings that the antidepressant-like effect of tofisopam is mediated by μ- and δ-opioid receptors. CONCLUSION This study is the first to show that tofisopam has antidepressant-like activity mediated by the serotonergic, catecholaminergic and opioidergic systems.
Collapse
Affiliation(s)
- Nazlı Turan Yücel
- Department of Pharmacology, Faculty of Pharmacy, Anadolu University, Eskişehir, Turkey
| | - Asaf Evrim Evren
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Anadolu University, Eskişehir, Turkey.,Pharmacy Services, Vocational School of Health Services, Bilecik Şeyh Edebali University, Bilecik, Turkey
| | - Ümmühan Kandemir
- Department of Pharmacology, Institute of Health Sciences, Anadolu University, Eskişehir, Turkey
| | - Özgür Devrim Can
- Department of Pharmacology, Faculty of Pharmacy, Anadolu University, Eskişehir, Turkey
| |
Collapse
|
4
|
Wang H, Cao D, Gillespie JC, Mendez RE, Selley DE, Liu-Chen LY, Zhang Y. Exploring the putative mechanism of allosteric modulations by mixed-action kappa/mu opioid receptor bitopic modulators. Future Med Chem 2021; 13:551-573. [PMID: 33590767 PMCID: PMC8027703 DOI: 10.4155/fmc-2020-0308] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 01/14/2021] [Indexed: 12/26/2022] Open
Abstract
The modulation and selectivity mechanisms of seven mixed-action kappa opioid receptor (KOR)/mu opioid receptor (MOR) bitopic modulators were explored. Molecular modeling results indicated that the 'message' moiety of seven bitopic modulators shared the same binding mode with the orthosteric site of the KOR and MOR, whereas the 'address' moiety bound with different subdomains of the allosteric site of the KOR and MOR. The 'address' moiety of seven bitopic modulators bound to different subdomains of the allosteric site of the KOR and MOR may exhibit distinguishable allosteric modulations to the binding affinity and/or efficacy of the 'message' moiety. Moreover, the 3-hydroxy group on the phenolic moiety of the seven bitopic modulators induced selectivity to the KOR over the MOR.
Collapse
Affiliation(s)
- Huiqun Wang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Danni Cao
- Center for Substance Abuse Research, Temple University Lewis Katz School of Medicine, Philadelphia, PA 19140, USA
| | - James C Gillespie
- Department of Pharmacology & Toxicology, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Rolando E Mendez
- Department of Pharmacology & Toxicology, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Dana E Selley
- Department of Pharmacology & Toxicology, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Lee-Yuan Liu-Chen
- Center for Substance Abuse Research, Temple University Lewis Katz School of Medicine, Philadelphia, PA 19140, USA
| | - Yan Zhang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USA
| |
Collapse
|
5
|
Park PSH. Supramolecular organization of rhodopsin in rod photoreceptor cell membranes. Pflugers Arch 2021; 473:1361-1376. [PMID: 33591421 DOI: 10.1007/s00424-021-02522-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Revised: 01/12/2021] [Accepted: 01/18/2021] [Indexed: 11/30/2022]
Abstract
Rhodopsin is the light receptor in rod photoreceptor cells that initiates scotopic vision. Studies on the light receptor span well over a century, yet questions about the organization of rhodopsin within the photoreceptor cell membrane still persist and a consensus view on the topic is still elusive. Rhodopsin has been intensely studied for quite some time, and there is a wealth of information to draw from to formulate an organizational picture of the receptor in native membranes. Early experimental evidence in apparent support for a monomeric arrangement of rhodopsin in rod photoreceptor cell membranes is contrasted and reconciled with more recent visual evidence in support of a supramolecular organization of rhodopsin. What is known so far about the determinants of forming a supramolecular structure and possible functional roles for such an organization are also discussed. Many details are still missing on the structural and functional properties of the supramolecular organization of rhodopsin in rod photoreceptor cell membranes. The emerging picture presented here can serve as a springboard towards a more in-depth understanding of the topic.
Collapse
Affiliation(s)
- Paul S-H Park
- Department of Ophthalmology and Visual Sciences, Case Western Reserve University, Cleveland, OH, 44106, USA.
| |
Collapse
|
6
|
Lei T, Hu Z, Ding R, Chen J, Li S, Zhang F, Pu X, Zhao N. Exploring the Activation Mechanism of a Metabotropic Glutamate Receptor Homodimer via Molecular Dynamics Simulation. ACS Chem Neurosci 2020; 11:133-145. [PMID: 31815422 DOI: 10.1021/acschemneuro.9b00425] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Metabotropic glutamate receptors of class C GPCRs exist as constitutive dimers, which play important roles in activating excitatory synapses of the central nervous system. However, the activation mechanism induced by agonists has not been clarified in experiments. To address the problem, we used microsecond all-atom molecular dynamics (MD) simulation couple with protein structure network (PSN) to explore the glutamate-induced activation for the mGluR1 homodimer. The results indicate that glutamate binding stabilizes not only the closure of Venus flytrap domains but also the polar interaction of LB2-LB2, in turn keeping the extracelluar domain in the active state. The activation of the extracelluar domain drives transmembrane domains (TMDs) of the two protomers closer and induces asymmetric activation for the TMD domains of the two protomers. One protomer with lower binding affinity to the agonist is activated, while the other protomer with higher binding energy is still in the inactive state. The PSN analysis identifies the allosteric regulation pathway from the ligand-binding pocket in the extracellular domain to the G-protein binding site in the intracellular TMD region and further reveals that the asymmetric activation is attributed to a combination of trans-pathway and cis-pathway regulations from two glumatates, rather than a single activation pathway. These observations could provide valuable molecular information for understanding of the structure and the implications in drug efficacy for the class C GPCR dimers.
Collapse
Affiliation(s)
- Ting Lei
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Zhenxin Hu
- College of Computer Science, Sichuan University, Chengdu 610064, China
| | - Ruolin Ding
- West China School of Stomatology, Sichuan University, Chengdu 610041, China
| | - Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Shiqi Li
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Fuhui Zhang
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Nanrong Zhao
- College of Chemistry, Sichuan University, Chengdu 610064, China
| |
Collapse
|
7
|
Barreto CAV, Baptista SJ, Preto AJ, Matos-Filipe P, Mourão J, Melo R, Moreira I. Prediction and targeting of GPCR oligomer interfaces. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 169:105-149. [PMID: 31952684 DOI: 10.1016/bs.pmbts.2019.11.007] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.
Collapse
Affiliation(s)
- Carlos A V Barreto
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal
| | - Salete J Baptista
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal; Centro de Ciências e Tecnologias Nucleares, Instituto Superior Técnico, Universidade de Lisboa, CTN, LRS, Portugal
| | - António José Preto
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal
| | - Pedro Matos-Filipe
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal
| | - Joana Mourão
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal; Institute for Interdisciplinary Research, University of Coimbra, Coimbra, Portugal
| | - Rita Melo
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal; Centro de Ciências e Tecnologias Nucleares, Instituto Superior Técnico, Universidade de Lisboa, CTN, LRS, Portugal
| | - Irina Moreira
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal; Science and Technology Faculty, University of Coimbra, Coimbra, Portugal.
| |
Collapse
|
8
|
Shen L, Yuan Y, Guo Y, Li M, Li C, Pu X. Probing the Druggablility on the Interface of the Protein-Protein Interaction and Its Allosteric Regulation Mechanism on the Drug Screening for the CXCR4 Homodimer. Front Pharmacol 2019; 10:1310. [PMID: 31787895 PMCID: PMC6855241 DOI: 10.3389/fphar.2019.01310] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 10/15/2019] [Indexed: 12/19/2022] Open
Abstract
Modulating protein–protein interactions (PPIs) with small drug-like molecules targeting it exhibits great promise in modern drug discovery. G protein-coupled receptors (GPCRs) are the largest family of targeted proteins and could form dimers in living biological cells through PPIs. However, compared to drug development of the orthosteric site, there has been lack of investigations on the druggability of the PPI interface for GPCRs and its functional implication on experiments. Thus, in order to address these issues, we constructed a novel computational strategy, which involved in molecular dynamics simulation, virtual screening and protein structure network (PSN), to study one representative GPCR homodimer (CXCR4). One druggable pocket was identified in the PPI interface and one small molecule targeting it was screened, which could strengthen PPI mainly through hydrophobic interaction between the benzene rings of the PPI molecule and TM4 of the receptor. The PSN results further reveals that the PPI molecule could increase the number of the allosteric regulation pathways between the druggable pocket of the dimer interface to the orthostatic site for the subunit A but only play minor role for the other subunit B, leading to the asymmetric change in the volume of the binding pockets for the two subunits (increase for the subunit A and minor change for the subunit B). Consequently, the screening performance of the subunit A to the antagonists is enhanced while the subunit B is unchanged nearly, implying that the PPI molecule may be beneficial to enhance the drug efficacies of the antagonists. In addition, one main regulation pathway with the highest frequency was identified for the subunit A, which consists of Trp1955.34–Tyr190ECL2–Val1965.35–Gln2005.39–Asp2626.58–Cys28N-term, revealing their importance in the allosteric regulation from the PPI molecule. The observations from the work could provide valuable information for the development of the PPI drug-like molecule for GPCRs.
Collapse
Affiliation(s)
- Liting Shen
- College of Chemistry, Sichuan University, Chengdu, China
| | - Yuan Yuan
- College of Management, Southwest University for Nationalities, Chengdu, China
| | - Yanzhi Guo
- College of Chemistry, Sichuan University, Chengdu, China
| | - Menglong Li
- College of Chemistry, Sichuan University, Chengdu, China
| | - Chuan Li
- College of Computer Science, Sichuan University, Chengdu, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu, China
| |
Collapse
|
9
|
Townsend-Nicholson A, Altwaijry N, Potterton A, Morao I, Heifetz A. Computational prediction of GPCR oligomerization. Curr Opin Struct Biol 2019; 55:178-184. [PMID: 31170578 DOI: 10.1016/j.sbi.2019.04.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2019] [Revised: 04/16/2019] [Accepted: 04/17/2019] [Indexed: 01/08/2023]
Abstract
There has been a recent and prolific expansion in the number of GPCR crystal structures being solved: in both active and inactive forms and in complex with ligand, with G protein and with each other. Despite this, there is relatively little experimental information about the precise configuration of GPCR oligomers during these different biologically relevant states. While it may be possible to identify the experimental conditions necessary to crystallize a GPCR preferentially in a specific structural conformation, computational approaches afford a potentially more tractable means of describing the probability of formation of receptor dimers and higher order oligomers. Ensemble-based computational methods based on structurally determined dimers, coupled with a computational workflow that uses quantum mechanical methods to analyze the chemical nature of the molecular interactions at a GPCR dimer interface, will generate the reproducible and accurate predictions needed to predict previously unidentified GPCR dimers and to inform future advances in our ability to understand and begin to precisely manipulate GPCR oligomers in biological systems. It may also provide information needed to achieve an increase in the number of experimentally determined oligomeric GPCR structures.
Collapse
Affiliation(s)
- Andrea Townsend-Nicholson
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom.
| | - Nojood Altwaijry
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom
| | - Andrew Potterton
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom; Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| |
Collapse
|
10
|
Rozo-Lugo C, Cuca-Suárez LE, Schmidt TJ, Coy-Barrera E. Tetrahydrobenzofuran-6(2 H)-one Neolignans from Ocotea heterochroma: Their Platelet Activating Factor (PAF) Antagonistic Activity and in Silico Insights into the PAF Receptor Binding Mode. JOURNAL OF NATURAL PRODUCTS 2018; 81:1968-1975. [PMID: 30188730 DOI: 10.1021/acs.jnatprod.8b00189] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Three new tetrahydrobenzofuran-6(2 H)-one-type neolignans, heterochromins A-C (1-3), along with a bicyclo[3.2.1]octane neolignan, cinerin C (4), were isolated from an ethanol extract from the leaves of Ocotea heterochroma, a native plant growing in the Colombo-Ecuadorian region of the Andes. The chemical structures of 1-3 were elucidated by spectroscopic methods. The platelet activating factor (PAF) antagonistic activity was tested in vitro for these compounds. Additionally, their binding mode to the PAF receptor was studied by molecular docking and molecular dynamics simulations in order to rationalize such activity. Heterochromin A (1) was found to be a potent PAF antagonist with a favorable molecular profile for interacting with the PAF receptor binding site.
Collapse
Affiliation(s)
- Catalina Rozo-Lugo
- Laboratorio de Investigación en Productos Naturales Vegetales, Departamento de Química , Universidad Nacional de Colombia, Ciudad Universitaria , 111321 , Bogotá D.C. , Colombia
- Laboratorio de Química Bioorgánica, Departamento de Química, Facultad de Ciencias Básicas y Aplicadas , Universidad Militar Nueva Granada , 250247 , Cajicá , Colombia
| | - Luis Enrique Cuca-Suárez
- Laboratorio de Investigación en Productos Naturales Vegetales, Departamento de Química , Universidad Nacional de Colombia, Ciudad Universitaria , 111321 , Bogotá D.C. , Colombia
| | - Thomas J Schmidt
- Institut für Pharmazeutische Biologie und Phytochemie (IPBP) , University of Münster, PharmaCampus , Corrensstraße 48 , D-48149 , Münster , Germany
| | - Ericsson Coy-Barrera
- Laboratorio de Química Bioorgánica, Departamento de Química, Facultad de Ciencias Básicas y Aplicadas , Universidad Militar Nueva Granada , 250247 , Cajicá , Colombia
| |
Collapse
|
11
|
Durdagi S, Erol I, Salmas RE, Aksoydan B, Kantarcioglu I. Oligomerization and cooperativity in GPCRs from the perspective of the angiotensin AT1 and dopamine D2 receptors. Neurosci Lett 2018; 700:30-37. [PMID: 29684528 DOI: 10.1016/j.neulet.2018.04.028] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2018] [Revised: 04/13/2018] [Accepted: 04/16/2018] [Indexed: 12/22/2022]
Abstract
G Protein-Coupled Receptors (GPCRs) can form homo- and heterodimers or constitute higher oligomeric clusters with other heptahelical GPCRs. In this article, multiscale molecular modeling approaches as well as experimental techniques which are used to study oligomerization of GPCRs are reviewed. In particular, the effect of dimerization/oligomerization to the ligand binding affinity of individual protomers and also on the efficacy of the oligomer are discussed by including diverse examples from the literature. In addition, possible allosteric effects that may emerge upon interaction of GPCRs with membrane components, like cholesterol, is also discussed. Investigation of these above-mentioned interactions may greatly contribute to the candidate molecule screening studies and development of novel therapeutics with fewer adverse effects.
Collapse
Affiliation(s)
- Serdar Durdagi
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University (BAU), Istanbul, Turkey; Neuroscience Program, Graduate School of Health Sciences, Bahcesehir University, Istanbul, Turkey.
| | - Ismail Erol
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University (BAU), Istanbul, Turkey; Department of Chemistry, Gebze Technical University, Kocaeli, Turkey
| | - Ramin Ekhteiari Salmas
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University (BAU), Istanbul, Turkey
| | - Busecan Aksoydan
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University (BAU), Istanbul, Turkey; Neuroscience Program, Graduate School of Health Sciences, Bahcesehir University, Istanbul, Turkey
| | - Isik Kantarcioglu
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University (BAU), Istanbul, Turkey; Bioengineering Program, Graduate School of Natural and Applied Sciences, Bahcesehir University, Istanbul, Turkey
| |
Collapse
|
12
|
Guidolin D, Marcoli M, Tortorella C, Maura G, Agnati LF. G protein-coupled receptor-receptor interactions give integrative dynamics to intercellular communication. Rev Neurosci 2018; 29:703-726. [DOI: 10.1515/revneuro-2017-0087] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Accepted: 01/01/2018] [Indexed: 01/14/2023]
Abstract
Abstract
The proposal of receptor-receptor interactions (RRIs) in the early 1980s broadened the view on the role of G protein-coupled receptors (GPCR) in the dynamics of the intercellular communication. RRIs, indeed, allow GPCR to operate not only as monomers but also as receptor complexes, in which the integration of the incoming signals depends on the number, spatial arrangement, and order of activation of the protomers forming the complex. The main biochemical mechanisms controlling the functional interplay of GPCR in the receptor complexes are direct allosteric interactions between protomer domains. The formation of these macromolecular assemblies has several physiologic implications in terms of the modulation of the signaling pathways and interaction with other membrane proteins. It also impacts on the emerging field of connectomics, as it contributes to set and tune the synaptic strength. Furthermore, recent evidence suggests that the transfer of GPCR and GPCR complexes between cells via the exosome pathway could enable the target cells to recognize/decode transmitters and/or modulators for which they did not express the pertinent receptors. Thus, this process may also open the possibility of a new type of redeployment of neural circuits. The fundamental aspects of GPCR complex formation and function are the focus of the present review article.
Collapse
Affiliation(s)
- Diego Guidolin
- Department of Neuroscience , University of Padova, via Gabelli 65 , I-35121 Padova , Italy
| | - Manuela Marcoli
- Department of Pharmacy and Center of Excellence for Biomedical Research , University of Genova , I-16126 Genova , Italy
| | - Cinzia Tortorella
- Department of Neuroscience , University of Padova, via Gabelli 65 , I-35121 Padova , Italy
| | - Guido Maura
- Department of Pharmacy and Center of Excellence for Biomedical Research , University of Genova , I-16126 Genova , Italy
| | - Luigi F. Agnati
- Department of Biomedical Sciences , University of Modena and Reggio Emilia , I-41121 Modena , Italy
- Department of Neuroscience , Karolinska Institutet , S-17177 Stockholm , Sweden
| |
Collapse
|
13
|
Liang T, Yuan Y, Wang R, Guo Y, Li M, Pu X, Li C. Structural Features and Ligand Selectivity for 10 Intermediates in the Activation Process of β 2-Adrenergic Receptor. ACS OMEGA 2017; 2:8557-8567. [PMID: 30023586 PMCID: PMC6045391 DOI: 10.1021/acsomega.7b01031] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/19/2017] [Accepted: 09/14/2017] [Indexed: 06/08/2023]
Abstract
It has already been suggested by researchers that there should be multiple intermediate states in the activation process for G-protein-coupled receptors (GPCRs). However, the intermediate states are very short-lived and hardly captured by the experiments, leading to very limited understanding of their structural features and drug efficacies. In this work, a novel joint strategy of targeted molecular dynamics simulation, conventional molecular dynamics simulation, and virtual screening is developed to address the problems. The results from 10 intermediate conformations obtained from the work reveal that the ligand pocket is very unstable and fluctuates between the inactive state and the active one in the case of ligand-free, in particular for ECL2 as a gate-keeper of the ligand-binding. The ligand-binding site could be stable in the active state with a small volume and a completely closed ECL2, only when the G-protein-binding region is fully activated. In addition, the activations of the ligand-binding pocket and G-protein-binding site are relatively independent and exhibit a loose allosteric coupling, which contributes to the existence of multiple intermediate conformations. Interestingly, the screening performance of the agonists does not increase on increasing the overall activity of the intermediate state, but is dependent on the activated extent of the ligand pocket. The receptor is prone to bind the agonist when closing ECL2 and reducing the ligand-binding pocket volume, whereas it is more favorable for binding the antagonist when opening ECL2 and increasing the pocket volume. These observations added to previous studies could help us better understand the activation mechanism of GPCRs and provide valuable information for drug design.
Collapse
Affiliation(s)
- Tao Liang
- College
of Chemistry and College of Computer Science, Sichuan University, No. 29 Jiuyanqiao Wangjiang Road, Chengdu 610064, People’s Republic
of China
| | - Yuan Yuan
- College
of Management, Southwest University for
Nationalities, No. 16 South Section 4, Yihuan Road, Chengdu 610041, People’s Republic
of China
| | - Ran Wang
- College
of Chemistry and College of Computer Science, Sichuan University, No. 29 Jiuyanqiao Wangjiang Road, Chengdu 610064, People’s Republic
of China
| | - Yanzhi Guo
- College
of Chemistry and College of Computer Science, Sichuan University, No. 29 Jiuyanqiao Wangjiang Road, Chengdu 610064, People’s Republic
of China
| | - Menglong Li
- College
of Chemistry and College of Computer Science, Sichuan University, No. 29 Jiuyanqiao Wangjiang Road, Chengdu 610064, People’s Republic
of China
| | - Xuemei Pu
- College
of Chemistry and College of Computer Science, Sichuan University, No. 29 Jiuyanqiao Wangjiang Road, Chengdu 610064, People’s Republic
of China
| | - Chuan Li
- College
of Chemistry and College of Computer Science, Sichuan University, No. 29 Jiuyanqiao Wangjiang Road, Chengdu 610064, People’s Republic
of China
| |
Collapse
|
14
|
Ramírez-Anguita JM, Rodríguez-Espigares I, Guixà-González R, Bruno A, Torrens-Fontanals M, Varela-Rial A, Selent J. Membrane cholesterol effect on the 5-HT 2A receptor: Insights into the lipid-induced modulation of an antipsychotic drug target. Biotechnol Appl Biochem 2017; 65:29-37. [PMID: 28877377 DOI: 10.1002/bab.1608] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The serotonin 5-hydroxytryptamine 2A (5-HT2A ) receptor is a G-protein-coupled receptor (GPCR) relevant for the treatment of CNS disorders. In this regard, neuronal membrane composition in the brain plays a crucial role in the modulation of the receptor functioning. Since cholesterol is an essential component of neuronal membranes, we have studied its effect on the 5-HT2A receptor dynamics through all-atom MD simulations. We find that the presence of cholesterol in the membrane increases receptor conformational variability in most receptor segments. Importantly, detailed structural analysis indicates that conformational variability goes along with the destabilization of hydrogen bonding networks not only within the receptor but also between receptor and lipids. In addition to increased conformational variability, we also find receptor segments with reduced variability. Our analysis suggests that this increased stabilization is the result of stabilizing effects of tightly bound cholesterol molecules to the receptor surface. Our finding contributes to a better understanding of membrane-induced alterations of receptor dynamics and points to cholesterol-induced stabilizing and destabilizing effects on the conformational variability of GPCRs.
Collapse
Affiliation(s)
- Juan Manuel Ramírez-Anguita
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF)-Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| | - Ismael Rodríguez-Espigares
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF)-Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| | - Ramon Guixà-González
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF)-Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain.,Institut für Medizinische Physik und Biophysik, Charité, Berlin, Germany
| | - Agostino Bruno
- Department of Pharmacy, Università degli studi di Parma, Italy, Parma
| | - Mariona Torrens-Fontanals
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF)-Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| | - Alejandro Varela-Rial
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF)-Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| | - Jana Selent
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF)-Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| |
Collapse
|
15
|
Wang YT, Chan YH. Understanding the molecular basis of agonist/antagonist mechanism of human mu opioid receptor through gaussian accelerated molecular dynamics method. Sci Rep 2017; 7:7828. [PMID: 28798303 PMCID: PMC5552784 DOI: 10.1038/s41598-017-08224-2] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2017] [Accepted: 07/10/2017] [Indexed: 01/01/2023] Open
Abstract
The most powerful analgesic and addictive properties of opiate alkaloids are mediated by the μ opioid receptor (MOR). The MOR has been extensively investigated as a drug target in the twentieth century, with numerous compounds of varying efficacy being identified. We employed molecular dynamics and Gaussian accelerated molecular dynamics techniques to identify the binding mechanisms of MORs to BU72 (agonist) and β-funaltrexamine (antagonist). Our approach theoretically suggests that the 34 residues (Lys209–Phe221 and Ile301–Cys321) of the MORs were the key regions enabling the two compounds to bind to the active site of the MORs. When the MORs were in the holo form, the key region was in the open conformation. When the MORs were in the apo form, the key region was in the closed conformation. The key region might be responsible for the selectivity of new MOR agonists and antagonists.
Collapse
Affiliation(s)
- Yeng-Tseng Wang
- Department of Biochemistry, College of Medicine, Kaohsiung Medical University, Kaohsiung, Taiwan. .,Center for Biomarkers and Biotech Drugs, Kaohsiung Medical University, Kaohsiung, Taiwan. .,Graduate Institute of Medicine, Kaohsiung Medical University, Kaohsiung, Taiwan. .,Department of Medical Research, Kaohsiung Medical University Hospital, Kaohsiung, Taiwan.
| | - Yang-Hsiang Chan
- Department of Chemistry, National Sun Yat-sen University, 70 Lien Hai Road, Kaohsiung, Taiwan
| |
Collapse
|
16
|
Altwaijry NA, Baron M, Wright DW, Coveney PV, Townsend-Nicholson A. An Ensemble-Based Protocol for the Computational Prediction of Helix-Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics. J Chem Theory Comput 2017; 13:2254-2270. [PMID: 28383913 PMCID: PMC5557214 DOI: 10.1021/acs.jctc.6b01246] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
![]()
The
accurate identification of the specific points of interaction
between G protein-coupled receptor (GPCR) oligomers is essential for
the design of receptor ligands targeting oligomeric receptor targets.
A coarse-grained molecular dynamics computer simulation approach would
provide a compelling means of identifying these specific protein–protein
interactions and could be applied both for known oligomers of interest
and as a high-throughput screen to identify novel oligomeric targets.
However, to be effective, this in silico modeling must provide accurate,
precise, and reproducible information. This has been achieved recently
in numerous biological systems using an ensemble-based all-atom molecular
dynamics approach. In this study, we describe an equivalent methodology
for ensemble-based coarse-grained simulations. We report the performance
of this method when applied to four different GPCRs known to oligomerize
using error analysis to determine the ensemble size and individual
replica simulation time required. Our measurements of distance between
residues shown to be involved in oligomerization of the fifth transmembrane
domain from the adenosine A2A receptor are in very good
agreement with the existing biophysical data and provide information
about the nature of the contact interface that cannot be determined
experimentally. Calculations of distance between rhodopsin, CXCR4,
and β1AR transmembrane domains reported to form contact
points in homodimers correlate well with the corresponding measurements
obtained from experimental structural data, providing an ability to predict
contact interfaces computationally. Interestingly, error analysis
enables identification of noninteracting regions. Our results confirm
that GPCR interactions can be reliably predicted using this novel
methodology.
Collapse
Affiliation(s)
- Nojood A Altwaijry
- Institute of Structural and Molecular Biology, Research Department of Structural and Molecular Biology, Division of Biosciences, University College London , London, WC1E 6BT, United Kingdom.,King Saud University , Riyadh, Kingdom of Saudi Arabia
| | - Michael Baron
- Institute of Structural and Molecular Biology, Research Department of Structural and Molecular Biology, Division of Biosciences, University College London , London, WC1E 6BT, United Kingdom
| | - David W Wright
- Centre for Computational Science, Department of Chemistry, University College London , London WC1H 0AJ, United Kingdom
| | - Peter V Coveney
- Centre for Computational Science, Department of Chemistry, University College London , London WC1H 0AJ, United Kingdom
| | - Andrea Townsend-Nicholson
- Institute of Structural and Molecular Biology, Research Department of Structural and Molecular Biology, Division of Biosciences, University College London , London, WC1E 6BT, United Kingdom
| |
Collapse
|
17
|
Methods used to study the oligomeric structure of G-protein-coupled receptors. Biosci Rep 2017; 37:BSR20160547. [PMID: 28062602 PMCID: PMC5398257 DOI: 10.1042/bsr20160547] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Revised: 01/05/2017] [Accepted: 01/06/2017] [Indexed: 02/02/2023] Open
Abstract
G-protein-coupled receptors (GPCRs), which constitute the largest family of cell surface receptors, were originally thought to function as monomers, but are now recognized as being able to act in a wide range of oligomeric states and indeed, it is known that the oligomerization state of a GPCR can modulate its pharmacology and function. A number of experimental techniques have been devised to study GPCR oligomerization including those based upon traditional biochemistry such as blue-native PAGE (BN-PAGE), co-immunoprecipitation (Co-IP) and protein-fragment complementation assays (PCAs), those based upon resonance energy transfer, FRET, time-resolved FRET (TR-FRET), FRET spectrometry and bioluminescence resonance energy transfer (BRET). Those based upon microscopy such as FRAP, total internal reflection fluorescence microscopy (TIRFM), spatial intensity distribution analysis (SpIDA) and various single molecule imaging techniques. Finally with the solution of a growing number of crystal structures, X-ray crystallography must be acknowledged as an important source of discovery in this field. A different, but in many ways complementary approach to the use of more traditional experimental techniques, are those involving computational methods that possess obvious merit in the study of the dynamics of oligomer formation and function. Here, we summarize the latest developments that have been made in the methods used to study GPCR oligomerization and give an overview of their application.
Collapse
|
18
|
Mukherjee S, Karolak A, Debant M, Buscaglia P, Renaudineau Y, Mignen O, Guida WC, Brooks WH. Molecular Dynamics Simulations of Membrane-Bound STIM1 to Investigate Conformational Changes during STIM1 Activation upon Calcium Release. J Chem Inf Model 2017; 57:335-344. [PMID: 28151650 DOI: 10.1021/acs.jcim.6b00475] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Calcium is involved in important intracellular processes, such as intracellular signaling from cell membrane receptors to the nucleus. Typically, calcium levels are kept at less than 100 nM in the nucleus and cytosol, but some calcium is stored in the endoplasmic reticulum (ER) lumen for rapid release to activate intracellular calcium-dependent functions. Stromal interacting molecule 1 (STIM1) plays a critical role in early sensing of changes in the ER's calcium level, especially when there is a sudden release of stored calcium from the ER. Inactive STIM1, which has a bound calcium ion, is activated upon ion release. Following activation of STIM1, there is STIM1-assisted initiation of extracellular calcium entry through channels in the cell membrane. This extracellular calcium entering the cell then amplifies intracellular calcium-dependent actions. At the end of the process, ER levels of stored calcium are reestablished. The main focus of this work was to study the conformational changes accompanying homo- or heterodimerization of STIM1. For this purpose, the ER luminal portion of STIM1 (residues 58-236), which includes the sterile alpha motif (SAM) domain plus the calcium-binding EF-hand domains 1 and 2 attached to the STIM1 transmembrane region (TM), was modeled and embedded in a virtual membrane. Next, molecular dynamics simulations were performed to study the conformational changes that take place during STIM1 activation and subsequent protein-protein interactions. Indeed, the simulations revealed exposure of residues in the EF-hand domains, which may be important for dimerization steps. Altogether, understanding conformational changes in STIM1 can help in drug discovery when targeting this key protein in intracellular calcium functions.
Collapse
Affiliation(s)
- Sreya Mukherjee
- Department of Chemistry, University of South Florida , Tampa, Florida 33620, United States
| | - Aleksandra Karolak
- Department of Chemistry, University of South Florida , Tampa, Florida 33620, United States
| | - Marjolaine Debant
- INSERM ESPRI, ERI29/EA2216 Laboratory of Immunotherapy and B Cell Pathologies, Laboratory of Immunology and Immunotherapy, CHRU Morvan, European University of Brittany , F29609 Brest, France.,Network "Ion channels and cancer-Cancéropole Grand Ouest (IC-CGO)" , F29609 Brest, France.,INSERM U1078, Brest University Medical School , F29609 Brest, France
| | - Paul Buscaglia
- Network "Ion channels and cancer-Cancéropole Grand Ouest (IC-CGO)" , F29609 Brest, France.,INSERM U1078, Brest University Medical School , F29609 Brest, France
| | - Yves Renaudineau
- INSERM ESPRI, ERI29/EA2216 Laboratory of Immunotherapy and B Cell Pathologies, Laboratory of Immunology and Immunotherapy, CHRU Morvan, European University of Brittany , F29609 Brest, France.,Network "Ion channels and cancer-Cancéropole Grand Ouest (IC-CGO)" , F29609 Brest, France
| | - Olivier Mignen
- Network "Ion channels and cancer-Cancéropole Grand Ouest (IC-CGO)" , F29609 Brest, France.,INSERM U1078, Brest University Medical School , F29609 Brest, France
| | - Wayne C Guida
- Department of Chemistry, University of South Florida , Tampa, Florida 33620, United States
| | - Wesley H Brooks
- Department of Chemistry, University of South Florida , Tampa, Florida 33620, United States
| |
Collapse
|
19
|
Maximova T, Moffatt R, Ma B, Nussinov R, Shehu A. Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics. PLoS Comput Biol 2016; 12:e1004619. [PMID: 27124275 PMCID: PMC4849799 DOI: 10.1371/journal.pcbi.1004619] [Citation(s) in RCA: 132] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.
Collapse
Affiliation(s)
- Tatiana Maximova
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Ryan Moffatt
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
- Department of Biongineering, George Mason University, Fairfax, Virginia, United States of America
- School of Systems Biology, George Mason University, Manassas, Virginia, United States of America
| |
Collapse
|
20
|
Kaczor AA, Guixà-González R, Carrió P, Poso A, Dove S, Pastor M, Selent J. Multi-Component Protein - Protein Docking Based Protocol with External Scoring for Modeling Dimers of G Protein-Coupled Receptors. Mol Inform 2015; 34:246-55. [PMID: 27490170 DOI: 10.1002/minf.201400088] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2014] [Accepted: 02/19/2015] [Indexed: 01/08/2023]
Abstract
In order to apply structure-based drug design techniques to GPCR complexes, it is essential to model their 3D structure. For this purpose, a multi-component protocol was derived based on protein-protein docking which generates populations of dimers compatible with membrane integration, considering all reasonable interfaces. At the next stage, we applied a scoring procedure based on up to eleven different parameters including shape or electrostatics complementarity. Two methods of consensus scoring were performed: (i) average scores of 100 best scored dimers with respect to each interface, and (ii) frequencies of interfaces among 100 best scored dimers. In general, our multi-component protocol gives correct indications for dimer interfaces that have been observed in X-ray crystal structures of GPCR dimers (opsin dimer, chemokine CXCR4 and CCR5 dimers, κ opioid receptor dimer, β1 adrenergic receptor dimer and smoothened receptor dimer) but also suggests alternative dimerization interfaces. Interestingly, at times these alternative interfaces are scored higher than the experimentally observed ones suggesting them to be also relevant in the life cycle of studied GPCR dimers. Further results indicate that GPCR dimer and higher-order oligomer formation may involve transmembrane helices (TMs) TM1-TM2-TM7, TM3-TM4-TM5 or TM4-TM5-TM6 but not TM1-TM2-TM3 or TM2-TM3-TM4 which is in general agreement with available experimental and computational data.
Collapse
Affiliation(s)
- Agnieszka A Kaczor
- School of Pharmacy, University of Eastern Finland, Yliopistonranta 1, P.O. Box 1627, FI-70211 Kuopio, Finland phone: +48815357365; fax: +48815357366. .,Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modeling Lab, Faculty of Pharmacy with Division for Medical Analytics, 4A Chodźki St., PL-20059 Lublin, Poland. .,Research Unit on Biomedical Informatics (GRIB), PRBB, Dr Aiguader 88, E-08003 Barcelona, Spain. .,Department of Pharmaceutical Chemistry II, Institute of Pharmacy, University of Regensburg, DE-93040 Regensburg, Germany.
| | - Ramon Guixà-González
- Research Unit on Biomedical Informatics (GRIB), PRBB, Dr Aiguader 88, E-08003 Barcelona, Spain
| | - Pau Carrió
- Research Unit on Biomedical Informatics (GRIB), PRBB, Dr Aiguader 88, E-08003 Barcelona, Spain
| | - Antti Poso
- School of Pharmacy, University of Eastern Finland, Yliopistonranta 1, P.O. Box 1627, FI-70211 Kuopio, Finland phone: +48815357365; fax: +48815357366
| | - Stefan Dove
- Department of Pharmaceutical Chemistry II, Institute of Pharmacy, University of Regensburg, DE-93040 Regensburg, Germany
| | - Manuel Pastor
- Research Unit on Biomedical Informatics (GRIB), PRBB, Dr Aiguader 88, E-08003 Barcelona, Spain
| | - Jana Selent
- Research Unit on Biomedical Informatics (GRIB), PRBB, Dr Aiguader 88, E-08003 Barcelona, Spain.
| |
Collapse
|
21
|
Pierdominici-Sottile G, Palma J. New insights into the meaning and usefulness of principal component analysis of concatenated trajectories. J Comput Chem 2014; 36:424-32. [PMID: 25516482 DOI: 10.1002/jcc.23811] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2014] [Revised: 11/12/2014] [Accepted: 11/23/2014] [Indexed: 11/12/2022]
Abstract
A comparison between different conformations of a given protein, relating both structure and dynamics, can be performed in terms of combined principal component analysis (combined-PCA). To that end, a trajectory is obtained by concatenating molecular dynamics trajectories of the individual conformations under comparison. Then, the principal components are calculated by diagonalizing the correlation matrix of the concatenated trajectory. Since the introduction of this approach in 1995 it has had a large number of applications. However, the interpretation of the eigenvectors and eigenvalues so obtained is based on intuitive foundations, because analytical expressions relating the concatenated correlation matrix with those of the individual trajectories under consideration have not been provided yet. In this article, we present such expressions for the cases of two, three, and an arbitrary number of concatenated trajectories. The formulas are simple and show what is to be expected and what is not to be expected from a combined-PCA. Their correctness and usefulness is demonstrated by discussing some representative examples. The results can be summarized in a simple sentence: the correlation matrix of a concatenated trajectory is given by the average of the individual correlation matrices plus the correlation matrix of the individual averages. From this it follows that the combined-PCA of trajectories belonging to different free energy basins provides information that could also be obtained by alternative and more straightforward means.
Collapse
Affiliation(s)
- Gustavo Pierdominici-Sottile
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, Sáenz Peña 352, Bernal, B1876BXD, Argentina
| | | |
Collapse
|
22
|
Understanding the effects on constitutive activation and drug binding of a D130N mutation in the β2 adrenergic receptor via molecular dynamics simulation. J Mol Model 2014; 20:2491. [PMID: 25342155 DOI: 10.1007/s00894-014-2491-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2014] [Accepted: 10/06/2014] [Indexed: 01/21/2023]
Abstract
G-protein-coupled receptors (GPCRs) are currently one of the largest families of drug targets. The constitutive activation induced by mutation of key GPCR residues is associated closely with various diseases. However, the structural basis underlying such activation and its role in drug binding has remained unclear. Herein, we used all-atom molecular dynamics simulations and free energy calculations to study the effects of a D130N mutation on the structure of β2 adrenergic receptor (β2AR) and its binding of the agonist salbutamol. The results indicate that the mutation caused significant changes in some key helices. In particular, the mutation leads to the departure of transmembrane 3 (TM3) from transmembrane 6 (TM6) and marked changes in the NPxxY region as well as the complete disruption of a key ionic lock, all of which contribute to the observed constitutive activation. In addition, the D130N mutation weakens some important H-bonds, leading to structural changes in these regions. Binding free energy calculations indicate that van der Waals and electrostatic interactions are the main driving forces in binding salbutamol; however, binding strength in the mutant β2AR is significantly enhanced mainly through modifying electrostatic interactions. Further analysis revealed that the increase in binding energy upon mutation stems mainly from the H-bonds formed between the hydroxyl group of salbutamol and the serine residues of TM5. This observation suggests that modifications of the H-bond groups of this drug could significantly influence drug efficacy in the treatment of diseases associated with this mutation.
Collapse
|
23
|
Ng HW, Laughton CA, Doughty SW. Molecular dynamics simulations of the adenosine A2a receptor in POPC and POPE lipid bilayers: effects of membrane on protein behavior. J Chem Inf Model 2014; 54:573-81. [PMID: 24460123 DOI: 10.1021/ci400463z] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Analysis of 300 ns (ns) molecular dynamics (MD) simulations of an adenosine A2a receptor (A2a AR) model, conducted in triplicate, in 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) bilayers reveals significantly different protein dynamical behavior. Principal component analysis (PCA) shows that the dissimilarities stem from interhelical rather than intrahelical motions. The difference in the hydrophobic thicknesses of these simulated lipid bilayers is potentially a significant reason for the observed difference in results. The distinct lipid headgroups might also lead to different molecular interactions and hence different protein loop motions. Overall, the A2a AR shows higher mobility and flexibility in POPC as compared to POPE.
Collapse
Affiliation(s)
- Hui Wen Ng
- School of Pharmacy, University of Nottingham Malaysia Campus , Jalan Broga, 43500 Semenyih, Selangor, Malaysia
| | | | | |
Collapse
|
24
|
Heifetz A, Barker O, Morris GB, Law R, Slack M, Biggin PC. Toward an understanding of agonist binding to human Orexin-1 and Orexin-2 receptors with G-protein-coupled receptor modeling and site-directed mutagenesis. Biochemistry 2013; 52:8246-60. [PMID: 24144388 PMCID: PMC3880013 DOI: 10.1021/bi401119m] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2013] [Revised: 10/21/2013] [Indexed: 02/06/2023]
Abstract
The class A G-protein-coupled receptors (GPCRs) Orexin-1 (OX1) and Orexin-2 (OX2) are located predominantly in the brain and are linked to a range of different physiological functions, including the control of feeding, energy metabolism, modulation of neuro-endocrine function, and regulation of the sleep-wake cycle. The natural agonists for OX1 and OX2 are two neuropeptides, Orexin-A and Orexin-B, which have activity at both receptors. Site-directed mutagenesis (SDM) has been reported on both the receptors and the peptides and has provided important insight into key features responsible for agonist activity. However, the structural interpretation of how these data are linked together is still lacking. In this work, we produced and used SDM data, homology modeling followed by MD simulation, and ensemble-flexible docking to generate binding poses of the Orexin peptides in the OX receptors to rationalize the SDM data. We also developed a protein pairwise similarity comparing method (ProS) and a GPCR-likeness assessment score (GLAS) to explore the structural data generated within a molecular dynamics simulation and to help distinguish between different GPCR substates. The results demonstrate how these newly developed methods of structural assessment for GPCRs can be used to provide a working model of neuropeptide-Orexin receptor interaction.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec
(U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, U.K.
| | - Oliver Barker
- Evotec
(U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, U.K.
| | - G. Benjamin Morris
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| | - Richard
J. Law
- Evotec
(U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, U.K.
| | - Mark Slack
- Evotec
AG, Manfred Eigen Campus,
Essener Bogen 7, 22419 Hamburg, Germany
| | - Philip C. Biggin
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
25
|
Moreira IS. Structural features of the G-protein/GPCR interactions. Biochim Biophys Acta Gen Subj 2013; 1840:16-33. [PMID: 24016604 DOI: 10.1016/j.bbagen.2013.08.027] [Citation(s) in RCA: 75] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2013] [Revised: 08/27/2013] [Accepted: 08/28/2013] [Indexed: 01/07/2023]
Abstract
BACKGROUND The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies. SCOPE OF REVIEW The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined. MAJOR CONCLUSIONS Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered. GENERAL SIGNIFICANCE In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.
Collapse
Affiliation(s)
- Irina S Moreira
- REQUIMTE/Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto, Portugal.
| |
Collapse
|
26
|
Bruno A, Amori L, Costantino G. Computational Insights into the Mechanism of Inhibition of OASS-A by a Small Molecule Inhibitor. Mol Inform 2013; 32:447-57. [PMID: 27481665 DOI: 10.1002/minf.201200174] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2012] [Accepted: 02/09/2012] [Indexed: 11/09/2022]
Abstract
O-Acetylserine sulfhydrylase (isoform A, OASS-A) is a PLP-dependent enzyme involved in the last step of cysteine biosynthesis in many pathogens. Many microorganisms use cysteine as the main building block for sulfur-containing antioxidants, and cysteine depletion in several pathogens resulted in a reduced antibiotic resistance, thus leading to the identification of OASS as novel suitable molecular targets to overcome antimicrobial resistances. The precise molecular mechanism of OASS-A inhibition by small peptides or by small molecule inhibitors is still unclear. To shed more lights on the structural basis underlying the inhibition mechanism for OASS, we engaged ourselves in studying the dynamic properties of this enzyme. In this paper, we describe a computational study involving unbiased MD simulations of OASS-A from Haemophilus influenzae (HiOASS) in its inhibitor free, PLP-bound form, and in complex with a pentapeptide inhibitor and with UPAR40, a small molecule which we have recently reported as a potent OASS-A inhibitors. We proposed that UPAR40 inhibits HiOASS-A through the stabilization of a closed conformation. Moreover, preliminary docking studies and sequence analysis allow us to speculate about the non-specificity of UPAR40 toward a particular OASS enzyme species or isoforms.
Collapse
Affiliation(s)
- Agostino Bruno
- Department of Pharmacy, Università degli Studi di Parma, Parco Area delle Scienze, Viale G. P. Usberti 27/A, 43124, Parma, Italy
| | - Laura Amori
- Department of Pharmacy, Università degli Studi di Parma, Parco Area delle Scienze, Viale G. P. Usberti 27/A, 43124, Parma, Italy
| | - Gabriele Costantino
- Department of Pharmacy, Università degli Studi di Parma, Parco Area delle Scienze, Viale G. P. Usberti 27/A, 43124, Parma, Italy.
| |
Collapse
|
27
|
Ng HW, Laughton CA, Doughty SW. Molecular dynamics simulations of the adenosine A2a receptor: structural stability, sampling, and convergence. J Chem Inf Model 2013; 53:1168-78. [PMID: 23514445 DOI: 10.1021/ci300610w] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.
Collapse
Affiliation(s)
- Hui Wen Ng
- School of Pharmacy, University of Nottingham Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor, Malaysia
| | | | | |
Collapse
|
28
|
Jastrzebska B, Orban T, Golczak M, Engel A, Palczewski K. Asymmetry of the rhodopsin dimer in complex with transducin. FASEB J 2013; 27:1572-84. [PMID: 23303210 DOI: 10.1096/fj.12-225383] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A large body of evidence for G-protein-coupled receptor (GPCR) oligomerization has accumulated over the past 2 decades. The smallest of these oligomers in vivo most likely is a dimer that buries 1000-Å(2) intramolecular surfaces and on stimulation forms a complex with heterotrimeric G protein in 2:1 stoichiometry. However, it is unclear whether each of the monomers adopts the same or a different conformation and function after activation of this dimer. With bovine rhodopsin (Rho) and its cognate bovine G-protein transducin (Gt) as a model system, we used the retinoid chromophores 11-cis-retinal and 9-cis-retinal to monitor each monomer of the dimeric GPCR within a stable complex with nucleotide-free Gt. We found that only 50% of Rho* in the Rho*-Gt complex is trapped in a Meta II conformation, while 50% evolves toward an opsin conformation and can be regenerated with 9-cis-retinal. We also found that all-trans-retinal can regenerate chromophore-depleted Rho*e complexed with Gt and FAK*TSA peptide containing Lys(296) with the attached all-trans retinoid (m/z of 934.5[MH](+)) was identified by mass spectrometry. Thus, our study shows that each of the monomers contributes unequally to the pentameric (2:1:1:1) complex of Rho dimer and Gt heterotrimer, validating the oligomeric structure of the complex and the asymmetry of the GPCR dimer, and revealing its structural/functional signature. This study provides a clear functional distinction between monomers of family A GPCRs in their oligomeric form.
Collapse
Affiliation(s)
- Beata Jastrzebska
- Department of Pharmacology, School of Medicine, Case Western Reserve University, Cleveland, Ohio 44106-4965, USA
| | | | | | | | | |
Collapse
|
29
|
Giraldo J. Modeling cooperativity effects in dimeric G protein-coupled receptors. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2013; 115:349-73. [PMID: 23415098 DOI: 10.1016/b978-0-12-394587-7.00008-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
G protein-coupled receptors organize into oligomeric arrangements to exert their function. In this chapter, three models of dimeric receptors, the two-state dimer receptor model, the metabotropic glutamate receptor model, and the asymmetric/symmetric three-state dimer receptor model are revisited focusing on the cooperative effects between their binding sites and the subunits they are composed of. The mathematical analysis reveals the complexity of the intra-receptor interactions providing insights on the mechanistic aspects of receptor function.
Collapse
Affiliation(s)
- Jesús Giraldo
- Laboratory of Systems Pharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Universitat Autònoma de Barcelona, Bellaterra, Spain
| |
Collapse
|
30
|
Cordomí A, Perez JJ. Structural Rearrangements of Rhodopsin Subunits in a Dimer Complex: a Molecular Dynamics Simulation Study. J Biomol Struct Dyn 2012; 27:127-47. [DOI: 10.1080/07391102.2009.10507303] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Arnau Cordomí
- a Dept d'Enginyeria Química , Technical University of Catalonia (UPC), ETS d'Enginyeria Industrial , Av. Diagonal 647, 08028 , Barcelona , Spain
| | - Juan J. Perez
- a Dept d'Enginyeria Química , Technical University of Catalonia (UPC), ETS d'Enginyeria Industrial , Av. Diagonal 647, 08028 , Barcelona , Spain
| |
Collapse
|
31
|
Light induced transmembrane proton gradient in artificial lipid vesicles reconstituted with photosynthetic reaction centers. J Bioenerg Biomembr 2012; 44:373-84. [DOI: 10.1007/s10863-012-9435-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2011] [Accepted: 03/19/2012] [Indexed: 10/28/2022]
|
32
|
Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties. PLoS Comput Biol 2012; 8:e1002473. [PMID: 22532793 PMCID: PMC3330085 DOI: 10.1371/journal.pcbi.1002473] [Citation(s) in RCA: 87] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2011] [Accepted: 02/26/2012] [Indexed: 11/19/2022] Open
Abstract
From computational simulations of a serotonin 2A receptor (5-HT2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization. The 5-HT2A receptor for the neurotransmitter serotonin (5-HT) belongs to family A (rhodopsin-like) G-protein coupled receptors (GPCRs), one of the most important classes of membrane proteins that are targeted by an extensive and diverse collection of external stimuli. Recently we learned that different ligands targeting the same GPCR can elicit different biological responses, but the mechanisms remain unknown. We address this fundamental question for the serotonin 5-HT2A receptor, because it is known to respond to the binding of structurally diverse ligands by producing similar stimuli in the cell, and to the binding of quite similar ligands with dramatically different responses. Molecular dynamics simulations of molecular models of the serotonin 5-HT2A receptor in complex with pharmacologically distinct ligands show the dynamic rearrangements of the receptor molecule to be different for these ligands, and the nature and extents of the rearrangements reflect the known pharmacological properties of the ligands as full, partial or inverse activators of the receptor. The different rearrangements of the receptor molecule are shown to produce different rearrangements of the surrounding membrane, a remodeling of the environment that can have differential ligand-determined effects on receptor function and association in the cell's membrane.
Collapse
|
33
|
Heifetz A, Morris GB, Biggin PC, Barker O, Fryatt T, Bentley J, Hallett D, Manikowski D, Pal S, Reifegerste R, Slack M, Law R. Study of Human Orexin-1 and -2 G-Protein-Coupled Receptors with Novel and Published Antagonists by Modeling, Molecular Dynamics Simulations, and Site-Directed Mutagenesis. Biochemistry 2012; 51:3178-97. [DOI: 10.1021/bi300136h] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Alexander Heifetz
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - G. Benjamin Morris
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Philip C. Biggin
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Oliver Barker
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - Tara Fryatt
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - Jonathan Bentley
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - David Hallett
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | | | - Sandeep Pal
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - Rita Reifegerste
- Evotec AG, Manfred Eigen Campus, Essener Bogen 7, 22419 Hamburg, Germany
| | - Mark Slack
- Evotec AG, Manfred Eigen Campus, Essener Bogen 7, 22419 Hamburg, Germany
| | - Richard Law
- Evotec (U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| |
Collapse
|
34
|
Bruno A, Costantino G. Molecular Dynamics Simulations of G Protein-Coupled Receptors. Mol Inform 2012; 31:222-30. [DOI: 10.1002/minf.201100138] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2011] [Accepted: 12/18/2011] [Indexed: 12/14/2022]
|
35
|
Fanelli F, De Benedetti PG. Update 1 of: computational modeling approaches to structure-function analysis of G protein-coupled receptors. Chem Rev 2011; 111:PR438-535. [PMID: 22165845 DOI: 10.1021/cr100437t] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Francesca Fanelli
- Dulbecco Telethon Institute, University of Modena and Reggio Emilia, via Campi 183, 41125 Modena, Italy.
| | | |
Collapse
|
36
|
Umanah GKE, Huang LY, Maccarone JM, Naider F, Becker JM. Changes in conformation at the cytoplasmic ends of the fifth and sixth transmembrane helices of a yeast G protein-coupled receptor in response to ligand binding. Biochemistry 2011; 50:6841-54. [PMID: 21728340 DOI: 10.1021/bi200254h] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The third intracellular loop (IL3) of G protein-coupled receptors (GPCRs) is an important contact domain between GPCRs and their G proteins. Previously, the IL3 of Ste2p, a Saccharomyces cerevisiae GPCR, was suggested to undergo a conformational change upon activation as detected by differential protease susceptibility in the presence and absence of ligand. In this study using disulfide cross-linking experiments we show that the Ste2p cytoplasmic ends of helix 5 (TM5) and helix 6 (TM6) that flank the amino and carboxyl sides of IL3 undergo conformational changes upon ligand binding, whereas the center of the IL3 loop does not. Single Cys substitution of residues in the middle of IL3 led to receptors that formed high levels of cross-linked Ste2p, whereas Cys substitution at the interface of IL3 and the contiguous cytoplasmic ends of TM5 and TM6 resulted in minimal disulfide-mediated cross-linked receptor. The alternating pattern of residues involved in cross-linking suggested the presence of a 3(10) helix in the middle of IL3. Agonist (WHWLQLKPGQPNleY) induced Ste2p activation reduced cross-linking mediated by Cys substitutions at the cytoplasmic ends of TM5 and TM6 but not by residues in the middle of IL3. Thus, the cytoplasmic ends of TM5 and TM6 undergo conformational change upon ligand binding. An α-factor antagonist (des-Trp, des-His-α-factor) did not influence disulfide-mediated Ste2p cross-linking, suggesting that the interaction of the N-terminus of α-factor with Ste2p is critical for inducing conformational changes at TM5 and TM6. We propose that the changes in conformation revealed for residues at the ends of TM5 and TM6 are affected by the presence of G protein but not G protein activation. This study provides new information about role of specific residues of a GPCR in signal transduction and how peptide ligand binding activates the receptor.
Collapse
Affiliation(s)
- George K E Umanah
- Department of Microbiology, University of Tennessee, Knoxville, Tennessee 37996, United States
| | | | | | | | | |
Collapse
|
37
|
Grossfield A. Recent progress in the study of G protein-coupled receptors with molecular dynamics computer simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2011; 1808:1868-78. [DOI: 10.1016/j.bbamem.2011.03.010] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/11/2011] [Revised: 02/23/2011] [Accepted: 03/21/2011] [Indexed: 01/28/2023]
|
38
|
Mobarec JC, Sanchez R, Filizola M. Modern homology modeling of G-protein coupled receptors: which structural template to use? J Med Chem 2010; 52:5207-16. [PMID: 19627087 DOI: 10.1021/jm9005252] [Citation(s) in RCA: 126] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The recent availability in the literature of new crystal structures of inactive G-protein coupled receptors (GPCRs) prompted us to study the extent to which these crystal structures constitute an advantage over the former prototypic rhodopsin template for homology modeling of the transmembrane (TM) region of human class A GPCRs. Our results suggest that better templates than those currently available are required by the majority of these GPCRs to generate homology models that are accurate enough for simple virtual screening aimed at computer-aided drug discovery. Thus, we investigated: (1) which class A GPCRs would have the highest impact as potential templates for homology modeling of other GPCRs, if their structures were solved, and (2) the extent to which multiple-template homology modeling (using all currently available GPCR crystal structures) provides an improvement over single-template homology modeling, as evaluated by the accuracy of rigid protein-flexible ligand docking on these models.
Collapse
Affiliation(s)
- Juan Carlos Mobarec
- Department of Structural and Chemical Biology, Mount Sinai School of Medicine, 1425 Madison Avenue, Box 1677, New York, NY 10029-6574, USA
| | | | | |
Collapse
|
39
|
Ge X, Roux B. Calculation of the standard binding free energy of sparsomycin to the ribosomal peptidyl-transferase P-site using molecular dynamics simulations with restraining potentials. J Mol Recognit 2010; 23:128-41. [PMID: 20151411 DOI: 10.1002/jmr.996] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The standard (absolute) binding free energy of the antibiotic sparsomycin with the 50S bacteria ribosomal subunit is calculated using molecular dynamics (MD) free energy perturbation (FEP) simulations with restraining potentials developed by Wang et al. [Biophys. J. 91:2798-2814 (2006)]. In the simulation protocol, restraining potentials are activated for the orientational and translational movements of the ligand relative to the binding site when it is decoupled from the binding pocket, and then released once the ligand fully interacts with the rest of the system. A reduced system is simulated to decrease the computational cost of the FEP/MD calculations and the effects of the surrounding atoms are incorporated using the generalized solvent boundary potential (GSBP) method. The loss of conformational freedom of the ligand upon binding is characterized using the potential of mean force (PMF) as a function of the root-mean-square deviation (RMSD) relative to the bound conformation. The number of water molecules in the binding pocket is allowed to fluctuate dynamically in response to the ligand during the calculations by combining FEP/MD with grand canonical Monte Carlo (GCMC) simulations. The calculated binding free energy is about -6 kcal/mol, which is in reasonable agreement with the experimental value. The information gleaned from this study provides new insight on the recognition of ribosome by sparsomycin and highlights the challenges in calculations of absolute binding free energies in these systems.
Collapse
Affiliation(s)
- Xiaoxia Ge
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, New York, NY, USA
| | | |
Collapse
|
40
|
Bortolato A, Mobarec JC, Provasi D, Filizola M. Progress in elucidating the structural and dynamic character of G Protein-Coupled Receptor oligomers for use in drug discovery. Curr Pharm Des 2010; 15:4017-25. [PMID: 20028319 DOI: 10.2174/138161209789824768] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
G Protein-Coupled Receptors (GPCRs) are the most targeted group of proteins for the development of therapeutic drugs. Until the last decade, structural information about this family of membrane proteins was relatively scarce, and their mechanisms of ligand binding and signal transduction were modeled on the assumption that GPCRs existed and functioned as monomeric entities. New crystal structures of native and engineered GPCRs, together with important biochemical and biophysical data that reveal structural details of the activation mechanism(s) of this receptor family, provide a valuable framework to improve dynamic molecular models of GPCRs with the ultimate goal of elucidating their allostery and functional selectivity. Since the dynamic movements of single GPCR protomers are likely to be affected by the presence of neighboring interacting subunits, oligomeric arrangements should be taken into account to improve the predictive ability of computer-assisted structural models of GPCRs for effective use in drug design.
Collapse
Affiliation(s)
- A Bortolato
- Department of Structural and Chemical Biology, Mount Sinai School of Medicine, New York, NY 10029-6574, USA
| | | | | | | |
Collapse
|
41
|
Bruno A, Guadix AE, Costantino G. Molecular dynamics simulation of the heterodimeric mGluR2/5HT(2A) complex. An atomistic resolution study of a potential new target in psychiatric conditions. J Chem Inf Model 2009; 49:1602-16. [PMID: 19422244 DOI: 10.1021/ci900067g] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Homo- and heterodimerization is becoming an assessed concept in G-protein coupled receptor (GPCR) pharmacology, and the notion that GPCRs may dimerize or oligomerize is allowing for a reinterpretation of some inconsistencies or anomalies and is providing medicinal chemists with potentially relevant novel molecular targets for a variety of therapeutic conditions. Recently, it has been reported that two unrelated GPCRs, namely class C metabotropic glutamate receptor type-2 (mGluR2) and class A 5HT(2A) serotoninergic receptor, can heterodimerize at the transmembrane domain level. We performed a 40 ns molecular dynamics simulation of the mGluR2/5HT(2A) heterocomplex constructed around a TM4/TM5 interface and embedded in an explicit phospholipidic bilayer surrounded by water molecules. In a separate experiment, the monomeric 5HT(2A) receptor was simulated for additional 40 ns under the same conditions. The analysis and the comparison of the two simulations allowed us to clearly identify a cross-talk between the two protomers and to put forward an effect of the heterodimerization on the shape of the binding pocket of 5HT(2A). This result provides the first molecular explanation for the reported allosteric effect of mGluR2 on 5HT(2A)-mediated response and suggests that the heterocomplex can be a more suitable target for in silico screening than the monomeric protomers.
Collapse
Affiliation(s)
- Agostino Bruno
- Dipartimento Farmaceutico, Via G. P. Usberti 27/A- Campus Universitario, Universita degli Studi di Parma, 43124 Parma, Italy
| | | | | |
Collapse
|
42
|
Coarse-grained modeling of allosteric regulation in protein receptors. Proc Natl Acad Sci U S A 2009; 106:14253-8. [PMID: 19706508 DOI: 10.1073/pnas.0901811106] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Allosteric regulation provides highly specific ligand recognition and signaling by transmembrane protein receptors. Unlike functions of protein molecular machines that rely on large-scale conformational transitions, signal transduction in receptors appears to be mediated by more subtle structural motions that are difficult to identify. We describe a theoretical model for allosteric regulation in receptors that addresses a fundamental riddle of signaling: What are the structural origins of the receptor agonism (specific signaling response to ligand binding)? The model suggests that different signaling pathways in bovine rhodopsin or human beta(2)-adrenergic receptor can be mediated by specific structural motions in the receptors. We discuss implications for understanding the receptor agonism, particularly the recently observed "biased agonism" (selected activation of specific signaling pathways), and for developing rational structure-based drug-design strategies.
Collapse
|
43
|
Khelashvili G, Grossfield A, Feller SE, Pitman MC, Weinstein H. Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations. Proteins 2009; 76:403-17. [PMID: 19173312 PMCID: PMC4101808 DOI: 10.1002/prot.22355] [Citation(s) in RCA: 100] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1-TM2-TM7-helix 8(H8) functional network that comprises the highly conserved NPxxY(x)(5,6)F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high-density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X-ray crystallography data on beta(2)-adrenergic GPCR; and the intracellular ends of TM2-TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid-body motions of the TMs in the TM1-TM2-TM7-H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function.
Collapse
Affiliation(s)
- George Khelashvili
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, New York, New York 10021, USA.
| | | | | | | | | |
Collapse
|
44
|
Parker MS, Sah R, Balasubramaniam A, Sallee FR, Sweatman T, Park EA, Parker SL. Dimers of the neuropeptide Y (NPY) Y2 receptor show asymmetry in agonist affinity and association with G proteins. J Recept Signal Transduct Res 2009; 28:437-51. [PMID: 18946765 DOI: 10.1080/10799890802447423] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
In conditions precluding activation of G proteins, the binding of agonists to dimers of the neuropeptide Y (NPY) Y2 receptor shows two components of similar size, but differing in affinity. The dimers of all NPY receptors are solubilized as approximately 180-kDa complexes containing one G protein alpha beta gamma trimer. These heteropentamers are stable to excess agonists, chelators, and alkylators. However, dispersion in the weak surfactant cholate releases approximately 300-kDa complexes. These findings indicate that both protomers in the Y2 dimer are associated with G protein heterotrimers, but the extent of interaction depends on affinity for the agonist peptide. The G protein in contact with the first-liganded, higher-affinity protomer should have a stronger interaction with the receptor and a larger probability of activation.
Collapse
Affiliation(s)
- M S Parker
- Department of Molecular Cell Sciences, University of Memphis, Memphis, Tennessee, USA.
| | | | | | | | | | | | | |
Collapse
|
45
|
Chang CW, Poteet E, Schetz JA, Gümüş ZH, Weinstein H. Towards a quantitative representation of the cell signaling mechanisms of hallucinogens: measurement and mathematical modeling of 5-HT1A and 5-HT2A receptor-mediated ERK1/2 activation. Neuropharmacology 2008; 56 Suppl 1:213-25. [PMID: 18762202 PMCID: PMC2635340 DOI: 10.1016/j.neuropharm.2008.07.049] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2008] [Revised: 07/23/2008] [Accepted: 07/25/2008] [Indexed: 11/27/2022]
Abstract
Through a multidisciplinary approach involving experimental and computational studies, we address quantitative aspects of signaling mechanisms triggered in the cell by the receptor targets of hallucinogenic drugs, the serotonin 5-HT2A receptors. To reveal the properties of the signaling pathways, and the way in which responses elicited through these receptors alone and in combination with other serotonin receptors' subtypes (the 5-HT1AR), we developed a detailed mathematical model of receptor-mediated ERK1/2 activation in cells expressing the 5-HT1A and 5-HT2A subtypes individually, and together. In parallel, we measured experimentally the activation of ERK1/2 by the action of selective agonists on these receptors expressed in HEK293 cells. We show here that the 5-HT1AR agonist Xaliproden HCl elicited transient activation of ERK1/2 by phosphorylation, whereas 5-HT2AR activation by TCB-2 led to higher, and more sustained responses. The 5-HT2AR response dominated the MAPK signaling pathway when co-expressed with 5-HT1AR, and diminution of the response by the 5-HT2AR antagonist Ketanserin could not be rescued by the 5-HT1AR agonist. Computational simulations produced qualitative results in good agreement with these experimental data, and parameter optimization made this agreement quantitative. In silico simulation experiments suggest that the deletion of the positive regulators PKC in the 5-HT2AR pathway, or PLA2 in the combined 5-HT1A/2AR model greatly decreased the basal level of active ERK1/2. Deletion of negative regulators of MKP and PP2A in 5-HT1AR and 5-HT2AR models was found to have even stronger effects. Under various parameter sets, simulation results implied that the extent of constitutive activity in a particular tissue and the specific drug efficacy properties may determine the distinct dynamics of the 5-HT receptor-mediated ERK1/2 activation pathways. Thus, the mathematical models are useful exploratory tools in the ongoing efforts to establish a mechanistic understanding and an experimentally testable representation of hallucinogen-specific signaling in the cellular machinery, and can be refined with quantitative, function-related information.
Collapse
MESH Headings
- Cell Line, Transformed
- Computer Simulation
- Dose-Response Relationship, Drug
- Extracellular Signal-Regulated MAP Kinases/metabolism
- Hallucinogens/pharmacology
- Humans
- Models, Biological
- Protein Binding/drug effects
- Radioligand Assay/methods
- Receptor, Serotonin, 5-HT1A/genetics
- Receptor, Serotonin, 5-HT1A/metabolism
- Receptor, Serotonin, 5-HT2A/genetics
- Receptor, Serotonin, 5-HT2A/metabolism
- Signal Transduction/drug effects
- Signal Transduction/physiology
- Time Factors
- Transfection
Collapse
Affiliation(s)
- Chiung-wen Chang
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, 1300 York Ave, New York, NY 10021 USA
| | - Ethan Poteet
- Department of Pharmacology and Neuroscience, University of North Texas Health Science Center, 3500 Camp Bowie Blvd. Fort Worth, TX 76107
| | - John A. Schetz
- Department of Pharmacology and Neuroscience, University of North Texas Health Science Center, 3500 Camp Bowie Blvd. Fort Worth, TX 76107
| | - Zeynep H. Gümüş
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, 1300 York Ave, New York, NY 10021 USA
- The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Medical College of Cornell University, 1300 York Ave, New York, NY 10021 USA
| | - Harel Weinstein
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, 1300 York Ave, New York, NY 10021 USA
- The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Medical College of Cornell University, 1300 York Ave, New York, NY 10021 USA
| |
Collapse
|
46
|
Han DS, Wang SX, Weinstein H. Active state-like conformational elements in the beta2-AR and a photoactivated intermediate of rhodopsin identified by dynamic properties of GPCRs. Biochemistry 2008; 47:7317-21. [PMID: 18558776 PMCID: PMC2664832 DOI: 10.1021/bi800442g] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
G-Protein-coupled receptors (GPCRs) adopt various functionally relevant conformational states in cell signaling processes. Recently determined crystal structures of rhodopsin and the beta 2-adrenergic receptor (beta 2-AR) offer insight into previously uncharacterized active conformations, but the molecular states of these GPCRs are likely to contain both inactive and active-like conformational elements. We have identified conformational rearrangements in the dynamics of the TM7-HX8 segment that relate to the properties of the conserved NPxxY(x)5,6F motif and show that they can be used to identify active state-like conformational elements in the corresponding regions of the new structures of rhodopsin and the beta 2-AR.
Collapse
Affiliation(s)
- Daniel S Han
- Department of Physiology and Biophysics, Weill Medical College, Cornell University, 1300 York Avenue, New York, New York 10021, USA
| | | | | |
Collapse
|
47
|
Niv MY, Filizola M. Influence of oligomerization on the dynamics of G-protein coupled receptors as assessed by normal mode analysis. Proteins 2008; 71:575-86. [PMID: 17963239 DOI: 10.1002/prot.21787] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The recently discovered impact of oligomerization on G-protein coupled receptor (GPCR) function further complicates the already challenging goal of unraveling the molecular and dynamic mechanisms of these receptors. To help understand the effect of oligomerization on the dynamics of GPCRs, we have compared the motion of monomeric, dimeric, and tetrameric arrangements of the prototypic GPCR rhodopsin, using an approximate-yet powerful-normal mode analysis (NMA) technique termed elastic network model (ENM). Moreover, we have used ENM to discriminate between putative dynamic mechanisms likely to account for the recently observed conformational rearrangement of the TM4,5-TM4,5 dimerization interface of GPCRs that occurs upon activation. Our results indicate: (1) significant perturbation of the normal modes (NMs) of the rhodopsin monomer upon oligomerization, which is mainly manifested at interfacial regions; (2) increased positive correlation among the transmembrane domains (TMs) and between the extracellular loop (EL) and TM regions of the rhodopsin protomer; (3) highest interresidue positive correlation at the interfaces between protomers; and (4) experimentally testable hypotheses of differential motional changes within different putative oligomeric arrangements.
Collapse
Affiliation(s)
- Masha Y Niv
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, New York, NY, USA
| | | |
Collapse
|
48
|
Dalrymple MB, Pfleger KDG, Eidne KA. G protein-coupled receptor dimers: functional consequences, disease states and drug targets. Pharmacol Ther 2008; 118:359-71. [PMID: 18486226 DOI: 10.1016/j.pharmthera.2008.03.004] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2008] [Accepted: 03/14/2008] [Indexed: 10/22/2022]
Abstract
With an ever-expanding need for reliable therapeutic agents that are highly effective and exhibit minimal deleterious side effects, a greater understanding of the mechanisms underlying G protein-coupled receptor (GPCR) regulation is fundamental. GPCRs comprise more than 30% of all therapeutic drug targets and it is likely that this will only increase as more orphan GPCRs are identified. The past decade has seen a dramatic shift in the prevailing concept of how GPCRs function, in particular the growing acceptance that GPCRs are capable of interacting with one another at a molecular level to form complexes, with significantly different pharmacological properties to their monomeric selves. While the ability of like-receptors to associate and form homodimers raises some interesting mechanistic issues, the possibility that unlike-receptors could heterodimerise in certain tissue types, producing a functionally unique signalling complex that binds specific ligands, provides an invaluable opportunity to refine and redefine pharmacological interventions with greater specificity and efficacy.
Collapse
Affiliation(s)
- Matthew B Dalrymple
- Laboratory for Molecular Endocrinology - GPCRs, Western Australian Institute for Medical Research and Centre for Medical Research, University of Western Australia, Nedlands, Perth, WA 6009, Australia
| | | | | |
Collapse
|
49
|
Mobarec JC, Filizola M. Advances in the Development and Application of Computational Methodologies for Structural Modeling of G-Protein Coupled Receptors. Expert Opin Drug Discov 2008; 3:343-355. [PMID: 19672320 DOI: 10.1517/17460441.3.3.343] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
BACKGROUND: Despite the large amount of experimental data accumulated in the past decade on G-protein coupled receptor (GPCR) structure and function, understanding of the molecular mechanisms underlying GPCR signaling is still far from being complete, thus impairing the design of effective and selective pharmaceuticals. OBJECTIVE: Understanding of GPCR function has been challenged even further by more recent experimental evidence that several of these receptors are organized in the cell membrane as homo- or hetero-oligomers, and that they may exhibit unique pharmacological properties. Given the complexity of these new signaling systems, researcher's efforts are turning increasingly to molecular modeling, bioinformatics and computational simulations for mechanistic insights of GPCR functional plasticity. METHODS: We review here current advances in the development and application of computational approaches to improve prediction of GPCR structure and dynamics, thus enhancing current understanding of GPCR signaling. RESULTS/CONCLUSIONS: Models resulting from use of these computational approaches further supported by experiments are expected to help elucidate the complex allosterism that propagates through GPCR complexes, ultimately aiming at successful structure-based rational drug design.
Collapse
Affiliation(s)
- Juan Carlos Mobarec
- Department of Structural and Chemical Biology, Mount Sinai School of Medicine, Icahn Medical Institute Building, 1425 Madison Avenue, Box 1677, New York, NY 10029-6574, Tel: 212-241-8634
| | | |
Collapse
|
50
|
Gonzalez A, Duran LS, Araya-Secchi R, Garate JA, Pessoa-Mahana CD, Lagos CF, Perez-Acle T. Computational modeling study of functional microdomains in cannabinoid receptor type 1. Bioorg Med Chem 2008; 16:4378-89. [PMID: 18342519 DOI: 10.1016/j.bmc.2008.02.070] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2007] [Revised: 02/18/2008] [Accepted: 02/21/2008] [Indexed: 10/22/2022]
Abstract
The seven transmembrane helices (TMH) G-protein-coupled receptors (GPCRs) constitute one of the largest superfamily of signaling proteins found in mammals. Some of its members, in which the cannabinoid (CB) receptors are included, stand out because their functional states can be modulated by a broad spectrum of effector molecules. The relative ligand promiscuity exhibited by these receptors could be related with particular attributes conferred by their molecular architecture and represents a motivating issue to be explored. In this regard, this study represents an effort to investigate the cannabinoid receptor type 1 (CB1) ligand recognition plasticity, using comparative modeling, molecular dynamics (MD) simulations and docking. Our results suggest that a cooperative set of subtle structural rearrangements within the TMHs provide to the CB1 protein the plasticity to reach alternate configurations. These changes include the relaxation of intramolecular constraints, the rotations, translations and kinks of the majority of TMHs and the reorganization of the ligand binding cavities.
Collapse
Affiliation(s)
- Angel Gonzalez
- Centre for Bioinformatics CBUC, Faculty of Biological Sciences, Pontificia Universidad Católica de Chile, Portugal 49-6513492, Santiago, Chile.
| | | | | | | | | | | | | |
Collapse
|