Beyond structure: emerging approaches to study GPCR dynamics

Orphan GPCRs in Neurodegenerative Disorders: Integrating Structural Biology and Drug Discovery Approaches

Neurodegenerative disorders, particularly Alzheimer's and Parkinson's diseases, continue to challenge modern medicine despite therapeutic advances. Orphan G-protein-coupled receptors (GPCRs) have emerged as promising targets in the central nervous system, offering new avenues for drug development. This review focuses on the structural biology of orphan GPCRs implicated in these disorders, providing a comprehensive analysis of their molecular architecture and functional mechanisms. We examine recent breakthroughs in structural determination techniques, such as cryo-electron microscopy and X-ray crystallography, which have elucidated the intricate conformations of these receptors. The review highlights how structural insights inform our understanding of orphan GPCR activation, ligand binding and signaling pathways. By integrating structural data with molecular pharmacology, we explore the potential of structure-guided approaches in developing targeted therapeutics toward orphan GPCRs. This structural-biology-centered perspective aims to deepen our comprehension of orphan GPCRs and guide future drug discovery efforts in neurodegenerative disorders.

Structural basis of psychedelic LSD recognition at dopamine D1 receptor

Understanding the kinetics of LSD in receptors and subsequent induced signaling is crucial for comprehending both the psychoactive and therapeutic effects of LSD. Despite extensive research on LSD's interactions with serotonin 2A and 2B receptors, its behavior on other targets, including dopamine receptors, has remained elusive. Here, we present cryo-EM structures of LSD/PF6142-bound dopamine D1 receptor (DRD1)-legobody complexes, accompanied by a β-arrestin-mimicking nanobody, NBA3, shedding light on the determinants of G protein coupling versus β-arrestin coupling. Structural analysis unveils a distinctive binding mode of LSD in DRD1, particularly with the ergoline moiety oriented toward TM4. Kinetic investigations uncover an exceptionally rapid dissociation rate of LSD in DRD1, attributed to the flexibility of extracellular loop 2 (ECL2). Moreover, G protein can stabilize ECL2 conformation, leading to a significant slowdown in ligand's dissociation rate. These findings establish a solid foundation for further exploration of G protein-coupled receptor (GPCR) dynamics and their relevance to signal transduction.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Structural Mass Spectrometry Captures Residue-Resolved Comprehensive Conformational Rearrangements of a G Protein-Coupled Receptor

G protein-coupled receptor (GPCR) structural studies with in-solution spectroscopic approaches have offered distinctive insights into GPCR activation and signaling that highly complement those yielded from structural snapshots by crystallography or cryo-EM. While most current spectroscopic approaches allow for probing structural changes at selected residues or loop regions, they are not suitable for capturing a holistic view of GPCR conformational rearrangements across multiple domains. Herein, we develop an approach based on limited proteolysis mass spectrometry (LiP-MS) to simultaneously monitor conformational alterations of a large number of residues spanning both flexible loops and structured transmembrane domains for a given GPCR. To benchmark LiP-MS for GPCR conformational profiling, we studied the adenosine 2A receptor (A2AR) in response to different ligand binding (agonist/antagonist/allosteric modulators) and G protein coupling. Systematic and residue-resolved profiling of A2AR conformational rearrangements by LiP-MS precisely captures structural mechanisms in multiple domains underlying ligand engagement, receptor activation, and allostery, and may also reflect local conformational flexibility. Furthermore, these residue-resolution structural fingerprints of the A2AR protein allow us to readily classify ligands of different pharmacology and distinguish the G protein-coupled state. Thus, our study provides a new structural MS approach that would be generalizable to characterizing conformational transition and plasticity for challenging integral membrane proteins.

Conformational dynamics underlying atypical chemokine receptor 3 activation

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon β-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of β-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

The cannabinoid receptor CB1 is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 states in apo conditions. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganization of key salt bridge and hydrogen bond networks known to control CB1 activation between states. For instance, a conserved D213-Y224 hydrogen bond and D184-K192 salt bridge interactions showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

"Photo-Adrenalines": Photoswitchable β2 -Adrenergic Receptor Agonists as Molecular Probes for the Study of Spatiotemporal Adrenergic Signaling

β2 -adrenergic receptor (β2 -AR) agonists are used for the treatment of asthma and chronic obstructive pulmonary disease, but also play a role in other complex disorders including cancer, diabetes and heart diseases. As the cellular and molecular mechanisms in various cells and tissues of the β2 -AR remain vastly elusive, we developed tools for this investigation with high temporal and spatial resolution. Several photoswitchable β2 -AR agonists with nanomolar activity were synthesized. The most potent agonist for β2 -AR with reasonable switching is a one-digit nanomolar active, trans-on arylazopyrazole-based adrenaline derivative and comprises valuable photopharmacological properties for further biological studies with high structural accordance to the native ligand adrenaline.

