The mechanism for ligand activation of the GPCR-G protein complex

Growth inhibitory peptides: a potential novel therapeutic approach to cancer treatment

Cancer remains a major global public health concern, with considerable interest in exploring biological molecules for cancer treatment and prevention. Growth inhibitory peptide (GIP), a promising new class of biological therapeutics, has drawn attention for its distinct anti-tumor properties. Derived from human alpha-fetoprotein (HAFP), this synthetic 34-amino-acid peptide has demonstrated substantial anti-tumor effects across various cancer cell lines, effectively inhibiting tumor cell proliferation, migration, and metastasis. Studies reveal that GIP mediates its effects through a range of mechanisms, including interactions with G protein-coupled receptors (GPCRs), anti-cell adhesion activities, inhibition of cell spreading and metastatic processes, morphological alterations, platelet aggregation inhibition, immune enhancement, cell membrane disruption, ion channel blockade, and cell cycle arrest. While GIP has exhibited promising anti-tumor activity in both in vitro and in vivo models, further investigation is essential to advance its development as a therapeutic drug, particularly regarding pharmacokinetics, safety profiles, storage stability, and clinical efficacy.

Atomistic Mechanism Underlying the Regulation of the GPA1 G Protein Signaling Pathway Mediated by the Gibberellin A1 Phytohormone Binding to the GCR1 Plant G-Protein-Coupled Receptor

We propose an atomistic mechanism suggesting that fundamental plant processes, including seed germination, root elongation, and flower and fruit production, may be regulated by phytohormones such as Gibberellin A1 (GA1) binding to the GCR1 plant G-protein-coupled receptor. This parallels the central roles of GPCRs in animals for vision, taste, smell, pain, depression, and nerve signaling, among others. Validating GCR1 as a genuine GPCR in plants, particularly its interaction with GPA1, G-protein, would mark a groundbreaking advancement in understanding plant processes, both biologically and agronomically. However, experimental confirmation of this interaction and evidence supporting the idea that binding of GA1 to GCR1 would regulate GPA1 activation are lacking. Indeed, the design of experiments to explore these hypotheses is impeded by the absence of structural information relating to interactions of GPA1 with the GA1-GCR1 complex. To address this gap, we employ molecular dynamics and metadynamics simulations to demonstrate that binding GPA1 to the GA1-GCR1 complex induces conformational changes that open up the Ras and Helical domains of GPA1 to release GDP for exchange with GTP, thereby enabling signaling. Our results suggest numerous mutations involving GA1 binding at the GCR1 site and the coupling of GCR1 to the GPA1 G-protein that could be used to validate (or not validate) our predicted mechanism. Such validation would serve as a foundation for devising strategies to design novel agonists and inverse agonists to provide precise control of crucial plant processes.

The regulatory functions of G protein-coupled receptors signaling pathways in B cell differentiation and development contributing to autoimmune diseases

Autoimmune diseases are characterized by a dysfunction of the immune system. Disruptions in the balance of B-cell dynamics and the increase in auto-antibody levels are pivotal in the triggering of several autoimmune disorders. All of this is inextricably linked to the differentiation, development, migration, and functional regulation of B cells in the human immune response. G protein-coupled receptors (GPCR) are recognized as crucial targets in drug development and play pivotal roles in both B cell differentiation and the underlying mechanisms of autoimmune diseases. However, there has been an inadequate comprehension of how GPCR intricately modulate B cell development and impact the pathogenesis of autoimmune diseases. Ligands and functions of GPCR-chemokine receptors including CXCR3, CXCR4, CXCR5 and CCR7, lipid receptors including S1PR1-5, cannabinoid receptor CB2 as well as orphan GPCR including GPR132, GPR183, GPR174, and P2RY8 in B cell differentiation and development, will be elaborated in this review. The roles these GPCR play in mediating B cells in several autoimmune diseases will also be discussed. The elucidation of the multifaceted mechanisms controlled by GPCR not only enriches our comprehension of immune responses but also provides a promising avenue for therapeutic interventions in the domain of autoimmune disorders.

Identification of allosteric sites and ligand-induced modulation in the dopamine receptor through large-scale alchemical mutation scan

G protein coupled receptors, particularly class A GPCRs are arguably the most important class of membrane receptors and preferred targets for drug development. Despite extensive research on how ligands modulate the receptor response, discovering new, highly specific ligands remains challenging. However, finding residues outside the conserved microswitches that affect the active-inactive state equilibrium and are specific for a certain receptor, can be beneficial for the design of ligands with higher receptor selectivity. Focusing on the human dopamine receptor 2 (DRD2), we uncover crucial residues for the activation modulation using alchemical non-equilibrium free energy calculations. Our findings match with literature on activation microswitches and experimental studies, while also uncovering novel important residues. Further, we analyzed mutation-induced changes in residue contact networks and found that modulating these networks can lead to a stabilization of the respective opposite state, an effect that could as well be achieved by well-engineered (small) ligands. This way we provide insights into the mechanism of action of the well-known drugs risperidone and bromocriptine and showcase on these two examples how our data can be used for the design of new ligands.

Influence of expression and purification protocols on Gα biochemical activity: kinetics of plant and mammalian G protein cycles

Heterotrimeric G proteins, composed of Gα, Gβ, and Gγ subunits, are a class of signal transduction complexes with broad roles in human health and agriculturally relevant plant physiological and developmental traits. In the classic paradigm, guanine nucleotide binding to the Gα subunit regulates the activation status of the complex. We sought to develop improved methods for heterologous expression and rapid purification of Gα subunits, initially targeting GPA1, the sole canonical Gα subunit of the model plant species, Arabidopsis thaliana. Compared to conventional methods, our expression methodology and rapid StrepII-tag mediated purification facilitates substantially higher yield, and isolation of protein with increased GTP binding and hydrolysis activities. Human GNAI1 purified using our approach displayed the expected binding and hydrolysis activities, indicating our protocol is applicable to mammalian Gα subunits, potentially including those for which purification of enzymatically active protein has been historically problematic. We subsequently utilized domain swaps of GPA1 and human GNAO1 to demonstrate that the inherent instability of GPA1 is a function of the interaction between the Ras and helical domains. Additionally, we found that GPA1-GNAO1 domain swaps partially uncouple the instability from the rapid nucleotide binding kinetics displayed by GPA1. In summary, our work provides insights into methods to optimally study heterotrimeric G proteins, and reveals roles of the helical domain in Gα kinetics and stability.

Functional dynamics of G protein-coupled receptors reveal new routes for drug discovery

