Exploring the free-energy landscape of GPCR activation

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.

Exploring the Activation Mechanism of the GPR183 Receptor

The G protein-coupled receptors (GPCRs) play a pivotal role in numerous biological processes as crucial cell membrane receptors. However, the dynamic mechanisms underlying the activation of GPR183, a specific GPCR, remain largely elusive. To address this, we employed computational simulation techniques to elucidate the activation process and key events associated with GPR183, including conformational changes from inactive to active state, binding interactions with the Gi protein complex, and GDP release. Our findings demonstrate that the association between GPR183 and the Gi protein involves the formation of receptor-specific conformations, the gradual proximity of the Gi protein to the binding pocket, and fine adjustments of the protein conformation, ultimately leading to a stable GPR183-Gi complex characterized by a high energy barrier. The presence of Gi protein partially promotes GPR183 activation, which is consistent with the observation of GPCR constitutive activity test experiments, thus illustrating the reliability of our calculations. Moreover, our study suggests the existence of a stable partially activated state preceding complete activation, providing novel avenues for future investigations. In addition, the relevance of GPR183 for various diseases, such as colitis, the response of eosinophils to Mycobacterium tuberculosis infection, antiviral properties, and pulmonary inflammation, has been emphasized, underscoring its therapeutic potential. Consequently, understanding the activation process of GPR183 through molecular dynamic simulations offers valuable kinetic insights that can aid in the development of targeted therapies.

Circulating tumor cells with metastasis-initiating competence survive fluid shear stress during hematogenous dissemination through CXCR4-PI3K/AKT signaling

To seed lethal secondary lesions, circulating tumor cells (CTCs) must survive all rate-limiting factors during hematogenous dissemination, including fluid shear stress (FSS) that poses a grand challenge to their survival. We thus hypothesized that CTCs with the ability to survive FSS in vasculature might hold metastasis-initiating competence. This study reported that FSS of physiologic magnitude selected a small subpopulation of suspended tumor cells in vitro with the traits of metastasis-initiating cells, including stemness, migration/invasion potential, cellular plasticity, and biophysical properties. These shear-selected cells generated local and metastatic tumors at the primary and distal sites efficiently, implicating their metastasis competence. Mechanistically, FSS activated the mechanosensitive protein CXCR4 and the downstream PI3K/AKT signaling, which were essential in shear-mediated selection of metastasis-competent CTCs. In summary, these findings conclude that CTCs with metastasis-initiating competence survive FSS during hematogenous dissemination through CXCR4-PI3K/AKT signaling, which may provide new therapeutic targets for the early prevention of tumor metastasis.

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.

The role of bile acid receptor TGR5 in regulating inflammatory signalling

Takeda G protein-coupled receptor 5 (TGR5) is a bile acid receptor, and its role in regulating metabolism after binding with bile acids has been established. Since the immune response depends on metabolism to provide biomolecules and energy to cope with challenging conditions, emerging evidence reveals the regulatory effects of TGR5 on the immune response. An in-depth understanding of the effect of TGR5 on immune regulation can help us disentangle the interaction of metabolism and immune response, accelerating the development of TGR5 as a therapeutic target. Herein, we reviewed more than 200 articles published in the last 20 years in PubMed, to discuss the roles of TGR5 in regulating inflammatory response, the molecular mechanism, as well as existing problems. Particularly, its anti-inflammation effect is emphasized.

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

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

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

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

Nutrient Deficiencies Impact on the Cellular and Metabolic Responses of Saxitoxin Producing Alexandrium minutum: A Transcriptomic Perspective

Dinoflagellate Alexandrium minutum Halim is commonly associated with harmful algal blooms (HABs) in tropical marine waters due to its saxitoxin production. However, limited information is available regarding the cellular and metabolic changes of A. minutum in nutrient-deficient environments. To fill this gap, our study aimed to investigate the transcriptomic responses of A. minutum under nitrogen and phosphorus deficiency. The induction of nitrogen and phosphorus deficiency resulted in the identification of 1049 and 763 differently expressed genes (DEGs), respectively. Further analysis using gene set enrichment analysis (GSEA) revealed 702 and 1251 enriched gene ontology (GO) terms associated with nitrogen and phosphorus deficiency, respectively. Our results indicate that in laboratory cultures, nitrogen deficiency primarily affects meiosis, carbohydrate catabolism, ammonium assimilation, ion homeostasis, and protein kinase activity. On the other hand, phosphorus deficiency primarily affects the carbon metabolic response, cellular ion transfer, actin-dependent cell movement, signalling pathways, and protein recycling. Our study provides valuable insights into biological processes and genes regulating A. minutum's response to nutrient deficiencies, furthering our understanding of the ecophysiological response of HABs to environmental change.

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

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

