An official website of the United States government
The long‐anticipated high‐resolution structures of the human melatonin G protein‐coupled receptors MT1and MT2, involved in establishing and maintaining circadian rhythm, were obtained in complex with two melatonin analogs and two approved anti‐insomnia and antidepression drugs using X‐ray free‐electron laser serial femtosecond crystallography. The structures shed light on the overall conformation and unusual structural features of melatonin receptors, as well as their ligand binding sites and the melatonergic pharmacophore, thereby providing insights into receptor subtype selectivity. The structures revealed an occluded orthosteric ligand binding site with a membrane‐buried channel for ligand entry in both receptors, and an additional putative ligand entry path in MT2from the extracellular side. This unexpected ligand entry mode contributes to facilitating the high specificity with which melatonin receptors bind their cognate ligand and exclude structurally similar molecules such as serotonin, the biosynthetic precursor of melatonin. Finally, the MT2structure allowed accurate mapping of type 2 diabetes‐related single‐nucleotide polymorphisms, where a clustering of residues in helices I and II on the protein–membrane interface was observed which could potentially influence receptor oligomerization. The role of receptor oligomerization is further discussed in light of the differential interaction of MT1and MT2with GPR50, a regulatory melatonin coreceptor. The melatonin receptor structures will facilitate design of selective tool compounds to further dissect the specific physiological function of each receptor subtype as well as provide a structural basis for next‐generation sleeping aids and other drugs targeting these receptors with higher specificity and fewer side effects.
more »
« less
Fine-tuning activation specificity of G-protein-coupled receptors via automated path searching
Physics-based simulation methods can grant atomistic insights into the molecular origin of the function of biomolecules. However, the potential of such approaches has been hindered by their low efficiency, including in the design of selective agonists where simulations of myriad protein–ligand combinations are necessary. Here, we describe an automated input-free path searching protocol that offers (within 14 d using Graphics Processing Unit servers) a minimum free energy path (MFEP) defined in high-dimension configurational space for activating sphingosine-1-phosphate receptors (S1PRs) by arbitrary ligands. The free energy distributions along the MFEP for four distinct ligands and three S1PRs reached a remarkable agreement with Bioluminescence Resonance Energy Transfer (BRET) measurements of G-protein dissociation. In particular, the revealed transition state structures pointed out toward two S1PR3 residues F263/I284, that dictate the preference of existing agonists CBP307 and BAF312 on S1PR1/5. Swapping these residues between S1PR1 and S1PR3 reversed their response to the two agonists in BRET assays. These results inspired us to design improved agonists with both strong polar head and bulky hydrophobic tail for higher selectivity on S1PR1. Through merely three in silico iterations, our tool predicted a unique compound scaffold. BRET assays confirmed that both chiral forms activate S1PR1 at nanomolar concentration, 1 to 2 orders of magnitude less than those for S1PR3/5. Collectively, these results signify the promise of our approach in fine agonist design for G-protein-coupled receptors.
more »
« less
- Award ID(s):
- 2142727
- PAR ID:
- 10586811
- Editor(s):
- Andricioaei, Ioan; Liu, Xiangyu
- Publisher / Repository:
- PNAS
- Date Published:
- Journal Name:
- PNAS nexus
- ISSN:
- 2752-6542
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Homomeric dimerization of metabotropic glutamate receptors (mGlus) is essential for the modulation of their functions and represents a promising avenue for the development of novel therapeutic approaches to address central nervous system diseases. Yet, the scarcity of detailed molecular and energetic data on mGlu2 impedes our in-depth comprehension of their activation process. Here, we employ computational simulation methods to elucidate the activation process and key events associated with the mGlu2, including a detailed analysis of its conformational transitions, the binding of agonists, Giprotein coupling, and the guanosine diphosphate (GDP) release. Our results demonstrate that the activation of mGlu2 is a stepwise process and several energy barriers need to be overcome. Moreover, we also identify the rate-determining step of the mGlu2’s transition from the agonist-bound state to its active state. From the perspective of free-energy analysis, we find that the conformational dynamics of mGlu2’s subunit follow coupled rather than discrete, independent actions. Asymmetric dimerization is critical for receptor activation. Our calculation results are consistent with the observation of cross-linking and fluorescent-labeled blot experiments, thus illustrating the reliability of our calculations. Besides, we also identify potential key residues in the Giprotein binding position on mGlu2, mGlu2 dimer’s TM6–TM6 interface, and Gi α5 helix by the change of energy barriers after mutation. The implications of our findings could lead to a more comprehensive grasp of class C G protein-coupled receptor activation.more » « less
-
Abstract G‐protein‐coupled receptors (GPCRs) are among the most important receptors in human physiology and pathology. They serve as master regulators of numerous key processes and are involved in as well as cause debilitating diseases. Consequently, GPCRs are among the most attractive targets for drug design and pharmaceutical interventions (>30% of drugs on the market). The glucagon‐like peptide 1 (GLP‐1) hormone receptor GLP1R is closely involved in insulin secretion by pancreatic β‐cells and constitutes a major druggable target for the development of anti‐diabetes and obesity agents. GLP1R structure was recently solved, with ligands, allosteric modulators and as part of a complex with its cognate G protein. However, the translation of this structural data into structure/function understanding remains limited. The current study functionally characterizes GLP1R with special emphasis on ligand and cellular partner binding interactions and presents a free‐energy landscape as well as a functional model of the activation cycle of GLP1R. Our results should facilitate a deeper understanding of the molecular mechanism underlying GLP1R activation, forming a basis for improved development of targeted therapeutics for diabetes and related disorders.more » « less
-
Abstract Angiogenesis, the formation of new vessels from existing vessels, is mediated by vascular endothelial growth factor (VEGF) and platelet-derived growth factor (PDGF). Despite discoveries supporting the cross-family interactions between VEGF and PDGF families, sharing the binding partners between them makes it challenging to identify growth factors that predominantly affect angiogenesis. Systems biology offers promises to untangle this complexity. Thus, in this study, we developed a mass-action kinetics-based computational model for cross-family interactions between VEGFs (VEGF-A, VEGF-B, and PlGF) and PDGFs (PDGF-AA, PDGF-AB, and PDGF-BB) with their receptors (VEGFR1, VEGFR2, NRP1, PDGFRα, and PDGFRβ). The model, parametrized with our literature mining and surface resonance plasmon assays, was validated by comparing the concentration of VEGFR1 complexes with a previously constructed angiogenesis model. The model predictions include five outcomes: 1) the percentage of free or bound ligands and 2) receptors, 3) the concentration of free ligands, 4) the percentage of ligands occupying each receptor, and 5) the concentration of ligands that is bound to each receptor. We found that at equimolar ligand concentrations (1 nM), PlGF and VEGF-A were the main binding partners of VEGFR1 and VEGFR2, respectively. Varying the density of receptors resulted in the following five outcomes: 1) Increasing VEGFR1 density depletes the free PlGF concentration, 2) increasing VEGFR2 density decreases PDGF:PDGFRα complexes, 3) increased NRP1 density generates a biphasic concentration of the free PlGF, 4) increased PDGFRα density increases PDGFs:PDGFRα binding, and 5) increasing PDGFRβ density increases VEGF-A:PDGFRβ. Our model offers a reproducible, fundamental framework for exploring cross-family interactions that can be extended to the tissue level or intracellular molecular level. Also, our model may help develop therapeutic strategies in pathological angiogenesis by identifying the dominant complex in the cell signaling. Author summaryNew blood vessel formation from existing ones is essential for growth, healing, and reproduction. However, when this process is disrupted—either too much or too little—it can contribute to diseases such as cancer and peripheral arterial disease. Two key families of proteins, vascular endothelial growth factors (VEGFs) and platelet-derived growth factors (PDGFs), regulate this process. Traditionally, scientists believed that VEGFs only bind to VEGF receptors and PDGFs to PDGF receptors. However, recent findings show that these proteins can interact with each other’s receptors, making it more challenging to understand and control blood vessel formation. To clarify these complex interactions, we combined computer modeling with biological data to map out which proteins bind to which receptors and to what extent. Our findings show that when VEGFs and PDGFs are present in equal amounts, VEGFs are the primary binding partners for VEGF receptors. We also explored how changes in receptor levels affect these interactions in disease-like conditions. This work provides a foundational computational model for studying cross-family interactions, which can be expanded to investigate tissue-level effects and processes inside cells. Ultimately, our model may help develop better treatments for diseases linked to abnormal blood vessel growth by identifying key protein-receptor interactions.more » « less
-
Understanding the activation mechanism of the μ-opioid receptor (μ-OR) and its selective coupling to the inhibitory G protein (Gi) is vital for pharmaceutical research aimed at finding treatments for the opioid overdose crisis. Many attempts have been made to understand the mechanism of the μ-OR activation, following the elucidation of new crystal structures such as the antagonist- and agonist-bound μ-OR. However, the focus has not been placed on the underlying energetics and specificity of the activation process. An energy-based picture would not only help to explain this coupling but also help to explore why other possible options are not common. For example, one would like to understand why μ-OR is more selective to Githan a stimulatory G protein (Gs). Our study used homology modeling and a coarse-grained model to generate all of the possible “end states” of the thermodynamic cycle of the activation of μ-OR. The end points were further used to generate reasonable intermediate structures of the receptor and the Gito calculate two-dimensional free energy landscapes. The results of the landscape calculations helped to propose a plausible sequence of conformational changes in the μ-OR and Gisystem and for exploring the path that leads to its activation. Furthermore, in silico alanine scanning calculations of the last 21 residues of the C terminals of Giand Gswere performed to shed light on the selective binding of Gito μ-OR. Overall, the present work appears to demonstrate the potential of multiscale modeling in exploring the action of G protein-coupled receptors.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