Understanding the molecular mechanisms of odorant binding and activation of the human OR52 family

Structural and mechanistic studies on human odorant receptors (ORs), key in olfactory signaling, are challenging because of their low surface expression in heterologous cells. The recent structure of OR51E2 bound to propionate provided molecular insight into odorant recognition, but the lack of an inactive OR structure limited understanding of the activation mechanism of ORs upon odorant binding. Here, we determined the cryo-electron microscopy structures of consensus OR52 (OR52cs), a representative of the OR52 family, in the ligand-free (apo) and octanoate-bound states. The apo structure of OR52cs reveals a large opening between transmembrane helices (TMs) 5 and 6. A comparison between the apo and active structures of OR52cs demonstrates the inward and outward movements of the extracellular and intracellular segments of TM6, respectively. These results, combined with molecular dynamics simulations and signaling assays, shed light on the molecular mechanisms of odorant binding and activation of the OR52 family.

Functional GPCR Expression in Eukaryotic LEXSY System

G protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins in the human genome, and represent one of the most important classes of drug targets. Their structural studies facilitate rational drug discovery. However, atomic structures of only about 20% of human GPCRs have been solved to date. Recombinant production of GPCRs for structural studies at a large scale is challenging due to their low expression levels and stability. Therefore, in this study, we explored the efficacy of the eukaryotic system LEXSY (Leishmania tarentolae) for GPCR production. We selected the human A2A adenosine receptor (A2AAR), as a model protein, expressed it in LEXSY, purified it, and compared with the same receptor produced in insect cells, which is the most popular expression system for structural studies of GPCRs. The A2AAR purified from both expression systems showed similar purity, stability, ligand-induced conformational changes and structural dynamics, with a remarkably higher protein yield in the case of LEXSY expression. Overall, our results suggest that LEXSY is a promising platform for large-scale production of GPCRs for structural studies.

Structural diversity of leukotriene G-protein coupled receptors

Dihydroxy acid leukotriene (LTB4) and cysteinyl leukotrienes (LTC4, LTD4, and LTE4) are inflammatory mediators derived from arachidonic acid via the 5-lipoxygenase pathway. While structurally similar, these two types of leukotrienes (LTs) exert their functions through interactions with two distinct G protein-coupled receptor (GPCR) families, BLT and CysLT receptors, which share low sequence similarity and belong to phylogenetically divergent GPCR groups. Selective antagonism of LT receptors has been proposed as a promising strategy for the treatment of many inflammation-related diseases including asthma and chronic obstructive pulmonary disease, rheumatoid arthritis, cystic fibrosis, diabetes, and several types of cancer. Selective CysLT1R antagonists are currently used as antiasthmatic drugs, however, there are no approved drugs targeting CysLT2 and BLT receptors. In this review, we highlight recently published structures of BLT1R and CysLTRs revealing unique structural features of the two receptor families. X-ray and cryo-EM data shed light on their overall conformations, differences in functional motifs involved in receptor activation, and details of the ligand-binding pockets. An unexpected binding mode of the selective antagonist BIIL260 in the BLT1R structure makes it the first example of a compound targeting the sodium-binding site of GPCRs and suggests a novel strategy for the receptor activity modulation. Taken together, these recent structural data reveal dramatic differences in the molecular architecture of the two LT receptor families and pave the way to new therapeutic strategies of selective targeting individual receptors with novel tool compounds obtained by the structure-based drug design approach.

Illuminating GPCR signaling mechanisms by NMR spectroscopy with stable-isotope labeled receptors

G protein-coupled receptors (GPCRs) exhibit remarkable structural plasticity, which underlies their capacity to recognize a wide range of extracellular molecules and interact with intracellular partner proteins. Nuclear magnetic resonance (NMR) spectroscopy is uniquely well-suited to investigate GPCR structural plasticity, enabled by stable-isotope "probes" incorporated into receptors that inform on structure and dynamics. Progress with stable-isotope labeling methods in Eukaryotic expression systems has enabled production of native or nearly-native human receptors with varied and complementary distributions of NMR probes. These advances have opened up new avenues for investigating the roles of conformational dynamics in signaling processes, including by mapping allosteric communication networks, understanding the specificity of GPCR interactions with partner proteins and exploring the impact of membrane environments on GPCR function.

Recovery of the histamine H3 receptor activity lost in yeast cells through error-prone PCR and in vivo selection

