Activation mechanism of the β2-adrenergic receptor

A shared mechanism for TNP-ATP recognition by members of the P2X receptor family

P2X receptors (P2X1-7) are non-selective cation channels involved in many physiological activities such as synaptic transmission, immunological modulation, and cardiovascular function. These receptors share a conserved mechanism to sense extracellular ATP. TNP-ATP is an ATP derivative acting as a nonselective competitive P2X antagonist. Understanding how it occupies the orthosteric site in the absence of agonism may help reveal the key allostery during P2X gating. However, TNP-ATP/P2X complexes (TNP-ATP/human P2X3 (hP2X3) and TNP-ATP/chicken P2X7 (ckP2X7)) with distinct conformations and different mechanisms of action have been proposed. Whether these represent species and subtype variations or experimental differences remains unclear. Here, we show that a common mechanism of TNP-ATP recognition exists for the P2X family members by combining enhanced conformation sampling, engineered disulfide bond analysis, and covalent occupancy. In this model, the polar triphosphate moiety of TNP-ATP interacts with the orthosteric site, while its TNP-moiety is deeply embedded in the head and dorsal fin (DF) interface, creating a restrictive allostery in these two domains that results in a partly enlarged yet ion-impermeable pore. Similar results were obtained from multiple P2X subtypes of different species, including ckP2X7, hP2X3, rat P2X2 (rP2X2), and human P2X1 (hP2X1). Thus, TNP-ATP uses a common mechanism for P2X recognition and modulation by restricting the movements of the head and DF domains which are essential for P2X activation. This knowledge is applicable to the development of new P2X inhibitors.

pKa Calculations of GPCRs: Understanding Protonation States in Receptor Activation

The increase in the available G protein-coupled receptor (GPCR) structures has been pivotal in helping to understand their activation process. However, the role of protonation-conformation coupling in GPCR activation still needs to be clarified. We studied the protonation behavior of the highly conserved Asp2.50 residue in five different class A GPCRs (active and inactive conformations) using a linear response approximation (LRA) pKa calculation protocol. We observed consistent differences (1.3 pK units) for the macroscopic pKa values between the inactive and active states of the A2AR and B2AR receptors, indicating the protonation of Asp2.50 during GPCR activation. This process seems to be specific and not conserved, as no differences were observed in the pKa values of the remaining receptors (CB1R, NT1R, and GHSR).

Quantifying Induced Dipole Effects in Small Molecule Permeation in a Model Phospholipid Bilayer

The cell membrane functions as a semipermeable barrier that governs the transport of materials into and out of cells. The bilayer features a distinct dielectric gradient due to the amphiphilic nature of its lipid components. This gradient influences various aspects of small molecule permeation and the folding and functioning of membrane proteins. Here, we employ polarizable molecular dynamics simulations to elucidate the impact of the electronic environment on the permeation process. We simulated eight distinct amino-acid side chain analogs within a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer using the Drude polarizable force field (FF). Our approach includes both unbiased and umbrella sampling simulations. By using a polarizable FF, we sought to investigate explicit dipole responses in relation to local electric fields along the membrane normal. We evaluate molecular dipole moments, which exhibit variation based on their localization within the membrane, and compare the outcomes with analogous simulations using the nonpolarizable CHARMM36 FF. This comparative analysis aims to discern characteristic differences in the free energy surfaces of permeation for the various amino-acid analogs. Our results provide the first systematic quantification of the impact of employing an explicitly polarizable FF in this context compared to the fixed-charge convention inherent to nonpolarizable FFs, which may not fully capture the influence of the membrane dielectric gradient.

The complex nature of CXCR4 mutations in WHIM syndrome

Heterozygous autosomal dominant mutations in the CXCR4 gene cause WHIM syndrome, a severe combined immunodeficiency disorder. The mutations primarily affect the C-terminal region of the CXCR4 chemokine receptor, specifically several potential phosphorylation sites critical for agonist (CXCL12)-mediated receptor internalization and desensitization. Mutant receptors have a prolonged residence time on the cell surface, leading to hyperactive signaling that is responsible for some of the symptoms of WHIM syndrome. Recent studies have shown that the situation is more complex than originally thought, as mutant WHIM receptors and CXCR4 exhibit different dynamics at the cell membrane, which also influences their respective cellular functions. This review examines the functional mechanisms of CXCR4 and the impact of WHIM mutations in both physiological and pathological conditions.

TRPV1 analgesics disturb core body temperature via a biased allosteric mechanism involving conformations distinct from that for nociception

Efforts on developing transient receptor potential vanilloid 1 (TRPV1) drugs for pain management have been hampered by deleterious hypo- or hyperthermia caused by TRPV1 agonists/antagonists. Here, we compared the effects of four antagonists on TRPV1 polymodal gating and core body temperature (CBT) in Trpv1+/+, Trpv1-/-, and Trpv1T634A/T634A. Neither the effect on proton gating nor drug administration route, hair coverage, CBT rhythmic fluctuations, or inflammation had any influence on the differential actions of TRPV1 drugs on CBT. We identified the S4-S5 linker region exposed to the vanilloid pocket of TRPV1 to be critical for hyperthermia associated with certain TRPV1 antagonists. PSFL2874, a TRPV1 antagonist we discovered, is effective against inflammatory pain but devoid of binding to the S4-S5 linker and inducing CBT changes. These findings implicate that biased allosteric mechanisms exist for TRPV1 coupling to nociception and CBT regulation, opening avenues for the development of non-opioid analgesics without affecting CBT.

Exploring GPCR conformational dynamics using single-molecule fluorescence

G protein-coupled receptors (GPCRs) are membrane proteins that transmit specific external stimuli into cells by changing their conformation. This conformational change allows them to couple and activate G-proteins to initiate signal transduction. A critical challenge in studying and inferring these structural dynamics arises from the complexity of the cellular environment, including the presence of various endogenous factors. Due to the recent advances in cell-expression systems, membrane-protein purification techniques, and labeling approaches, it is now possible to study the structural dynamics of GPCRs at a single-molecule level both in vitro and in live cells. In this review, we discuss state-of-the-art techniques and strategies for expressing, purifying, and labeling GPCRs in the context of single-molecule research. We also highlight four recent studies that demonstrate the applications of single-molecule microscopy in revealing the dynamics of GPCRs. These techniques are also useful as complementary methods to verify the results obtained from other structural biology tools like cryo-electron microscopy and x-ray crystallography.