G protein-coupled receptors (GPCRs) are the largest human membrane protein family that transduce extracellular signals into cellular responses. They are major pharmacological targets, with approximately 26% of marketed drugs targeting GPCRs, primarily at their orthosteric binding site. Despite their prominence, predicting the pharmacological effects of novel GPCR-targeting drugs remains challenging due to the complex functional dynamics of these receptors. Recent advances in X-ray crystallography, cryo-electron microscopy, spectroscopic techniques and molecular simulations have enhanced our understanding of receptor conformational dynamics and ligand interactions with GPCRs. These developments have revealed novel ligand-binding modes, mechanisms of action and druggable pockets. In this Review, we highlight such aspects for recently discovered small-molecule drugs and drug candidates targeting GPCRs, focusing on three categories: allosteric modulators, biased ligands, and bivalent and bitopic compounds. Although studies so far have largely been retrospective, integrating structural data on ligand-induced receptor functional dynamics into the drug discovery pipeline has the potential to guide the identification of drug candidates with specific abilities to modulate GPCR interactions with intracellular effector proteins such as G proteins and β-arrestins, enabling more tailored selectivity and efficacy profiles.

AlphaFold3 versus experimental structures: assessment of the accuracy in ligand-bound G protein-coupled receptors

G protein-coupled receptors (GPCRs) are critical drug targets involved in numerous physiological processes, yet many of their structures remain unresolved due to inherent flexibility and diverse ligand interactions. This study systematically evaluates the accuracy of AlphaFold3-predicted GPCR structures compared to experimentally determined structures, with a primary focus on ligand-bound states. Our analysis reveals that while AlphaFold3 shows improved performance over AlphaFold2 in predicting overall GPCR backbone architecture, significant discrepancies persist in ligand-binding poses, particularly for ions, peptides, and proteins. Despite advancements, these limitations constrain the utility of AlphaFold3 models in functional studies and structure-based drug design, where high-resolution details of ligand interactions are crucial. We assess the accuracy of predicted structures across various ligand types, quantifying deviations in binding pocket geometries and ligand orientations. Our findings highlight specific challenges in the computational prediction of ligand-bound GPCR structures, emphasizing areas where further refinement is needed. This study provides valuable insights for researchers using AlphaFold3 in GPCR studies, underscores the ongoing necessity for experimental structure determination, and offers direction for improving protein-ligand interaction predictions in future computational models.

CXC Chemokine Ligand 12 Facilitates Gi Protein Binding to CXC Chemokine Receptor 4 by Stabilizing Packing of the Proline-Isoleucine-Phenylalanine Motif: Insights from Automated Path Searching

The activation of G protein-coupled receptors (GPCRs) is a complex multibody multievent process involving agonist binding, receptor activation, G protein coupling, and subsequent G protein activation. The order and energetics of these events, though crucial for the rational design of selective GPCR drugs, are challenging to characterize and remain largely underexplored. Here, we employed molecular dynamics simulations and our recently developed traveling salesman-based automated path searching (TAPS) algorithm to efficiently locate the minimum free-energy paths for the coupling events of the CXC chemokine receptor 4 (CXCR4) with its endogenous ligand CXC chemokine ligand 12 (CXCL12) and Gi protein. We show that, after overcoming three low energy barriers (3.24-6.89 kcal/mol), Gi alone can precouple with CXCR4 even without CXCL12, consistent with previous reports on the existence of the apo CXCR4-Gi complex and our NanoBiT experiments. The highest barrier of 6.89 kcal/mol in this process corresponds to the packing of the proline-isoleucine-phenylalanine (PIF) motif of CXCR4. Interestingly, without Gi, CXCL12 alone cannot activate CXCR4 (high barrier of 18.89 kcal/mol). Instead, it can enhance Gi coupling by circumventing the energy barrier of PIF packing. Shedding new light on the activation mechanism of CXCR4, these results present TAPS as a promising tool for uncovering complete activation pathways of GPCRs and the corresponding agonist design.

Ligand-Induced Biased Activation of GPCRs: Recent Advances and New Directions from In Silico Approaches

G-protein coupled receptors (GPCRs) are the largest family of membrane proteins engaged in transducing signals from the extracellular environment into the cell. GPCR-biased signaling occurs when two different ligands, sharing the same binding site, induce distinct signaling pathways. This selective signaling offers significant potential for the design of safer and more effective drugs. Although its molecular mechanism remains elusive, big efforts are made to try to explain this mechanism using a wide range of methods. Recent advances in computational techniques and AI technology have introduced a variety of simulations and machine learning tools that facilitate the modeling of GPCR signal transmission and the analysis of ligand-induced biased signaling. In this review, we present the current state of in silico approaches to elucidate the structural mechanism of GPCR-biased signaling. This includes molecular dynamics simulations that capture the main interactions causing the bias. We also highlight the major contributions and impacts of transmembrane domains, loops, and mutations in mediating biased signaling. Moreover, we discuss the impact of machine learning models on bias prediction and diffusion-based generative AI to design biased ligands. Ultimately, this review addresses the future directions for studying the biased signaling problem through AI approaches.

Structural Insight into the Inactive/Active States of 5-HT1AR and Molecular Mechanisms of Electric Fields in Modulating 5-HT1AR

Probing the differences between inactive/active states of the serotonin 1A receptor (5-HT1AR) and the dynamic receptor conformations is vital for understanding signaling transduction pathways and diverse physiological responses. Here, we compared the conformational features between the inactive and active states of 5-HT1AR and explored the role of serotonin in the activation process of 5-HT1AR by using molecular dynamics (MD) simulations. The results show that the position of TM6 and the arrangements of key motifs exhibit distinctions in the inactive and active states of 5-HT1AR. The binding of serotonin to 5-HT1AR is mostly driven by hydrophobic, aromatic stacking, anion-π, and H-bonding interactions. We also performed additional MD simulations with electric fields (EFs) of 0.01 and 0.03 V/nm to investigate the effects of EFs on the conformation of the 5-HT1AR-serotonin complex. The conformational change of 5-HT1AR and the inward movement of TM6 are increased with the field strength, indicative of a dependence on the strength of the EF. The EF of 0.03 V/nm affects the binding behaviors of serotonin with 5-HT1AR and further disturbs the activation of 5-HT1AR by serotonin. This study first reveals atomic-level information about the distinct features between inactive and active states of 5-HT1AR and demonstrates the pivotal role of EF in modulating the 5-HT1AR-ligand complex.

Ligand-Independent Spontaneous Activation of Purinergic P2Y6 Receptor Under Cell Culture Soft Substrate

G protein-coupled receptors (GPCRs) exist in the conformational equilibrium between inactive state and active state, where the proportion of active state in the absence of a ligand determines the basal activity of GPCRs. Although many GPCRs have different basal activity, it is still unclear whether physiological stresses such as substrate stiffness affect the basal activity of GPCRs. In this study, we identified that purinergic P2Y6 receptor (P2Y6R) induced spontaneous Ca2+ oscillation without a nucleotide ligand when cells were cultured in a silicon chamber. This P2Y6R-dependent Ca2+ oscillation was absent in cells cultured in glass dishes. Coating substrates, including collagen, laminin, and fibronectin, did not affect the P2Y6R spontaneous activity. Mutation of the extracellular Arg-Gly-Asp (RGD) motif of P2Y6R inhibited spontaneous activity. Additionally, extracellular Ca2+ was required for P2Y6R-dependent spontaneous Ca2+ oscillation. The GPCR screening assay identified cells expressing 10 GPCRs, including purinergic P2Y1R, P2Y2R, and P2Y6R, that exhibited spontaneous Ca2+ oscillation under cell culture soft substrate. Our results suggest that stiffness of the cell adhesion surface modulates spontaneous activities of several GPCRs, including P2Y6R, through a ligand-independent mechanism.

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.

