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

Adrenoceptor Desensitization: Current Understanding of Mechanisms

G-protein coupled receptors (GPCRs) transduce a wide range of extracellular signals. They are key players in the majority of biologic functions including vision, olfaction, chemotaxis, and immunity. However, as essential as most of them are to body function and homeostasis, overactivation of GPCRs has been implicated in many pathologic diseases such as cancer, asthma, and heart failure (HF). Therefore, an important feature of G protein signaling systems is the ability to control GPCR responsiveness, and one key process to control overstimulation involves initiating receptor desensitization. A number of steps are appreciated in the desensitization process, including cell surface receptor phosphorylation, internalization, and downregulation. Rapid or short-term desensitization occurs within minutes and involves receptor phosphorylation via the action of intracellular protein kinases, the binding of β-arrestins, and the consequent uncoupling of GPCRs from their cognate heterotrimeric G proteins. On the other hand, long-term desensitization occurs over hours to days and involves receptor downregulation or a decrease in cell surface receptor protein level. Of the proteins involved in this biologic phenomenon, β-arrestins play a particularly significant role in both short- and long-term desensitization mechanisms. In addition, β-arrestins are involved in the phenomenon of biased agonism, where the biased ligand preferentially activates one of several downstream signaling pathways, leading to altered cellular responses. In this context, this review discusses the different patterns of desensitization of the α 1-, α 2- and the β adrenoceptors and highlights the role of β-arrestins in regulating physiologic responsiveness through desensitization and biased agonism. SIGNIFICANCE STATEMENT: A sophisticated network of proteins orchestrates the molecular regulation of GPCR activity. Adrenoceptors are GPCRs that play vast roles in many physiological processes. Without tightly controlled desensitization of these receptors, homeostatic imbalance may ensue, thus precipitating various diseases. Here, we critically appraise the mechanisms implicated in adrenoceptor desensitization. A better understanding of these mechanisms helps identify new druggable targets within the GPCR desensitization machinery and opens exciting therapeutic fronts in the treatment of several pathologies.

Computational Insights into Prostaglandin E2 Ligand Binding and Activation of G-Protein-Coupled Receptors

G-protein coupled receptors (GPCRs) are eukaryotic integral membrane proteins that regulate signal transduction cascade pathways implicated in a variety of human diseases and are consequently of interest as drug targets. For this reason, it is of interest to investigate the way in which specific ligands bind and trigger conformational changes in the receptor during activation and how this in turn modulates intracellular signaling. In the present study, we investigate the way in which the ligand Prostaglandin E2 interacts with three GPCRs in the E-prostanoid family: EP1, EP2, and EP3. We examine information transfer pathways based on long-time scale molecular dynamics simulations using transfer entropy and betweenness centrality to measure the physical transfer of information among residues in the system. We monitor specific residues involved in binding to the ligand and investigate how the information transfer behavior of these residues changes upon ligand binding. Our results provide key insights that enable a deeper understanding of EP activation and signal transduction functioning pathways at the molecular level, as well as enabling us to make some predictions about the activation pathway for the EP1 receptor, for which little structural information is currently available. Our results should advance ongoing efforts in the development of potential therapeutics targeting these receptors.

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.

Coevolution-Driven Method for Efficiently Simulating Conformational Changes in Proteins Reveals Molecular Details of Ligand Effects in the β2AR Receptor

With the advent of AI-powered structure prediction, the scientific community is inching closer to solving protein folding. An unresolved enigma, however, is to accurately, reliably, and deterministically predict alternative conformational states that are crucial for the function of, e.g., transporters, receptors, or ion channels where conformational cycling is innately coupled to protein function. Accurately discovering and exploring all conformational states of membrane proteins has been challenging due to the need to retain atomistic detail while enhancing the sampling along interesting degrees of freedom. The challenges include but are not limited to finding which degrees of freedom are relevant, how to accelerate the sampling along them, and then quantifying the populations of each micro- and macrostate. In this work, we present a methodology that finds relevant degrees of freedom by combining evolution and physics through machine learning and apply it to the conformational sampling of the β2 adrenergic receptor. In addition to predicting new conformations that are beyond the training set, we have computed free energy surfaces associated with the protein's conformational landscape. We then show that the methodology is able to quantitatively predict the effect of an array of ligands on the β2 adrenergic receptor activation through the discovery of new metastable states not present in the training set. Lastly, we also stake out the structural determinants of activation and inactivation pathway signaling through different ligands and compare them to functional experiments to validate our methodology and potentially gain further insights into the activation mechanism of the β2 adrenergic receptor.

Dynamic Insights into the Self-Activation Pathway and Allosteric Regulation of the Orphan G-Protein-Coupled Receptor GPR52