G protein-coupled receptors (GPCRs) are the largest protein family in humans and are important drug targets. Yeast, especially Saccharomyces cerevisiae, is a useful host for modifying the function and stability of GPCRs through protein engineering, which is advantageous for mammalian cells. When GPCRs are expressed in yeast, their function is often impaired. In this study, we performed random mutagenesis using error-prone PCR and then an in vivo screening to obtain mutants that recovered the activity of the human histamine H3 receptor (H3R), which loses its signaling function when expressed in yeast. Four mutations with recovered activity were identified after screening. Three of the mutations were identified near the DRY and NPxxY motifs of H3R, which are important for activation and are commonly found in class A GPCRs. The mutants responded exclusively to the yeast YB1 strain harboring Gi-chimera proteins, showing retention of G protein specificity. Analysis of one of the mutants with recovered activity, C415R, revealed that it maintained its ligand-binding characteristics. The strategy used in this study may enable the recovery of the activity of other GPCRs that do not function in S. cerevisiae and may be useful in creating GPCRs mutants stabilized in their active conformations.

Conformational Dynamics of the Activated GLP-1 Receptor-Gs Complex Revealed by Cross-Linking Mass Spectrometry and Integrative Structure Modeling

Despite advances in characterizing the structures and functions of G protein-coupled receptors (GPCRs), our understanding of GPCR activation and signaling is still limited by the lack of information on conformational dynamics. It is particularly challenging to study the dynamics of GPCR complexes with their signaling partners because of their transient nature and low stability. Here, by combining cross-linking mass spectrometry (CLMS) with integrative structure modeling, we map the conformational ensemble of an activated GPCR-G protein complex at near-atomic resolution. The integrative structures describe heterogeneous conformations for a high number of potential alternative active states of the GLP-1 receptor-G_s complex. These structures show marked differences from the previously determined cryo-EM structure, especially at the receptor-G_s interface and in the interior of the G_s heterotrimer. Alanine-scanning mutagenesis coupled with pharmacological assays validates the functional significance of 24 interface residue contacts only observed in the integrative structures, yet absent in the cryo-EM structure. Through the integration of spatial connectivity data from CLMS with structure modeling, our study provides a new approach that is generalizable to characterizing the conformational dynamics of GPCR signaling complexes.

Sub-millisecond conformational dynamics of the A2A adenosine receptor revealed by single-molecule FRET

The complex pharmacology of G-protein-coupled receptors (GPCRs) is defined by their multi-state conformational dynamics. Single-molecule Förster Resonance Energy Transfer (smFRET) is well suited to quantify dynamics for individual protein molecules; however, its application to GPCRs is challenging. Therefore, smFRET has been limited to studies of inter-receptor interactions in cellular membranes and receptors in detergent environments. Here, we performed smFRET experiments on functionally active human A_2A adenosine receptor (A_2AAR) molecules embedded in freely diffusing lipid nanodiscs to study their intramolecular conformational dynamics. We propose a dynamic model of A_2AAR activation that involves a slow (>2 ms) exchange between the active-like and inactive-like conformations in both apo and antagonist-bound A_2AAR, explaining the receptor's constitutive activity. For the agonist-bound A_2AAR, we detected faster (390 ± 80 µs) ligand efficacy-dependent dynamics. Our work establishes a general smFRET platform for GPCR investigations that can potentially be used for drug screening and/or mechanism-of-action studies.

Allosteric Modulators of Dopamine D2 Receptors for Fine-Tuning of Dopaminergic Neurotransmission in CNS Diseases: Overview, Pharmacology, Structural Aspects and Synthesis

Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays a hot topic in medicinal chemistry. Allosteric modulators, i.e., compounds which bind in a receptor site topologically distinct from orthosteric sites, exhibit a number of advantages. They are more selective, safer and display a ceiling effect which prevents overdosing. Allosteric modulators of dopamine D₂ receptor are potential drugs against a number of psychiatric and neurological diseases, such as schizophrenia and Parkinson's disease. In this review, an insightful summary of current research on D₂ receptor modulators is presented, ranging from their pharmacology and structural aspects of ligand-receptor interactions to their synthesis.

The Nature of Functional Features of Different Classes of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are a critical family in the human proteome and are involved in various physiological processes. They are also the most important drug target, with approximately 30% of approved drugs acting on such receptors. The members of the family are divided into six classes based on their structural and functional characteristics. Understanding their structural-functional relationships will benefit us in future drug development. In this article, we investigate the features of protein function, structure, and energy that describe the dynamics of the GPCR activation process between different families. GPCRs straddle the cell membrane and transduce signals from outside the membrane into the cell. During the process, the conformational change in GPCRs that is activated by the binding of signal molecules is essential. During the binding process, different types of signal molecules result in different signal transfer efficiencies. Therefore, the GPCR classes show a variety of structures and activation processes. Based on the experimental crystal structures, we modeled the activation process of the β2 adrenergic receptor (β2AR), glucagon receptor (GCGR), and metabotropic glutamate receptor 2 (mGluR2), which represent class A, B, and C GPCRs, respectively. We calculated their activation free-energy landscapes and analyzed the structure-energy-function relationship. The results show a consistent picture of the activation mechanisms between different types of GPCRs. This could also provide us a way to understand other signal transduction proteins.