Computational insights into diverse binding modes of the allosteric modulator and their regulation on dopamine D1 receptor

Allosteric drugs hold the promise of addressing many challenges in the current drug development of GPCRs. However, the molecular mechanism underlying their allosteric modulations remain largely elusive. The dopamine D1 receptor (DRD1), a member of Class A GPCRs, is critical for treating psychiatric disorders, and LY3154207 serves as its promising positive allosteric modulator (PAM). In the work, we utilized extensive Gaussian-accelerated molecular dynamics simulations (a total of 41μs) for the first time probe the diverse binding modes of the allosteric modulator and their regulation effects, based on the DRD1 and LY3154207 as representative. Our simulations identify four binding modes of LY3154207 (one boat mode, two metastable vertical modes and a novel cleft-anchored mode), in which the boat mode is the most stable while there three modes are similar in the stability. However, it is interesting to observed that the most stable boat mode inversely exhibits the weakest positive allosteric effect on influencing the orthosteric ligand binding and maintaining the activity of the transducer binding site. It should result from its induced weaker correlation between the allosteric site and the orthosteric site, and between the orthosteric site and the transducer binding site than the other three binding modes, as well as its weakened interaction between a crucial activation-related residue (S2025.46) and the orthosteric ligand (dopamine). Overall, the work offers atomic-level information to advance our understanding of the complex allosteric regulation on GPCRs, which is beneficial to the allosteric modulator design and development.

Computational Insights into the Conformational Dynamics of HIV-1 Vpr in a Lipid Bilayer for Ion Channel Modeling

HIV-1 Vpr is a multifunctional accessory protein consisting of 96 amino acids that play a critical role in viral pathogenesis. Among its diverse range of activities, Vpr can create a cation-selective ion channel within the plasma membrane. However, the oligomeric state of this channel has not yet been elucidated. In this study, we investigated the conformational dynamics of Vpr helices to model the ion channel topology. First, we employed a series of multiscale simulations to investigate the specific structure of monomeric Vpr in a membrane model. During the lipid bilayer self-assembly coarse grain simulation, the C-terminal helix (residues 56-77) effectively formed the transmembrane region, while the N-terminal helix exhibited an amphipathic nature by associating horizontally with a single leaflet. All-atom molecular dynamics (MD) simulations of full-length Vpr inside a phospholipid bilayer show that the C-terminal helix remains very stable inside the bilayer core in a vertical orientation. Subsequently, using the predicted C-terminal helix orientation and conformation, various oligomeric states (ranging from tetramer to heptamer) possibly forming the Vpr ion channel were built and further evaluated. Among these models, the pentameric form exhibited consistent stability in MD simulations and displayed a compatible conformation for a water-assisted ion transport mechanism. This study provides structural insights into the ion channel activity of the Vpr protein and the foundation for developing therapeutics against HIV-1 Vpr-related conditions.

Calcium-Driven In Silico Inactivation of a Human Olfactory Receptor

Conformational changes as well as molecular determinants related to the activation and inactivation of olfactory receptors are still poorly understood due to the intrinsic difficulties in the structural determination of this GPCR family. Here, we perform, for the first time, the in silico inactivation of human olfactory receptor OR51E2, highlighting the possible role of calcium in this receptor state transition. Using molecular dynamics simulations, we show that a divalent ion in the ion binding site, coordinated by two acidic residues at positions 2.50 and 3.39 conserved across most ORs, stabilizes the receptor in its inactive state. In contrast, protonation of the same two acidic residues is not sufficient to drive inactivation within the microsecond timescale of our simulations. Our findings suggest a novel molecular mechanism for OR inactivation, potentially guiding experimental validation and offering insights into the possible broader role of divalent ions in GPCR signaling.

Structural insights into the orthosteric inhibition of P2X receptors by non-ATP analog antagonists

P2X receptors are extracellular ATP-gated ion channels that form homo- or heterotrimers and consist of seven subtypes. They are expressed in various tissues, including neuronal and nonneuronal cells, and play critical roles in physiological processes such as neurotransmission, inflammation, pain, and cancer. As a result, P2X receptors have attracted considerable interest as drug targets, and various competitive inhibitors have been developed. However, although several P2X receptor structures from different subtypes have been reported, the limited structural information of P2X receptors in complex with competitive antagonists hampers the understanding of orthosteric inhibition, hindering the further design and optimization of those antagonists for drug discovery. We determined the cryogenic electron microscopy (cryo-EM) structures of the mammalian P2X7 receptor in complex with two classical competitive antagonists of pyridoxal-5'-phosphate derivatives, pyridoxal-5'-phosphate-6-(2'-naphthylazo-6'-nitro-4',8'-disulfonate) (PPNDS) and pyridoxal phosphate-6-azophenyl-2',5'-disulfonic acid (PPADS), and performed structure-based mutational analysis by patch-clamp recording as well as molecular dynamics (MD) simulations. Our structures revealed the orthosteric site for PPADS/PPNDS, and structural comparison with the previously reported apo- and ATP-bound structures showed how PPADS/PPNDS binding inhibits the conformational changes associated with channel activation. In addition, structure-based mutational analysis identified key residues involved in the PPNDS sensitivity of P2X1 and P2X3, which are known to have higher affinity for PPADS/PPNDS than other P2X subtypes.

Switching Go̅-Martini for Investigating Protein Conformational Transitions and Associated Protein-Lipid Interactions

Proteins are dynamic biomolecules that can transform between different conformational states when exerting physiological functions, which is difficult to simulate using all-atom methods. Coarse-grained (CG) Go̅-like models are widely used to investigate large-scale conformational transitions, which usually adopt implicit solvent models and therefore cannot explicitly capture the interaction between proteins and surrounding molecules, such as water and lipid molecules. Here, we present a new method, named Switching Go̅-Martini, to simulate large-scale protein conformational transitions between different states, based on the switching Go̅ method and the CG Martini 3 force field. The method is straightforward and efficient, as demonstrated by the benchmarking applications for multiple protein systems, including glutamine binding protein (GlnBP), adenylate kinase (AdK), and β2-adrenergic receptor (β2AR). Moreover, by employing the Switching Go̅-Martini method, we can not only unveil the conformational transition from the E2Pi-PL state to E1 state of the type 4 P-type ATPase (P4-ATPase) flippase ATP8A1-CDC50 but also provide insights into the intricate details of lipid transport.

