Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β₂-adrenergic receptor (β₂AR) and the cholecystokinin (CCK) receptor subfamily. The β₂AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β₂AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK₁R shows functional sensitivity to membrane cholesterol while CCK₂R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β₂AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-T_m) shows that mutation of a conserved CRAC sequence on TM7 of the β₂AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.     

Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties

From computational simulations of a serotonin 2A receptor (5-HT(2A)R) 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-HT(2A)R, 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.     

A mechanistic role of Helix 8 in GPCRs: Computational modeling of the dopamine D2 receptor interaction with the GIPC1–PDZ-domain

Helix-8 (Hx8) is a structurally conserved amphipathic helical motif in class-A GPCRs, adjacent to the C-terminal sequence that is responsible for PDZ-domain-recognition. The Hx8 segment in the dopamine D2 receptor (D2R) constitutes the C-terminal segment and we investigate its role in the function of D2R by studying the interaction with the PDZ-containing GIPC1 using homology models based on the X-ray structures of very closely related analogs: the D3R for the D2R model, and the PDZ domain of GIPC2 for GIPC1–PDZ. The mechanism of this interaction was investigated with all-atom unbiased molecular dynamics (MD) simulations that reveal the role of the membrane in maintaining the helical fold of Hx8, and with biased MD simulations to elucidate the energy drive for the interaction with the GIPC1–PDZ. We found that it becomes more favorable energetically for Hx8 to adopt the extended conformation observed in all PDZ–ligand complexes when it moves away from the membrane, and that C-terminus palmitoylation of D2R enhanced membrane penetration by the Hx8 backbone. De-palmitoylation enables Hx8 to move out into the aqueous environment for interaction with the PDZ domain. All-atom unbiased MD simulations of the full D2R–GIPC1-PDZ complex in sphingolipid/cholesterol membranes show that the D2R carboxyl C-terminus samples the region of the conserved GFGL motif located on the carboxylate-binding loop of the GIPC1–PDZ, and the entire complex distances itself from the membrane interface. Together, these results outline a likely mechanism of Hx8 involvement in the interaction of the GPCR with PDZ-domains in the course of signaling.     

Cholesterol Modulates the Dimer Interface of the β2-Adrenergic Receptor via Cholesterol Occupancy Sites

The β₂-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β₂-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β₂-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.     

Activation mechanisms of the first sphingosine‐1‐phosphate receptor

Activation of the first sphingosine‐1‐phosphate receptor (S1PR₁) promotes permeability of the blood brain barrier, astrocyte and neuronal protection, and lymphocyte egress from secondary lymphoid tissues. Although an agonist often activates the S1PR₁, the receptor exhibits high levels of basal activity. In this study, we performed long‐timescale molecular dynamics and accelerated molecular dynamics (aMD) simulations to investigate activation mechanisms of the ligand‐free (apo) S1PR₁. In the aMD enhanced sampling simulations, we observed four independent events of activation, which is characterized by close interaction between Y311^7.53 and Y221^5.58 and increased distance between the intracellular ends of transmembrane (TM) helices 3 and 6. Although TM helices TM3, TM6, TM5 and, TM7 are associated with GPCR activation, we discovered that their movements are not necessarily correlated during activation. Instead, TM5 showed a decreased correlation with each of these regions during activation. During activation of the apo receptor, Y221^5.58 and Y311^7.53 became more solvated, because a water channel formed in the intracellular pocket. Additionally, a lipid molecule repeatedly entered the receptor between the extracellular ends of TM1 and TM7, providing important insights into the pathway of ligand entry into the S1PR₁.     

Differential dynamics of the serotonin1A receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

Abstract

The serotonin_1A receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin_1A receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3?μs, to analyze the effect of cholesterol on the structure and dynamics of the serotonin_1A receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin_1A receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.     

Cholesterol Modulates the Dimer Interface of the β2-Adrenergic Receptor via Cholesterol Occupancy Sites

The β₂-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β₂-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β₂-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.     

The Role of Water in Activation Mechanism of Human N-Formyl Peptide Receptor 1 (FPR1) Based on Molecular Dynamics Simulations