Metabotropic GABAB Receptor Activation Induced by G Protein Coupling

G protein-coupled receptors (GPCRs) play central roles in regulating cellular responses through heterotrimeric G proteins (GP). Extensive studies have elucidated the complex cellular signaling mediated by GPCRs that accompany dynamic conformational changes upon activation. However, there has been less focus on the role of the GP on the activation process, particularly for class C GPCRs that function as obligate dimers. Herein, we report the pivotal role of GP coupling on the dynamic activation process for the metabotropic γ-aminobutyric acid receptor (GABABR) based on extensive atomistic simulations. We find that GP coupling triggers drastic conformational changes in the GABABR transmembrane domain (TMD), while an agonist alone is insufficient to shift the equilibrium state from the inactive to the active states. These conformational changes induced by GP coupling destabilize the inactive TM5/TM5 interface, shifting the equilibrium toward the activated TM6/TM6 interface. This active role of the GP in activation provides fresh insights into the activation mechanism of GABABR and perhaps other class C GPCRs. These insights should aid in the development of more potent and selective drugs.

Fusarium graminearum Ste2 and Ste3 receptors undergo peroxidase-induced heterodimerization when expressed heterologously in Saccharomyces cerevisiae

Fusarium graminearum FgSte2 and FgSte3 are G-protein-coupled receptors (GPCRs) shown to play roles in hyphal chemotropism and fungal plant pathogenesis in response to activity arising from host-secreted peroxidases. Here, we follow up on the observation that chemotropism is dependent on both FgSte2 and FgSte3 being present; testing the possibility that this might be due to formation of an FgSte2-FgSte3 heterodimer. Bioluminescence resonance energy transfer (BRET) analyses were conducted in Saccharomyces cerevisiae, where the addition of horse radish peroxidase (HRP) was found to increase the transfer of energy from the inducibly expressed FgSte3-Nano luciferase donor, to the constitutively expressed FgSte2-yellow fluorescent protein (YFP) acceptor, compared to controls. A partial response was also detected when an HRP-derived ligand-containing extract was enriched from F. graminearum spores and applied instead of HRP. In contrast, substitution with pheromones or an unrelated bovine GPCR, rhodopsin-YFP used as acceptor, eliminated all BRET responses. Interaction results were validated by affinity pulldown and receptor expression was validated by confocal immunofluorescence microscopy. Taken together these findings demonstrate the formation of HRP and HRP-derived ligand stimulated heterodimers between FgSte2 and FgSte3. Outcomes are discussed from the context of the roles of ligands and reactive oxygen species in GPCR dimerization.

Charge Movements and Conformational Changes: Biophysical Properties and Physiology of Voltage-Dependent GPCRs

G protein-coupled receptors (GPCRs) regulate multiple cellular functions and represent important drug targets. More than 20 years ago, it was noted that GPCR activation (agonist binding) and signaling (G protein activation) are dependent on the membrane potential (VM). While it is now proven that many GPCRs display an intrinsic voltage dependence, the molecular mechanisms of how GPCRs sense depolarization of the plasma membrane are less well defined. This review summarizes the current knowledge of voltage-dependent signaling in GPCRs. We describe how voltage dependence was discovered in muscarinic receptors, present an overview of GPCRs that are regulated by voltage, and show how biophysical properties of GPCRs led to the discovery of voltage-sensing mechanisms in those receptors. Furthermore, we summarize physiological functions that have been shown to be regulated by voltage-dependent GPCR signaling of endogenous receptors in excitable tissues, such as the nervous system or the heart. Finally, we discuss challenges that remain in analyzing voltage-dependent signaling of GPCRs in vivo and present an outlook on experimental applications of the interesting concept of GPCR signaling.

Minute-timescale free-energy calculations reveal a pseudo-active state in the adenosine A2A receptor activation mechanism

G protein-coupled receptors (GPCRs) are membrane proteins targeted by over one-third of marketed drugs. Understanding their activation mechanism is essential for precise regulation of drug pharmacological response. In this work, we elucidate the conformational landscape of the adenosine A2A receptor (A2AR) activation mechanism in its basal apo form and under different ligand-bound conditions through minute-timescale free-energy calculations. We identified a pseudo-active state (pAs) of the A2AR apo form, stabilized by specific "microswitch" residues, including a salt bridge established between the conserved residues R5.66 and E6.30. The pAs enables A2AR to couple with Gs protein upon rearrangement of the intracellular end of transmembrane helix 6, providing unprecedented structural insights into receptor function and signaling dynamics. Our simulation protocol is versatile and can be adapted to study the activation of any GPCRs, potentially making it a valuable tool for drug design and "biased signaling" studies.

Agonist activation to open the Gα subunit of the GPCR-G protein precoupled complex defines functional agonist activation of TAS2R5

G protein-coupled receptors (GPCRs) regulate multiple cellular responses and represent highly successful therapeutic targets. The mechanisms by which agonists activate the G protein are unclear for many GPCR families, including the bitter taste receptors (TAS2Rs). We ascertained TAS2R5 properties by live cell-based functional assays, direct binding affinity measurements using optical resonators, and atomistic molecular dynamics simulations. We focus on three agonists that exhibit a wide range of signal transduction in cells despite comparable ligand-receptor binding energies derived from direct experiment and computation. Metadynamics simulations revealed that the critical barrier to activation is ligand-induced opening of the G protein between the α-helical (AH) and Ras-like domains of Gα subunit from a precoupled TAS2R5-G protein state to the fully activated state. A moderate agonist opens the AH-Ras cleft from 22 Å to 31 Å with an energy gain of -4.8 kcal mol-1, making GDP water-exposed for signaling. A high-potency agonist had an energy gain of -11.1 kcal mol-1. The low-potency agonist is also exothermic for Gα opening, but with an energy gain of only -1.4 kcal mol-1. This demonstrates that TAS2R5 agonist-bound functional potencies are derived from energy gains in the transition from a precoupled complex at the level of Gα opening. Our experimental and computational study provides insights into the activation mechanism of signal transduction that provide a basis for rational design of new drugs.

GPCR-BSD: a database of binding sites of human G-protein coupled receptors under diverse states