Unraveling Activation-Related Rearrangements and Intrinsic Divergence from Ligand-Specific Conformational Changes of the Dopamine D3 and D2 Receptors

Effective rational drug discovery hinges on understanding the functional states of the target protein and distinguishing it from homologues. However, for the G protein coupled receptors, both activation-related conformational changes (ACCs) and intrinsic divergence among receptors can be misled or obscured by ligand-specific conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics (MD) simulation results of the receptors bound with various ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors, including the extracellular portion of TM5 (TM5e) and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found more outward tilting of TM6e in the D2R compared to the D3R in both the experimental structures and simulations bound with ligands in different scaffolds. However, this difference was drastically reduced in the simulations bound with nonselective agonist quinpirole, suggesting a misleading effect of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may also obscure intrinsic divergence. Importantly, our MD simulations revealed divergence in the dynamics of these receptors. Specifically, the D2R exhibited heightened flexibility compared to the D3R in the extracellular loops and TMs 5e, 6e, and 7e, associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting the D2R and D3R with more precise pharmacological profiles.

Determinants of Neutral Antagonism and Inverse Agonism in the β2-Adrenergic Receptor

Free-energy profiles for the activation/deactivation of the β2-adrenergic receptor (ADRB2) with neutral antagonist and inverse agonist ligands have been determined with well-tempered multiple-walker (MW) metadynamics simulations. The inverse agonists carazolol and ICI118551 clearly favor single inactive conformational minima in both the binary and ternary ligand-receptor-G-protein complexes, in accord with the inverse-agonist activity of the ligands. The behavior of neutral antagonists is more complex, as they seem also to affect the recruitment of the G-protein. The results are analyzed in terms of the conformational states of the well-known microswitches that have been proposed as indicators of receptor activity.

State-dependent dynamics of extramembrane domains in the β2 -adrenergic receptor

G protein-coupled receptors (GPCRs) are membrane-bound signaling proteins that play an essential role in cellular signaling processes. Due to their intrinsic function of transmitting internal signals in response to external cues, these receptors are adapted to be highly dynamic in nature. The β2 -adrenergic receptor (β2 AR) is a representative member of the family that has been extensively analyzed in terms of its structure and activation. Although the structure of the transmembrane domain has been characterized in the different functional states of the receptor, the conformational dynamics of the extramembrane domains, especially the intrinsically disordered regions are still emerging. In this study, we analyze the state-dependent dynamics of extramembrane domains of β2 AR using atomistic molecular dynamics simulations. We introduce a parameter, the residue excess dynamics that allows us to better quantify receptor dynamics. Using this measure, we show that the dynamics of the extramembrane domains are sensitive to the receptor state. Interestingly, the ligand-bound intermediateR ' state shows the maximal dynamics compared to either the active R*G or inactive R states. Ligand binding appears to be correlated with high residue excess dynamics that are dampened upon G protein coupling. The intracellular loop-3 (ICL3) domain has a tendency to flip towards the membrane upon ligand binding, which could contribute to receptor "priming." We highlight an important ICL1-helix-8 interplay that is broken in the ligand-bound state but is retained in the active state. Overall, our study highlights the importance of characterizing the functional dynamics of the GPCR loop domains.

Fine-tuning activation specificity of G-protein-coupled receptors via automated path searching

Physics-based simulation methods can grant atomistic insights into the molecular origin of the function of biomolecules. However, the potential of such approaches has been hindered by their low efficiency, including in the design of selective agonists where simulations of myriad protein-ligand combinations are necessary. Here, we describe an automated input-free path searching protocol that offers (within 14 d using Graphics Processing Unit servers) a minimum free energy path (MFEP) defined in high-dimension configurational space for activating sphingosine-1-phosphate receptors (S1PRs) by arbitrary ligands. The free energy distributions along the MFEP for four distinct ligands and three S1PRs reached a remarkable agreement with Bioluminescence Resonance Energy Transfer (BRET) measurements of G-protein dissociation. In particular, the revealed transition state structures pointed out toward two S1PR3 residues F263/I284, that dictate the preference of existing agonists CBP307 and BAF312 on S1PR1/5. Swapping these residues between S1PR1 and S1PR3 reversed their response to the two agonists in BRET assays. These results inspired us to design improved agonists with both strong polar head and bulky hydrophobic tail for higher selectivity on S1PR1. Through merely three in silico iterations, our tool predicted a unique compound scaffold. BRET assays confirmed that both chiral forms activate S1PR1 at nanomolar concentration, 1 to 2 orders of magnitude less than those for S1PR3/5. Collectively, these results signify the promise of our approach in fine agonist design for G-protein-coupled receptors.

Engineering CREB-activated promoters for adenosine-induced gene expression

Accumulation of the ribonucleoside, adenosine (ADO), triggers a cAMP response element binding protein (CREB)-mediated signaling pathway to suppress the function of immune cells in tumors. Here, we describe a collection of CREB-activated promoters that allow for strong and tunable ADO-induced gene expression in human cells. By optimizing number of CREB transcription factor binding sites and altering the core promoter region of CREB-based hybrid promoters, we created synthetic constructs that drive gene expression to higher levels than strong, endogenous mammalian promoters in the presence of ADO. These synthetic promoters are induced up to 47-fold by ADO, with minimal expression in their "off" state. We further determine that our CREB-based promoters are activated by other compounds that act as signaling analogs, and that combinatorial addition of ADO and these compounds has a synergistic impact on gene expression. Surprisingly, we also detail how background ADO degradation caused by the common cell culture media additive, fetal bovine serum (FBS), confounds experiments designed to determine ADO dose-responsiveness. We show that only after long-term heat deactivation of FBS can our synthetic promoters enable gene expression induction at physiologically relevant levels of ADO. Finally, we demonstrate that the strength of a CREB-based promoter is enhanced by incorporating other transcription factor binding sites.

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.

Characterizing conformational states in GPCR structures using machine learning