Electrically stimulated droplet injector for reduced sample consumption in serial crystallography

With advances in X-ray free-electron lasers (XFELs), serial femtosecond crystallography (SFX) has enabled the static and dynamic structure determination for challenging proteins such as membrane protein complexes. In SFX with XFELs, the crystals are typically destroyed after interacting with a single XFEL pulse. Therefore, thousands of new crystals must be sequentially introduced into the X-ray beam to collect full data sets. Because of the serial nature of any SFX experiment, up to 99% of the sample delivered to the X-ray beam during its "off-time" between X-ray pulses is wasted due to the intrinsic pulsed nature of all current XFELs. To solve this major problem of large and often limiting sample consumption, we report on improvements of a revolutionary sample-saving method that is compatible with all current XFELs. We previously reported 3D-printed injection devices coupled with gas dynamic virtual nozzles (GDVNs) capable of generating samples containing droplets segmented by an immiscible oil phase for jetting crystal-laden droplets into the path of an XFEL. Here, we have further improved the device design by including metal electrodes inducing electrowetting effects for improved control over droplet generation frequency to stimulate the droplet release to matching the XFEL repetition rate by employing an electrical feedback mechanism. We report the improvements in this electrically triggered segmented flow approach for sample conservation in comparison with a continuous GDVN injection using the microcrystals of lysozyme and 3-deoxy-D-manno-octulosonate 8-phosphate synthase and report the segmented flow approach for sample injection applied at the Macromolecular Femtosecond Crystallography instrument at the Linear Coherent Light Source for the first time.

A Statistical Journey through the Topological Determinants of the β2 Adrenergic Receptor Dynamics

Activation of G-protein-coupled receptors (GPCRs) is mediated by molecular switches throughout the transmembrane region of the receptor. In this work, we continued along the path of a previous computational study wherein energy transport in the β2 Adrenergic Receptor (β2-AR) was examined and allosteric switches were identified in the molecular structure through the reorganization of energy transport networks during activation. In this work, we further investigated the allosteric properties of β2-AR, using Protein Contact Networks (PCNs). In this paper, we report an extensive statistical analysis of the topological and structural properties of β2-AR along its molecular dynamics trajectory to identify the activation pattern of this molecular system. The results show a distinct character to the activation that both helps to understand the allosteric switching previously identified and confirms the relevance of the network formalism to uncover relevant functional features of protein molecules.

The Relaxin-3 Receptor, RXFP3, Is a Modulator of Aging-Related Disease

During the aging process our body becomes less well equipped to deal with cellular stress, resulting in an increase in unrepaired damage. This causes varying degrees of impaired functionality and an increased risk of mortality. One of the most effective anti-aging strategies involves interventions that combine simultaneous glucometabolic support with augmented DNA damage protection/repair. Thus, it seems prudent to develop therapeutic strategies that target this combinatorial approach. Studies have shown that the ADP-ribosylation factor (ARF) GTPase activating protein GIT2 (GIT2) acts as a keystone protein in the aging process. GIT2 can control both DNA repair and glucose metabolism. Through in vivo co-regulation analyses it was found that GIT2 forms a close coexpression-based relationship with the relaxin-3 receptor (RXFP3). Cellular RXFP3 expression is directly affected by DNA damage and oxidative stress. Overexpression or stimulation of this receptor, by its endogenous ligand relaxin 3 (RLN3), can regulate the DNA damage response and repair processes. Interestingly, RLN3 is an insulin-like peptide and has been shown to control multiple disease processes linked to aging mechanisms, e.g., anxiety, depression, memory dysfunction, appetite, and anti-apoptotic mechanisms. Here we discuss the molecular mechanisms underlying the various roles of RXFP3/RLN3 signaling in aging and age-related disorders.

Structural basis of GPCR coupling to distinct signal transducers: implications for biased signaling

Three classes of G-protein-coupled receptor (GPCR) partners - G proteins, GPCR kinases, and arrestins - preferentially bind active GPCRs. Our analysis suggests that the structures of GPCRs bound to these interaction partners available today do not reveal a clear conformational basis for signaling bias, which would have enabled the rational design of biased GRCR ligands. In view of this, three possibilities are conceivable: (i) there are no generalizable GPCR conformations conducive to binding a particular type of partner; (ii) subtle differences in the orientation of individual residues and/or their interactions not easily detectable in the receptor-transducer structures determine partner preference; or (iii) the dynamics of GPCR binding to different types of partners rather than the structures of the final complexes might underlie transducer bias.

Sampling alternative conformational states of transporters and receptors with AlphaFold2

