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 β
-adrenergic receptors (β
AR), 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 β
AR 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 R131
and Y326
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 β
AR 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.
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.
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 β2-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.
Activation of heterotrimeric G-proteins (Gαβγ) downstream to receptor tyrosine kinases (RTKs) is a well-established crosstalk between the signaling pathways mediated by G-protein coupled receptors ...(GPCRs) and RTKs. While GPCR serves as a guanine exchange factor (GEF) in the canonical activation of Gα that facilitates the exchange of GDP for GTP, the mechanism through which RTK phosphorylations induce Gα activation remains unclear. Recent experimental studies revealed that the epidermal growth factor receptor (EGFR), a well-known RTK, phosphorylates the helical domain tyrosine residues Y154 and Y155 and accelerates the GDP release from the Gαi3, a subtype of Gα-protein. Using well-tempered metadynamics and extensive unbiased molecular dynamics simulations, we captured the GDP release event and identified the intermediates between bound and unbound states through Markov state models. In addition to weakened salt bridges at the domain interface, phosphorylations induced the unfolding of helix αF, which contributed to increased flexibility near the hinge region, facilitating a greater distance between domains in the phosphorylated Gαi3. Although the larger domain separation in the phosphorylated system provided an unobstructed path for the nucleotide, the accelerated release of GDP was attributed to increased fluctuations in several conserved regions like P-loop, switch 1, and switch 2. Overall, this study provides atomistic insights into the activation of G-proteins induced by RTK phosphorylations and identifies the specific structural motifs involved in the process. The knowledge gained from the study could establish a foundation for targeting non-canonical signaling pathways and developing therapeutic strategies against the ailments associated with dysregulated G-protein signaling.
Display omitted
•Well-tempered metadynamics captured the GDP dissociation for unphosphorylated and phosphorylated G-proteins.•Markov state model reveals metastable states of GDP release.•Unfolding of helix αF in phosphorylated G-protein facilitates large interdomain separation.•Phosphorylation increases fluctuations in GDP binding regions, promoting faster GDP release.
G-protein-coupled receptors (GPCRs) constitute the largest family of membrane-bound proteins involved in translation of extracellular signals into intracellular responses. They regulate diverse ...physiological and pathophysiological processes, and hence, they are prime drug targets for therapeutic intervention. In spite of the recent advancements in membrane protein crystallography, limited information is available on the molecular signatures of activation of GPCRs. Although few studies have been reported for class A GPCRs, the activation mechanism of class B GPCRs remains unexplored. Corticotropin-releasing factor 1 receptor (CRF1R), a class B GPCR, is associated with various disease conditions including stress, anxiety, and irritable bowel syndrome. Here, we report the activation of CRF1R using accelerated molecular dynamics simulations of the apo receptor. The breakage of His1552.50–Glu2093.50 and Glu2093.50–Thr3166.42 interactions is found to be crucial in transition of the receptor to its active conformation. Compared to the inactive crystal structure, major structural rearrangements occurred in the intracellular region of the transmembrane (TM) domain upon activation: TM3 twisted away from TM2, and an opening of the G-protein binding site occurred as a result of the outward movements of TM5 and TM6 from the helical bundle. Further, an inward tilt of TM7 toward the helical core is observed at the extracellular side, in agreement with recent findings (Coin et al. Cell 2013, 155, 1258–1269), where it is proposed that this movement helps in establishing favorable interactions with peptide agonist. Moreover, different allosteric pathways in the inactive and active states are identified using the correlations in torsion angle space. The inactive state is found to be less dynamic as compared to the putative active state of the receptor. Results from the current study could present a model for class B GPCRs activation and aid in the design of CRF1R modulators against brain and metabolic disorders.
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.
In order to understand the role of osmolytes in regulating physicochemical behavior of proteins, we investigated the influence of protein destabilizing (urea and guanidinium chloride) and stabilizing ...osmolytes (TMAO, glycerol, and betaine) on a model salt-bridge (SB) formed between structural analogues of arginine and glutamate/aspartate sidechains in a solvent continuum using first-principles quantum chemical calculations based on DFT and MP2 methods. The binding strength of the osmolyte with the SB is found to be in the order of betaine > TMAO > Gdm
+
> glycerol > urea. The osmolytes (TMAO and betaine) that preferentially bind to the SB cation have a marginal influence on SB stability. Also, pure π-π stacking interaction between Gdm
+
and the SB cation plays an insignificant role in destabilizing the SB. In fact, the interaction strength of osmolytes with the SB anion mainly determines the stability of SB. For instance, a competition between Gdm
+
and the SB cation to bind with the SB anion is responsible for instability and subsequent dissociation of the SB. The competition provided by other osmolytes is too weak to break the SB. Exploiting this information, we designed three structural derivatives of Gdm
+
, all having a stronger interaction with SB anion, and thereby show a stronger SB dissociation potential. Furthermore, we find an excellent linear anti-correlation between SB interaction energy and the energy of interaction between osmolyte and the SB anion, which suggests that by knowing only the strength of osmolyte acetate interaction, one can predict the influence of osmolytes on the salt-bridge instability. This information is useful in fine-tuning the SB dissociation power of Gdm
+
, which has a practical significance in obtaining the mechanistic insight into the influence of GdmCl on protein stability. Our results also provide a basis for understanding the chemistry of other ion-pairs formed between a cationic hydrogen donor and an anionic acceptor.
The salt-bridge (SB) interaction energy and the energy of interaction between osmolyte and SB anion are anti-linearly correlated, suggesting that by merely knowing osmolyte acetate interaction, one might predict the influence of osmolytes on a SB.
There is ample evidence that both native functions and pathogenic aggregation of α-synuclein are intimately dependent on lipid interactions and fatty acid type; the regulatory mechanism however ...remains unclear. In the present work, using extensive atomistic molecular dynamics simulations and enhanced-sampling, we have focused on exploring the mechanism of fatty acid dependent regulation of monomeric α-Syn100 in a native synaptic vesicle-like membrane. Our results show that α-Syn100 spontaneously binds to the membrane through its N-terminal region (residues 1–34), where the depth of membrane insertion, the structure, and orientation of the membrane-bound α-Syn100 and its impact on membrane structure are modulated by docosahexaenoic acid (DHA). DHA is a polyunsaturated fatty acid abundantly found in the brain and known to promote the oligomerization of α-synuclein. We found that DHA exhibits marked propensity to interact with monomeric α-Syn100 and modulates the microenvironment of the protein by preferentially sorting DHA-containing phospholipids, depleting other phospholipids and cholesterol as well as increasing the proportion of anionic to neutral lipids in the immediate vicinity of the protein. Owing to the unique conformational flexibility, DHA chains form more lipid-packing defects in the membrane and efficiently coat the membrane-embedded surface of the protein, compared to the saturated and monounsaturated fatty acids. DHA thus makes the bilayer more amiable to protein adsorption and less prone to α-synuclein-induced perturbation associated with cytotoxicity. Indeed, in the absence of DHA, we observed significant thinning of the local bilayer membrane induced by α-Syn100. Though α-Syn100 is predominantly α-helical in membranes studied here, in the presence of DHA we observe formation of β-sheet/β-strands in the C-terminal region (residues 35–100) of α-Syn100, which is extended out from the membrane surface. Notably, DHA induces β structure in the NAC domain of α-Syn100 and promotes extended conformations as well as large solvent exposure of this hydrophobic domain, properties that are known to facilitate self-assembly of α-synuclein. To the best of our knowledge, this study for the first time provides the atomistic insights into DHA-induced regulatory mechanism of monomeric α-synuclein, having implications in protein structure and its physiological/pathological functions.
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. ...In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
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.