G protein-coupled receptors (GPCRs) play a pivotal role in signal transduction and represent attractive targets for drug development. Recent advances in structural biology have provided insights into GPCR conformational states, which are critical for understanding their signaling pathways and facilitating structure-based drug discovery. In this study, we introduce a machine learning approach for conformational state annotation of GPCRs. We represent GPCR conformations as high-dimensional feature vectors, incorporating information about amino acid residue pairs involved in the activation pathway. Using a dataset of GPCR conformations in inactive and active states obtained through molecular dynamics simulations, we trained machine learning models to distinguish between inactive-like and active-like conformations. The developed model provides interpretable predictions and can be used for the large-scale analysis of molecular dynamics trajectories of GPCRs.

The Pivotal Distinction between Antagonists' and Agonists' Binding into Dopamine D4 Receptor-MD and FMO/PIEDA Studies

The dopamine D4 receptor (D4R) is a promising therapeutic target in widespread diseases, and the search for novel agonists and antagonists appears to be clinically relevant. The mechanism of binding to the receptor (R) for antagonists and agonists varies. In the present study, we conducted an in-depth computational study, teasing out key similarities and differences in binding modes, complex dynamics, and binding energies for D4R agonists and antagonists. The dynamic network method was applied to investigate the communication paths between the ligand (L) and G-protein binding site (GBS) of human D4R. Finally, the fragment molecular orbitals with pair interaction energy decomposition analysis (FMO/PIEDA) scheme was used to estimate the binding energies of L-R complexes. We found that a strong salt bridge with D3.32 initiates the inhibition of the dopamine D4 receptor. This interaction also occurs in the binding of agonists, but the change in the receptor conformation to the active state starts with interaction with cysteine C3.36. Such a mechanism may arise in the case of agonists unable to form a hydrogen bond with the serine S5.46, considered, so far, to be crucial in the activation of GPCRs. The energy calculations using the FMO/PIEDA method indicate that antagonists show higher residue occupancy of the receptor binding site than agonists, suggesting they could form relatively more stable complexes. Additionally, antagonists were characterized by repulsive interactions with S5.46 distinguishing them from agonists.

Allosteric modulation of G protein-coupled receptors as a novel therapeutic strategy in neuropathic pain

Neuropathic pain is a debilitating pathological condition that presents significant therapeutic challenges in clinical practice. Unfortunately, current pharmacological treatments for neuropathic pain lack clinical efficacy and often lead to harmful adverse reactions. As G protein-coupled receptors (GPCRs) are widely distributed throughout the body, including the pain transmission pathway and descending inhibition pathway, the development of novel neuropathic pain treatments based on GPCRs allosteric modulation theory is gaining momentum. Extensive research has shown that allosteric modulators targeting GPCRs on the pain pathway can effectively alleviate symptoms of neuropathic pain while reducing or eliminating adverse effects. This review aims to provide a comprehensive summary of the progress made in GPCRs allosteric modulators in the treatment of neuropathic pain, and discuss the potential benefits and adverse factors of this treatment. We will also concentrate on the development of biased agonists of GPCRs, and based on important examples of biased agonist development in recent years, we will describe universal strategies for designing structure-based biased agonists. It is foreseeable that, with the continuous improvement of GPCRs allosteric modulation and biased agonist theory, effective GPCRs allosteric drugs will eventually be available for the treatment of neuropathic pain with acceptable safety.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor

Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.

β-adrenergic receptor regulates embryonic epithelial extensibility through actomyosin inhibition

During morphogenesis, epithelial tissues reshape and expand to cover the body and organs. The molecular mechanisms of this deformability remain elusive. Here, we investigate the role of the β-adrenergic receptor (ADRB) in orchestrating actomyosin contractility, pivotal for epithelial extensibility. Chemical screens on Xenopus laevis embryos pinpointed ADRB2 as a principal regulator. ADRB2 promotes actomyosin relaxation, facilitating apical cell area expansion during body elongation. In contrast, ADRB2 knockdown results in heightened cell contraction, marked by synchronous oscillation of F-actin and myosin, impeding body elongation. ADRB2 mutants with reduced affinity for ligand binding lack the function to induce cellular relaxation, highlighting the ligand's essential roles even in the developing epidermis. Our findings unveil ADRB2's critical contribution to extensibility of the epidermis and subsequent body elongation during development. This study also offers insights into the physiology of mature epithelial organs deformed by the smooth muscle response to the adrenergic autonomic nervous system.

Is the Functional Response of a Receptor Determined by the Thermodynamics of Ligand Binding?

For an effective drug, strong binding to the target protein is a prerequisite, but it is not enough. To produce a particular functional response, drugs need to either block the proteins' functions or modulate their activities by changing their conformational equilibrium. The binding free energy of a compound to its target is routinely calculated, but the timescales for the protein conformational changes are prohibitively long to be efficiently modeled via physics-based simulations. Thermodynamic principles suggest that the binding free energies of the ligands with different receptor conformations may infer their efficacy. However, this hypothesis has not been thoroughly validated. We present an actionable protocol and a comprehensive study to show that binding thermodynamics provides a strong predictor of the efficacy of a ligand. We apply the absolute binding free energy perturbation method to ligands bound to active and inactive states of eight G protein-coupled receptors and a nuclear receptor and then compare the resulting binding free energies. We find that carefully designed restraints are often necessary to efficiently model the corresponding conformational ensembles for each state. Our method achieves unprecedented performance in classifying ligands as agonists or antagonists across the various investigated receptors, all of which are important drug targets.

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.

Engineered odorant receptors illuminate structural principles of odor discrimination

A central challenge in olfaction is understanding how the olfactory system detects and distinguishes odorants with diverse physicochemical properties and molecular configurations. Vertebrate animals perceive odors via G protein-coupled odorant receptors (ORs). In humans, ~400 ORs enable the sense of smell. The OR family is composed of two major classes: Class I ORs are tuned to carboxylic acids while Class II ORs, representing the vast majority of the human repertoire, respond to a wide variety of odorants. How ORs recognize chemically diverse odorants remains poorly understood. A fundamental bottleneck is the inability to visualize odorant binding to ORs. Here, we uncover fundamental molecular properties of odorant-OR interactions by employing engineered ORs crafted using a consensus protein design strategy. Because such consensus ORs (consORs) are derived from the 17 major subfamilies of human ORs, they provide a template for modeling individual native ORs with high sequence and structural homology. The biochemical tractability of consORs enabled four cryoEM structures of distinct consORs with unique ligand recognition properties. The structure of a Class I consOR, consOR51, showed high structural similarity to the native human receptor OR51E2 and yielded a homology model of a related member of the human OR51 family with high predictive power. Structures of three Class II consORs revealed distinct modes of odorant-binding and activation mechanisms between Class I and Class II ORs. Thus, the structures of consORs lay the groundwork for understanding molecular recognition of odorants by the OR superfamily.