Equilibrium fluctuations and triggered conformational changes often underlie the functional cycles of membrane proteins. For example, transporters mediate the passage of molecules across cell membranes by alternating between inward- and outward-facing states, while receptors undergo intracellular structural rearrangements that initiate signaling cascades. Although the conformational plasticity of these proteins has historically posed a challenge for traditional de novo protein structure prediction pipelines, the recent success of AlphaFold2 (AF2) in CASP14 culminated in the modeling of a transporter in multiple conformations to high accuracy. Given that AF2 was designed to predict static structures of proteins, it remains unclear if this result represents an underexplored capability to accurately predict multiple conformations and/or structural heterogeneity. Here, we present an approach to drive AF2 to sample alternative conformations of topologically diverse transporters and G-protein-coupled receptors that are absent from the AF2 training set. Whereas models of most proteins generated using the default AF2 pipeline are conformationally homogeneous and nearly identical to one another, reducing the depth of the input multiple sequence alignments by stochastic subsampling led to the generation of accurate models in multiple conformations. In our benchmark, these conformations spanned the range between two experimental structures of interest, with models at the extremes of these conformational distributions observed to be among the most accurate (average template modeling score of 0.94). These results suggest a straightforward approach to identifying native-like alternative states, while also highlighting the need for the next generation of deep learning algorithms to be designed to predict ensembles of biophysically relevant states.

Co-flow injection for serial crystallography at X-ray free-electron lasers

Serial femtosecond crystallography (SFX) is a powerful technique that exploits X-ray free-electron lasers to determine the structure of macro-molecules at room temperature. Despite the impressive exposition of structural details with this novel crystallographic approach, the methods currently available to introduce crystals into the path of the X-ray beam sometimes exhibit serious drawbacks. Samples requiring liquid injection of crystal slurries consume large quantities of crystals (at times up to a gram of protein per data set), may not be compatible with vacuum configurations on beamlines or provide a high background due to additional sheathing liquids present during the injection. Proposed and characterized here is the use of an immiscible inert oil phase to supplement the flow of sample in a hybrid microfluidic 3D-printed co-flow device. Co-flow generation is reported with sample and oil phases flowing in parallel, resulting in stable injection conditions for two different resin materials experimentally. A numerical model is presented that adequately predicts these flow-rate conditions. The co-flow generating devices reduce crystal clogging effects, have the potential to conserve protein crystal samples up to 95% and will allow degradation-free light-induced time-resolved SFX.

Development of enhanced conformational sampling methods to probe the activation landscape of GPCRs

G protein-coupled receptors (GPCRs) make up the largest superfamily of integral membrane proteins and play critical signal transduction roles in many physiological processes. Developments in molecular biology, biophysical, biochemical, pharmacological, and computational techniques aimed at these important therapeutic targets are beginning to provide unprecedented details on the structural as well as functional basis of their pleiotropic signaling mediated by G proteins, β arrestins, and other transducers. This pleiotropy presents a pharmacological challenge as the same ligand-receptor interaction can cause a therapeutic effect as well as an undesirable on-target side-effect through different downstream pathways. GPCRs don't function as simple binary on-off switches but as finely tuned shape-shifting machines described by conformational ensembles, where unique subsets of conformations may be responsible for specific signaling cascades. X-ray crystallography and more recently cryo-electron microscopy are providing snapshots of some of these functionally-important receptor conformations bound to ligands and/or transducers, which are being utilized by computational methods to describe the dynamic conformational energy landscape of GPCRs. In this chapter, we review the progress in computational conformational sampling methods based on molecular dynamics and discrete sampling approaches that have been successful in complementing biophysical and biochemical studies on these receptors in terms of their activation mechanisms, allosteric effects, actions of biased ligands, and effects of pathological mutations. Some of the sampled simulation time scales are beginning to approach receptor activation time scales. The list of conformational sampling methods and example uses discussed is not exhaustive but includes representative examples that have pushed the limits of classical molecular dynamics and discrete sampling methods to describe the activation energy landscape of GPCRs.

Molecular interaction of HIC, an agonist of P2Y1 receptor, and its role in prostate cancer apoptosis

Prostate cancer is a heterogeneous, slow growing asymptomatic cancer that predominantly affects man. A purinergic G-protein coupled receptor, P2Y1R, is targeted for its therapeutic value since it plays a crucial role in many key molecular events of cancer progression and invasion. Our previous study demonstrated that indoline derivative, 1 ((1-(2-Hydroxy-5-nitrophenyl) (4-hydroxyphenyl) methyl)indoline-4‑carbonitrile; HIC), stimulates prostate cancer cell (PCa) growth inhibition via P2Y1R. However, the mode of interaction of P2Y1R with HIC involved in this process remains unclear. Here, we have reported the molecular interactions of HIC with P2Y1R. Molecular dynamics simulation was performed that revealed the stable specific binding of the protein-ligand complex. In vitro analysis has shown increased apoptosis of PCa-cells, PC3, and DU145, upon specific interaction of P2Y1R-HIC. This was further validated using siRNA analysis that showed a higher percentage of apoptotic cells in PCa-cells transfected with P2Y-siRNA-MRS2365 than P2Y-siRNA-HIC treatment. Decreased mitochondrial membrane potential (MMP) activity and reduced glutathione (GSH) level show their role in P2Y1R-HIC mediated apoptosis. These in silico and in vitro results confirmed that HIC could induce mitochondrial apoptotic signaling through the P2Y1R activation. Thus, HIC being a potential ligand upon interaction with P2Y1R might have therapeutic value for the treatment of prostate cancer.