Prediction of dynamic allostery for the transmembrane domain of the sweet taste receptor subunit, TAS1R3

The sweet taste receptor plays an essential role as an energy sensor by detecting carbohydrates. However, the dynamic mechanisms of receptor activation remain unclear. Here, we describe the interactions between the transmembrane domain of the G protein-coupled sweet receptor subunit, TAS1R3, and allosteric modulators. Molecular dynamics simulations reproduced species-specific sensitivity to ligands. We found that a human-specific sweetener, cyclamate, interacted with the mouse receptor as a negative allosteric modulator. Agonist-induced allostery during receptor activation was found to destabilize the intracellular part of the receptor, which potentially interfaces with the Gα subunit, through ionic lock opening. A common human variant (R757C) of the TAS1R3 exhibited a reduced response to sweet taste, in support of our predictions. Furthermore, histidine residues in the binding site acted as pH-sensitive microswitches to modulate the sensitivity to saccharin. This study provides important insights that may facilitate the prediction of dynamic activation mechanisms for other G protein-coupled receptors.

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

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

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.

Computational study of the conformational ensemble of CX3C chemokine receptor 1 (CX3CR1) and its interactions with antagonist and agonist ligands

The CX3C chemokine receptor 1 (CX3CR1), a member of the class A of G Protein-Coupled Receptors (GPCR) superfamily, and its ligand fractalkine constitute an important biochemical axis that influence many cellular pathways involving homeostatic and inflammatory processes. They participate in the activation, chemotaxis and recruitment of multiple immunological cells such as microglia, macrophages and monocytes, and play a critical role in neuroinflammatory conditions such as Alzheimer's disease and multiple sclerosis, in the recovery from central nervous system injuries, in several chronic, peripheral inflammatory entities and in some infective processes including HIV-AIDS. In this work we present the study of the CX3CR1 receptor employing extensive atomistic Molecular Dynamics (MD) simulations with the aim to characterize the conformational ensemble of the receptor in the presence of its antagonist and agonist ligands. We analyzed the receptor conformational changes and described interactions within its key regions and the bounded ligands to identify their notable differences. Finally, we classify the features that would allow the identification of patterns that characterize a functional state to contribute to the understanding of the complexity of the GPCR superfamily.

Multi-state modeling of G-protein coupled receptors at experimental accuracy

The family of G-protein coupled receptors (GPCRs) is one of the largest protein families in the human genome. GPCRs transduct chemical signals from extracellular to intracellular regions via a conformational switch between active and inactive states upon ligand binding. While experimental structures of GPCRs remain limited, high-accuracy computational predictions are now possible with AlphaFold2. However, AlphaFold2 only predicts one state and is biased toward either the active or inactive conformation depending on the GPCR class. Here, a multi-state prediction protocol is introduced that extends AlphaFold2 to predict either active or inactive states at very high accuracy using state-annotated templated GPCR databases. The predicted models accurately capture the main structural changes upon activation of the GPCR at the atomic level. For most of the benchmarked GPCRs (10 out of 15), models in the active and inactive states were closer to their corresponding activation state structures. Median RMSDs of the transmembrane regions were 1.12 Å and 1.41 Å for the active and inactive state models, respectively. The models were more suitable for protein-ligand docking than the original AlphaFold2 models and template-based models. Finally, our prediction protocol predicted accurate GPCR structures and GPCR-peptide complex structures in GPCR Dock 2021, a blind GPCR-ligand complex modeling competition. We expect that high accuracy GPCR models in both activation states will promote understanding in GPCR activation mechanisms and drug discovery for GPCRs. At the time, the new protocol paves the way towards capturing the dynamics of proteins at high-accuracy via machine-learning methods.

Class A and C GPCR Dimers in Neurodegenerative Diseases

