Structures of β1-adrenergic receptor in complex with Gs and ligands of different efficacies

FFParam-v2.0: A Comprehensive Tool for CHARMM Additive and Drude Polarizable Force-Field Parameter Optimization and Validation

Developing production quality CHARMM force-field (FF) parameters is a very detailed process involving a variety of calculations, many of which are specific for the molecule of interest. The first version of FFParam was developed as a standalone Python package designed for the optimization of electrostatic and bonded parameters of the CHARMM additive and polarizable Drude FFs by using quantum mechanical (QM) target data. The new version of FFParam has multiple new capabilities for FF parameter optimization and validation, with an emphasis on the ability to use condensed-phase target data in optimization. FFParam-v2 allows optimization of Lennard-Jones (LJ) parameters using potential energy scans of interactions between selected atoms in a molecule and noble gases, viz., He and Ne, and through condensed-phase calculations, from which experimental observables such as heats of vaporization and free energies of solvation may be obtained. This functionality serves as a gold standard for both optimizing parameters and validating the performance of the final parameters. A new bonded parameter optimization algorithm has been introduced to account for simultaneously optimizing multiple molecules sharing parameters. FFParam-v2 also supports the comparison of normal modes and the potential energy distribution of internal coordinates towards each normal mode obtained from QM and molecular mechanics calculations. Such comparison capability is vital to validate the balance among various bonded parameters that contribute to the complex normal modes of molecules. User interaction has been extended beyond the original graphical user interface to include command-line interface capabilities that allow for integration of FFParam in workflows, thereby facilitating the automation of parameter optimization. With these new functionalities, FFParam is a more comprehensive parameter optimization tool for both beginners and advanced users.

Real-Time Determination of Intracellular cAMP Reveals Functional Coupling of Gs Protein to the Melatonin MT1 Receptor

Melatonin is a neuroendocrine hormone that regulates the circadian rhythm and many other physiological processes. Its functions are primarily exerted through two subtypes of human melatonin receptors, termed melatonin type-1 (MT1) and type-2 (MT2) receptors. Both MT1 and MT2 receptors are generally classified as Gi-coupled receptors owing to their well-recognized ability to inhibit cAMP accumulation in cells. However, it remains an enigma as to why melatonin stimulates cAMP production in a number of cell types that express the MT1 receptor. To address if MT1 can dually couple to Gs and Gi proteins, we employed a highly sensitive luminescent biosensor (GloSensorTM) to monitor the real-time changes in the intracellular cAMP level in intact live HEK293 cells that express MT1 and/or MT2. Our results demonstrate that the activation of MT1, but not MT2, leads to a robust enhancement on the forskolin-stimulated cAMP formation. In contrast, the activation of either MT1 or MT2 inhibited cAMP synthesis driven by the activation of the Gs-coupled β2-adrenergic receptor, which is consistent with a typical Gi-mediated response. The co-expression of MT1 with Gs enabled melatonin itself to stimulate cAMP production, indicating a productive coupling between MT1 and Gs. The possible existence of a MT1-Gs complex was supported through molecular modeling as the predicted complex exhibited structural and thermodynamic characteristics that are comparable to that of MT1-Gi. Taken together, our data reveal that MT1, but not MT2, can dually couple to Gs and Gi proteins, thereby enabling the bi-directional regulation of adenylyl cyclase to differentially modulate cAMP levels in cells that express different complements of MT1, MT2, and G proteins.

Binding kinetics drive G protein subtype selectivity at the β1-adrenergic receptor

G protein-coupled receptors (GPCRs) bind to different G protein α-subtypes with varying degrees of selectivity. The mechanism by which GPCRs achieve this selectivity is still unclear. Using 13C methyl methionine and 19F NMR, we investigate the agonist-bound active state of β1AR and its ternary complexes with different G proteins in solution. We find the receptor in the ternary complexes adopts very similar conformations. In contrast, the full agonist-bound receptor active state assumes a conformation differing from previously characterised activation intermediates or from β1AR in ternary complexes. Assessing the kinetics of binding for the agonist-bound receptor with different G proteins, we find the increased affinity of β1AR for Gs results from its much faster association with the receptor. Consequently, we suggest a kinetic-driven selectivity gate between canonical and secondary coupling which arises from differential favourability of G protein binding to the agonist-bound receptor active state.

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.

Recognition of methamphetamine and other amines by trace amine receptor TAAR1