Unraveling activation-related rearrangements and intrinsic divergence from ligand-induced conformational changes of the dopamine D3 and D2 receptors

Effective rational drug discovery targeting a specific protein hinges on understanding their functional states and distinguishing it from homologs. However, for the G protein coupled receptors, both the activation-related conformational changes (ACCs) and the intrinsic divergence among receptors can be misled or obscured by ligand-induced conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics simulation results of the receptors bound with different ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors including TM5e and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found pronounced outward tilting of TM6e in the D2R compared to the D3R in both experimental structures and simulations with ligands in different scaffolds. This tilting was drastically reduced in the simulations of the receptors bound with nonselective full agonist quinpirole, suggesting a misleading impact of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may obscure intrinsic divergence. In addition, our analysis showed that the impact of the nonconserved TM1 propagated to conserved Trp7.40 and Glu2.65, both are ligand binding residues. We also found that the D2R exhibited heightened flexibility compared to the D3R in the extracellular portions of TMs 5, 6, and 7, potentially associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting D2R or D3R with more precise pharmacological profiles.

Collective and Coordinated Conformational Changes Determine Agonism or Antagonism at the Human Trace Amine-Associated Receptor 1

The human trace amine-associated receptor (hTAAR1), a G protein-coupled receptor, has been postulated as a new target in the treatment of neuropsychiatric conditions. The mechanism associated with activation or inactivation by agonists or antagonists in hTAAR1 and other GPCRs has not yet been fully elucidated. In this study, we combined computational methods including homology modeling, docking, and molecular dynamic simulations to reveal novel conformational changes associated with agonist and antagonist interactions in hTAAR1. Our findings suggest a differential cascade of coordinated movements based on the presence of either an agonist or antagonist and primarily involving the second extracellular loop, transmembrane domain 5, and the third intracellular domains of hTAAR1. Our study provides an opportunity to predict the effects on new ligands with agonistic or antagonistic activity at hTAAR1 based on the reported conformational changes.

Subtle Structural Modification of a Synthetic Cannabinoid Receptor Agonist Drastically Increases its Efficacy at the CB1 Receptor

The emergence of synthetic cannabinoid receptor agonists (SCRAs) as illicit psychoactive substances has posed considerable public health risks, including fatalities. Many SCRAs exhibit much higher efficacy and potency compared with the phytocannabinoid Δ9-tetrahydrocannabinol (THC) at the cannabinoid receptor 1 (CB1R), leading to dramatic differences in signaling levels that can be toxic. In this study, we investigated the structure-activity relationships of aminoalkylindole SCRAs at CB1Rs, focusing on 5F-pentylindoles containing an amide linker attached to different head moieties. Using in vitro bioluminescence resonance energy transfer assays, we identified a few SCRAs exhibiting significantly higher efficacy in engaging the Gi protein and recruiting β-arrestin than the reference CB1R full agonist CP55940. Importantly, the extra methyl group on the head moiety of 5F-MDMB-PICA, as compared to that of 5F-MMB-PICA, led to a large increase in efficacy and potency at the CB1R. This pharmacological observation was supported by the functional effects of these SCRAs on glutamate field potentials recorded in hippocampal slices. Molecular modeling and simulations of the CB1R models bound with both of the SCRAs revealed critical structural determinants contributing to the higher efficacy of 5F-MDMB-PICA and how these subtle differences propagated to the receptor-G protein interface. Thus, we find that apparently minor structural changes in the head moiety of SCRAs can cause major changes in efficacy. Our results highlight the need for close monitoring of the structural modifications of newly emerging SCRAs and their potential for toxic drug responses in humans.

Activation/Deactivation Free-Energy Profiles for the β2-Adrenergic Receptor: Ligand Modes of Action

We use enhanced-sampling simulations with an effective collective variable to study the activation of the β2-adrenergic receptor in the presence of ligands with different efficacy. The free-energy profiles are computed for the ligand-free (apo) receptor and binary (apo-receptor + G-protein α-subunit and receptor + ligand) and ternary complexes. The results are not only compatible with available experiments but also allow unprecedented structural insight into the nature of GPCR conformations along the activation pathway and their role in the activation mechanism. In particular, the simulations reveal an unexpected mode of action of partial agonists such as salmeterol and salbutamol that arises already in the binary complex without the G-protein. Specific differences in the polar interactions with residues in TM5, which are required to stabilize an optimal TM6 conformation that facilitates G-protein binding and receptor activation, play a major role in differentiating them from full agonists.

Interplay of thermodynamics and evolution within the ternary ligand-GPCR-G protein complex

Recent studies on G protein-coupled receptors (GPCRs) dynamics report that GPCRs adopt a wide range of conformations that coexist in equilibrium, with the apo state of a GPCR having a high entropy. The formation of a ligand-GPCR-transducer complex comes with a reduction of conformational space and therefore with an entropic cost. We hypothesize that the availability of binding partners, their binding affinity and the rigidity of the respective binding sites are reflected in a distinct degree of sequence conservation to balance the energetic cost of intra- and extracellular binding events. Here, we outline the current findings in delineating the conformational space and include sequential conservation of many-to-many ligand-receptor systems to discuss the entropic cost that comes with GPCR signal transduction.

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.

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.

The World of GPCR dimers - Mapping dopamine receptor D2 homodimers in different activation states and configuration arrangements

G protein-coupled receptors (GPCRs) are known to dimerize, but the molecular and structural basis of GPCR dimers is not well understood. In this study, we developed a computational framework to generate models of symmetric and asymmetric GPCR dimers using different monomer activation states and identified their most likely interfaces with molecular details. We chose the dopamine receptor D2 (D2R) homodimer as a case study because of its biological relevance and the availability of structural information. Our results showed that transmembrane domains 4 and 5 (TM4 and TM5) are mostly found at the dimer interface of the D2R dimer and that these interfaces have a subset of key residues that are mostly nonpolar from TM4 and TM5, which was in line with experimental studies. In addition, TM2 and TM3 appear to be relevant for D2R dimers. In some cases, the inactive configuration is unaffected by the partnered protomer, whereas in others, the active protomer adopts the properties of an inactive receptor. Additionally, the β-arrestin configuration displayed the properties of an active receptor in the absence of an agonist, suggesting that a switch to another meta-state during dimerization occurred. Our findings are consistent with the experimental data, and this method can be adapted to study heterodimers and potentially extended to include additional proteins such as G proteins or β-arrestins. In summary, this approach provides insight into the impact of the conformational status of partnered protomers on the overall quaternary GPCR macromolecular structure and dynamics.