Neurodegenerative diseases affect over 30 million people worldwide with an ascending trend. Most individuals suffering from these irreversible brain damages belong to the elderly population, with onset between 50 and 60 years. Although the pathophysiology of such diseases is partially known, it remains unclear upon which point a disease turns degenerative. Moreover, current therapeutics can treat some of the symptoms but often have severe side effects and become less effective in long-term treatment. For many neurodegenerative diseases, the involvement of G proteincoupled receptors (GPCRs), which are key players of neuronal transmission and plasticity, has become clearer and holds great promise in elucidating their biological mechanism. With this review, we introduce and summarize class A and class C GPCRs, known to form heterodimers or oligomers to increase their signalling repertoire. Additionally, the examples discussed here were shown to display relevant alterations in brain signalling and had already been associated with the pathophysiology of certain neurodegenerative diseases. Lastly, we classified the heterodimers into two categories of crosstalk, positive or negative, for which there is known evidence.

Computational and experimental approaches to probe GPCR activation and signaling

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

Mechanism Study of Proteins under Membrane Environment

Membrane proteins play crucial roles in various physiological processes, including molecule transport across membranes, cell communication, and signal transduction. Approximately 60% of known drug targets are membrane proteins. There is a significant need to deeply understand the working mechanism of membrane proteins in detail, which is a challenging work due to the lack of available membrane structures and their large spatial scale. Membrane proteins carry out vital physiological functions through conformational changes. In the current study, we utilized a coarse-grained (CG) model to investigate three representative membrane protein systems: the TMEM16A channel, the family C GPCRs mGlu2 receptor, and the P4-ATPase phospholipid transporter. We constructed the reaction pathway of conformational changes between the two-end structures. Energy profiles and energy barriers were calculated. These data could provide reasonable explanations for TMEM16A activation, the mGlu2 receptor activation process, and P4-ATPase phospholipid transport. Although they all belong to the members of membrane proteins, they behave differently in terms of energy. Our work investigated the working mechanism of membrane proteins and could give novel insights into other membrane protein systems of interest.

Hydration-mediated G-protein-coupled receptor activation

Although G-protein–coupled receptors (GPCRs) control vast physiological pathways, their activation remains chemically and physically enigmatic. Our osmotic stress studies of the visual receptor rhodopsin have redefined the standard model of GPCR signaling by revealing the essential role of bulk water. We show results consistent with a large number of water molecules flooding the rhodopsin interior during activation to stabilize the effector binding conformation. These results suggest a model of GPCR activation in which the receptor becomes solvent-swollen upon formation of the active state. We thus demonstrate the mechanism whereby water acts as a powerful allosteric modulator of a pharmacologically important membrane protein family.

The Rhodopsin family of G-protein–coupled receptors (GPCRs) comprises the targets of nearly a third of all pharmaceuticals. Despite structural water present in GPCR X-ray structures, the physiological relevance of these solvent molecules to rhodopsin signaling remains unknown. Here, we show experimental results consistent with the idea that rhodopsin activation in lipid membranes is coupled to bulk water movements into the protein. To quantify hydration changes, we measured reversible shifting of the metarhodopsin equilibrium due to osmotic stress using an extensive series of polyethylene glycol (PEG) osmolytes. We discovered clear evidence that light activation entails a large influx of bulk water (∼80–100 molecules) into the protein, giving insight into GPCR activation mechanisms. Various size polymer osmolytes directly control rhodopsin activation, in which large solutes are excluded from rhodopsin and dehydrate the protein, favoring the inactive state. In contrast, small osmolytes initially forward shift the activation equilibrium until a quantifiable saturation point is reached, similar to gain-of-function protein mutations. For the limit of increasing osmolyte size, a universal response of rhodopsin to osmotic stress is observed, suggesting it adopts a dynamic, hydrated sponge-like state upon photoactivation. Our results demand a rethinking of the role of water dynamics in modulating various intermediates in the GPCR energy landscape. We propose that besides bound water, an influx of bulk water plays a necessary role in establishing the active GPCR conformation that mediates signaling.

The Impact of Membrane Protein Diffusion on GPCR Signaling

