Activation of corticotropin-releasing factor 1 receptor: insights from molecular dynamics simulations

Biased Signaling in Mutated Variants of β2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations

The molecular basis of receptor bias in G protein-coupled receptors (GPCRs) caused by mutations that preferentially activate specific intracellular transducers over others remains poorly understood. Two experimentally identified biased variants of β2-adrenergic receptors (β2AR), a prototypical GPCR, are a triple mutant (T68F, Y132A, and Y219A) and a single mutant (Y219A); the former bias the receptor toward the β-arrestin pathway by disfavoring G protein engagement, while the latter induces G protein signaling explicitly due to selection against GPCR kinases (GRKs) that phosphorylate the receptor as a prerequisite of β-arrestin binding. Though rigorous characterizations have revealed functional implications of these mutations, the atomistic origin of the observed transducer selectivity is not clear. In this study, we investigated the allosteric mechanism of receptor bias in β2AR using microseconds of all-atom Gaussian accelerated molecular dynamics (GaMD) simulations. Our observations reveal distinct rearrangements in transmembrane helices, intracellular loop 3, and critical residues R1313.50 and Y3267.53 in the conserved motifs D(E)RY and NPxxY for the mutant receptors, leading to their specific transducer interactions. Moreover, partial dissociation of G protein from the receptor core is observed in the simulations of the triple mutant in contrast to the single mutant and wild-type receptor. The reorganization of allosteric communications from the extracellular agonist BI-167107 to the intracellular receptor-transducer interfaces drives the conformational rearrangements responsible for receptor bias in the single and triple mutants. The molecular insights into receptor bias of β2AR presented here could improve the understanding of biased signaling in GPCRs, potentially opening new avenues for designing novel therapeutics with fewer side-effects and superior efficacy.

Structural Validation by the G-Factor Properly Regulates Boost Potentials Imposed in Conformational Sampling of Proteins

Free energy landscapes (FELs) of proteins are indispensable for evaluating thermodynamic properties. Molecular dynamics (MD) simulation is a computational method for calculating FELs; however, conventional MD simulation frequently fails to search a broad conformational subspace due to its accessible timescale, which results in the calculation of an unreliable FEL. To search a broad subspace, an external bias can be imposed on a protein system, and biased sampling tends to cause a strong perturbation that might collapse the protein structures, indicating that the strength of the external bias should be properly regulated. This regulation can be challenging, and empirical parameters are frequently employed to impose an optimal bias. To address this issue, several methods regulate the external bias by referring to system energies. Herein, we focused on protein structural information for this regulation. In this study, a well-established structural indicator (the G-factor) was used to obtain structural information. Based on the G-factor, we proposed a scheme for regulating biased sampling, which is referred to as a G-factor-based external bias limiter (GERBIL). With GERBIL, the configurations were structurally validated by the G-factor during biased sampling. As an example of biased sampling, an accelerated MD (aMD) simulation was adopted in GERBIL (aMD-GERBIL), whereby the aMD simulation was repeatedly performed by increasing the strength of the boost potential. Furthermore, the configurations sampled by the aMD simulation were structurally validated by their G-factor values, and aMD-GERBIL stopped increasing the strength of the boost potential when the sampled configurations were regarded as low-quality (collapsed) structures. This structural validation is regarded as a "Brake" of the boost potential. For demonstrations, aMD-GERBIL was applied to globular proteins (ribose binding and maltose-binding proteins) to promote their large-amplitude open-closed transitions and successfully identify their domain motions.

Molecular Insights into Phosphorylation-Induced Allosteric Conformational Changes in a β2-Adrenergic Receptor