The Molecular Mechanism of Positive Allosteric Modulation at the Dopamine D1 Receptor

The dopamine D1 receptor (D1R) is a promising target for treating various psychiatric disorders. While upregulation of D1R activity has shown potential in alleviating motor and cognitive symptoms, orthosteric agonists have limitations, restricting their clinical applications. However, the discovery of several allosteric compounds specifically targeting the D1R, such as LY3154207, has opened new therapeutic avenues. Based on the cryo-EM structures of the D1R, we conducted molecular dynamics simulations to investigate the binding and allosteric mechanisms of LY3154207. Our simulations revealed that LY3154207 preferred the horizontal orientation above intracellular loop 2 (IL2) and stabilized the helical conformation of IL2. Moreover, LY3154207 binding induced subtle yet significant changes in key structural motifs and their neighboring residues. Notably, a cluster of residues centered around the Na+-binding site became more compact, while interactions involving the PIF motif and its neighboring residues were loosened upon LY3154207 binding, consistent with their role in opening the intracellular crevice for receptor activation. Additionally, we identified an allosteric pathway likely responsible for the positive allosteric effect of LY3154207 in enhancing Gs protein coupling. This mechanistic understanding of LY3154207's allosteric action at the D1R paves the way for the rational design of more potent and effective allosteric modulators.

A non-canonical mechanism of GPCR activation

The goal of designing safer, more effective drugs has led to tremendous interest in molecular mechanisms through which ligands can precisely manipulate signaling of G-protein-coupled receptors (GPCRs), the largest class of drug targets. Decades of research have led to the widely accepted view that all agonists-ligands that trigger GPCR activation-function by causing rearrangement of the GPCR's transmembrane helices, opening an intracellular pocket for binding of transducer proteins. Here we demonstrate that certain agonists instead trigger activation of free fatty acid receptor 1 by directly rearranging an intracellular loop that interacts with transducers. We validate the predictions of our atomic-level simulations by targeted mutagenesis; specific mutations which disrupt interactions with the intracellular loop convert these agonists into inverse agonists. Further analysis suggests that allosteric ligands could regulate signaling of many other GPCRs via a similar mechanism, offering rich possibilities for precise control of pharmaceutically important targets.

Membrane mediated mechanical stimuli produces distinct active-like states in the AT1 receptor

The Angiotensin II Type 1 (AT1) receptor is one of the most widely studied GPCRs within the context of biased signaling. While the AT1 receptor is activated by agonists such as the peptide AngII, it can also be activated by mechanical stimuli such as membrane stretch or shear in the absence of a ligand. Despite the importance of mechanical activation of the AT1 receptor in biological processes such as vasoconstriction, little is known about the structural changes induced by external physical stimuli mediated by the surrounding lipid membrane. Here, we present a systematic simulation study that characterizes the activation of the AT1 receptor under various membrane environments and mechanical stimuli. We show that stability of the active state is highly sensitive to membrane thickness and tension. Structural comparison of membrane-mediated vs. agonist-induced activation shows that the AT1 receptor has distinct active conformations. This is supported by multi-microsecond free energy calculations that show unique landscapes for the inactive and various active states. Our modeling results provide structural insights into the mechanical activation of the AT1 receptor and how it may produce different functional outcomes within the framework of biased agonism.

Identification of a β-arrestin-biased negative allosteric modulator for the β2-adrenergic receptor

Catecholamine-stimulated β2-adrenergic receptor (β2AR) signaling via the canonical Gs-adenylyl cyclase-cAMP-PKA pathway regulates numerous physiological functions, including the therapeutic effects of exogenous β-agonists in the treatment of airway disease. β2AR signaling is tightly regulated by GRKs and β-arrestins, which together promote β2AR desensitization and internalization as well as downstream signaling, often antithetical to the canonical pathway. Thus, the ability to bias β2AR signaling toward the Gs pathway while avoiding β-arrestin-mediated effects may provide a strategy to improve the functional consequences of β2AR activation. Since attempts to develop Gs-biased agonists and allosteric modulators for the β2AR have been largely unsuccessful, here we screened small molecule libraries for allosteric modulators that selectively inhibit β-arrestin recruitment to the receptor. This screen identified several compounds that met this profile, and, of these, a difluorophenyl quinazoline (DFPQ) derivative was found to be a selective negative allosteric modulator of β-arrestin recruitment to the β2AR while having no effect on β2AR coupling to Gs. DFPQ effectively inhibits agonist-promoted phosphorylation and internalization of the β2AR and protects against the functional desensitization of β-agonist mediated regulation in cell and tissue models. The effects of DFPQ were also specific to the β2AR with minimal effects on the β1AR. Modeling, mutagenesis, and medicinal chemistry studies support DFPQ derivatives binding to an intracellular membrane-facing region of the β2AR, including residues within transmembrane domains 3 and 4 and intracellular loop 2. DFPQ thus represents a class of biased allosteric modulators that targets an allosteric site of the β2AR.

The molecular mechanism of positive allosteric modulation at the dopamine D1 receptor

The dopamine D1 receptor (D1R) is a promising target for treating various psychiatric disorders. While upregulation of D1R activity has shown potential in alleviating motor and cognitive symptoms, orthosteric agonists have limitations, restricting their clinical applications. However, the discovery of several allosteric compounds specifically targeting the D1R, such as LY3154207, has opened new therapeutic avenues. Based on the cryo-EM structures of the D1R, we conducted molecular dynamics simulations to investigate the binding and allosteric mechanisms of LY3154207. Our simulations revealed that LY3154207 preferred the horizontal orientation above intracellular loop 2 (IL2) and stabilized the helical conformation of IL2. Moreover, LY3154207 binding induced subtle yet significant changes in key structural motifs and their neighboring residues. Notably, a cluster of residues centered around the Na + binding site became more compact, while interactions involving the PIF motif and its neighboring residues were loosened upon LY3154207 binding, consistent with their role in opening the intracellular crevice for receptor activation. Additionally, we identified an allosteric pathway likely responsible for the positive allosteric effect of LY3154207 in enhancing Gs protein coupling. This mechanistic understanding of LY3154207's allosteric action at the D1R pave the way for the rational design of more potent and effective allosteric modulators.