The Formyl Peptide Receptor 1 (FPR1) is an important chemotaxis receptor involved in various aspects of host defense and inflammatory processes. We constructed a model of FPR1 using as a novel template the chemokine receptor CXCR4 from the same branch of the phylogenetic tree of G-protein-coupled receptors. The previously employed template of rhodopsin contained a bulge at the extracellular part of TM2 which directly influenced binding of ligands. We also conducted molecular dynamics (MD) simulations of FPR1 in the apo form as well as in a form complexed with the agonist fMLF and the antagonist tBocMLF in the model membrane. During all MD simulation of the fMLF-FPR1 complex a water molecule transiently bridged the hydrogen bond between W254^6.48 and N108^3.35 in the middle of the receptor. We also observed a change in the cytoplasmic part of FPR1 of a rotamer of the Y301^7.53 residue (tyrosine rotamer switch). This effect facilitated movement of more water molecules toward the receptor center. Such rotamer of Y301^7.53 was not observed in any crystal structures of GPCRs which can suggest that this state is temporarily formed to pass the water molecules during the activation process. The presence of a distance between agonist and residues R201^5.38 and R205^5.42 on helix TM5 may suggest that the activation of FPR1 is similar to the activation of β-adrenergic receptors since their agonists are separated from serine residues on helix TM5. The removal of water molecules bridging these interactions in FPR1 can result in shrinking of the binding site during activation similarly to the shrinking observed in β-ARs. The number of GPCR crystal structures with agonists is still scarce so the designing of new ligands with agonistic properties is hampered, therefore homology modeling and docking can provide suitable models. Additionally, the MD simulations can be beneficial to outline the mechanisms of receptor activation and the agonist/antagonist sensing.     

Solid-state NMR and molecular dynamics simulations reveal the oligomeric ion-channels of TM2-GABAA stabilized by intermolecular hydrogen bonding

The second transmembrane (TM2) domain of GABA_A receptor forms the inner-lining surface of chloride ion-channel and plays important roles in the function of the receptor protein. In this study, we report the first structure of TM2 in lipid bilayers determined using solid-state NMR and MD simulations. The interatomic ¹³C-¹⁵N distances measured from REDOR magic angle spinning experiments on multilamellar vesicles, containing a TM2 peptide site specifically labeled with ¹³C′ and ¹⁵N isotopes, were used to determine the secondary structure of the peptide. The ¹⁵N chemical shift and ¹H-¹⁵N dipolar coupling parameters measured from PISEMA experiments on mechanically aligned phospholipid bilayers, containing a TM2 peptide site specifically labeled with ¹⁵N isotopes, under static conditions were used to determine the membrane orientation of the peptide. Our results reveal that the TM2 peptide forms an alpha helical conformation with a tilted transmembrane orientation, which is unstable as a monomer but stable as pentameric oligomers as indicated by MD simulations. Even though the peptide consists of a number of hydrophilic residues, the transmembrane folding of the peptide is stabilized by intermolecular hydrogen bondings between the side chains of Ser and Thr residues as revealed by MD simulations. The results also suggest that peptide-peptide interactions in the tilted transmembrane orientation overcome the hydrophobic mismatch between the peptide and bilayer thickness.     

Elastic network normal mode dynamics reveal the GPCR activation mechanism

G‐protein‐coupled receptors (GPCR) are a family of membrane‐embedded metabotropic receptors which translate extracellular ligand binding into an intracellular response. Here, we calculate the motion of several GPCR family members such as the M2 and M3 muscarinic acetylcholine receptors, the A_2A adenosine receptor, the β₂‐adrenergic receptor, and the CXCR4 chemokine receptor using elastic network normal modes. The normal modes reveal a dilation and a contraction of the GPCR vestibule associated with ligand passage, and activation, respectively. Contraction of the vestibule on the extracellular side is correlated with cavity formation of the G‐protein binding pocket on the intracellular side, which initiates intracellular signaling. Interestingly, the normal modes of rhodopsin do not correlate well with the motion of other GPCR family members. Electrostatic potential calculation of the GPCRs reveal a negatively charged field around the ligand binding site acting as a siphon to draw‐in positively charged ligands on the membrane surface. Altogether, these results expose the GPCR activation mechanism and show how conformational changes on the cell surface side of the receptor are allosterically translated into structural changes on the inside. Proteins 2014; 82:579–586. © 2013 Wiley Periodicals, Inc.     

Cholesterol‐β1AR interaction versus cholesterol‐β2AR interaction