Biased and allosteric modulation of bone cell-expressing G protein-coupled receptors as a novel approach to osteoporosis therapy

On the cellular level, osteoporosis (OP) is a result of imbalanced bone remodeling, in which osteoclastic bone resorption outcompetes osteoblastic bone formation. Currently available OP medications include both antiresorptive and bone-forming drugs. However, their long-term use in OP patients, mainly in postmenopausal women, is accompanied by severe side effects. Notably, the fundamental coupling between bone resorption and formation processes underlies the existence of an undesirable secondary outcome that bone anabolic or anti-resorptive drugs also reduce bone formation. This drawback requires the development of anti-OP drugs capable of selectively stimulating osteoblastogenesis and concomitantly reducing osteoclastogenesis. We propose that the application of small synthetic biased and allosteric modulators of bone cell receptors, which belong to the G-protein coupled receptors (GPCR) family, could be the key to resolving the undesired anti-OP drug selectivity. This approach is based on the capacity of these GPCR modulators, unlike the natural ligands, to trigger signaling pathways that promote beneficial effects on bone remodeling while blocking potentially deleterious effects. Under the settings of OP, an optimal anti-OP drug should provide fine-tuned regulation of downstream effects, for example, intermittent cyclic AMP (cAMP) elevation, preservation of Ca²⁺ balance, stimulation of osteoprotegerin (OPG) and estrogen production, suppression of sclerostin secretion, and/or preserved/enhanced canonical β-catenin/Wnt signaling pathway. As such, selective modulation of GPCRs involved in bone remodeling presents a promising approach in OP treatment. This review focuses on the evidence for the validity of our hypothesis.

Electron paramagnetic resonance spectroscopy on G-protein-coupled receptors: Adopting strategies from related model systems

Membrane proteins, including ion channels, transporters and G-protein-coupled receptors (GPCRs), play a significant role in various physiological processes. Many of these proteins are difficult to express in large quantities, imposing crucial experimental restrictions. Nevertheless, there is now a wide variety of studies available utilizing electron paramagnetic resonance (EPR) spectroscopic techniques that expand experimental accessibility by using relatively small quantities of protein. Here, we give an overview starting from basic strategies in EPR on membrane proteins with a focus on GPCRs, while emphasizing several applications from recent years. We highlight how the arsenal of EPR-based techniques may provide significant further contributions to understanding the complex molecular machinery and energetic phenomena responsible for seamless workflow in essential biological processes.

Structural insights into ligand recognition and activation of angiotensin receptors

G protein-coupled angiotensin II receptors, AT₁R and AT₂R, are integral components of the renin-angiotensin system (RAS) that regulates blood pressure and fluid balance in humans. While AT₁R is a well-established target of angiotensin receptor blockers (ARBs) for managing hypertension and a prime system for studying biased signaling, AT₂R has been recognized as a promising target against neuropathic pain and lung fibrosis. In this review, we discuss how recent structural advances illuminate ligand-binding modes and subtype selectivity, shared and distinct features of the receptors, their transducer-coupling patterns, and downstream signaling responses. We also underscore the key ATR aspects that require further studies to fully appreciate the mechanistic framework that fine-tunes their cellular and physiological functions, providing untapped potential for drug discovery.

Molecular mechanisms of metabotropic GABAB receptor function

Metabotropic γ-aminobutyric acid G protein-coupled receptors (GABA_B) represent one of the two main types of inhibitory neurotransmitter receptors in the brain. These receptors act both pre- and postsynaptically by modulating the transmission of neuronal signals and are involved in a range of neurological diseases, from alcohol addiction to epilepsy. A series of recent cryo-EM studies revealed critical details of the activation mechanism of GABA_B Structures are now available for the receptor bound to ligands with different modes of action, including antagonists, agonists, and positive allosteric modulators, and captured in different conformational states from the inactive apo to the fully active state bound to a G protein. These discoveries provide comprehensive insights into the activation of the GABA_B receptor, which not only broaden our understanding of its structure, pharmacology, and physiological effects but also will ultimately facilitate the discovery of new therapeutic drugs and neuromodulators.

Cryo-electron microscopy structure of the antidiuretic hormone arginine-vasopressin V2 receptor signaling complex