Spatiotemporal signal shaping in G protein-coupled receptor (GPCR) signaling is now a well-established and accepted notion to explain how signaling specificity can be achieved by a superfamily sharing only a handful of downstream second messengers. Dozens of Gs-coupled GPCR signals ultimately converge on the production of cAMP, a ubiquitous second messenger. This idea is almost always framed in terms of local concentrations, the differences in which are maintained by means of spatial separation. However, given the dynamic nature of the reaction-diffusion processes at hand, the dynamics, in particular the local diffusional properties of the receptors and their cognate G proteins, are also important. By combining some first principle considerations, simulated data, and experimental data of the receptors diffusing on the membranes of living cells, we offer a short perspective on the modulatory role of local membrane diffusion in regulating GPCR-mediated cell signaling. Our analysis points to a diffusion-limited regime where the effective production rate of activated G protein scales linearly with the receptor–G protein complex’s relative diffusion rate and to an interesting role played by the membrane geometry in modulating the efficiency of coupling.

Application of Coarse-Grained (CG) Models to Explore Conformational Pathway of Large-Scale Protein Machines

Protein machines are clusters of protein assemblies that function in order to control the transfer of matter and energy in cells. For a specific protein machine, its working mechanisms are not only determined by the static crystal structures, but also related to the conformational transition dynamics and the corresponding energy profiles. With the rapid development of crystallographic techniques, the spatial scale of resolved structures is reaching up to thousands of residues, and the concomitant conformational changes become more and more complicated, posing a great challenge for computational biology research. Previously, a coarse-grained (CG) model aiming at conformational free energy evaluation was developed and showed excellent ability to reproduce the energy profiles by accurate electrostatic interaction calculations. In this study, we extended the application of the CG model to a series of large-scale protein machine systems. The spike protein trimer of SARS-CoV-2, ATP citrate lyase (ACLY) tetramer, and P4-ATPases systems were carefully studied and discussed as examples. It is indicated that the CG model is effective to depict the energy profiles of the conformational pathway between two endpoint structures, especially for large-scale systems. Both the energy change and energy barrier between endpoint structures provide reasonable mechanism explanations for the associated biological processes, including the opening of receptor binding domain (RBD) of spike protein, the phospholipid transportation of P4-ATPase, and the loop translocation of ACLY. Taken together, the CG model provides a suitable alternative in mechanistic studies related to conformational change in large-scale protein machines.

The Nature of Noradrenergic Volume Transmission From Locus Coeruleus to Brainstem Mesencephalic Trigeminal Sensory Neurons

Noradrenergic neurons in the locus coeruleus (LC) release noradrenaline (NA) that acts via volume transmission to activate extrasynaptic G-protein coupled receptors (GPCRs) in target cells throughout the brain. As the closest projection, the dorsal LC laterally adjoins the mesencephalic trigeminal nucleus (MTN), in which proprioceptive primary sensory neurons innervating muscle spindles of jaw-closing muscles are exceptionally located. MTN neurons express α₂-adrenergic receptors (α₂-ARs) and display hyperpolarization-activated cyclic nucleotide-gated (HCN) currents (Ihs), which is downregulated by α₂-AR activation. To quantify the activity-dependent outcome of volume transmission of NA from LC to MTN, we investigated how direct LC activation inhibits Ih in MTN neurons by performing dual whole-cell recordings from LC and MTN neurons. Repetition of 20 Hz spike-train evoked with 1-s current-pulse in LC neurons every 30 s resulted in a gradual decrease in Ih evoked every 30 s, revealing a Hill-type relationship between the number of spike-trains in LC neurons and the degree of Ih inhibition in MTN neurons. On the other hand, when microstimulation was applied in LC every 30 s, an LC neuron repeatedly displayed a transient higher-frequency firing followed by a tonic firing at 5–10 Hz for 30 s. This subsequently caused a similar Hill-type inhibition of Ih in the simultaneously recorded MTN neuron, but with a smaller Hill coefficient, suggesting a lower signal transduction efficacy. In contrast, 20 Hz activity induced by a 1-s pulse applied every 5–10 s caused only a transient facilitation of Ih inhibition followed by a forced termination of Ih inhibition. Thus, the three modes of LC activities modulated the volume transmission to activate α₂-adrenergic GPCR to differentially inhibit Ih in MTN neurons.

The Many Faces of G Protein-Coupled Receptor 143, an Atypical Intracellular Receptor

