Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein

Exploring the constitutive activation mechanism of the class A orphan GPR20

GPR20, an orphan G protein-coupled receptor (GPCR), shows significant expression in intestinal tissue and represents a potential therapeutic target to treat gastrointestinal stromal tumors. GPR20 performs high constitutive activity when coupling with Gi. Despite the pharmacological importance of GPCR constitutive activation, determining the mechanism has long remained unclear. In this study, we explored the constitutive activation mechanism of GPR20 through large-scale unbiased molecular dynamics simulations. Our results unveil the allosteric nature of constitutively activated GPCR signal transduction involving extracellular and intracellular domains. Moreover, the constitutively active state of the GPR20 requires both the N-terminal cap and Gi protein. The N-terminal cap of GPR20 functions like an agonist and mediates long-range activated conformational shift. Together with the previous study, this study enhances our knowledge of the self-activation mechanism of the orphan receptor, facilitates the drug discovery efforts that target GPR20.

Early Events in β2AR Dimer Dynamics Mediated by Activation-Related Microswitches

G-Protein-Coupled Receptors (GPCRs) make up around 3-4% of the human genome and are the targets of one-third of FDA-approved drugs. GPCRs typically exist as monomers but also aggregate to form higher-order oligomers, including dimers. β2AR, a pharmacologically relevant GPCR, is known to be targeted for the treatment of asthma and cardiovascular diseases. The activation of β2AR at the dimer level remains under-explored. In the current study, molecular dynamics (MD) simulations have been performed to understand activation-related structural changes in β2AR at the dimer level. The transition from inactive to active and vice versa has been studied by starting the simulations in the apo, agonist-bound, and inverse agonist-bound β2AR dimers for PDB ID: 2RH1 and PDB ID: 3P0G, respectively. A cumulative total of around 21-μs simulations were performed. Residue-based distances, RMSD, and PCA calculations suggested that either of the one monomer attained activation-related features for the apo and agonist-bound β2AR dimers. The TM5 and TM6 helices within the two monomers were observed to be in significant variation in all the simulations. TM5 bulge and proximity of TM2 and TM7 helices may be contributing to one of the early events in activation. The dimeric interface between TM1 and helix 8 were observed to be well maintained in the apo and agonist-bound simulations. The presence of inverse agonists favored inactive features in both the monomers. These key features of activation known for monomers were observed to have an impact on β2AR dimers, thereby providing an insight into the oligomerization mechanism of GPCRs.

Using the sympathetic system, beta blockers and alpha-2 agonists, to address acute respiratory distress syndrome

Acute Respiratory Distress Syndrome (ARDS) manifests as an acute inflammatory lung injury characterized by persistent hypoxemia, featuring a swift onset, high mortality, and predominantly supportive care as the current therapeutic approach, while effective treatments remain an area of active investigation. Adrenergic receptors (AR) play a pivotal role as stress hormone receptors, extensively participating in various inflammatory processes by initiating downstream signaling pathways. Advancements in molecular biology and pharmacology continually unveil the physiological significance of distinct AR subtypes. Interventions targeting these subtypes have the potential to induce specific alterations in cellular and organismal functions, presenting a promising avenue as a therapeutic target for managing ARDS. This article elucidates the pathogenesis of ARDS and the basic structure and function of AR. It also explores the relationship between AR and ARDS from the perspective of different AR subtypes, aiming to provide new insights for the improvement of ARDS.

Ligand-Receptor Interactions and Structure-Function Relationships in Off-Target Binding of the β3-Adrenergic Agonist Mirabegron to α1A-Adrenergic Receptors

The β3-adrenoceptor agonist mirabegron is available for the treatment of storage symptoms of overactive bladder, including frequency, urgency, and incontinence. The off-target effects of mirabegron include binding to α1-adrenoceptors, which are central in the treatment of voiding symptoms. Here, we examined the structure-function relationships in the binding of mirabegron to a cryo-electron microscopy structure of α1A. The binding was simulated by docking mirabegron to a 3D structure of a human α1A-adrenoceptor (7YMH) using Autodock Vina. The simulations identified two binding states: slope orientation involving 10 positions and horizontal binding to the receptor surface involving 4 positions. No interactions occurred with positions constituting the α1A binding pocket, including Asp-106, Ser-188, or Phe-312, despite the positioning of the phenylethanolamine moiety in transmembrane regions close to the binding pocket by contact with Phe-288, -289, and Val-107. Contact with the unique positions of α1A included the transmembrane Met-292 during slope binding and exosite Phe-86 during horizontal binding. Exosite binding in slope orientation involved contact of the anilino part, rather than the aminothiazol end, to Ile-178, Ala-103, and Asn-179. In conclusion, contact with Met-292 and Phe-86, which are unique positions of α1A, accounts for mirabegron binding to α1A. Because of its lack of interactions with the binding pocket, mirabegron has lower affinity compared to α1A-blockers and no effects on voiding symptoms.

Elucidating molecular mechanisms of protoxin-II state-specific binding to the human NaV1.7 channel

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

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 basis of agonist specificity of α1A-adrenergic receptor

α1-adrenergic receptors (α1-ARs) play critical roles in the cardiovascular and nervous systems where they regulate blood pressure, cognition, and metabolism. However, the lack of specific agonists for all α1 subtypes has limited our understanding of the physiological roles of different α1-AR subtypes, and led to the stagnancy in agonist-based drug development for these receptors. Here we report cryo-EM structures of α1A-AR in complex with heterotrimeric G-proteins and either the endogenous common agonist epinephrine or the α1A-AR-specific synthetic agonist A61603. These structures provide molecular insights into the mechanisms underlying the discrimination between α1A-AR and α1B-AR by A61603. Guided by the structures and corresponding molecular dynamics simulations, we engineer α1A-AR mutants that are not responsive to A61603, and α1B-AR mutants that can be potently activated by A61603. Together, these findings advance our understanding of the agonist specificity for α1-ARs at the molecular level, opening the possibility of rational design of subtype-specific agonists.