Two 8‐µs all‐atom molecular dynamics simulations have been performed on the two highly homologous G protein‐coupled receptor (GPCR) subtypes, β₁‐ and β₂‐adrenergic receptors, which were embedded in a lipid bilayer with randomly dispersed cholesterol molecules. During the simulations, cholesterol molecules accumulate to different surface regions of the two receptors, suggesting the subtype specificity of cholesterol–β‐adrenergic receptor interaction and providing some clues to the physiological difference of the two subtypes. Meanwhile, comparison between the two receptors in interacting with cholesterols shed some new light on general determinants of cholesterol binding to GPCRs. Our results indicate that although the concave surface, charged residues and aromatic residues are important, neither of these stabilizing factors is indispensable for a cholesterol interaction site. Different combinations of these factors lead to the diversified binding modes of cholesterol binding to the receptors. Our long‐time simulations, for the first time, revealed the pathway of a cholesterol molecule entering the consensus cholesterol motif (CCM) site, and the binding process of cholesterol to CCM is accompanied by a side chain flipping of the conserved Trp4.50. Moreover, the simulation results suggest that the I‐/V‐/L‐rich region on the extracellular parts of helix 6 might be an alternatively conserved cholesterol‐binding site for the class‐A GPCRs. Proteins 2014; 82:760–770. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.     

Structural dynamics of Smoothened (SMO) in the ciliary membrane and its interaction with membrane lipids

The Smoothened receptor (SMO, a 7 pass transmembrane domain, Class F GPCR family protein) plays a crucial role in the Hedgehog (HH) signaling pathway, which is involved in embryonic development and is implicated in various types of cancer throughout the animal kingdom. In the absence of HH signaling, SMO is inhibited by Patched 1 (PTC1; a 12 pass transmembrane domain protein), which is localized in the primary cilia. HH binding leads to the dislocation of PTC1 from the cilia, thus making way for SMO to localize in the primary cilia, as an essential prerequisite for its activation. We have carried out MARTINI coarse-grained molecular dynamics simulations of SMO in POPC and in ciliary membrane models, respectively, to study the interactions of SMO with cholesterol and other lipid molecules in the ciliary membrane, and to gain molecular-level insights into the role of the primary cilia in shaping the functional dynamics of SMO. We are able to identify the interaction of membrane cholesterols with definite sites and domains within SMO and relate them with known cholesterol-binding sequence and structure motifs. We show that cholesterol interactions with the transmembrane domain TMD, unlike those with the cysteine-rich domain (CRD) and the intracellular domain (ICD), are through residues belonging to known cholesterol-binding motifs. Notably, a few persistent interactions of cholesterol with lower TM cholesterol-binding domains are governed by the presence of multiple cholesterol-binding motifs. These analyses have also helped to identify and define a strict cholesterol consensus motif (CCM), which may well steer cholesterol into the hitherto identified binding sites within the TMD of SMO. We have also reported the interaction of phosphatidylinositol 4-phosphate with the intracellular region of transmembrane (TM) helices (TM1, TM3, TM4, and TM5), intracellular loop1, helix8, and Arg/Lys clusters of the ICD. Structural analysis of SMO domains shows significant changes in the CRD and ICD, during the course of the simulation. Further detailed analysis of the dynamics of the TMD reveals the movements of TM5, TM6, and TM7, linked with the helix8, which are possibly involved in shaping the conformational disposition of the ICD. The movement of these TM helices could possibly be a consequence of interactions involving the extracellular domain and extracellular loops. In addition, our analysis also shows that phosphatidylinositol-4-phosphate (PI4P), along with some ICD cholesterols, are implicated in anchoring SMO in the membrane.     

The conformation of the cytoplasmic helix 8 of the CB1 cannabinoid receptor using NMR and circular dichroism