G-protein coupled receptors (GPCRs), the largest family of membrane proteins in human body, involve a great variety of biological processes and thus have become highly valuable drug targets. By binding with ligands (e.g., drugs), GPCRs switch between active and inactive conformational states, thereby performing functions such as signal transmission. The changes in binding pockets under different states are important for a better understanding of drug-target interactions. Therefore it is critical, as well as a practical need, to obtain binding sites in human GPCR structures. We report a database (called GPCR-BSD) that collects 127,990 predicted binding sites of 803 GPCRs under active and inactive states (thus 1,606 structures in total). The binding sites were identified from the predicted GPCR structures by executing three geometric-based pocket prediction methods, fpocket, CavityPlus and GHECOM. The server provides query, visualization, and comparison of the predicted binding sites for both GPCR predicted and experimentally determined structures recorded in PDB. We evaluated the identified pockets of 132 experimentally determined human GPCR structures in terms of pocket residue coverage, pocket center distance and redocking accuracy. The evaluation showed that fpocket and CavityPlus methods performed better and successfully predicted orthosteric binding sites in over 60% of the 132 experimentally determined structures. The GPCR Binding Site database is freely accessible at https://gpcrbs.bigdata.jcmsc.cn . This study not only provides a systematic evaluation of the commonly-used fpocket and CavityPlus methods for the first time but also meets the need for binding site information in GPCR studies.

Steviol rebaudiosides bind to four different sites of the human sweet taste receptor (T1R2/T1R3) complex explaining confusing experiments

Sucrose provides both sweetness and energy by binding to both Venus flytrap domains (VFD) of the heterodimeric sweet taste receptor (T1R2/T1R3). In contrast, non-caloric sweeteners such as sucralose and aspartame only bind to one specific domain (VFD2) of T1R2, resulting in high-intensity sweetness. In this study, we investigate the binding mechanism of various steviol glycosides, artificial sweeteners, and a negative allosteric modulator (lactisole) at four distinct binding sites: VFD2, VFD3, transmembrane domain 2 (TMD2), and TMD3 through binding experiments and computational docking studies. Our docking results reveal multiple binding sites for the tested ligands, including the radiolabeled ligands. Our experimental evidence demonstrates that the C20 carboxy terminus of the Gα protein can bind to the intracellular region of either TMD2 or TMD3, altering GPCR affinity to the high-affinity state for steviol glycosides. These findings provide a mechanistic understanding of the structure and function of this heterodimeric sweet taste receptor.

The G Protein-First Mechanism for Activation of the Class B Glucagon-like Peptide 1 Receptor Coupled to N-Terminal Domain-Mediated Conformational Progression

Recently, there has been a great deal of excitement about new glucagon-like peptide 1 receptor (GLP-1R) agonists (e.g., semaglutide and tirzepatide) that have received FDA approval for type 2 diabetes and obesity. Although effective, these drugs come with side effects that limit their use. While research efforts continue to focus intensively on long-lasting, orally administered GLP-1R medications with fewer side effects, a major impediment to developing improved GLP-1R medications is that the mechanism by which an agonist activates GLP-1R to imitate signaling is not known. Here we present and validate the G protein (GP)-first mechanism for the GLP-1R supported by extensive atomistic simulations. We propose that GLP-1R is preactivated through the formation of a GLP-1R-GP precoupled complex at the cell membrane prior to ligand binding. Despite a transmembrane helix 6 (TM6)-bentout conformation characteristic of activated GLP-1R, this precoupled complex remains unactivated until an agonist binds to elicit signaling. Notably, this new hypothesis offers a unified and predictive model for the activities of a series of full and partial agonists, including the peptides ExP5, GLP-1(7-36), and GLP-1(9-36). Most surprisingly, our simulations reveal an N-terminus domain (NTD)-swing/agonist-insertion mechanism wherein the long extracellular NTD of GLP-1R tightly holds the C-terminal half of the peptide agonist and progressively shifts the N-terminal head of the peptide to facilitate insertion into the orthosteric pocket. Our findings provide novel mechanistic insights into the activation and function of class B GPCRs and should provide a realistic basis for structure-based ligand design.

Effects of an external static EF on the conformational transition of 5-HT1A receptor: A molecular dynamics simulation study

The serotonin receptor subtype 1A (5-HT1AR), one of the G-protein-coupled receptor (GPCR) family, has been implicated in several neurological conditions. Understanding the activation and inactivation mechanism of 5-HT1AR at the molecular level is critical for discovering novel therapeutics in many diseases. Recently there has been a growing appreciation for the role of external electric fields (EFs) in influencing the structure and activity of biomolecules. In this study, we used molecular dynamics (MD) simulations to examine conformational features of active states of 5-HT1AR and investigate the effect of an external static EF with 0.02 V/nm applied on the active state of 5-HT1AR. Our results showed that the active state of 5-HT1AR maintained the native structure, while the EF led to structural modifications in 5-HT1AR, particularly inducing the inward movement of transmembrane helix 6 (TM6). Furthermore, it disturbed the conformational switches associated with activation in the CWxP, DRY, PIF, and NPxxY motifs, consequently predisposing an inclination towards the inactive-like conformation. We also found that the EF led to an overall increase in the dipole moment of 5-HT1AR, encompassing TM6 and pivotal amino acids. The analyses of conformational properties of TM6 showed that the changed secondary structure and decreased solvent exposure occurred upon the EF condition. The interaction of 5-HT1AR with the membrane lipid bilayer was also altered under the EF. Our findings reveal the molecular mechanism underlying the transition of 5-HT1AR conformation induced by external EFs, which offer potential novel insights into the prospect of employing structure-based EF applications for GPCRs.

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.

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.

Descriptive molecular pharmacology of the δ opioid receptor (DOR): A computational study with structural approach

This work focuses on the δ receptor (DOR), a G protein-coupled receptor (GPCR) belonging to the opioid receptor group. DOR is expressed in numerous tissues, particularly within the nervous system. Our study explores computationally the receptor's interactions with various ligands, including opiates and opioid peptides. It elucidates how these interactions influence the δ receptor response, relevant in a wide range of health and pathological processes. Thus, our investigation aims to explore the significance of DOR as an incoming drug target for pain relief and neurodegenerative diseases and as a source for novel opioid non-narcotic analgesic alternatives. We analyze the receptor's structural properties and interactions using Molecular Dynamics (MD) simulations and Gaussian-accelerated MD across different functional states. To thoroughly assess the primary differences in the structural and conformational ensembles across our different simulated systems, we initiated our study with 1 μs of conventional Molecular Dynamics. The strategy was chosen to encompass the full activation cycle of GPCRs, as activation processes typically occur within this microsecond range. Following the cMD, we extended our study with an additional 100 ns of Gaussian accelerated Molecular Dynamics (GaMD) to enhance the sampling of conformational states. This simulation approach allowed us to capture a comprehensive range of dynamic interactions and conformational changes that are crucial for GPCR activation as influenced by different ligands. Our study includes comparing agonist and antagonist complexes to uncover the collective patterns of their functional states, regarding activation, blocking, and inactivation of DOR, starting from experimental data. In addition, we also explored interactions between agonist and antagonist molecules from opiate and opioid classifications to establish robust structure-activity relationships. These interactions have been systematically quantified using a Quantitative Structure-Activity Relationships (QSAR) model. This research significantly contributes to our understanding of this significant pharmacological target, which is emerging as an attractive subject for drug development.