The antidiuretic hormone arginine-vasopressin (AVP) forms a signaling complex with the V2 receptor (V2R) and the G_s protein, promoting kidney water reabsorption. Molecular mechanisms underlying activation of this critical G protein-coupled receptor (GPCR) signaling system are still unknown. To fill this gap of knowledge, we report here the cryo-electron microscopy structure of the AVP-V2R-G_s complex. Single-particle analysis revealed the presence of three different states. The two best maps were combined with computational and nuclear magnetic resonance spectroscopy constraints to reconstruct two structures of the ternary complex. These structures differ in AVP and G_s binding modes. They reveal an original receptor-G_s interface in which the Gα_s subunit penetrates deep into the active V2R. The structures help to explain how V2R R137H or R137L/C variants can lead to two severe genetic diseases. Our study provides important structural insights into the function of this clinically relevant GPCR signaling complex.

Allosteric communication regulates ligand-specific GPCR activity

G protein-coupled receptors (GPCRs) are membrane-bound proteins that are ubiquitously expressed in many cell types and take part in mediating multiple signaling pathways. GPCRs are dynamic proteins and exist in an equilibrium between an ensemble of conformational states such as inactive and fully active states. This dynamic nature of GPCRs is one of the factors that confers their basal activity even in the absence of any ligand-mediated activation. Ligands selectively bind and stabilize a subset of the conformations from the ensemble leading to a shift in the equilibrium toward the inactive or the active state depending on the nature of the ligand. This ligand-selective effect is achieved through allosteric communication between the ligand binding site and G protein or β-arrestin coupling site. Similarly, the G protein coupling to the receptor exerts the allosteric effect on the ligand binding region leading to increased binding affinity for agonists and decreased affinity for antagonists or inverse agonists. In this review, we enumerate the current state of our understanding of the mechanism of allosteric communication in GPCRs with a specific focus on the critical role of computational methods in delineating the residues involved in allosteric communication. Analyzing allosteric communication mechanism using molecular dynamics simulations has revealed (a) a structurally conserved mechanism of allosteric communication that regulates the G protein coupling, (b) a rational structure-based approach to designing selective ligands, and (c) an approach to designing allosteric GPCR mutants that are either ligand and G protein or β-arrestin selective.

GPCR oligomerization as a target for antidepressants: Focus on GPR39

At present most of the evidence for the relevance of oligomerization for the pharmacology of depression comes from in vitro studies which identified oligomers, and from neuropsychopharmacological studies of receptors which participate in oligomerization. For example, behavioural and biochemical studies in knockout animals suggest that GPR39 may mediate the antidepressant action of monoaminergic antidepressants. We have recently found long-lasting antidepressant-like effects of GPR39 agonist, thus suggesting GPR39 as a target for the development of novel antidepressant drugs. In vitro studies have shown that GPR39 oligomerizes with other GPCRs. Oligomerization of GPR39 should thus be considered in relation to the development of new antidepressants targeting this receptor as well as antidepressants targeting other receptors that may form complexes with GPR39. Here, we summarize recent data suggestive of the importance of oligomerization for the pharmacology of depression and discuss approaches for validation of this phenomenon.

GPR52 accelerates fatty acid biosynthesis in a ligand-dependent manner in hepatocytes and in response to excessive fat intake in mice

Gpr52 is an orphan G-protein-coupled receptor of unknown physiological function. We found that Gpr52-deficient (Gpr52^−/−) mice exhibit leanness associated with reduced liver weight, decreased hepatic de novo lipogenesis, and enhanced insulin sensitivity. Treatment of the hepatoma cell line HepG2 cells with c11, the synthetic GPR52 agonist, increased fatty acid biosynthesis, and GPR52 knockdown (KD) abolished the lipogenic action of c11. In addition, c11 induced the expressions of lipogenic enzymes (SCD1 and ELOVL6), whereas these inductions were attenuated by GPR52-KD. In contrast, cholesterol biosynthesis was not increased by c11, but its basal level was significantly suppressed by GPR52-KD. High-fat diet (HFD)-induced increase in hepatic expression of Pparg2 and its targets (Scd1 and Elovl6) was absent in Gpr52^−/− mice with alleviated hepatosteatosis. Our present study showed that hepatic GPR52 promotes the biosynthesis of fatty acid and cholesterol in a ligand-dependent and a constitutive manner, respectively, and Gpr52 participates in HFD-induced fatty acid synthesis in liver.

•

Hepatosteatosis is inherently an adaptive response to overnutrition to store energy

•

On the other hand, it can be a pathological condition causing insulin resistance

•

High-fat diet increases PPARγ2 expression and lipogenesis in liver via GPR52

•

Gpr52 ablation protects mice from developing hepatosteatosis and insulin resistance

A set of common movements within GPCR-G-protein complexes from variability analysis of cryo-EM datasets