The cytoplasmic helix domain (fourth cytoplasmic loop, helix 8) of numerous GPCRs such as rhodopsin and the β-adrenergic receptor exhibits unique structural and functional characteristics. Computational models also predict the existence of such a structural motif within the CB1 cannabinoid receptor, another member of the G-protein coupled receptor superfamily. To gain insights into the conformational properties of this GPCR component, a peptide corresponding to helix 8 of the CB1 receptor with a small contiguous segment from transmembrane helix 7 (TM7) was chemically synthesized and its secondary structure determined by circular dichroism (CD) and solution NMR spectroscopy. Our studies in DPC and SDS micelles revealed significant α-helical structure while in an aqueous medium, the peptide exhibited a random coil configuration. The relative orientation of helix 8 within the CB1 receptor was obtained from intermolecular ³¹P-¹H and ¹H-¹H NOE measurements. Our results suggest that in the presence of an amphipathic membrane environment, helix 8 assumes an alpha helical structure with an orientation parallel to the phospholipid membrane surface and perpendicular to TM7. In this model, positively charged side chains interact with the lipid headgroups while the other polar side chains face the aqueous region. The above observations may be relevant to the activation/deactivation of the CB1 receptor.     

A model for how Gβγ couples Gα to GPCR

Representing ∼5% of the human genome, G-protein-coupled receptors (GPCRs) are a primary target for drug discovery; however, the molecular details of how they couple to heterotrimeric G protein subunits are incompletely understood. Here, I propose a hypothetical initial docking model for the encounter between GPCR and Gβγ that is defined by transient interactions between the cytosolic surface of the GPCR and the prenyl moiety and the tripeptide motif, asparagine–proline–phenylalanine (NPF), in the C-terminus of the Gγ subunit. Analysis of class A GPCRs reveals a conserved NPF binding site formed by the interaction of the TM1 and H8. Functional studies using differentially prenylated proteins and peptides further suggest that the intracellular hydrophobic core of the GPCR is a prenyl binding site. Upon binding TM1 and H8 of GPCRs, the propensity of the C-terminal region of Gγ to convert into an α helix allows it to extend into the hydrophobic core of the GPCR, facilitating the GPCR active state. Conservation of the NPF motif in Gγ isoforms and interacting residues in TM1 and H8 suggest that this is a general mechanism of GPCR–G protein signaling. Analysis of the rhodopsin dimer also suggests that Gγ–rhodopsin interactions may facilitate GPCR dimer transactivation.     

Role of rhodopsin N-terminus in structure and function of rhodopsin-bitter taste receptor chimeras

The bitter taste receptors (T2Rs) belong to the G protein-coupled receptor (GPCR) superfamily. In humans, bitter taste sensation is mediated by 25 T2Rs. Structure–function studies on T2Rs are impeded by the low-level expression of these receptors. Different lengths of rhodopsin N-terminal sequence inserted at the N-terminal region of T2Rs are commonly used to express these receptors in heterologous systems. While the additional sequences were reported, to enhance the expression of the T2Rs, the local structural perturbations caused by these sequences and its effect on receptor function or allosteric ligand binding were not characterized. In this study, we elucidated how different lengths of rhodopsin N-terminal sequence effect the structure and function of the bitter taste receptor, T2R4. Guided by molecular models of T2R4 built using a rhodopsin crystal structure as template, we constructed chimeric T2R4 receptors containing the rhodopsin N-terminal 33 and 38 amino acids. The chimeras were functionally characterized using calcium imaging, and receptor expression was determined by flow cytometry. Our results show that rhodopsin N-terminal 33 amino acids enhance expression of T2R4 by 2.5-fold and do not cause perturbations in the receptor structure.     

FoldGPCR: Structure prediction protocol for the transmembrane domain of G protein‐coupled receptors from class A

Building reliable structural models of G protein‐coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template‐based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1–3 Å Cα RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter‐residue contacts between the template and target, and exploits an all‐atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor‐ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human β₂‐adrenergic receptor (β₂AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X‐ray crystal structure of the β₂AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix‐helix packing interactions that distinguish rhodopsin from β₂AR. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Internal hydration increases during activation of the G-protein-coupled receptor rhodopsin

Rhodopsin, the membrane protein responsible for dim-light vision, until recently was the only G-protein-coupled receptor (GPCR) with a known crystal structure. As a result, there is enormous interest in studying its structure, dynamics, and function. Here we report the results of three all-atom molecular dynamics simulations, each at least 1.5 μs, which predict that substantial changes in internal hydration play a functional role in rhodopsin activation. We confirm with ¹H magic angle spinning NMR that the increased hydration is specific to the metarhodopsin-I intermediate. The internal water molecules interact with several conserved residues, suggesting that changes in internal hydration may be important during the activation of other GPCRs. The results serve to illustrate the synergism of long-time-scale molecular dynamics simulations and NMR in enhancing our understanding of GPCR function.