The large conformational flexibility of G protein-coupled receptors (GPCRs) has been a puzzle in structural and pharmacological studies for the past few decades. Apart from structural rearrangements induced by ligands, enzymatic phosphorylations by GPCR kinases (GRKs) at the carboxy-terminal tail (C-tail) of a GPCR also make conformational alterations to the transmembrane helices and facilitates the binding of one of its transducer proteins named β-arrestin. The phosphorylation-induced conformational transition of the receptor that causes specific binding to β-arrestin but prevents the association of other transducers such as G proteins lacks atomistic understanding and is elusive to experimental studies. Using microseconds of all-atom conventional and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigate the allosteric mechanism of phosphorylation induced-conformational changes in β₂-adrenergic receptor, a well-characterized GPCR model system. Free energy profiles reveal that the phosphorylated receptor samples a new conformational state in addition to the canonical active state corroborating with recent nuclear magnetic resonance experimental findings. The new state has a smaller intracellular cavity that is likely to accommodate β-arrestin better than G protein. Using contact map and inter-residue interaction energy calculations, we found the phosphorylated C-tail adheres to the cytosolic surface of the transmembrane domain of the receptor. Transfer entropy calculations show that the C-tail residues drive the correlated motions of TM residues, and the allosteric signal is relayed via several residues at the cytosolic surface. Our results also illustrate how the redistribution of inter-residue nonbonding interaction couples with the allosteric communication from the phosphorylated C-tail to the transmembrane. Atomistic insight into phosphorylation-induced β-arrestin specific conformation is therapeutically important to design drugs with higher efficacy and fewer side effects. Our results, therefore, open novel opportunities to fine-tune β-arrestin bias in GPCR signaling.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

Understanding the Mechanism of Activation/Deactivation of GLP-1R via Accelerated Molecular Dynamics Simulation

Glucagon-like peptide-1 receptor (GLP-1R), as a member of the class B G protein-coupled receptors (GPCRs), plays a crucial role in regulating blood glucose level signal recognition through its activation. The conformation changes during the activation pathway are of particular importance for its function. To investigate the activation mechanism of GLP-1R, the crystal structures of active and inactive forms are chosen to perform a total of 2 μs of accelerated molecular dynamics (aMD) simulations and 400ns of conventional molecular dynamics (cMD) simulations. With the aid of structural analysis and potential of mean force (PMF) calculations, we reveal the role of different helices in the activation and deactivation process and obtain the intermediate states during activation and deactivation that are difficult to capture in experiments. Protein structure network (PSN) was utilised to clarify the allosteric communication pathways of activation and deactivation and reveal the mechanisms of its activation and deactivation. The results could advance our understanding of the activation mechanism of GLP-1R and the related drug design.

Functional effect of nobiletin as a food-derived allosteric modulator of mouse CRFR2β in skeletal muscle

Activation of corticotropin-releasing factor receptor 2β (CRFR2β) results in increased skeletal muscle mass and the prevention of muscle atrophy. Using a luciferase reporter assay, we screened 357 functional food factors that activate CRFR2β and, subsequently, confirmed that nobiletin (NBT) increases CRFR2β activity. Additionally, we found that NBT augments the activity of the endogenous peptide ligand urocortin 2 (Ucn2) in a concentration-dependent manner. Computational simulation of CRFR2β confirmed that transmembrane domains (TMs) 1 and 2 are important for the synergistic activity of NBT and also identified important amino acids in these domains. Finally, we demonstrated that a co-administration of Ucn2 and NBT increases the hypertrophic signal in mouse skeletal muscle. These observations demonstrate that NBT can activate CRFR2β and amplify the agonistic activity of Ucn2 and that such food-derived molecules have the potential to enhance endogenous G protein-coupled receptor ligand activities and contribute to the maintenance of skeletal muscle mass and function.

Proton Sensitivity of Corticotropin-Releasing Hormone Receptor 1 Signaling to Proopiomelanocortin in Male Mice

Because an acidic cellular microenvironment is engendered by inflammation and may determine cell differentiation, we elucidated the impact of acidic conditions on induction of proopiomelanocortin (POMC) expression. Here, we demonstrate mechanisms for proton sensitivity of CRH receptor 1 (CRHR1) signaling to POMC and ACTH production. Low pH (6.8) resulted in doubling of POMC expression and ACTH production in pituitary cell line AtT-20 and in primary mouse pituitary cells. Using CRISPR knockout, we show that CRHR1 is necessary for acid-induced POMC expression, and this induction is mediated by CRHR1 histidine residues and calmodulin-dependent protein kinase II in both pituitary corticotroph cells and in nonpituitary cell lines expressing ectopic ACTH. In contrast, CRH ligand binding affinity to CRHR1 was decreased with acidic pH, implying that proton-induced POMC expression prevails in acidic conditions independently of CRH ligand binding. The results indicate that proton-induced CRHR1 signaling regulates ACTH production in response to an acidic microenvironment.