Relevance of G protein-coupled receptor (GPCR) dynamics for receptor activation, signalling bias and allosteric modulation

G protein-coupled receptors (GPCRs) are one of the major drug targets. In recent years, computational drug design for GPCRs has mainly focused on static structures obtained through X-ray crystallography, cryogenic electron microscopy (cryo-EM) or in silico modelling as a starting point for virtual screening campaigns. However, GPCRs are highly flexible entities with the ability to adopt different conformational states that elicit different physiological responses. Including this knowledge in the drug discovery pipeline can help to tailor novel conformation-specific drugs with an improved therapeutic profile. In this review, we outline our current knowledge about GPCR dynamics that is relevant for receptor activation, signalling bias and allosteric modulation. Ultimately, we highlight new technological implementations such as time-resolved X-ray crystallography and cryo-EM as well as computational algorithms that can contribute to a more comprehensive understanding of receptor dynamics and its relevance for GPCR functionality.

TMC6 functions as a GPCR-like receptor to sense noxious heat via Gαq signaling

Thermosensation is vital for the survival, propagation, and adaption of all organisms, but its mechanism is not fully understood yet. Here, we find that TMC6, a membrane protein of unknown function, is highly expressed in dorsal root ganglion (DRG) neurons and functions as a Gαq-coupled G protein-coupled receptor (GPCR)-like receptor to sense noxious heat. TMC6-deficient mice display a substantial impairment in noxious heat sensation while maintaining normal perception of cold, warmth, touch, and mechanical pain. Further studies show that TMC6 interacts with Gαq via its intracellular C-terminal region spanning Ser780 to Pro810. Specifically disrupting such interaction using polypeptide in DRG neurons, genetically ablating Gαq, or pharmacologically blocking Gαq-coupled GPCR signaling can replicate the phenotype of TMC6 deficient mice regarding noxious heat sensation. Noxious heat stimulation triggers intracellular calcium release from the endoplasmic reticulum (ER) of TMC6- but not control vector-transfected HEK293T cell, which can be significantly inhibited by blocking PLC or IP3R. Consistently, noxious heat-induced intracellular Ca2+ release from ER and action potentials of DRG neurons largely reduced when ablating TMC6 or blocking Gαq/PLC/IP3R signaling pathway as well. In summary, our findings indicate that TMC6 can directly function as a Gαq-coupled GPCR-like receptor sensing noxious heat.

Influence of expression and purification protocols on Gα biochemical activity: kinetics of plant and mammalian G protein cycles

Heterotrimeric G proteins are a class of signal transduction complexes with broad roles in human health and agriculturally important plant traits. In the classic paradigm, guanine nucleotide binding to the Gα subunit regulates the activation status of the complex. Using the Arabidopsis thaliana Gα subunit, GPA1, we developed a rapid StrepII-tag mediated purification method that facilitates isolation of protein with increased enzymatic activities as compared to conventional methods, and is demonstrably also applicable to mammalian Gα subunits. We subsequently utilized domain swaps of GPA1 and human GNAO1 to demonstrate the instability of recombinant GPA1 is a function of the interaction between the Ras and helical domains, and can be partially uncoupled from the rapid nucleotide binding kinetics displayed by GPA1.

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.

GPR41 and GPR43: From development to metabolic regulation

G-protein-coupled receptors are a diverse class of cell surface receptors that orchestrate numerous physiological functions. The G-protein-coupled receptors, GPR41 and GPR43, sense short-chain fatty acids (SCFAs), which are metabolites of dietary fermentation by the host's intestinal bacteria. These receptors have gained attention as potential therapeutic targets against various diseases because of their SCFA-mediated beneficial effects on the host's intestinal health. Mounting evidence has associated the activity of these receptors with chronic metabolic diseases, including obesity, diabetes, inflammation, and cardiovascular disease. However, despite intensive research using various strategies, including gene knockout (KO) mouse models, evidence about the precise roles of GPR41 and GPR43 in disease treatment remains inconsistent. Here, we comprehensively review the latest findings from functional studies of the signaling mechanisms that underlie the activities of GPR41 and GPR43, as well as highlight their multifaceted roles in health and disease. We anticipate that this knowledge will guide future research priorities and the development of effective therapeutic interventions.

A pattern recognition artificial olfactory system based on human olfactory receptors and organic synaptic devices

Neuromorphic sensors, designed to emulate natural sensory systems, hold the promise of revolutionizing data extraction by facilitating rapid and energy-efficient analysis of extensive datasets. However, a challenge lies in accurately distinguishing specific analytes within mixtures of chemically similar compounds using existing neuromorphic chemical sensors. In this study, we present an artificial olfactory system (AOS), developed through the integration of human olfactory receptors (hORs) and artificial synapses. This AOS is engineered by interfacing an hOR-functionalized extended gate with an organic synaptic device. The AOS generates distinct patterns for odorants and mixtures thereof, at the molecular chain length level, attributed to specific hOR-odorant binding affinities. This approach enables precise pattern recognition via training and inference simulations. These findings establish a foundation for the development of high-performance sensor platforms and artificial sensory systems, which are ideal for applications in wearable and implantable devices.

CX3CL1 (Fractalkine)-CX3CR1 Axis in Inflammation-Induced Angiogenesis and Tumorigenesis

The chemotactic cytokine fractalkine (FKN, chemokine CX3CL1) has unique properties resulting from the combination of chemoattractants and adhesion molecules. The soluble form (sFKN) has chemotactic properties and strongly attracts T cells and monocytes. The membrane-bound form (mFKN) facilitates diapedesis and is responsible for cell-to-cell adhesion, especially by promoting the strong adhesion of leukocytes (monocytes) to activated endothelial cells with the subsequent formation of an extracellular matrix and angiogenesis. FKN signaling occurs via CX3CR1, which is the only known member of the CX3C chemokine receptor subfamily. Signaling within the FKN-CX3CR1 axis plays an important role in many processes related to inflammation and the immune response, which often occur simultaneously and overlap. FKN is strongly upregulated by hypoxia and/or inflammation-induced inflammatory cytokine release, and it may act locally as a key angiogenic factor in the highly hypoxic tumor microenvironment. The importance of the FKN/CX3CR1 signaling pathway in tumorigenesis and cancer metastasis results from its influence on cell adhesion, apoptosis, and cell migration. This review presents the role of the FKN signaling pathway in the context of angiogenesis in inflammation and cancer. The mechanisms determining the pro- or anti-tumor effects are presented, which are the cause of the seemingly contradictory results that create confusion regarding the therapeutic goals.

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Structural and pharmacological characterization of a medium-chain fatty acid receptor GPR84 in common carp (Cyprinus carpio)