GPCRs transform extracellular stimuli into a physiological response by activating an intracellular signaling cascade initiated via binding to G proteins. Orphan G protein-coupled receptors (GPCRs) hold the potential to pave the way for development of new, innovative therapeutic strategies. In this review we will introduce G protein-coupled receptor 143 (GPR143), an enigmatic receptor in terms of classification within the GPCR superfamily and localization. GPR143 has not been assigned to any of the GPCR families due to the lack of common structural motifs. Hence we will describe the most important motifs of classes A and B and compare them to the protein sequence of GPR143. While a precise function for the receptor has yet to be determined, the protein is expressed abundantly in pigment producing cells. Many GPR143 mutations cause X-linked Ocular Albinism Type 1 (OA1, Nettleship-Falls OA), which results in hypopigmentation of the eyes and loss of visual acuity due to disrupted visual system development and function. In pigment cells of the skin, loss of functional GPR143 results in abnormally large melanosomes (organelles in which pigment is produced). Studies have shown that the receptor is localized internally, including at the melanosomal membrane, where it may function to regulate melanosome size and/or facilitate protein trafficking to the melanosome through the endolysosomal system. Numerous additional roles have been proposed for GPR143 in determining cancer predisposition, regulation of blood pressure, development of macular degeneration and signaling in the brain, which we will briefly describe as well as potential ligands that have been identified. Furthermore, GPR143 is a promiscuous receptor that has been shown to interact with multiple other melanosomal proteins and GPCRs, which strongly suggests that this orphan receptor is likely involved in many different physiological actions.

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.

Exploring the Activation Process of the β2AR-Gs Complex