Within over 800 members of G-protein-coupled receptors, there are numerous orphan receptors whose endogenous ligands are largely unknown, providing many opportunities for novel drug discovery. However, the lack of an in-depth understanding of the intrinsic working mechanism for orphan receptors severely limits the related rational drug design. The G-protein-coupled receptor 52 (GPR52) is a unique orphan receptor that constitutively increases cellular 5'-cyclic adenosine monophosphate (cAMP) levels without binding any exogenous agonists and has been identified as a promising therapeutic target for central nervous system disorders. Although recent structural biology studies have provided snapshots of both active and inactive states of GPR52, the mechanism of the conformational transition between these states remains unclear. Here, an acceptable self-activation pathway for GPR52 was proposed through 6 μs Gaussian accelerated molecular dynamics (GaMD) simulations, in which the receptor spontaneously transitions from the active state to that matching the inactive crystal structure. According to the three intermediate states of the receptor obtained by constructing a reweighted potential of mean force, how the allosteric regulation occurs between the extracellular orthosteric binding pocket and the intracellular G-protein-binding site is revealed. Combined with the independent gradient model, several important microswitch residues and the allosteric communication pathway that directly links the two regions are both identified. Transfer entropy calculations not only reveal the complex allosteric signaling within GPR52 but also confirm the unique role of ECL2 in allosteric regulation, which is mutually validated with the results of GaMD simulations. Overall, this work elucidates the allosteric mechanism of GPR52 at the atomic level, providing the most detailed information to date on the self-activation of the orphan receptor.

Steric Communication between Dynamic Components on DNA Nanodevices

Biomolecular nanotechnology has helped emulate basic robotic capabilities such as defined motion, sensing, and actuation in synthetic nanoscale systems. DNA origami is an attractive approach for nanorobotics, as it enables creation of devices with complex geometry, programmed motion, rapid actuation, force application, and various kinds of sensing modalities. Advanced robotic functions like feedback control, autonomy, or programmed routines also require the ability to transmit signals among subcomponents. Prior work in DNA nanotechnology has established approaches for signal transmission, for example through diffusing strands or structurally coupled motions. However, soluble communication is often slow and structural coupling of motions can limit the function of individual components, for example to respond to the environment. Here, we introduce an approach inspired by protein allostery to transmit signals between two distal dynamic components through steric interactions. These components undergo separate thermal fluctuations where certain conformations of one arm will sterically occlude conformations of the distal arm. We implement this approach in a DNA origami device consisting of two stiff arms each connected to a base platform via a flexible hinge joint. We demonstrate the ability for one arm to sterically regulate both the range of motion and the conformational state (latched or freely fluctuating) of the distal arm, results that are quantitatively captured by mesoscopic simulations using experimentally informed energy landscapes for hinge-angle fluctuations. We further demonstrate the ability to modulate signal transmission by mechanically tuning the range of thermal fluctuations and controlling the conformational states of the arms. Our results establish a communication mechanism well-suited to transmit signals between thermally fluctuating dynamic components and provide a path to transmitting signals where the input is a dynamic response to parameters like force or solution conditions.

Biased Activation Mechanism Induced by GPCR Heterodimerization: Observations from μOR/δOR Dimers

GPCRs regulate multiple intracellular signaling cascades. Biasedly activating one signaling pathway over the others provides additional clinical utility to optimize GPCR-based therapies. GPCR heterodimers possess different functions from their monomeric states, including their selectivity to different transducers. However, the biased signaling mechanism induced by the heterodimerization remains unclear. Motivated by the issue, we select an important GPCR heterodimer (μOR/δOR heterodimer) as a case and use microsecond Gaussian accelerated molecular dynamics simulation coupled with potential of mean force and protein structure network (PSN) to probe mechanisms regarding the heterodimerization-induced constitutive β-arrestin activity and efficacy change of the agonist DAMGO. The results show that only the lowest energy state of the μOR/δOR heterodimer, which adopts a slightly outward shift of TM6 and an ICL2 conformation close to the receptor core, can selectively accommodate β-arrestins. PSN further reveals important roles of H8, ICL1, and ICL2 in regulating the constitutive β-arrestin-biased activity for the apo μOR/δOR heterodimer. In addition, the heterodimerization can allosterically alter the binding mode of DAMGO mainly by means of W7.35. Consequently, DAMGO transmits the structural signal mainly through TM6 and TM7 in the dimer, rather than TM3 similar to the μOR monomer, thus changing the efficacy of DAMGO from a balanced agonist to the β-arrestin-biased one. On the other side, the binding of DAMGO to the heterodimer can stabilize μOR/δOR heterodimers through a stronger interaction of TM1/TM1 and H8/H8, accordingly enhancing the interaction of μOR with δOR and the binding affinity of the dimer to the β-arrestin. The agonist DAMGO does not change main compositions of the regulation network from the dimer interface to the transducer binding pocket of the μOR protomer, but induces an increase in the structural communication of the network, which should contribute to the enhanced β-arrestin coupling. Our observations, for the first time, reveal the molecular mechanism of the biased signaling induced by the heterodimerization for GPCRs, which should be beneficial to more comprehensively understand the GPCR bias signaling.

Protein Conformational Space at the Edge of Allostery: Turning a Non-allosteric Malate Dehydrogenase into an "Allosterized" Enzyme using Evolution Guided Punctual Mutations

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.