The Structure of the CGRP and Related Receptors

The canonical CGRP receptor is a complex between calcitonin receptor-like receptor (CLR), a family B G-protein-coupled receptor (GPCR) and receptor activity-modifying protein 1 (RAMP1). A third protein, receptor component protein (RCP) is needed for coupling to Gs. CGRP can interact with other RAMP-receptor complexes, particularly the AMY1 receptor formed between the calcitonin receptor (CTR) and RAMP1. Crystal structures are available for the binding of CGRP_27-37 [D³¹,P³⁴,F³⁵] to the extracellular domain (ECD) of CLR and RAMP1; these show that extreme C-terminal amide of CGRP interacts with W84 of RAMP1 but the rest of the analogue interacts with CLR. Comparison with the crystal structure of a fragment of the allied peptide adrenomedullin bound to the ECD of CLR/RAMP2 confirms the importance of the interaction of the ligand C-terminus and the RAMP in determining pharmacology specificity, although the RAMPs probably also have allosteric actions. A cryo-electron microscope structure of calcitonin bound to the full-length CTR associated with Gs gives important clues as to the structure of the complete receptor and suggests that the N-terminus of CGRP makes contact with His^5.40b, high on TM5 of CLR. However, it is currently not known how the RAMPs interact with the TM bundle of any GPCR. Major challenges remain in understanding how the ECD and TM domains work together to determine ligand specificity, and how G-proteins influence this and the role of RCP. It seems likely that allosteric mechanisms are particularly important as are the dynamics of the receptors.

Assessment of Conformational State Transitions of Class B GPCRs Using Molecular Dynamics

Class B G protein-coupled receptors (GPCRs) comprise a family of 15 peptide-binding members, which are crucial targets for endocrine, metabolic, and stress-related disorders. While their protein structures and dynamics remain largely unclear, computer modeling and simulations represent a promising means to help solve such puzzles. Herein, we present a basic introduction to the methodology of molecular dynamics (MD) simulations and two analytical methods to assess the conformational ensembles and transitions of Class B GPCRs, using our recent studies of the human pituitary adenylate cyclase activating polypeptide (PAC₁) receptor as an example. From long MD simulations, conformational ensembles with different roles in ligand binding and receptor activation are sampled to establish four states identified as either "open" or "closed" for the PAC₁ receptor. Next, the dynamical network can be applied to analyze the simulations and identify key features within each conformational ensemble, which help distinguish the ligand-bound states of the PAC₁ receptor from the ligand-free one. Further, the Markov State Model has emerged as a key approach to construct the transition network and connect the GPCR ensembles, providing detailed information for the transition pathways and kinetics. For the ligand-free PAC₁ receptor, the transitions within the closed states are near 10-30 times faster than the open-closed transitions, which is likely related to the activation mechanism of the receptor. Overall, long MD simulations and analyses are useful to assess conformational transitions for the Class B GPCRs and to gain mechanistic insight, which is difficult to obtain using other methods.

Receptor activity-modifying protein dependent and independent activation mechanisms in the coupling of calcitonin gene-related peptide and adrenomedullin receptors to Gs

