1
|
Ghosh S, de March CA, Branciamore S, Kaleem S, Matsunami H, Vaidehi N. Sequence coevolution and structure stabilization modulate olfactory receptor expression. Biophys J 2022; 121:830-840. [PMID: 35065915 PMCID: PMC8947990 DOI: 10.1016/j.bpj.2022.01.015] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Revised: 12/14/2021] [Accepted: 01/19/2022] [Indexed: 11/29/2022] Open
Abstract
Olfactory receptors (ORs) belong to class A G-protein coupled receptors (GPCRs) and are activated by a variety of odorants. To date, there is no three-dimensional structure of an OR available. One of the major bottlenecks in obtaining purified protein for structural studies of ORs is their poor expression in heterologous cells. To design mutants that enhance expression and thereby enable protein purification, we first identified computable physical properties that recapitulate OR and class A GPCR expression and further conducted an iterative computational prediction-experimental test cycle and generated human OR mutants that express as high as biogenic amine receptors for which structures have been solved. In the process of developing the computational method to recapitulate the expression of ORs in membranes, we identified properties, such as amino acid sequence coevolution, and the strength of the interactions between intracellular loop 1 (ICL1) and the helix 8 region of ORs, to enhance their heterologous expression. We identified mutations that are directly located in these regions as well as other mutations not located in these regions but allosterically strengthen the ICL1-helix 8 enhance expression. These mutants also showed functional responses to known odorants. This method to enhance heterologous expression of mammalian ORs will facilitate high-throughput "deorphanization" of ORs, and enable OR purification for biochemical and structural studies to understand odorant-OR interactions.
Collapse
Affiliation(s)
- Soumadwip Ghosh
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA, USA
| | - Claire A. de March
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA
| | - Sergio Branciamore
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA, USA
| | - Sahar Kaleem
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA
| | - Hiroaki Matsunami
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA; Department of Neurobiology, Duke Institute for Brain Sciences, Duke University School of Medicine, Durham, NC, USA.
| | - Nagarajan Vaidehi
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA, USA.
| |
Collapse
|
2
|
Suresh R, Subramaniam V. Molecular dynamics simulation involved in expounding the activation of adrenoceptors by sympathetic nervous system signaling. Struct Chem 2020. [DOI: 10.1007/s11224-020-01553-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
3
|
Muk S, Ghosh S, Achuthan S, Chen X, Yao X, Sandhu M, Griffor MC, Fennell KF, Che Y, Shanmugasundaram V, Qiu X, Tate CG, Vaidehi N. Machine Learning for Prioritization of Thermostabilizing Mutations for G-Protein Coupled Receptors. Biophys J 2019; 117:2228-2239. [PMID: 31703801 DOI: 10.1016/j.bpj.2019.10.023] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2019] [Revised: 10/11/2019] [Accepted: 10/17/2019] [Indexed: 01/01/2023] Open
Abstract
Although the three-dimensional structures of G-protein coupled receptors (GPCRs), the largest superfamily of drug targets, have enabled structure-based drug design, there are no structures available for 87% of GPCRs. This is due to the stiff challenge in purifying the inherently flexible GPCRs. Identifying thermostabilized mutant GPCRs via systematic alanine scanning mutations has been a successful strategy in stabilizing GPCRs, but it remains a daunting task for each GPCR. We developed a computational method that combines sequence-, structure-, and dynamics-based molecular properties of GPCRs that recapitulate GPCR stability, with four different machine learning methods to predict thermostable mutations ahead of experiments. This method has been trained on thermostability data for 1231 mutants, the largest publicly available data set. A blind prediction for thermostable mutations of the complement factor C5a receptor 1 retrieved 36% of the thermostable mutants in the top 50 prioritized mutants compared to 3% in the first 50 attempts using systematic alanine scanning.
Collapse
Affiliation(s)
- Sanychen Muk
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, California
| | - Soumadwip Ghosh
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, California
| | - Srisairam Achuthan
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, California
| | | | - XiaoJie Yao
- Discovery Sciences, Pfizer, Groton, Connecticut
| | - Manbir Sandhu
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, California
| | | | | | - Ye Che
- Discovery Sciences, Pfizer, Groton, Connecticut
| | | | - Xiayang Qiu
- Discovery Sciences, Pfizer, Groton, Connecticut
| | | | - Nagarajan Vaidehi
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, California.
| |
Collapse
|
4
|
Jana S, Ghosh S, Muk S, Levy B, Vaidehi N. Prediction of Conformation Specific Thermostabilizing Mutations for Class A G Protein-Coupled Receptors. J Chem Inf Model 2019; 59:3744-3754. [PMID: 31408606 DOI: 10.1021/acs.jcim.9b00175] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
G protein-coupled receptors (GPCRs) are highly flexible and prone to denaturation during protein extraction in detergents and purification. This poses a huge challenge to purify a conformationally homogeneous solution of GPCRs. Thermostabilizing mutations have been used widely to purify and obtain crystal structures of several GPCRs. However, identifying thermostabilizing mutations for GPCRs remains a tedious and expensive task as they are not transferable even among closely related GPCRs. Additionally, the mutations stabilizing one conformational state of a GPCR do not always stabilize other conformational state(s) of the same GPCR. Previously we developed a computational method, LiticonDesign, for rapid prediction of thermostabilizing mutations for a specific GPCR conformation. In this study, we have used LiticonDesign to predict thermostabilizing mutations for the agonist bound active-intermediate state of the human adenosine receptor (A2AR) using the structure of the inactive state of the same GPCR and vice versa. Our study shows that the thermostable mutation predictions using LiticonDesign, for an active-intermediate state of a GPCR (A2AR in our case), requires a homology model that is derived from an active/active-intermediate state GPCR structure as a template. Similarly, the homology models derived from inactive state GPCR conformations are better in predicting the thermostable mutations for the inactive state of A2AR. Overall, LiticonDesign method is not only efficient in predicting thermostabilizing mutations for a given GPCR sequence but also can recover conformation specific mutations for a state of interest, if a suitable starting structure of desired conformation is chosen.
Collapse
Affiliation(s)
- Suvamay Jana
- Department of Computation and Quantitative Medicine , Beckman Research Institute of the City of Hope , 1500 E. Duarte Road , Duarte , California 91010 , United States
| | - Soumadwip Ghosh
- Department of Computation and Quantitative Medicine , Beckman Research Institute of the City of Hope , 1500 E. Duarte Road , Duarte , California 91010 , United States
| | - Sanychen Muk
- Department of Computation and Quantitative Medicine , Beckman Research Institute of the City of Hope , 1500 E. Duarte Road , Duarte , California 91010 , United States
| | - Benjamin Levy
- Department of Computation and Quantitative Medicine , Beckman Research Institute of the City of Hope , 1500 E. Duarte Road , Duarte , California 91010 , United States
| | - Nagarajan Vaidehi
- Department of Computation and Quantitative Medicine , Beckman Research Institute of the City of Hope , 1500 E. Duarte Road , Duarte , California 91010 , United States
| |
Collapse
|
5
|
Ghosh S, Bierig T, Lee S, Jana S, Löhle A, Schnapp G, Tautermann CS, Vaidehi N. Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors. J Chem Theory Comput 2018; 14:6574-6585. [PMID: 30359017 DOI: 10.1021/acs.jctc.8b00602] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D6.33 and K303E7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D6.33 and R307E7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D6.33/R307E7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.
Collapse
Affiliation(s)
- Soumadwip Ghosh
- Department of Molecular Imaging and Therapy , Beckman Research Institute of the City of Hope , 1500 East Duarte Road , Duarte , California 91010 , United States
| | - Tobias Bierig
- Department of Medicinal Chemistry , Boehringer Ingelheim Pharma GmbH & Company KG , Birkendorfer Strasse 65 , D-88397 Biberach an der Riss , Germany
| | | | | | | | - Gisela Schnapp
- Department of Medicinal Chemistry , Boehringer Ingelheim Pharma GmbH & Company KG , Birkendorfer Strasse 65 , D-88397 Biberach an der Riss , Germany
| | - Christofer S Tautermann
- Department of Medicinal Chemistry , Boehringer Ingelheim Pharma GmbH & Company KG , Birkendorfer Strasse 65 , D-88397 Biberach an der Riss , Germany
| | - Nagarajan Vaidehi
- Department of Molecular Imaging and Therapy , Beckman Research Institute of the City of Hope , 1500 East Duarte Road , Duarte , California 91010 , United States
| |
Collapse
|
6
|
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
|
7
|
How Can Mutations Thermostabilize G-Protein-Coupled Receptors? Trends Pharmacol Sci 2015; 37:37-46. [PMID: 26547284 DOI: 10.1016/j.tips.2015.09.005] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2015] [Revised: 09/14/2015] [Accepted: 09/21/2015] [Indexed: 01/04/2023]
Abstract
Structures of over 30 different G-protein-coupled receptors (GPCRs) have advanced our understanding of cell signaling and have provided a foundation for structure-guided drug design. This exciting progress has required the development of three complementary methods to facilitate GPCR crystallization, one of which is the thermostabilization of receptors by systematic mutagenesis. However, the reason why a particular mutation, or combination of mutations, stabilizes the receptor is not always evident from a static crystal structure. Molecular dynamics (MD) simulations have been used to identify and estimate the energetic factors that affect thermostability through comparing the dynamics of the thermostabilized receptors with structure-based models of the wild-type receptor. The data indicate that receptors are stabilized through a combination of factors, including an increase in receptor rigidity, a decrease in collective motion, reduced stress at specific residues, and the presence of ordered water molecules. Predicting thermostabilizing mutations computationally represents a major challenge for the field.
Collapse
|
8
|
Xiao X, Zeng X, Yuan Y, Gao N, Guo Y, Pu X, Li M. Understanding the conformation transition in the activation pathway of β2 adrenergic receptor via a targeted molecular dynamics simulation. Phys Chem Chem Phys 2015; 17:2512-22. [DOI: 10.1039/c4cp04528a] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The conformation transition in the activation pathway of β2 adrenergic receptor was explored mainly using a target molecular dynamics simulation.
Collapse
Affiliation(s)
- Xiuchan Xiao
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Xiaojun Zeng
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Yuan Yuan
- College of Management
- Southwest University for Nationalities
- Chengdu
- People's Republic of China
| | - Nan Gao
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Yanzhi Guo
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Xuemei Pu
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Menglong Li
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| |
Collapse
|
9
|
Bhattacharya S, Lee S, Grisshammer R, Tate CG, Vaidehi N. Rapid Computational Prediction of Thermostabilizing Mutations for G Protein-Coupled Receptors. J Chem Theory Comput 2014; 10:5149-5160. [PMID: 25400524 PMCID: PMC4230369 DOI: 10.1021/ct500616v] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2014] [Indexed: 01/22/2023]
Abstract
![]()
G protein-coupled
receptors (GPCRs) are highly dynamic and often
denature when extracted in detergents. Deriving thermostable mutants
has been a successful strategy to stabilize GPCRs in detergents, but
this process is experimentally tedious. We have developed a computational
method to predict the position of the thermostabilizing mutations
for a given GPCR sequence. We have validated the method against experimentally
measured thermostability data for single mutants of the β1-adrenergic receptor (β1AR), adenosine A2A receptor (A2AR) and neurotensin receptor 1 (NTSR1).
To make these predictions we started from homology models of these
receptors of varying accuracies and generated an ensemble of conformations
by sampling the rigid body degrees of freedom of transmembrane helices.
Then, an all-atom force field function was used to calculate the enthalpy
gain, known as the “stability score” upon mutation of
every residue, in these receptor structures, to alanine. For all three
receptors, β1AR, A2AR, and NTSR1, we observed
that mutations of hydrophobic residues in the transmembrane domain
to alanine that have high stability scores correlate with high experimental
thermostability. The prediction using the stability score improves
when using an ensemble of receptor conformations compared to a single
structure, showing that receptor flexibility is important. We also
find that our previously developed LITiCon method for generating conformation
ensembles is similar in performance to predictions using ensembles
obtained from microseconds of molecular dynamics simulations (which
is computationally hundred times slower than LITiCon). We improved
the thermostability prediction by including other properties such
as residue-based stress and the extent of allosteric communication
by each residue in the stability score. Our method is the first step
toward a computational method for rapid prediction of thermostable
mutants of GPCRs.
Collapse
Affiliation(s)
- Supriyo Bhattacharya
- Division of Immunology, Beckman Research Institute of the City of Hope , 1500 East Duarte Rd, Duarte, California 91010, United States
| | - Sangbae Lee
- Division of Immunology, Beckman Research Institute of the City of Hope , 1500 East Duarte Rd, Duarte, California 91010, United States
| | - Reinhard Grisshammer
- Membrane Protein Structure Function Unit, National Institute of Neurological Disorders and Stroke, National Institutes of Health , Department of Health and Human Services, Rockville, Maryland 20852, United States
| | - Christopher G Tate
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus , Francis Crick Avenue, Cambridge CB2 0QH, United Kingdom
| | - Nagarajan Vaidehi
- Division of Immunology, Beckman Research Institute of the City of Hope , 1500 East Duarte Rd, Duarte, California 91010, United States
| |
Collapse
|
10
|
Vaidehi N, Bhattacharya S, Larsen AB. Structure and dynamics of G-protein coupled receptors. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 796:37-54. [PMID: 24158800 DOI: 10.1007/978-94-007-7423-0_3] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
G-protein coupled receptors (GPCRs) are seven helical transmembrane proteins that mediate cell-to-cell communication. They also form the largest superfamily of drug targets. Hence detailed studies of the three dimensional structure and dynamics are critical to understanding the functional role of GPCRs in signal transduction pathways, and for drug design. In this chapter we compare the features of the crystal structures of various biogenic amine receptors, such as β1 and β2 adrenergic receptors, dopamine D3 receptor, M2 and M3 muscarinic acetylcholine receptors. This analysis revealed that conserved residues are located facing the inside of the transmembrane domain in these GPCRs improving the efficiency of packing of these structures. The NMR structure of the chemokine receptor CXCR1 without any ligand bound, shows significant dynamics of the transmembrane domain, especially the helical kink angle on the transmembrane helix6. The activation mechanism of the β2-adrenergic receptor has been studied using multiscale computational methods. The results of these studies showed that the receptor without any ligand bound, samples conformations that resemble some of the structural characteristics of the active state of the receptor. Ligand binding stabilizes some of the conformations already sampled by the apo receptor. This was later observed in the NMR study of the dynamics of human β2-adrenergic receptor. The dynamic nature of GPCRs leads to a challenge in obtaining purified receptors for biophysical studies. Deriving thermostable mutants of GPCRs has been a successful strategy to reduce the conformational heterogeneity and stabilize the receptors. This has lead to several crystal structures of GPCRs. However, the cause of how these mutations lead to thermostability is not clear. Computational studies are beginning to shed some insight into the possible structural basis for the thermostability. Molecular Dynamics simulations studying the conformational ensemble of thermostable mutants have shown that the stability could arise from both enthalpic and entropic factors. There are regions of high stress in the wild type GPCR that gets relieved upon mutation conferring thermostability.
Collapse
Affiliation(s)
- Nagarajan Vaidehi
- Division of Immunology, Beckman Research Institute of the City of Hope, 1500, E. Duarte Road, Duarte, CA, 91010, USA,
| | | | | |
Collapse
|
11
|
Niesen MJM, Bhattacharya S, Grisshammer R, Tate CG, Vaidehi N. Thermostabilization of the β1-adrenergic receptor correlates with increased entropy of the inactive state. J Phys Chem B 2013; 117:7283-91. [PMID: 23697892 DOI: 10.1021/jp403207c] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The dynamic nature of GPCRs is a major hurdle in their purification and crystallization. Thermostabilization can facilitate GPCR structure determination, as has been shown by the structure of the thermostabilized β1-adrenergic receptor (β1AR) mutant, m23-β1AR, which has been thermostabilized in the inactive state. However, it is unclear from the structure how the six thermostabilizing mutations in m23-β1AR affect receptor dynamics. We have used molecular dynamics simulations in explicit solvent to compare the conformational ensembles for both wild type β1AR (wt-β1AR) and m23-β1AR. Thermostabilization results in an increase in the number of accessible microscopic conformational states within the inactive state ensemble, effectively increasing the side chain entropy of the inactive state at room temperature, while suppressing large-scale main chain conformational changes that lead to activation. We identified several diverse mechanisms of thermostabilization upon mutation. These include decrease of long-range correlated movement between residues in the G-protein coupling site to the extracellular region (Y227A(5.58), F338M(7.48)), formation of new hydrogen bonds (R68S), and reduction of local stress (Y227(5.58), F327(7.37), and F338(7.48)). This study provides insights into microscopic mechanisms underlying thermostability that leads to an understanding of the effect of these mutations on the structure of the receptor.
Collapse
Affiliation(s)
- Michiel J M Niesen
- Division of Immunology, Beckman Research Institute of the City of Hope, 1500 E. Duarte Road, Duarte, California 91010, USA
| | | | | | | | | |
Collapse
|
12
|
Levit A, Barak D, Behrens M, Meyerhof W, Niv MY. Homology model-assisted elucidation of binding sites in GPCRs. Methods Mol Biol 2013; 914:179-205. [PMID: 22976029 DOI: 10.1007/978-1-62703-023-6_11] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
G protein-coupled receptors (GPCRs) are important mediators of cell signaling and a major family of drug targets. Despite recent breakthroughs, experimental elucidation of GPCR structures remains a formidable challenge. Homology modeling of 3D structures of GPCRs provides a practical tool for elucidating the structural determinants governing the interactions of these important receptors with their ligands. The working model of the binding site can then be used for virtual screening of additional ligands that may fit this site, for determining and comparing specificity profiles of related receptors, and for structure-based design of agonists and antagonists. The current review presents the protocol and enumerates the steps for modeling and validating the residues involved in ligand binding. The main stages include (a) modeling the receptor structure using an automated fragment-based approach, (b) predicting potential binding pockets, (c) docking known binders, (d) analyzing predicted interactions and comparing with positions that have been shown to bind ligands in other receptors, (e) validating the structural model by mutagenesis.
Collapse
Affiliation(s)
- Anat Levit
- Institute of Biochemistry, Food Science, and Nutrition, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University, Rehovot, Israel
| | | | | | | | | |
Collapse
|
13
|
LITiCon: a discrete conformational sampling computational method for mapping various functionally selective conformational states of transmembrane helical proteins. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2012; 914:167-78. [PMID: 22976028 DOI: 10.1007/978-1-62703-023-6_10] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/20/2023]
Abstract
G-Protein-coupled receptors (GPCRs) are seven helical transmembrane proteins that mediate cell signaling thereby controlling many important physiological and pathological functions. GPCRs get activated upon ligand binding and trigger the signal transduction process. GPCRs exist in multiple inactive and active conformations, and there is a finite population of the active and inactive states even in the ligand-free condition. An understanding of the nature of the conformational ensemble sampled by GPCRs and the atomic level mechanism of the conformational transitions require a combination of computational methods and experimental techniques. We have developed a coarse grained discrete conformational sampling computational method called "LITiCon" to map the conformational ensemble sampled by GPCRs in the presence and absence of ligands. The LITiCon method can also be used to predict functional selective conformational states starting from the inactive state of the receptor. LITiCon has been applied to map the conformational ensemble of β2-adrenergic receptor, a class A GPCR. We have shown that β2-adrenergic receptor samples a larger conformational space in the ligand-free state and that different ligands select and stabilize conformations from this ensemble. In this review we describe the LITiCon method in detail and elucidate the uses and pitfalls of this method.
Collapse
|
14
|
Trzaskowski B, Latek D, Yuan S, Ghoshdastider U, Debinski A, Filipek S. Action of molecular switches in GPCRs--theoretical and experimental studies. Curr Med Chem 2012; 19:1090-109. [PMID: 22300046 PMCID: PMC3343417 DOI: 10.2174/092986712799320556] [Citation(s) in RCA: 333] [Impact Index Per Article: 27.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2011] [Revised: 12/30/2011] [Accepted: 01/02/2012] [Indexed: 01/14/2023]
Abstract
G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.
Collapse
Affiliation(s)
- B Trzaskowski
- Faculty of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | | | | | | | | | | |
Collapse
|
15
|
A crystal clear solution for determining G-protein-coupled receptor structures. Trends Biochem Sci 2012; 37:343-52. [PMID: 22784935 DOI: 10.1016/j.tibs.2012.06.003] [Citation(s) in RCA: 106] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2012] [Revised: 05/31/2012] [Accepted: 06/05/2012] [Indexed: 11/21/2022]
Abstract
G-protein-coupled receptors (GPCRs) are medically important membrane proteins that are targeted by over 30% of small molecule drugs. At the time of writing, 15 unique GPCR structures have been determined, with 77 structures deposited in the PDB database, which offers new opportunities for drug development and for understanding the molecular mechanisms of GPCR activation. Many different factors have contributed to this success, but if there is one single factor that can be singled out as the foundation for producing well-diffracting GPCR crystals, it is the stabilisation of the detergent-solubilised receptor-ligand complex. This review will focus predominantly on one of the successful strategies for the stabilisation of GPCRs, namely the thermostabilisation of GPCRs using systematic mutagenesis coupled with thermostability assays. Structures of thermostabilised GPCRs bound to a wide variety of ligands have been determined, which has led to an understanding of ligand specificity; why some ligands act as agonists as opposed to partial or inverse agonists; and the structural basis for receptor activation.
Collapse
|
16
|
Wheatley M, Wootten D, Conner MT, Simms J, Kendrick R, Logan RT, Poyner DR, Barwell J. Lifting the lid on GPCRs: the role of extracellular loops. Br J Pharmacol 2012; 165:1688-1703. [PMID: 21864311 DOI: 10.1111/j.1476-5381.2011.01629.x] [Citation(s) in RCA: 224] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
Abstract
GPCRs exhibit a common architecture of seven transmembrane helices (TMs) linked by intracellular loops and extracellular loops (ECLs). Given their peripheral location to the site of G-protein interaction, it might be assumed that ECL segments merely link the important TMs within the helical bundle of the receptor. However, compelling evidence has emerged in recent years revealing a critical role for ECLs in many fundamental aspects of GPCR function, which supported by recent GPCR crystal structures has provided mechanistic insights. This review will present current understanding of the key roles of ECLs in ligand binding, activation and regulation of both family A and family B GPCRs.
Collapse
Affiliation(s)
- M Wheatley
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - D Wootten
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - M T Conner
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - J Simms
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - R Kendrick
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - R T Logan
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - D R Poyner
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| | - J Barwell
- School of Biosciences, University of Birmingham, Birmingham, UKDrug Discovery Biology Laboratory, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, AustraliaDepartment of Pharmacology, Monash University, Parkville, Victoria, AustraliaSchool of Life and Health Sciences, Aston University, Aston Triangle, Birmingham, UK
| |
Collapse
|
17
|
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
|
18
|
Baker JG, Proudman RGW, Tate CG. The pharmacological effects of the thermostabilising (m23) mutations and intra and extracellular (β36) deletions essential for crystallisation of the turkey β-adrenoceptor. NAUNYN-SCHMIEDEBERG'S ARCHIVES OF PHARMACOLOGY 2011; 384:71-91. [PMID: 21547538 PMCID: PMC3116118 DOI: 10.1007/s00210-011-0648-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/15/2011] [Accepted: 04/18/2011] [Indexed: 11/30/2022]
Abstract
The X-ray crystal structure of the turkey β-adrenoceptor has recently been determined. However, mutations were introduced into the native receptor that was essential for structure determination. These may cause alterations to the receptor pharmacology. It is therefore essential to understand the effects of these mutations on the pharmacological characteristics of the receptor. This study examined the pharmacological effects of both the m23 mutations and the β36 deletions, both alone and then in combination in the β36-m23 mutant used in the crystallisation and structure determination of the turkey β-adrenoceptor. Stable CHO-K1 cell lines were made of each of the receptor mutants and the affinity and efficacy of ligands assessed by (3)H-CGP 12177 whole cell ligand binding, (3)H-cAMP accumulation, and CRE-SPAP gene transcription assays. The m23 mutations reduced affinity for agonists, partial agonists and neutral antagonists by about tenfold whilst the β36 deletions alone had no effect on ligand affinity. Both sets of changes appeared to reduce the agonist activation of the receptor. Both the m23 and the β36 receptors retained two active agonist-induced receptor conformations similar to that of the original tβtrunc receptor. The combined β36-m23 receptor bound ligands with similar affinity to the m23 receptor; however, agonist activation was only observed with a few agonists including the catecholamines. Although the combination of mutations severely reduced the activation ability, the final crystallised receptor (β36-m23) was still a fully functional receptor capable of binding agonist and antagonist ligands and activating intracellular agonist responses.
Collapse
Affiliation(s)
- Jillian G Baker
- Institute of Cell Signalling, C Floor Medical School, University of Nottingham, Queen's Medical Centre, UK.
| | | | | |
Collapse
|
19
|
Warne T, Moukhametzianov R, Baker JG, Nehmé R, Edwards PC, Leslie AG, Schertler GF, Tate CG. The structural basis for agonist and partial agonist action on a β(1)-adrenergic receptor. Nature 2011; 469:241-4. [PMID: 21228877 PMCID: PMC3023143 DOI: 10.1038/nature09746] [Citation(s) in RCA: 488] [Impact Index Per Article: 37.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2010] [Accepted: 12/09/2010] [Indexed: 12/02/2022]
Abstract
β-adrenergic receptors (βARs) are G-protein-coupled receptors (GPCRs) that activate intracellular G proteins upon binding catecholamine agonist ligands such as adrenaline and noradrenaline. Synthetic ligands have been developed that either activate or inhibit βARs for the treatment of asthma, hypertension or cardiac dysfunction. These ligands are classified as either full agonists, partial agonists or antagonists, depending on whether the cellular response is similar to that of the native ligand, reduced or inhibited, respectively. However, the structural basis for these different ligand efficacies is unknown. Here we present four crystal structures of the thermostabilized turkey (Meleagris gallopavo) β(1)-adrenergic receptor (β(1)AR-m23) bound to the full agonists carmoterol and isoprenaline and the partial agonists salbutamol and dobutamine. In each case, agonist binding induces a 1 Å contraction of the catecholamine-binding pocket relative to the antagonist bound receptor. Full agonists can form hydrogen bonds with two conserved serine residues in transmembrane helix 5 (Ser(5.42) and Ser(5.46)), but partial agonists only interact with Ser(5.42) (superscripts refer to Ballesteros-Weinstein numbering). The structures provide an understanding of the pharmacological differences between different ligand classes, illuminating how GPCRs function and providing a solid foundation for the structure-based design of novel ligands with predictable efficacies.
Collapse
Affiliation(s)
- Tony Warne
- MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 0QH, UK
| | | | | | - Rony Nehmé
- MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 0QH, UK
| | | | | | - Gebhard F.X. Schertler
- Joint corresponding authors: MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 0QH, UK, , , Telephone +44-(0)1223-402266, Fax +44-(0)1223-213556
| | | |
Collapse
|
20
|
Multiscale computational methods for mapping conformational ensembles of G-protein-coupled receptors. COMPUTATIONAL CHEMISTRY METHODS IN STRUCTURAL BIOLOGY 2011; 85:253-80. [DOI: 10.1016/b978-0-12-386485-7.00007-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|