G-Protein-coupled receptors (GPCRs) belong to an important family of integral membrane receptor proteins that are essential for a variety of transmembrane signaling process, such as vision, olfaction, and hormone responses. They are also involved in many human diseases (Alzheimer's, heart diseases, etc.) and are therefore common drug targets. Thus, understanding the details of the GPCR activation process is a task of major importance. Various types of crystal structures of GPCRs have been solved either at stable end-point states or at possible intermediate states. However, the detailed mechanism of the activation process is still poorly understood. For example, it is not completely clear when the nucleotide release from the G protein occurs and how the key residues on α5 contribute to the coupling process and further affect the binding specificity. In this work we show by free energy analysis that the guanosine diphosphate (GDP) molecule could be released from the Gs protein when the binding cavity is half open. This occurs during the transition to the Gs open state, which is the rate-determining step in the system conformational change. We also account for the experimentally observed slow-down effects by the change of the reaction barriers after mutations. Furthermore, we identify potential key residues on α5 and validated their significance by site-directed mutagenesis, which illustrates that computational works have predictive value even for complex biophysical systems. The methodology of the current work may be applied to other biophysical systems of interest.

Dissecting the Conformational Free Energy of a Small Peptide in Solution

The free energy surface of a small peptide was analyzed based on an unbiased microsecond molecular dynamics simulation. The peptide sampled disordered conformational ensembles of distinct compactness, and its free energy was decomposed into separate contributions from the intramolecular potential energy, conformational entropy, and solvation free energy. The latter was further broken down into enthalpic and entropic contributions due to peptide-water and water-water interactions. This decomposition was enabled by a generalized linear response relation between the peptide-water interaction energy and the solvation free energy, which was empirically parametrized by explicit solvation free energy calculations for representative peptide conformations. This full dissection of the peptide free energy identifies individual contributions that stabilize and destabilize compact and extended peptide conformational ensembles and reveals the origin of a free energy barrier associated with transitions between them.

Coevolution underlies GPCR-G protein selectivity and functionality

G protein-coupled receptors (GPCRs) regulate diverse physiological events, which makes them as the major targets for many approved drugs. G proteins are downstream molecules that receive signals from GPCRs and trigger cell responses. The GPCR-G protein selectivity mechanism on how they properly and timely interact is still unclear. Here, we analyzed model GPCRs (i.e. HTR, DAR) and Gα proteins with a coevolutionary tool, statistical coupling analysis. The results suggested that 5-hydroxytryptamine receptors and dopamine receptors have common conserved and coevolved residues. The Gα protein also have conserved and coevolved residues. These coevolved residues were implicated in the molecular functions of the analyzed proteins. We also found specific coevolving pairs related to the selectivity between GPCR and G protein were identified. We propose that these results would contribute to better understandings of not only the functional residues of GPCRs and Gα proteins but also GPCR-G protein selectivity mechanisms.

The Knockout for G Protein-Coupled Receptor-Like PfSR25 Increases the Susceptibility of Malaria Parasites to the Antimalarials Lumefantrine and Piperaquine but Not to Medicine for Malaria Venture Compounds

Previously we have reported that the G protein-coupled receptor (GPCR)-like PfSR25 in Plasmodium falciparum is a potassium (K⁺) sensor linked to intracellular calcium signaling and that knockout parasites (PfSR25-) are more susceptible to oxidative stress and antimalarial compounds. Here, we explore the potential role of PfSR25 in susceptibility to the antimalarial compounds atovaquone, chloroquine, dihydroartemisinin, lumefantrine, mefloquine, piperaquine, primaquine, and pyrimethamine and the Medicine for Malaria Venture (MMV) compounds previously described to act on egress/invasion (MMV006429, MMV396715, MMV019127, MMV665874, MMV665878, MMV665785, and MMV66583) through comparative assays with PfSR25- and 3D7 parasite strains, using flow cytometry assays. The IC₅₀ and IC₉₀ results show that lumefantrine and piperaquine have greater activity on the PfSR25- parasite strain when compared to 3D7. For MMV compounds, we found no differences between the strains except for the compound MMV665831, which we used to investigate the store-operated calcium entry (SOCE) mechanism. The results suggest that PfSR25 may be involved in the mechanism of action of the antimalarials lumefantrine and piperaquine. Our data clearly show that MMV665831 does not affect calcium entry in parasites after we depleted their internal calcium pools with thapsigargin. The results demonstrated here shed light on new possibilities on the antimalarial mechanism, bringing evidence of the involvement of the GPCR-like PfSR25.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

Lysolipid Chain Length Switches Agonistic to Antagonistic G Protein-Coupled Receptor Modulation

Activation of lysolipid-sensitive G protein-coupled receptors (GPCR) depends not only on lysolipid class but also on the length and degree of saturation of their respective hydrophobic tails. Positive regulation of these signaling networks caused by the lipid chain length specificity of upstream phospholipases is firmly established. Nonagonistic lysolipid homologues, featuring incompatible lipid tails, have been suggested to indirectly modulate GPCR signaling by delaying agonist catabolism. Nonetheless, recent results seem inconsistent with this hypothesis. Utilizing a simplified lysolipid-GPCR signaling assay based on the established lysophosphatidylglucoside-GPR55 signaling axis in primary sensory neurons, we demonstrate that short-chain ligand homologues directly modulate receptor activation via a potent competitive antagonistic activity. Considering the well-documented tissue-specific concentration of lysolipid homologues, we propose that endogenous lysolipids with insufficient chain length for stable receptor activation exert an antagonistic activity, effectively representing a negative control mechanism for GPCR-associated lysolipid signaling.

Exploring the activation pathway and Gi-coupling specificity of the μ-opioid receptor

Understanding the activation mechanism of the μ-opioid receptor (μ-OR) and its selective coupling to the inhibitory G protein (Gi) is vital for pharmaceutical research aimed at finding treatments for the opioid overdose crisis. Many attempts have been made to understand the mechanism of the μ-OR activation, following the elucidation of new crystal structures such as the antagonist- and agonist-bound μ-OR. However, the focus has not been placed on the underlying energetics and specificity of the activation process. An energy-based picture would not only help to explain this coupling but also help to explore why other possible options are not common. For example, one would like to understand why μ-OR is more selective to Gi than a stimulatory G protein (Gs). Our study used homology modeling and a coarse-grained model to generate all of the possible "end states" of the thermodynamic cycle of the activation of μ-OR. The end points were further used to generate reasonable intermediate structures of the receptor and the Gi to calculate two-dimensional free energy landscapes. The results of the landscape calculations helped to propose a plausible sequence of conformational changes in the μ-OR and Gi system and for exploring the path that leads to its activation. Furthermore, in silico alanine scanning calculations of the last 21 residues of the C terminals of Gi and Gs were performed to shed light on the selective binding of Gi to μ-OR. Overall, the present work appears to demonstrate the potential of multiscale modeling in exploring the action of G protein-coupled receptors.

Structural insights into emergent signaling modes of G protein-coupled receptors

G protein-coupled receptors (GPCRs) represent the largest family of cell membrane proteins, with >800 GPCRs in humans alone, and recognize highly diverse ligands, ranging from photons to large protein molecules. Very important to human medicine, GPCRs are targeted by about 35% of prescription drugs. GPCRs are characterized by a seven-transmembrane α-helical structure, transmitting extracellular signals into cells to regulate major physiological processes via heterotrimeric G proteins and β-arrestins. Initially viewed as receptors whose signaling via G proteins is delimited to the plasma membrane, it is now recognized that GPCRs signal also at various intracellular locations, and the mechanisms and (patho)physiological relevance of such signaling modes are actively investigated. The propensity of GPCRs to adopt different signaling modes is largely encoded in the structural plasticity of the receptors themselves and of their signaling complexes. Here, we review emerging modes of GPCR signaling via endosomal membranes and the physiological implications of such signaling modes. We further summarize recent structural insights into mechanisms of GPCR activation and signaling. We particularly emphasize the structural mechanisms governing the continued GPCR signaling from endosomes and the structural aspects of the GPCR resensitization mechanism and discuss the recently uncovered and important roles of lipids in these processes.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

Role of Actin Cytoskeleton in Dynamics and Function of the Serotonin1A Receptor

G protein-coupled receptors (GPCRs) are important membrane proteins in higher eukaryotes that carry out a vast array of cellular signaling and act as major drug targets. The serotonin_1A receptor is a prototypical member of the GPCR family and is implicated in neuropsychiatric disorders such as anxiety and depression, besides serving as an important drug target. With an overall goal of exploring the functional consequence of altered receptor dynamics, in this work, we probed the role of the actin cytoskeleton in the dynamics, ligand binding, and signaling of the serotonin_1A receptor. We monitored receptor dynamics utilizing single particle tracking, which provides information on relative distribution of receptors in various diffusion modes in addition to diffusion coefficient. We show here that the short-term diffusion coefficient of the receptor increases upon actin destabilization by cytochalasin D. In addition, analysis of individual trajectories shows that there are changes in relative populations of receptors undergoing various types of diffusion upon actin destabilization. The release of dynamic constraint was evident by an increase in the radius of confinement of the receptor upon actin destabilization. The functional implication of such actin destabilization was manifested as an increase in specific agonist binding and downstream signaling, monitored by measuring reduction in cellular cAMP levels. These results bring out the interdependence of GPCR dynamics with cellular signaling.

Comprehensive structural modeling and preparation of human 5-HT2A G-protein coupled receptor in functionally active form

The serotonin 2A receptor (5-HT_2A R) is an important member of the G-protein coupled receptor (GPCR) family involved in an array of neuromodulatory functions. Although the high-resolution structures of truncated versions of GPCRs, captured in ligand-bound conformational states, are available, the structures lack several functional regions, which have crucial roles in receptor response. Here, in order to understand the structure and dynamics of the ligand-free form of the receptor, we have performed meticulous modeling of the 5-HT_2A R with the third intracellular loop (ICL3). Our analyses revealed that the ligand-free ground state structure of 5-HT_2A R has marked distinction with ligand-bound conformations of 5-HT₂ subfamily proteins and exhibits extensive backbone flexibility across the loop regions, suggesting the importance of purifying the receptor in its native form for further studies. Hence, we have standardized a strategy that efficiently increases the expression of 5-HT_2A R by infecting Sf9 cells with a very low multiplicity of infection of baculovirus in conjunction with production boost additive and subsequently, purify the full-length receptor. Furthermore, we have optimized the selective over-expression of glycosylated and nonglycosylated forms of the receptor merely by switching the postinfection growth time, a method that has not been reported earlier.

Kinetic asymmetry allows macromolecular catalysts to drive an information ratchet

Molecular machines carry out their function by equilibrium mechanical motions in environments that are far from thermodynamic equilibrium. The mechanically equilibrated character of the trajectories of the macromolecule has allowed development of a powerful theoretical description, reminiscent of Onsager’s trajectory thermodynamics, that is based on the principle of microscopic reversibility. Unlike the situation at thermodynamic equilibrium, kinetic parameters play a dominant role in determining steady-state concentrations away from thermodynamic equilibrium, and kinetic asymmetry provides a mechanism by which chemical free-energy released by catalysis can drive directed motion, molecular adaptation, and self-assembly. Several examples drawn from the recent literature, including a catenane-based chemically driven molecular rotor and a synthetic molecular assembler or pump, are discussed.

The mechanism by which macromolecular catalysts use energy from exergonic reactions to move, adapt, and assemble has been unclear. In this Perspective article, R. Dean Astumian shows that in addition to disequilibrium of the catalyzed reaction, kinetic asymmetry is the essential feature required to drive non-equilibrium response by an information ratchet mechanism.

A free-energy landscape for the glucagon-like peptide 1 receptor GLP1R

G-protein-coupled receptors (GPCRs) are among the most important receptors in human physiology and pathology. They serve as master regulators of numerous key processes and are involved in as well as cause debilitating diseases. Consequently, GPCRs are among the most attractive targets for drug design and pharmaceutical interventions (>30% of drugs on the market). The glucagon-like peptide 1 (GLP-1) hormone receptor GLP1R is closely involved in insulin secretion by pancreatic β-cells and constitutes a major druggable target for the development of anti-diabetes and obesity agents. GLP1R structure was recently solved, with ligands, allosteric modulators and as part of a complex with its cognate G protein. However, the translation of this structural data into structure/function understanding remains limited. The current study functionally characterizes GLP1R with special emphasis on ligand and cellular partner binding interactions and presents a free-energy landscape as well as a functional model of the activation cycle of GLP1R. Our results should facilitate a deeper understanding of the molecular mechanism underlying GLP1R activation, forming a basis for improved development of targeted therapeutics for diabetes and related disorders.

Studying GPCR conformational dynamics by single molecule fluorescence

Over the last decades, G protein coupled receptors (GPCRs) have experienced a tremendous amount of attention, which has led to a boost of structural and pharmacological insights on this large membrane protein superfamily involved in various essential physiological functions. Recently, evidence has emerged that, rather than being activated by ligands in an on/off manner switching from an inactive to an active state, GPCRs exhibit high structural flexibility in the absence and even in the presence of ligands. So far the physiological as well as pharmacological impact of this structural flexibility remains largely unexplored albeit its potential role in precisely fine-tuning receptor function and regulating the specificity of signal transduction into the cell. By complementing other biophysical approaches, single molecule fluorescence (SMF) offers the advantage of monitoring structural dynamics in biomolecules in real-time, with minimal structural invasiveness and in the context of complex biological environments. In this review a general introduction to GPCR structural dynamics is given followed by a presentation of SMF methods used to explore them. Particular attention is paid to single molecule Förster resonance energy transfer (smFRET), a key method to measure actual distance changes between two probes, and highlight conformational changes occurring at timescales relevant for protein conformational movements. The available literature reporting on GPCR structural dynamics by SMF is discussed with a focus on the newly gained biological insights on receptor activation and signaling, in particular for the β2 adrenergic and the metabotropic glutamate receptors.

Dynamics and Mechanisms in the Recruitment and Transference of Histone Chaperone CIA/ASF1

The recruitment and transference of proteins through protein-protein interactions is a general process involved in various biological functions in cells. Despite the importance of this general process, the dynamic mechanism of how proteins are recruited and transferred from one interacting partner to another remains unclear. In this study, we investigated the dynamic mechanisms of recruitment and translocation of histone chaperone CIA/ASF1 for nucleosome disassembly by exploring the conformational space and the free energy profile of unbound DBD(CCG1) and CIA/ASF1-bound DBD(CCG1) systems through extensive molecular dynamics simulations. It was found that there exists three metastable conformational states for DBD(CCG1), an unbound closed state, a CIA/ASF1-bound half-open state, and an open state. The free energy landscape shows that the closed state and the half-open state are separated by a high free energy barrier, while the half-open state and the open state are connected with a moderate free energy increase. The high free energy barrier between the closed and half-open states explains why DBD(CCG1) can recruit CIA/ASF1 and remain in the binding state during the transportation. In addition, the asymmetric binding of CIA/ASF1 on DBD(CCG1) allows DBD(CCG1) to adopt the open state by moving one of its two domains, such that the exposed domain of DBD(CCG1) is able to recognize the acetylated histone H4 tails. As such, CIA/ASF1 has a chance to translocate from DBD(CCG1) to histone, which is also facilitated by the moderate energy increase from the bound half-open state to the open state of DBD(CCG1). These findings suggest that the recruitment and transference of histone chaperone CIA/ASF1 is highly favored by its interaction with DBD(CCG1) via conformational selection and asymmetric binding, which may represent a general mechanism of similar biological processes.

Precision medicine review: rare driver mutations and their biophysical classification

How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.