Calcitonin gene-related peptide (CGRP) or adrenomedullin (AM) receptors are heteromers of the calcitonin receptor-like receptor (CLR), a class B G protein-coupled receptor, and one of three receptor activity-modifying proteins (RAMPs). How CGRP and AM activate CLR and how this process is modulated by RAMPs is unclear. We have defined how CGRP and AM induce Gs-coupling in CLR-RAMP heteromers by measuring the effect of targeted mutagenesis in the CLR transmembrane domain on cAMP production, modeling the active state conformations of CGRP and AM receptors in complex with the Gs C-terminus and conducting molecular dynamics simulations in an explicitly hydrated lipidic bilayer. The largest effects on receptor signaling were seen with H295A^5.40b, I298A^5.43b, L302A^5.47b, N305A^5.50b, L345A^6.49b and E348A^6.52b, F349A^6.53b and H374A^7.47b (class B numbering in superscript). Many of these residues are likely to form part of a group in close proximity to the peptide binding site and link to a network of hydrophilic and hydrophobic residues, which undergo rearrangements to facilitate Gs binding. Residues closer to the extracellular loops displayed more pronounced RAMP or ligand-dependent effects. Mutation of H374^7.47b to alanine increased AM potency 100-fold in the CGRP receptor. The molecular dynamics simulation showed that TM5 and TM6 pivoted around TM3. The data suggest that hydrophobic interactions are more important for CLR activation than other class B GPCRs, providing new insights into the mechanisms of activation of this class of receptor. Furthermore the data may aid in the understanding of how RAMPs modulate the signaling of other class B GPCRs.

Molecular Mechanism of Nucleotide-Dependent Allosteric Regulation in AMP-Activated Protein Kinase

The AMP-activated protein kinase (AMPK), a central enzyme in the regulation of energy homeostasis, is an important drug target for type 2 diabetes, obesity, and cancer. Binding of adenosine nucleotides to the regulatory γ-subunit tightly regulates the activity of this enzyme. Though recent crystal structures of AMPK have provided important insights into the allosteric activation of AMPK, molecular details of the regulatory mechanism of AMPK activation is still elusive. Here, we have performed extensive all-atom molecular dynamics (MD) simulations and shown that the kinase domain (KD) and γ-subunit come closer resulting in a more compact heterotrimeric AMPK complex in AMP-bound state compared to the ATP-bound state. The binding of ATP at site 3 of regulatory γ-subunit allosterically inhibits AMPK by destabilizing different regulatory regions of α-subunit: the autoinhibitory domain, the linker region, and the activation loop of the kinase core. The catalytically important residues experience a change in mechanical stress, and major rearrangements in community structure derived from residue-residue interaction energy-based network are observed in KD and α-linker region upon binding of different nucleotides. Our results also highlight the role of conserved charged residues forming an ionic network near the site 3 of γ-subunit in allosteric communications.

Receptor Activity-modifying Protein-directed G Protein Signaling Specificity for the Calcitonin Gene-related Peptide Family of Receptors

The calcitonin gene-related peptide (CGRP) family of G protein-coupled receptors (GPCRs) is formed through the association of the calcitonin receptor-like receptor (CLR) and one of three receptor activity-modifying proteins (RAMPs). Binding of one of the three peptide ligands, CGRP, adrenomedullin (AM), and intermedin/adrenomedullin 2 (AM2), is well known to result in a Gα_s-mediated increase in cAMP. Here we used modified yeast strains that couple receptor activation to cell growth, via chimeric yeast/Gα subunits, and HEK-293 cells to characterize the effect of different RAMP and ligand combinations on this pathway. We not only demonstrate functional couplings to both Gα_s and Gα_q but also identify a Gα_i component to CLR signaling in both yeast and HEK-293 cells, which is absent in HEK-293S cells. We show that the CGRP family of receptors displays both ligand- and RAMP-dependent signaling bias among the Gα_s, Gα_i, and Gα_q/11 pathways. The results are discussed in the context of RAMP interactions probed through molecular modeling and molecular dynamics simulations of the RAMP-GPCR-G protein complexes. This study further highlights the importance of RAMPs to CLR pharmacology and to bias in general, as well as identifying the importance of choosing an appropriate model system for the study of GPCR pharmacology.

In silico Exploration of the Conformational Universe of GPCRs