The medium-chain fatty acid receptor GPR84, a member of the G protein-coupled receptor family, is mainly expressed in macrophages and microglia, and is involved in the regulation of inflammatory responses and retinal development in mammals and amphibians. However, structure, tissue distribution, and pharmacology of this receptor have rarely been reported in fish. In this study, we cloned the coding sequence (CDS) of common carp GPR84 (ccGPR84), examined its tissue distribution, and explored its cellular signaling function. The results showed that the CDS of ccGPR84 is 1191 bp and encodes a putative protein with 396 amino acids. Phylogenetic and chromosomal synteny analyses revealed that ccGPR84 was evolutionarily conserved with Cyprinids. Real-time quantitative PCR (qPCR) indicated that ccGPR84 was predominantly expressed in the intestine and spleen. Luciferase reporter assay demonstrated that nonanoic acid, capric acid (decanoic acid), undecanoic acid and lauric acid could inhibit cAMP signaling pathway and activate MAPK/ERK signaling pathway, while the potencies of these four fatty acids on the two signaling pathways were different. Lauric acid has the highest inhibitory potency on cAMP signaling pathway, followed by undecanoic acid, nonanoic acid, and capric acid. While for MAPK/ERK signaling pathway, nonanoic acid has the highest activation potency, followed by undecanoic acid, capric acid, and lauric acid. These findings lay the foundation for revealing the roles of different medium-chain fatty acids in the inflammatory response of common carp.

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.

Fully activated structure of the sterol-bound Smoothened GPCR-Gi protein complex

Smoothened (SMO) is an oncoprotein and signal transducer in the Hedgehog signaling pathway that regulates cellular differentiation and embryogenesis. As a member of the Frizzled (Class F) family of G protein-coupled receptors (GPCRs), SMO biochemically and functionally interacts with Gi family proteins. However, key molecular features of fully activated, G protein-coupled SMO remain elusive. We present the atomistic structure of activated human SMO complexed with the heterotrimeric Gi protein and two sterol ligands, equilibrated at 310 K in a full lipid bilayer at physiological salt concentration and pH. In contrast to previous experimental structures, our equilibrated SMO complex exhibits complete breaking of the pi-cation interaction between R4516.32 and W5357.55, a hallmark of Class F receptor activation. The Gi protein couples to SMO at seven strong anchor points similar to those in Class A GPCRs: intracellular loop 1, intracellular loop 2, transmembrane helix 6, and helix 8. On the path to full activation, we find that the extracellular cysteine-rich domain (CRD) undergoes a dramatic tilt, following a trajectory suggested by positions of the CRD in active and inactive experimental SMO structures. Strikingly, a sterol ligand bound to a shallow transmembrane domain (TMD) site in the initial structure migrates to a deep TMD pocket found exclusively in activator-bound SMO complexes. Thus, our results indicate that SMO interacts with Gi prior to full activation to break the molecular lock, form anchors with Gi subunits, tilt the CRD, and facilitate migration of a sterol ligand in the TMD to an activated position.

Integrating unsupervised language model with multi-view multiple sequence alignments for high-accuracy inter-chain contact prediction

Accurate identification of inter-chain contacts in the protein complex is critical to determine the corresponding 3D structures and understand the biological functions. We proposed a new deep learning method, ICCPred, to deduce the inter-chain contacts from the amino acid sequences of the protein complex. This pipeline was built on the designed deep residual network architecture, integrating the pre-trained language model with three multiple sequence alignments (MSAs) from different biological views. Experimental results on 709 non-redundant benchmarking protein complexes showed that the proposed ICCPred significantly increased inter-chain contact prediction accuracy compared to the state-of-the-art approaches. Detailed data analyses showed that the significant advantage of ICCPred lies in the utilization of pre-trained transformer language models which can effectively extract the complementary co-evolution diversity from three MSAs. Meanwhile, the designed deep residual network enhances the correlation between the co-evolution diversity and the patterns of inter-chain contacts. These results demonstrated a new avenue for high-accuracy deep-learning inter-chain contact prediction that is applicable to large-scale protein-protein interaction annotations from sequence alone.

Interaction of Synthetic Cannabinoid Receptor Agonists with Cannabinoid Receptor I: Insights into Activation Molecular Mechanism

Synthetic cannabinoid receptor agonists (SCRAs) have become a wide group of new psychoactive substances since the 2010s. For the last few years, the X-ray structures of the complexes of cannabinoid receptor I (CB1) with SCRAs as well as the complexes of CB1 with its antagonist have been published. Based on those data, SCRA-CB1 interactions are analyzed in detail, using molecular modeling and molecular dynamics simulations. The molecular mechanism of the conformational transformation of the transmembrane domain of CB1 caused by its interaction with SCRA is studied. These conformational changes allosterically modulate the CB1-Gi complex, providing activation of the Gi protein. Based on the X-ray-determined structures of the CB1-ligand complexes, a stable apo conformation of inactive CB1 with a relatively low potential barrier of receptor activation was modeled. For that model, molecular dynamic simulations of SCRA binding to CB1 led to the active state of CB1, which allowed us to explore the key features of this activation and the molecular mechanism of the receptor's structural transformation. The simulated CB1 activation is in accordance with the previously published experimental data for the activation at protein mutations or structural changes of ligands. The key feature of the suggested activation mechanism is the determination of the stiff core of the CB1 transmembrane domain and the statement that the entire conformational transformation of the receptor to the active state is caused by a shift of alpha helix TM7 relative to this core. The shift itself is caused by protein-ligand interactions. It was verified via steered molecular dynamics simulations of the X-ray-determined structures of the inactive receptor, which resulted in the active conformation of CB1 irrespective of the placement of agonist ligand in the receptor's active site.

Molecular Cloning, Tissue Distribution, and Pharmacological Characterization of GPR84 in Grass Carp (Ctenopharyngodon Idella)

The G-protein-coupled receptor GPR84, activated by medium-chain fatty acids, primarily expressed in macrophages and microglia, is involved in inflammatory responses and retinal development in mammals and amphibians. However, our understanding of its structure, function, tissue expression, and signaling pathways in fish is limited. In this study, we cloned and characterized the coding sequence of GPR84 (ciGPR84) in grass carp. A phylogenetic analysis revealed its close relationship with bony fishes. High expression levels of GPR84 were observed in the liver and spleen. The transfection of HEK293T cells with ciGPR84 demonstrated its responsiveness to medium-chain fatty acids and diindolylmethane (DIM). Capric acid, undecanoic acid, and lauric acid activated ERK and inhibited cAMP signaling. Lauric acid showed the highest efficiency in activating the ERK pathway, while capric acid was the most effective in inhibiting cAMP signaling. Notably, DIM did not activate GPR84 in grass carp, unlike in mammals. These findings provide valuable insights for mitigating chronic inflammation in grass carp farming and warrant further exploration of the role of medium-chain fatty acids in inflammation regulation in this species.

Dynamic allosteric networks drive adenosine A1 receptor activation and G-protein coupling