G-protein coupled receptors (GPCRs) are among the most versatile signal transducers in the cell. Once activated, GPCRs sample a large conformational space and couple to G-proteins to initiate distinct signaling pathways. The dynamical behavior of GPCR-G-protein complexes is difficult characterize structurally, and it might hinder obtaining routine high-resolution density maps in single-particle reconstructions. Here, we used variability analysis on the rhodopsin-G_i-Fab16 complex cryo-EM dataset, and the results provide insights into the dynamic nature of the receptor-complex interaction. We compare the outcome of this analysis with recent results obtained on the cannabinoid-G_i- and secretin-G_s-receptor complexes. Despite differences related to the biochemical compositions of the three samples, a set of consensus movements emerges. We anticipate that systematic variability analysis on GPCR-G-protein complexes may provide useful information not only at the biological level, but also for improving the preparation of more stable samples for cryo-EM single-particle analysis.

Conformational rearrangement during activation of a metabotropic glutamate receptor

G protein-coupled receptors (GPCRs) relay information across cell membranes through conformational coupling between the ligand-binding domain and cytoplasmic signaling domain. In dimeric class C GPCRs, the mechanism of this process, which involves propagation of local ligand-induced conformational changes over 12 nm through three distinct structural domains, is unknown. Here, we used single-molecule FRET (smFRET) and live-cell imaging and found that metabotropic glutamate receptor 2 (mGluR2) interconverts between four conformational states, two of which were previously unknown, and activation proceeds through the conformational selection mechanism. Furthermore, the conformation of the ligand-binding domains and downstream domains are weakly coupled. We show that the intermediate states act as conformational checkpoints for activation and control allosteric modulation of signaling. Our results demonstrate a mechanism for activation of mGluRs where ligand binding controls the proximity of signaling domains, analogous to some receptor kinases. This design principle may be generalizable to other biological allosteric sensors.

Comprehensive gene expression profiling of human astrocytes treated with a hepatic encephalopathy-inducible factor, alpha 1-antichymotripsin

Astrocytes are major glial cells that play a critical role in brain homeostasis. Abnormalities in astrocytic function, such as hepatic encephalopathy (HE) during acute liver failure, can result in brain death following brain edema and the associated astrocyte swelling. Recently, we have identified alpha 1-antichymotripsin (ACT) to be a biomarker candidate for HE. ACT induces astrocyte swelling by upregulating aquaporin 4 (AQP4); however, the causal connection between these proteins is not clear yet. In this study, we utilized a microarray profile to screen the differentially expressed genes (DEGs) in astrocytes treated with ACT. We then performed Gene Ontology, REACTOME, and the comprehensive resource of mammalian protein complexes (CORUM) enrichment analyses of the identified DEGs. The results of these analyses indicated that the DEGs were enriched in pathways activating adenylate cyclase (AC)-coupled G protein-coupled receptors (GPCRs) and therefore were involved in the cyclic adenosine monophosphate (cAMP) signaling. These results indicate that ACT may act as a ligand of Gs-GPCRs and subsequently upregulate cAMP. As cAMP is known to upregulate AQP4 in astrocytes, these results suggest that ACT may upregulate AQP4 by activating AC GPCRs and therefore serve as a therapeutic target for acute HE.

•

The expression of AQP4 is induced by alpha 1-antichymotripsin dependent activation of adenylate cyclase.

•

Combination treatment of alpha 1-antichymotripsin and arginine vasopressin enhance the activation of adenylate cyclase.

•

Alpha 1-antichemotripsin might be a therapeutic target for acute hepatic encephalopathy.

Statistics for the analysis of molecular dynamics simulations: providing P values for agonist-dependent GPCR activation

Molecular dynamics (MD) is the common computational technique for assessing efficacy of GPCR-bound ligands. Agonist efficacy measures the capability of the ligand-bound receptor of reaching the active state in comparison with the free receptor. In this respect, agonists, neutral antagonists and inverse agonists can be considered. A collection of MD simulations of both the ligand-bound and the free receptor are needed to provide reliable conclusions. Variability in the trajectories needs quantification and proper statistical tools for meaningful and non-subjective conclusions. Multiple-factor (time, ligand, lipid) ANOVA with repeated measurements on the time factor is proposed as a suitable statistical method for the analysis of agonist-dependent GPCR activation MD simulations. Inclusion of time factor in the ANOVA model is consistent with the time-dependent nature of MD. Ligand and lipid factors measure agonist and lipid influence on receptor activation. Previously reported MD simulations of adenosine A2a receptor (A2aR) are reanalyzed with this statistical method. TM6–TM3 and TM7–TM3 distances are selected as dependent variables in the ANOVA model. The ligand factor includes the presence or absence of adenosine whereas the lipid factor considers DOPC or DOPG lipids. Statistical analysis of MD simulations shows the efficacy of adenosine and the effect of the membrane lipid composition. Subsequent application of the statistical methodology to NECA A2aR agonist, with resulting P values in consistency with its pharmacological profile, suggests that the method is useful for ligand comparison and potentially for dynamic structure-based virtual screening.