ALX/FPR2 Activation by Stereoisomers of D1 Resolvins Elucidating with Molecular Dynamics Simulation

Chronic inflammation contributes to several diseases, but its resolution is driven by specialized pro-resolving mediators (SPM) such as resolvin D1 (RvD1) and its epimer aspirin-triggered resolvin D1 (AT-RvD1), both biosynthesized from ω-3 fatty docosahexaenoic acid (DHA). RvD1 and AT-RvD1 have anti-inflammatory and pro-resolution potentials, and their effects could be mediated by formyl peptide receptor type 2 receptor ALX/FPR2, a G-protein-coupled receptor (GPCR). In this work, we performed 44 μs of molecular dynamics simulations with two complexes: FPR2@AT-RvD1 and FPR2@RvD1. Our results show the following: (i) in the AT-RvD1 simulations, the ALX/FPR2 receptor remained in the active state in 62% of the frames, while in the RVD1 simulations, the receptor remained in the active state in 74% of the frames; (ii) two residues, R201 and R205, of ALX/FPR2 appear, establishing interactions with both resolvins in all simulations (22 in total); (iii) RvD1 hydrogen bonds with R201 and R205 presented higher frequency than AT-RvD1; and (iv) residues R201 and R205 are the two receptor hotspots, demonstrated by the binding free calculations. Such results show that the ALX/FPR2 receptor remained in the active state for longer in the FPR2@RvD1 simulations than in the FPR2@AT-RvD1 simulations.

Mechanistic Elucidation of Activation/Deactivation Signal Transduction within Neurotensin Receptor 1 Triggered by 'Driver Chemical Groups' of Modulators: A Comparative Molecular Dynamics Simulation

Small-molecule modulators of neurotensin receptor 1 (NTSR1), a class A G-protein-coupled receptor (GPCR), has emerged as promising therapeutic agent for psychiatric disorders and cancer. Interestingly, a chemical group substitution in NTSR1 modulators can launch different types of downstream regulation, highlighting the significance of deciphering the internal fine-tuning mechanism. Here, we conducted a synergistic application of a Gaussian accelerated molecular dynamics simulation, a conventional molecular dynamics simulation, and Markov state models (MSM) to investigate the underlying mechanism of 'driver chemical groups' of modulators triggering inverse signaling. The results indicated that the flexibility of the leucine moiety in NTSR1 agonists contributes to the inward displacement of TM7 through a loosely coupled allosteric pathway, while the rigidity of the adamantane moiety in NTSR1 antagonists leads to unfavorable downward transduction of agonistic signaling. Furthermore, we found that R322^6.54, Y319^6.51, F353^7.42, R148^3.32, S356^7.45, and S357^7.46 may play a key role in inducing the activation of NTSR1. Together, our findings not only highlight the ingenious signal transduction within class A GPCRs but also lay a foundation for the development of targeted drugs harboring different regulatory functions of NTSR1.

The dynamics of agonist-β2-adrenergic receptor activation induced by binding of GDP-bound Gs protein

There is considerable uncertainty about the mechanism by which the β₂-adrenergic receptor (β₂AR) is activated. Here we use molecular metadynamics computations to predict the mechanism by which an agonist induces the activation of the β₂AR and its cognate Gs protein. We found that binding agonist alone to the inactive β₂AR does not break the ionic lock and hence does not drive the β₂AR towards the activated conformation. However, we found that attaching the inactive Gs protein to the agonist-bound inactive β₂AR (containing the ionic lock) leads to partial insertion of Gαs-α5 into the core of β₂AR, which breaks the ionic lock, leading to activation of the Gs protein coupled to β₂AR. Upon activation, the Gαs protein undergoes a remarkable opening of the GDP binding pocket, making the GDP available for exchange or release. Concomitantly, Gαs-α5 undergoes a remarkable expansion in the β₂AR cytoplasmic region after the ionic lock is broken, inducing TM6 to displace outward by ~5 Å from TM3.

A subtle structural modification of a synthetic cannabinoid receptor agonist drastically increases its efficacy at the CB1 receptor

The emergence of synthetic cannabinoid receptor agonists (SCRAs) as illicit psychoactive substances has posed considerable public health risks that include fatalities. Many SCRAs exhibit much higher efficacy and potency, compared with the phytocannabinoid Δ9-tetrahydrocannabinol (THC), at the cannabinoid receptor 1 (CB1R), a G protein-coupled receptor involved in modulating neurotransmitter release. In this study, we investigated structure activity relationships (SAR) of aminoalkylindole SCRAs at CB1Rs, focusing on 5F-pentylindoles containing an amide linker attached to different head moieties. Using in vitro bioluminescence resonance energy transfer (BRET) assays, we identified a few of SCRAs exhibiting significantly higher efficacy in engaging the Gi protein and recruiting β-arrestin than the reference CB1R full agonist CP55940. Importantly, adding a methyl group at the head moiety of 5F-MMB-PICA yielded 5F-MDMB-PICA, an agonist exhibiting a large increase in efficacy and potency at the CB1R. This pharmacological observation was supported by a functional assay of the effects of these SCRAs on glutamate field potentials recorded in hippocampal slices. Molecular modeling and simulations of the CB1R bound with either of the SCRAs revealed critical structural determinants contributing to the higher efficacy of 5F-MDMB-PICA, and how these subtle differences propagated to the receptor-G protein interface. Thus, we find that apparently minor structural changes in the head moiety of SCRAs can cause major changes in efficacy. Our results highlight the need for close monitoring of structural modifications of newly emerging SCRAs and their potential for toxic drug responses in humans.

Stabilization of pre-existing neurotensin receptor conformational states by β-arrestin-1 and the biased allosteric modulator ML314