G-protein coupled receptors (GPCRs) present specific activation pathways and signaling among receptor subtypes. Hence, an extensive knowledge of the structural dynamics of the receptor is critical for the development of therapeutics. Here, we target the adenosine A1 receptor (A1R), for which a negligible number of drugs have been approved. We combine molecular dynamics simulations, enhanced sampling techniques, network theory and pocket detection to decipher the activation pathway of A1R, decode the allosteric networks and identify transient pockets. The A1R activation pathway reveal hidden intermediate and pre-active states together with the inactive and fully-active states observed experimentally. The protein energy networks computed throughout these conformational states successfully unravel the extra and intracellular allosteric centers and the communication pathways that couples them. We observe that the allosteric networks are dynamic, being increased along activation and fine-tuned in presence of the trimeric G-proteins. Overlap of transient pockets and energy networks uncover how the allosteric coupling between pockets and distinct functional regions of the receptor is altered along activation. By an in-depth analysis of the bridge between activation pathway, energy networks and transient pockets, we provide a further understanding of A1R. This information can be useful to ease the design of allosteric modulators for A1R.

Dynamics of the Apo µ-Opioid Receptor in Complex with Gi Protein

Opioid receptors, particularly the µ-opioid receptor (μOR), play a pivotal role in mediating the analgesic and addictive effects of opioid drugs. G protein signaling is an important pathway of μOR function, usually associated with painkilling effects. However, the molecular mechanisms underlying the interaction between the μOR and G protein remain poorly understood. In this study, we employed classical all-atom molecular dynamics simulations to investigate the structural changes occurring with the μOR-G protein complex under two different conditions: with the G protein in the apo form (open) and with the GDP bound G protein (closed, holo form). The receptor was in the apo form and active conformation in both cases, and the simulation time comprised 1µs for each system. In order to assess the effect of the G protein coupling on the receptor activation state, three parameters were monitored: the correlation of the distance between TM3 and TM6 and the RMSD of the NPxxYA motif; the universal activation index (A100); and the χ2 dihedral distribution of residue W2936.48. When complexed with the open G protein, receptor conformations with intermediate activation state prevailed throughout the molecular dynamics, whereas in the condition with the closed G protein, mostly inactive conformations of the receptor were observed. The major effect of the G protein in the receptor conformation comes from a steric hindrance involving an intracellular loop of the receptor and a β-sheet region of the G protein. This suggests that G-protein precoupling is essential for receptor activation, but this fact is not sufficient for complete receptor activation.

Dissecting the novel abilities of aripiprazole: The generation of anti-colorectal cancer effects by targeting Gαq via HTR2B

Colorectal cancer (CRC) is a type of malignant tumor that seriously threatens human health and life, and its treatment has always been a difficulty and hotspot in research. Herein, this study for the first time reports that antipsychotic aripiprazole (Ari) against the proliferation of CRC cells both in vitro and in vivo, but with less damage in normal colon cells. Mechanistically, the results showed that 5-hydroxytryptamine 2B receptor (HTR2B) and its coupling protein G protein subunit alpha q (Gαq) were highly distributed in CRC cells. Ari had a strong affinity with HTR2B and inhibited HTR2B downstream signaling. Blockade of HTR2B signaling suppressed the growth of CRC cells, but HTR2B was not found to have independent anticancer activity. Interestingly, the binding of Gαq to HTR2B was decreased after Ari treatment. Knockdown of Gαq not only restricted CRC cell growth, but also directly affected the anti-CRC efficacy of Ari. Moreover, an interaction between Ari and Gαq was found in that the mutation at amino acid 190 of Gαq reduced the efficacy of Ari. Thus, these results confirm that Gαq coupled to HTR2B was a potential target of Ari in mediating CRC proliferation. Collectively, this study provides a novel effective strategy for CRC therapy and favorable evidence for promoting Ari as an anticancer agent.

Cholesterol Biases the Conformational Landscape of the Chemokine Receptor CCR3: A MAS SSNMR-Filtered Molecular Dynamics Study

Cholesterol directs the pathway of ligand-induced G protein-coupled receptor (GPCR) signal transduction. The GPCR C-C motif chemokine receptor 3 (CCR3) is the principal chemotactic receptor for eosinophils, with roles in cancer metastasis and autoinflammatory conditions. Recently, we discovered a direct correlation between bilayer cholesterol and increased agonist-triggered CCR3 signal transduction. However, the allosteric molecular mechanism escalating ligand affinity and G protein coupling is unknown. To study cholesterol-guided CCR3 conformational selection, we implement comparative, objective measurement of protein architectures by scoring shifts (COMPASS) to grade model structures from molecular dynamics simulations. In this workflow, we scored predicted chemical shifts against 2-dimensional solid-state NMR 13C-13C correlation spectra of U-15N,13C-CCR3 samples prepared with and without cholesterol. Our analysis of trajectory model structures uncovers that cholesterol induces site-specific conformational restraint of extracellular loop (ECL) 2 and conserved motion in transmembrane helices and ECL3 not observed in simulations of bilayers with only phosphatidylcholine lipids. PyLipID analysis implicates direct cholesterol agency in CCR3 conformational selection and dynamics. Residue-residue contact scoring shows that cholesterol biases the conformational selection of the orthosteric pocket involving Y411.39, Y1133.32, and E2877.39. Lastly, we observe contact remodeling in activation pathway residues centered on the initial transmission switch, Na+ pocket, and R3.50 in the DRY motif. Our observations have unique implications for understanding of CCR3 ligand recognition and specificity and provide mechanistic insight into how cholesterol functions as an allosteric regulator of CCR3 signal transduction.

Ligand and G-protein selectivity in the κ-opioid receptor

The κ-opioid receptor (KOR) represents a highly desirable therapeutic target for treating not only pain but also addiction and affective disorders1. However, the development of KOR analgesics has been hindered by the associated hallucinogenic side effects2. The initiation of KOR signalling requires the Gi/o-family proteins including the conventional (Gi1, Gi2, Gi3, GoA and GoB) and nonconventional (Gz and Gg) subtypes. How hallucinogens exert their actions through KOR and how KOR determines G-protein subtype selectivity are not well understood. Here we determined the active-state structures of KOR in a complex with multiple G-protein heterotrimers-Gi1, GoA, Gz and Gg-using cryo-electron microscopy. The KOR-G-protein complexes are bound to hallucinogenic salvinorins or highly selective KOR agonists. Comparisons of these structures reveal molecular determinants critical for KOR-G-protein interactions as well as key elements governing Gi/o-family subtype selectivity and KOR ligand selectivity. Furthermore, the four G-protein subtypes display an intrinsically different binding affinity and allosteric activity on agonist binding at KOR. These results provide insights into the actions of opioids and G-protein-coupling specificity at KOR and establish a foundation to examine the therapeutic potential of pathway-selective agonists of KOR.

ELMO1 Deficiency Reduces Neutrophil Chemotaxis in Murine Peritonitis