Trace amine-associated receptor 1 (TAAR1), the founding member of a nine-member family of trace amine receptors, is responsible for recognizing a range of biogenic amines in the brain, including the endogenous β-phenylethylamine (β-PEA)1 as well as methamphetamine2, an abused substance that has posed a severe threat to human health and society3. Given its unique physiological role in the brain, TAAR1 is also an emerging target for a range of neurological disorders including schizophrenia, depression and drug addiction2,4,5. Here we report structures of human TAAR1-G-protein complexes bound to methamphetamine and β-PEA as well as complexes bound to RO5256390, a TAAR1-selective agonist, and SEP-363856, a clinical-stage dual agonist for TAAR1 and serotonin receptor 5-HT1AR (refs. 6,7). Together with systematic mutagenesis and functional studies, the structures reveal the molecular basis of methamphetamine recognition and underlying mechanisms of ligand selectivity and polypharmacology between TAAR1 and other monoamine receptors. We identify a lid-like extracellular loop 2 helix/loop structure and a hydrogen-bonding network in the ligand-binding pockets, which may contribute to the ligand recognition in TAAR1. These findings shed light on the ligand recognition mode and activation mechanism for TAAR1 and should guide the development of next-generation therapeutics for drug addiction and various neurological disorders.

Catecholamine Derivatives: Natural Occurrence, Structural Diversity, and Biological Activity

Catecholamines (CAs) are aromatic amines containing a 3,4-dihydroxyphenyl nucleus and an amine side chain. Representative CAs included the endogenous neurotransmitters epinephrine, norepinephrine, and dopamine. CAs and their derivatives are good resources for the development of sympathomimetic or central nervous system drugs, while they also provide ligands important for G-protein coupled receptor (GPCR) research. CAs are of broad interest in the fields of chemical, biological, medical, and material sciences due to their high adhesive capacities, chemical reactivities, metal-chelating abilities, redox activities, excellent biocompatibilities, and ease of degradability. Herein, we summarize CAs derivatives isolated and identified from microorganisms, plants, insects, and marine invertebrates in recent decades, alongside their wide range of reported biological activities. The aim of this review is to provide an overview of the structural and biological diversities of CAs, the regularity of their natural occurrences, and insights toward future research and development pertinent to this important class of naturally occurring compounds.

Structural titration reveals Ca2+-dependent conformational landscape of the IP3 receptor

Inositol 1,4,5-trisphosphate receptors (IP3Rs) are endoplasmic reticulum Ca2+ channels whose biphasic dependence on cytosolic Ca2+ gives rise to Ca2+ oscillations that regulate fertilization, cell division and cell death. Despite the critical roles of IP3R-mediated Ca2+ responses, the structural underpinnings of the biphasic Ca2+ dependence that underlies Ca2+ oscillations are incompletely understood. Here, we collect cryo-EM images of an IP3R with Ca2+ concentrations spanning five orders of magnitude. Unbiased image analysis reveals that Ca2+ binding does not explicitly induce conformational changes but rather biases a complex conformational landscape consisting of resting, preactivated, activated, and inhibited states. Using particle counts as a proxy for relative conformational free energy, we demonstrate that Ca2+ binding at a high-affinity site allows IP3Rs to activate by escaping a low-energy resting state through an ensemble of preactivated states. At high Ca2+ concentrations, IP3Rs preferentially enter an inhibited state stabilized by a second, low-affinity Ca2+ binding site. Together, these studies provide a mechanistic basis for the biphasic Ca2+-dependence of IP3R channel activity.

Biased allosteric activation of ketone body receptor HCAR2 suppresses inflammation

Hydroxycarboxylic acid receptor 2 (HCAR2), modulated by endogenous ketone body β-hydroxybutyrate and exogenous niacin, is a promising therapeutic target for inflammation-related diseases. HCAR2 mediates distinct pathophysiological events by activating Gi/o protein or β-arrestin effectors. Here, we characterize compound 9n as a Gi-biased allosteric modulator (BAM) of HCAR2 and exhibit anti-inflammatory efficacy in RAW264.7 macrophages via a specific HCAR2-Gi pathway. Furthermore, four structures of HCAR2-Gi complex bound to orthosteric agonists (niacin or monomethyl fumarate), compound 9n, and niacin together with compound 9n simultaneously reveal a common orthosteric site and a unique allosteric site. Combined with functional studies, we decipher the action framework of biased allosteric modulation of compound 9n on the orthosteric site. Moreover, co-administration of compound 9n with orthosteric agonists could enhance anti-inflammatory effects in the mouse model of colitis. Together, our study provides insight to understand the molecular pharmacology of the BAM and facilitates exploring the therapeutic potential of the BAM with orthosteric drugs.

A Novel "Activation Switch" Motif Common to All Aminergic Receptors

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

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.