The structural plasticity of G protein coupled receptors (GPCRs) leads to a conformational universe going from inactive to active receptor states with several intermediate states. Many of them have not been captured yet and their role for GPCR activation is not well understood. The study of this conformational space and the transition dynamics between different receptor populations is a major challenge in molecular biophysics. The rational design of effector molecules that target such receptor populations allows fine-tuning receptor signalling with higher specificity to produce drugs with safer therapeutic profiles. In this minireview, we outline highly conserved receptor regions which are considered determinant for the establishment of distinct receptor states. We then discuss in-silico approaches such as dimensionality reduction methods and Markov State Models to explore the GPCR conformational universe and exploit the obtained conformations through structure-based drug design.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Transmembrane signal transduction by peptide hormones via family B G protein-coupled receptors

Although family B G protein-coupled receptors (GPCRs) contain only 15 members, they play key roles in transmembrane signal transduction of hormones. Family B GPCRs are drug targets for developing therapeutics for diseases ranging from metabolic to neurological disorders. Despite their importance, the molecular mechanism of activation of family B GPCRs remains largely unexplored due to the challenges in expression and purification of functional receptors to the quantity for biophysical characterization. Currently, there is no crystal structure available of a full-length family B GPCR. However, structures of key domains, including the extracellular ligand binding regions and seven-helical transmembrane regions, have been solved by X-ray crystallography and NMR, providing insights into the mechanisms of ligand recognition and selectivity, and helical arrangements within the cell membrane. Moreover, biophysical and biochemical methods have been used to explore functions, key residues for signaling, and the kinetics and dynamics of signaling processes. This review summarizes the current knowledge of the signal transduction mechanism of family B GPCRs at the molecular level and comments on the challenges and outlook for mechanistic studies of family B GPCRs.

Analysis of the Free Energy Landscapes for the Opening-Closing Dynamics of the Maltose Transporter ATPase MalK2 Using Enhanced-Sampling Molecular Dynamics Simulation

Protein dynamics are considered significant for many physiological processes, such as metabolism, biomolecular recognition, and the regulation of several vital cellular processes. Due to their flexibility, proteins may stay in different substates with or without the existence of the cognate substrates. To describe these phenomena, two models have been proposed: the "induced fit" and the "conformational selection" mechanisms. In this study, we used MalK2, the subunits that mainly include the nucleotide-binding domains (NBDs) of the maltose transporter from Escherichia coli, as a target to understand the NBD dimerization mechanism. Accelerated and conventional molecular dynamics have been performed. The results revealed that Mg-ATP binding to MalK2 led to a significant change in the free energy profile and thus stabilized the closed conformation. On the contrary, when Mg-ATP was removed, the open conformation would be favored. The fact that ligand binding induces a drastic free energy change leads to a significant inference: MalK2 dimerization would occur through the induced-fit mechanism rather than the conformational selection mechanism. This study sheds new light on the NBD dimerization mechanism and would be of wide applicability to other ABC transporters.

Linear Response Path Following: A Molecular Dynamics Method To Simulate Global Conformational Changes of Protein upon Ligand Binding

Molecular functions of proteins are often fulfilled by global conformational changes that couple with local events such as the binding of ligand molecules. High molecular complexity of proteins has, however, been an obstacle to obtain an atomistic view of the global conformational transitions, imposing a limitation on the mechanistic understanding of the functional processes. In this study, we developed a new method of molecular dynamics (MD) simulation called the linear response path following (LRPF) to simulate a protein's global conformational changes upon ligand binding. The method introduces a biasing force based on a linear response theory, which determines a local reaction coordinate in the configuration space that represents linear coupling between local events of ligand binding and global conformational changes and thus provides one with fully atomistic models undergoing large conformational changes without knowledge of a target structure. The overall transition process involving nonlinear conformational changes is simulated through iterative cycles consisting of a biased MD simulation with an updated linear response force and a following unbiased MD simulation for relaxation. We applied the method to the simulation of global conformational changes of the yeast calmodulin N-terminal domain and successfully searched out the end conformation. The atomistically detailed trajectories revealed a sequence of molecular events that properly lead to the global conformational changes and identified key steps of local-global coupling that induce the conformational transitions. The LRPF method provides one with a powerful means to model conformational changes of proteins such as motors and transporters where local-global coupling plays a pivotal role in their functional processes.