The neurotensin receptor 1 (NTS1) is a G protein-coupled receptor (GPCR) with promise as a drug target for the treatment of pain, schizophrenia, obesity, addiction, and various cancers. A detailed picture of the NTS1 structural landscape has been established by X-ray crystallography and cryo-EM and yet, the molecular determinants for why a receptor couples to G protein versus arrestin transducers remain poorly defined. We used 13CεH3-methionine NMR spectroscopy to show that binding of phosphatidylinositol-4,5-bisphosphate (PIP2) to the receptor's intracellular surface allosterically tunes the timescale of motions at the orthosteric pocket and conserved activation motifs - without dramatically altering the structural ensemble. β-arrestin-1 further remodels the receptor ensemble by reducing conformational exchange kinetics for a subset of resonances, whereas G protein coupling has little to no effect on exchange rates. A β-arrestin biased allosteric modulator transforms the NTS1:G protein complex into a concatenation of substates, without triggering transducer dissociation, suggesting that it may function by stabilizing signaling incompetent G protein conformations such as the non-canonical state. Together, our work demonstrates the importance of kinetic information to a complete picture of the GPCR activation landscape.

Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics

Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M₄ muscarinic acetylcholine receptor (M₄ mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M₄ mAChR. By determining the cryo-electron microscopy structures of the M₄ mAChR, bound to a cognate G_i1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand-receptor-transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.

General Metadynamics Protocol To Simulate Activation/Deactivation of Class A GPCRs: Proof of Principle for the Serotonin Receptor

We present a generally applicable metadynamics protocol for characterizing the activation free-energy profiles of class A G-protein coupled receptors and a proof-of-principle study for the 5HT_1A-receptor. The almost universal A¹⁰⁰ activation index, which depends on five inter-helix distances, is used as the single collective variable in well-tempered multiple-walker metadynamics simulations. Here, we show free-energy profiles for the serotonin receptor as binary (apo-receptor + G-protein-α-subunit and receptor + ligand) and ternary complexes with two prototypical orthosteric ligands: the full agonist serotonin and the partial agonist aripiprazole. Our results are not only compatible with previously reported experimental and computational data, but they also allow differences between active and inactive conformations to be determined in unprecedented atomic detail, and with respect to the so-called microswitches that have been suggested as determinants of activation, giving insight into their role in the activation mechanism.

Distinct activation mechanisms regulate subtype selectivity of Cannabinoid receptors

Design of cannabinergic subtype selective ligands is challenging because of high sequence and structural similarities of cannabinoid receptors (CB₁ and CB₂). We hypothesize that the subtype selectivity of designed selective ligands can be explained by the ligand binding to the conformationally distinct states between cannabinoid receptors. Analysis of ~ 700 μs of unbiased simulations using Markov state models and VAMPnets identifies the similarities and distinctions between the activation mechanism of both receptors. Structural and dynamic comparisons of metastable intermediate states allow us to observe the distinction in the binding pocket volume change during CB₁ and CB₂ activation. Docking analysis reveals that only a few of the intermediate metastable states of CB₁ show high affinity towards CB₂ selective agonists. In contrast, all the CB₂ metastable states show a similar affinity for these agonists. These results mechanistically explain the subtype selectivity of these agonists by deciphering the activation mechanism of cannabinoid receptors.

Computational insights into ligand-induced G protein and β-arrestin signaling of the dopamine D1 receptor

The dopamine D1 receptor (D1R), is a class A G protein coupled-receptor (GPCR) which has been a promising drug target for psychiatric and neurological disorders such as Parkinson's disease (PD). Previous studies have suggested that therapeutic effects can be realized by targeting the β-arrestin signaling pathway of dopamine receptors, while overactivation of the G protein-dependent pathways leads to side effects, such as dyskinesias. Therefore, it is highly desirable to develop a D1R ligand that selectively regulates the β-arrestin pathway. Currently, most D1R agonists are signaling-balanced and stimulate both G protein and β-arrestin pathways, with a few reports of G protein biased ligands. However, identification and characterization of β-arrestin biased D1R agonists has been a challenge thus far. In this study, we implemented Gaussian accelerated molecular dynamics (GaMD) simulations to provide valuable computational insights into the possible underlying molecular mechanism of the different signaling properties of two catechol and two non-catechol D1R agonists that are either G protein biased or signaling-balanced. Dynamic network analysis further identified critical residues in the allosteric signaling network of D1R for each ligand at different conformational or binding states. Some of these residues are crucial for G protein or arrestin signals of GPCRs based on previous studies. Finally, we provided a molecular design strategy which can be utilized by medicinal chemists to develop potential β-arrestin biased D1R ligands. The proposed hypotheses are experimentally testable and can guide the development of safer and more effective medications for a variety of CNS disorders.

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

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-β2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and β2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Molecular Dynamics Simulation and Pharmacoinformatic Integrated Analysis of Bioactive Phytochemicals from Azadirachta indica (Neem) to Treat Diabetes Mellitus

Diabetes mellitus is a chronic hormonal and metabolic disorder in which our body cannot generate necessary insulin or does not act in response to it, accordingly, ensuing in discordantly high blood sugar (glucose) levels. Diabetes mellitus can lead to systemic dysfunction in the multiorgan system, including cardiac dysfunction, severe kidney disease, lowered quality of life, and increased mortality risk from diabetic complications. To uncover possible therapeutic targets to treat diabetes mellitus, the in silico drug design technique is widely used, which connects the ligand molecules with target proteins to construct a protein-ligand network. To identify new therapeutic targets for type 2 diabetes mellitus, Azadirachta indica is subjected to phytochemical screening using in silico molecular docking, pharmacokinetic behavior analysis, and simulation-based molecular dynamic analysis. This study has analyzed around 63 phytochemical compounds, and the initial selection of the compounds was made by analyzing their pharmacokinetic properties by comparing them with Lipinski’s rule of 5. The selected compounds were subjected to molecular docking. The top four ligand compounds were reported along with the control drug nateglinide based on their highest negative molecular binding affinity. The protein-ligand interaction of selected compounds has been analyzed to understand better how compounds interact with the targeted protein structure. The results of the in silico analysis revealed that 7-Deacetyl-7-oxogedunin had the highest negative docking score of −8.9 Kcal/mol and also demonstrated standard stability in a 100 ns molecular dynamic simulation performed with insulin receptor ectodomain. It has been found that these substances may rank among the essential supplementary antidiabetic drugs for treating type 2 diabetes mellitus. It is suggested that more in vivo and in vitro research studies be carried out to support the conclusions drawn from this in silico research strategy.