Peritoneal inflammation remains a major cause of treatment failure in patients with kidney failure who receive peritoneal dialysis. Peritoneal inflammation is characterized by an increase in neutrophil infiltration. However, the molecular mechanisms that control neutrophil recruitment in peritonitis are not fully understood. ELMO and DOCK proteins form complexes which function as guanine nucleotide exchange factors to activate the small GTPase Rac to regulate F-actin dynamics during chemotaxis. In the current study, we found that deletion of the Elmo1 gene causes defects in chemotaxis and the adhesion of neutrophils. ELMO1 plays a role in the fMLP-induced activation of Rac1 in parallel with the PI3K and mTORC2 signaling pathways. Importantly, we also reveal that peritoneal inflammation is alleviated in Elmo1 knockout mice in the mouse model of thioglycollate-induced peritonitis. Our results suggest that ELMO1 functions as an evolutionarily conserved regulator for the activation of Rac to control the chemotaxis of neutrophils both in vitro and in vivo. Our results suggest that the targeted inhibition of ELMO1 may pave the way for the design of novel anti-inflammatory therapies for peritonitis.

Predicted Three-Dimensional Structure of the GCR1 Putative GPCR in Arabidopsis thaliana and Its Binding to Abscisic Acid and Gibberellin A1

GCR1 has been proposed as a plant analogue to animal G-protein-coupled receptors that can promote or regulate several physiological processes by binding different phytohormones. For instance, abscisic acid (ABA) and gibberellin A1 (GA1) have been shown to promote or regulate germination and flowering, root elongation, dormancy, and biotic and abiotic stresses, among others. They may act through binding to GCR1, which would put GCR1 at the heart of key signaling processes of agronomic importance. Unfortunately, this GPCR function has yet to be fully validated due to the lack of an X-ray or cryo-EM 3D atomistic structure for GCR1. Here, we used the primary sequence data from Arabidopsis thaliana and the GEnSeMBLE complete sampling method to examine 13 trillion possible packings of the 7 transmembrane helical domains corresponding to GCR1 to downselect an ensemble of 25 configurations likely to be accessible to the binding of ABA or GA1. We then predicted the best binding sites and energies for both phytohormones to the best GCR1 configurations. To provide the basis for the experimental validation of our predicted ligand-GCR1 structures, we identify several mutations that should improve or weaken the interactions. Such validations could help establish the physiological role of GCR1 in plants.

Allosteric Regulation of G-Protein-Coupled Receptors: From Diversity of Molecular Mechanisms to Multiple Allosteric Sites and Their Ligands

Allosteric regulation is critical for the functioning of G protein-coupled receptors (GPCRs) and their signaling pathways. Endogenous allosteric regulators of GPCRs are simple ions, various biomolecules, and protein components of GPCR signaling (G proteins and β-arrestins). The stability and functional activity of GPCR complexes is also due to multicenter allosteric interactions between protomers. The complexity of allosteric effects caused by numerous regulators differing in structure, availability, and mechanisms of action predetermines the multiplicity and different topology of allosteric sites in GPCRs. These sites can be localized in extracellular loops; inside the transmembrane tunnel and in its upper and lower vestibules; in cytoplasmic loops; and on the outer, membrane-contacting surface of the transmembrane domain. They are involved in the regulation of basal and orthosteric agonist-stimulated receptor activity, biased agonism, GPCR-complex formation, and endocytosis. They are targets for a large number of synthetic allosteric regulators and modulators, including those constructed using molecular docking. The review is devoted to the principles and mechanisms of GPCRs allosteric regulation, the multiplicity of allosteric sites and their topology, and the endogenous and synthetic allosteric regulators, including autoantibodies and pepducins. The allosteric regulation of chemokine receptors, proteinase-activated receptors, thyroid-stimulating and luteinizing hormone receptors, and beta-adrenergic receptors are described in more detail.

Intermediate-state-trapped mutants pinpoint G protein-coupled receptor conformational allostery

Understanding the roles of intermediate states in signaling is pivotal to unraveling the activation processes of G protein-coupled receptors (GPCRs). However, the field is still struggling to define these conformational states with sufficient resolution to study their individual functions. Here, we demonstrate the feasibility of enriching the populations of discrete states via conformation-biased mutants. These mutants adopt distinct distributions among five states that lie along the activation pathway of adenosine A2A receptor (A2AR), a class A GPCR. Our study reveals a structurally conserved cation-π lock between transmembrane helix VI (TM6) and Helix8 that regulates cytoplasmic cavity opening as a "gatekeeper" for G protein penetration. A GPCR activation process based on the well-discerned conformational states is thus proposed, allosterically micro-modulated by the cation-π lock and a previously well-defined ionic interaction between TM3 and TM6. Intermediate-state-trapped mutants will also provide useful information in relation to receptor-G protein signal transduction.

Metadynamics simulations leveraged by statistical analyses and artificial intelligence-based tools to inform the discovery of G protein-coupled receptor ligands

G Protein-Coupled Receptors (GPCRs) are a large family of membrane proteins with pluridimensional signaling profiles. They undergo ligand-specific conformational changes, which in turn lead to the differential activation of intracellular signaling proteins and the consequent triggering of a variety of biological responses. This conformational plasticity directly impacts our understanding of GPCR signaling and therapeutic implications, as do ligand-specific kinetic differences in GPCR-induced transducer activation/coupling or GPCR-transducer complex stability. High-resolution experimental structures of ligand-bound GPCRs in the presence or absence of interacting transducers provide important, yet limited, insights into the highly dynamic process of ligand-induced activation or inhibition of these receptors. We and others have complemented these studies with computational strategies aimed at characterizing increasingly accurate metastable conformations of GPCRs using a combination of metadynamics simulations, state-of-the-art algorithms for statistical analyses of simulation data, and artificial intelligence-based tools. This minireview provides an overview of these approaches as well as lessons learned from them towards the identification of conformational states that may be difficult or even impossible to characterize experimentally and yet important to discover new GPCR ligands.

Determination of the nanoscale electrical properties of olfactory receptor hOR1A1 and their dependence on ligand binding: Towards the development of capacitance-operated odorant biosensors

The transduction of odorant binding into cellular signaling by olfactory receptors (ORs) is not understood and knowing its mechanism would enable developing new pharmacology and biohybrid electronic detectors of volatile organic compounds bearing high sensitivity and selectivity. The electrical characterization of ORs in bulk experiments is subject to microscopic models and assumptions. We have directly determined the nanoscale electrical properties of ORs immobilized in a fixed orientation, and their change upon odorant binding, using electrochemical scanning tunneling microscopy (EC-STM) in near-physiological conditions. Recordings of current versus time, distance, and electrochemical potential allows determining the OR impedance parameters and their dependence with odorant binding. Our results allow validating OR structural-electrostatic models and their functional activation processes, and anticipating a novel macroscopic biosensor based on ORs.

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

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