Unifying family A GPCR theories of activation

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

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

Elucidating the Activation Mechanism of the Proton-sensing GPR68 Receptor

GPR68 is a proton-sensing G-protein Coupled Receptor (GPCR) involved in a variety of physiological processes and disorders including neoplastic pathologies. While GPR68 and few other GPCRs have been shown to be activated by a decrease in the extracellular pH, the molecular mechanism of their activation remains largely unknown. In this work, we used a combined computational and in vitro approach to provide new insight into the activation mechanism of the receptor. Molecular Dynamics simulations of GPR68 were used to model the changes in residue interactions and motions triggered by pH. Global and local rearrangements consistent with partial activation were observed upon protonation of the inactive state. Selected extracellular histidine and transmembrane acidic residues were found to have significantly upshifted pKa values during the simulations, consistently with their previously hypothesised role in activation through changes in protonation state. Moreover, a novel pairing between histidine and acidic residues in the extracellular region was highlighted by both sequence analyses and simulation data and tested through site-directed mutagenesis. At last, we identified a previously unknown hydrophobic lock in the extracellular region that might stabilise the inactive conformation and regulate the transition to the active state.

Balancing G protein selectivity and efficacy in the adenosine A2A receptor

The adenosine A2A receptor (A2AR) engages several G proteins, notably Go and its cognate Gs protein. This coupling promiscuity is facilitated by a dynamic ensemble, revealed by 19F nuclear magnetic resonance imaging of A2AR and G protein. Two transmembrane helix 6 (TM6) activation states, formerly associated with partial and full agonism, accommodate the differing volumes of Gs and Go. While nucleotide depletion biases TM7 toward a fully active state in A2AR-Gs, A2AR-Go is characterized by a dynamic inactive/intermediate fraction. Molecular dynamics simulations reveal that the NPxxY motif, a highly conserved switch, establishes a unique configuration in the A2AR-Go complex, failing to stabilize the helix-8 interface with Gs, and adoption of the active state. The resulting TM7 dynamics hamper G protein coupling, suggesting kinetic gating may be responsible for reduced efficacy in the noncognate G protein complex. Thus, dual TM6 activation states enable greater diversity of coupling partners while TM7 dynamics dictate coupling efficacy.

G protein-coupled receptors and traditional Chinese medicine: new thinks for the development of traditional Chinese medicine

G protein-coupled receptors (GPCRs) widely exist in vivo and participate in many physiological processes, thus emerging as important targets for drug development. Approximately 30% of the Food and Drug Administration (FDA)-approved drugs target GPCRs. To date, the 'one disease, one target, one molecule' strategy no longer meets the demands of drug development. Meanwhile, small-molecule drugs account for 60% of FDA-approved drugs. Traditional Chinese medicine (TCM) has garnered widespread attention for its unique theoretical system and treatment methods. TCM involves multiple components, targets and pathways. Centered on GPCRs and TCM, this paper discusses the similarities and differences between TCM and GPCRs from the perspectives of syndrome of TCM, the consistency of TCM's multi-component and multi-target approaches and the potential of GPCRs and TCM in the development of novel drugs. A novel strategy, 'simultaneous screening of drugs and targets', was proposed and applied to the study of GPCRs. We combine GPCRs with TCM to facilitate the modernisation of TCM, provide valuable insights into the rational application of TCM and facilitate the research and development of novel drugs. This study offers theoretical support for the modernisation of TCM and introduces novel ideas for development of safe and effective drugs.

Developing insulin-like peptide 5-based antagonists for the G protein-coupled receptor, RXFP4

Human insulin-like peptide 5 (INSL5) is a gut hormone produced by colonic L-cells, and its biological functions are mediated by Relaxin Family Peptide Receptor 4 (RXFP4). Our preliminary data indicated that RXFP4 agonists are potential drug leads for the treatment of constipation. More recently, we designed and developed a novel RXFP4 antagonist, A13-nR that was shown to block agonist-induced activity in cells and animal models. We showed that A13-nR was able to block agonist-induced increases in colon motility in mice of both genders that express the receptor, RXFP4. Our data also showed that colorectal propulsion induced by intracolonic administration of short-chain fatty acids was antagonized by A13-nR. Therefore, A13-nR is an important research tool and potential drug lead for the treatment of colon motility disorders, such as bacterial diarrhea. However, A13-nR acted as a partial agonist at high concentrations in vitro and demonstrated modest antagonist potency (∼35 nM). Consequently, the primary objective of this study is to pinpoint novel modifications to A13-nR that eliminate partial agonist effects while preserving or augmenting antagonist potency. In this work, we detail the creation of a series of A13-nR-modified analogues, among which analogues 3, 4, and 6 demonstrated significantly improved RXFP4 affinity (∼3 nM) with reduced partial agonist activity, enhanced antagonist potency (∼10 nM) and maximum agonist inhibition (∼80 %) when compared with A13-nR. These compounds have potential as candidates for further preclinical evaluations, marking a significant stride toward innovative therapeutics for colon motility disorders.

GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

Snapshot of the cannabinoid receptor 1-arrestin complex unravels the biased signaling mechanism

Cannabis activates the cannabinoid receptor 1 (CB1), which elicits analgesic and emotion regulation benefits, along with adverse effects, via Gi and β-arrestin signaling pathways. However, the lack of understanding of the mechanism of β-arrestin-1 (βarr1) coupling and signaling bias has hindered drug development targeting CB1. Here, we present the high-resolution cryo-electron microscopy structure of CB1-βarr1 complex bound to the synthetic cannabinoid MDMB-Fubinaca (FUB), revealing notable differences in the transducer pocket and ligand-binding site compared with the Gi protein complex. βarr1 occupies a wider transducer pocket promoting substantial outward movement of the TM6 and distinctive twin toggle switch rearrangements, whereas FUB adopts a different pose, inserting more deeply than the Gi-coupled state, suggesting the allosteric correlation between the orthosteric binding pocket and the partner protein site. Taken together, our findings unravel the molecular mechanism of signaling bias toward CB1, facilitating the development of CB1 agonists.

Single-molecule visualization of human A2A adenosine receptor activation by a G protein and constitutively activating mutations

Mutations that constitutively activate G protein-coupled receptors (GPCRs), known as constitutively activating mutations (CAMs), modify cell signaling and interfere with drugs, resulting in diseases with limited treatment options. We utilize fluorescence imaging at the single-molecule level to visualize the dynamic process of CAM-mediated activation of the human A2A adenosine receptor (A2AAR) in real time. We observe an active-state population for all CAMs without agonist stimulation. Importantly, activating mutations significantly increase the population of an intermediate state crucial for receptor activation, notably distinct from the addition of a partner G protein. Activation kinetics show that while CAMs increase the frequency of transitions to the intermediate state, mutations altering sodium sensitivity increase transitions away from it. These findings indicate changes in GPCR function caused by mutations may be predicted based on whether they favor or disfavor formation of an intermediate state, providing a framework for designing receptors with altered functions or therapies that target intermediate states.

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

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

Dynamic Insights into the Self-Activation Pathway and Allosteric Regulation of the Orphan G-Protein-Coupled Receptor GPR52

Within over 800 members of G-protein-coupled receptors, there are numerous orphan receptors whose endogenous ligands are largely unknown, providing many opportunities for novel drug discovery. However, the lack of an in-depth understanding of the intrinsic working mechanism for orphan receptors severely limits the related rational drug design. The G-protein-coupled receptor 52 (GPR52) is a unique orphan receptor that constitutively increases cellular 5'-cyclic adenosine monophosphate (cAMP) levels without binding any exogenous agonists and has been identified as a promising therapeutic target for central nervous system disorders. Although recent structural biology studies have provided snapshots of both active and inactive states of GPR52, the mechanism of the conformational transition between these states remains unclear. Here, an acceptable self-activation pathway for GPR52 was proposed through 6 μs Gaussian accelerated molecular dynamics (GaMD) simulations, in which the receptor spontaneously transitions from the active state to that matching the inactive crystal structure. According to the three intermediate states of the receptor obtained by constructing a reweighted potential of mean force, how the allosteric regulation occurs between the extracellular orthosteric binding pocket and the intracellular G-protein-binding site is revealed. Combined with the independent gradient model, several important microswitch residues and the allosteric communication pathway that directly links the two regions are both identified. Transfer entropy calculations not only reveal the complex allosteric signaling within GPR52 but also confirm the unique role of ECL2 in allosteric regulation, which is mutually validated with the results of GaMD simulations. Overall, this work elucidates the allosteric mechanism of GPR52 at the atomic level, providing the most detailed information to date on the self-activation of the orphan receptor.

GPCRana: A web server for quantitative analysis of GPCR structures

G protein-coupled receptors (GPCRs) attract tremendous attention from both industrial and academic researchers with currently over 900 released structures. Structural analysis is widely used to understand receptor functionality and pharmacology, but more user-friendly tools are needed. Residue-residue contact score (RRCS) is an atomic distance-based method that allows a quantitative description of GPCR structures. Here, we present GPCRana, a web server that provides a user-friendly interface to analyze GPCR structures. After uploading selected structures, GPCRana immediately generates a comprehensive report covering four aspects: (i) RRCS for all residue pairs incorporated with real-time 3D visualization; (ii) ligand-receptor interactions; (iii) activation pathway analysis; and (iv) RRCS_TMs that indicates the global movements of transmembrane helices. Moreover, conformational changes between two structures can be analyzed. Applying GPCRana on AlphaFold2-predicted models reveals differentiated inter-helical packing forms in a receptor-dependent manner. Our web server offers a fast and precise way to study GPCR structures and is freely available at http://gpcranalysis.com/#/.

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

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

CovET: A covariation-evolutionary trace method that identifies protein structure-function modules

Measuring the relative effect that any two sequence positions have on each other may improve protein design or help better interpret coding variants. Current approaches use statistics and machine learning but rarely consider phylogenetic divergences which, as shown by Evolutionary Trace studies, provide insight into the functional impact of sequence perturbations. Here, we reframe covariation analyses in the Evolutionary Trace framework to measure the relative tolerance to perturbation of each residue pair during evolution. This approach (CovET) systematically accounts for phylogenetic divergences: at each divergence event, we penalize covariation patterns that belie evolutionary coupling. We find that while CovET approximates the performance of existing methods to predict individual structural contacts, it performs significantly better at finding structural clusters of coupled residues and ligand binding sites. For example, CovET found more functionally critical residues when we examined the RNA recognition motif and WW domains. It correlates better with large-scale epistasis screen data. In the dopamine D2 receptor, top CovET residue pairs recovered accurately the allosteric activation pathway characterized for Class A G protein-coupled receptors. These data suggest that CovET ranks highest the sequence position pairs that play critical functional roles through epistatic and allosteric interactions in evolutionarily relevant structure-function motifs. CovET complements current methods and may shed light on fundamental molecular mechanisms of protein structure and function.

Structural Basis for Agonistic Activity and Selectivity toward Melatonin Receptors hMT1 and hMT2

Glaucoma, a major ocular neuropathy originating from a progressive degeneration of retinal ganglion cells, is often associated with increased intraocular pressure (IOP). Daily IOP fluctuations are physiologically influenced by the antioxidant and signaling activities of melatonin. This endogenous modulator has limited employment in treating altered IOP disorders due to its low stability and bioavailability. The search for low-toxic compounds as potential melatonin agonists with higher stability and bioavailability than melatonin itself could start only from knowing the molecular basis of melatonergic activity. Thus, using a computational approach, we studied the melatonin binding toward its natural macromolecular targets, namely melatonin receptors 1 (MT1) and 2 (MT2), both involved in IOP signaling regulation. Besides, agomelatine, a melatonin-derivative agonist and, at the same time, an atypical antidepressant, was also included in the study due to its powerful IOP-lowering effects. For both ligands, we evaluated both stability and ligand positioning inside the orthosteric site of MTs, mapping the main molecular interactions responsible for receptor activation. Affinity values in terms of free binding energy (ΔG_bind) were calculated for the selected poses of the chosen compounds after stabilization through a dynamic molecular docking protocol. The results were compared with experimental in vivo effects, showing a higher potency and more durable effect for agomelatine with respect to melatonin, which could be ascribed both to its higher affinity for hMT2 and to its additional activity as an antagonist for the serotonin receptor 5-HT2c, in agreement with the in silico results.

Homology Modeling of the G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are therapeutically important family of membrane proteins. Despite growing number of experimental structures available for GPCRs, homology modeling remains a relevant method for studying these receptors and for discovering new ligands for them. Here we describe the state-of-the-art methods for modeling GPCRs, starting from template selection, through fine-tuning sequence alignment to model refinement.

AlphaFold2 versus experimental structures: evaluation on G protein-coupled receptors

As important drug targets, G protein-coupled receptors (GPCRs) play pivotal roles in a wide range of physiological processes. Extensive efforts of structural biology have been made on the study of GPCRs. However, a large portion of GPCR structures remain unsolved due to structural instability. Recently, AlphaFold2 has been developed to predict structure models of many functionally important proteins including all members of the GPCR family. Herein we evaluated the accuracy of GPCR structure models predicted by AlphaFold2. We revealed that AlphaFold2 could capture the overall backbone features of the receptors. However, the predicted models and experimental structures were different in many aspects including the assembly of the extracellular and transmembrane domains, the shape of the ligand-binding pockets, and the conformation of the transducer-binding interfaces. These differences impeded the use of predicted structure models in the functional study and structure-based drug design of GPCRs, which required reliable high-resolution structural information.

Biased agonists differentially modulate the receptor conformation ensembles in Angiotensin II type 1 receptor

The structural features that contribute to the efficacy of biased agonists targeting G protein-coupled receptors (GPCRs) towards G proteins or β-arrestin (β-arr) signaling pathways is nebulous, although such knowledge is critical in designing biased ligands. The dynamics of the agonist-GPCR complex is one of the critical factors in determining agonist bias. Angiotensin II type I receptor (AT1R) is an ideal model system to study the molecular basis of bias since it has multiple β-arr2 and Gq protein biased agonists as well as experimentally solved three dimensional structures. Using Molecular Dynamics (MD) simulations for the Angiotensin II type I receptor (AT1R) bound to ten different agonists, we infer that the agonist bound receptor samples conformations with different relative weights, from both the inactive and active state ensembles of the receptor. This concept is perhaps extensible to other class A GPCRs. Such a weighted mixed ensemble recapitulates the inter-residue distance distributions measured for different agonists bound AT1R using DEER experiments. The ratio of the calculated relative strength of the allosteric communication to β-arr2 vs Gq coupling sites scale similarly to the experimentally measured bias factors. Analysis of the inter-residue distance distributions of the activation microswitches involved in class A GPCR activation suggests that β-arr2 biased agonists turn on different combination of microswitches with different relative strengths of activation. We put forth a model that activation microswitches behave like rheostats that tune the relative efficacy of the biased agonists toward the two signaling pathways. Finally, based on our data we propose that the agonist specific residue contacts in the binding site elicit a combinatorial response in the microswitches that in turn differentially modulate the receptor conformation ensembles resulting in differences in coupling to Gq and β-arrestin.

Structure of full-length TSH receptor in complex with antibody K1-70™

Determination of the full-length thyroid-stimulating hormone receptor (TSHR) structure by cryo-electron microscopy (cryo-EM) is described. The TSHR complexed with human monoclonal TSHR autoantibody K1-70™ (a powerful inhibitor of TSH action) was detergent solubilised, purified to homogeneity and analysed by cryo-EM. The structure (global resolution 3.3 Å) is a monomer with all three domains visible: leucine-rich domain (LRD), hinge region (HR) and transmembrane domain (TMD). The TSHR extracellular domain (ECD, composed of the LRD and HR) is positioned on top of the TMD extracellular surface. Extensive interactions between the TMD and ECD are observed in the structure, and their analysis provides an explanation of the effects of various TSHR mutations on TSHR constitutive activity and on ligand-induced activation. K1-70™ is seen to be well clear of the lipid bilayer. However, superimposition of M22™ (a human monoclonal TSHR autoantibody which is a powerful stimulator of the TSHR) on the cryo-EM structure shows that it would clash with the bilayer unless the TSHR HR rotates upwards as part of the M22™ binding process. This rotation could have an important role in TSHR stimulation by M22™ and as such provides an explanation as to why K1-70™ blocks the binding of TSH and M22™ without activating the receptor itself.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

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.

Engineered Allosteric Regulation of Protein Function

Allosteric regulation of proteins has been utilized to study various aspects of cell signaling, from unicellular events to organism-wide phenotypes. However, traditional methods of allosteric regulation, such as constitutively active mutants and inhibitors, lack tight spatiotemporal control. This often leads to unintended signaling consequences that interfere with data interpretation. To overcome these obstacles, researchers employed protein engineering approaches that enable tight control of protein function through allosteric mechanisms. These methods provide high specificity as well as spatial and temporal precision in regulation of protein activity in vitro and in vivo. In this review, we focus on the recent advancements in engineered allosteric regulation and discuss the various bioengineered allosteric techniques available now, from chimeric GPCRs to chemogenetic and optogenetic switches. We highlight the benefits and pitfalls of each of these techniques as well as areas in which future improvements can be made. Additionally, we provide a brief discussion on implementation of engineered allosteric regulation approaches, demonstrating that these tools can shed light on elusive biological events and have the potential to be utilized in precision medicine.

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.

An Interpretable Convolutional Neural Network Framework for Analyzing Molecular Dynamics Trajectories: a Case Study on Functional States for G-Protein-Coupled Receptors

Molecular dynamics (MD) simulations have made great contribution to revealing structural and functional mechanisms for many biomolecular systems. However, how to identify functional states and important residues from vast conformation space generated by MD remains challenging; thus an intelligent navigation is highly desired. Despite intelligent advantages of deep learning exhibited in analyzing MD trajectory, its black-box nature limits its application. To address this problem, we explore an interpretable convolutional neural network (CNN)-based deep learning framework to automatically identify diverse active states from the MD trajectory for G-protein-coupled receptors (GPCRs), named the ICNNMD model. To avoid the information loss in representing the conformation structure, the pixel representation is introduced, and then the CNN module is constructed to efficiently extract features followed by a fully connected neural network to realize the classification task. More importantly, we design a local interpretable model-agnostic explanation interpreter for the classification result by local approximation with a linear model, through which important residues underlying distinct active states can be quickly identified. Our model showcases higher than 99% classification accuracy for three important GPCR systems with diverse active states. Notably, some important residues in regulating different biased activities are successfully identified, which are beneficial to elucidating diverse activation mechanisms for GPCRs. Our model can also serve as a general tool to analyze MD trajectory for other biomolecular systems. All source codes are freely available at https://github.com/Jane-Liu97/ICNNMD for aiding MD studies.

Allosteric modulation of dopamine D2L receptor in complex with Gi1 and Gi2 proteins: the effect of subtle structural and stereochemical ligand modifications

Background

Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays one of the hot topics in drug discovery. In particular, allosteric modulators of D₂ receptor have been proposed as potential modern therapeutics to treat schizophrenia and Parkinson’s disease.

Methods

To address some subtle structural and stereochemical aspects of allosteric modulation of D₂ receptor, we performed extensive in silico studies of both enantiomers of two compounds (compound 1 and compound 2), and one of them (compound 2) was synthesized as a racemate in-house and studied in vitro.

Results

Our molecular dynamics simulations confirmed literature reports that the R enantiomer of compound 1 is a positive allosteric modulator of the D_2L receptor, while its S enantiomer is a negative allosteric modulator. Moreover, based on the principal component analysis (PCA), we hypothesized that both enantiomers of compound 2 behave as silent allosteric modulators, in line with our in vitro studies. PCA calculations suggest that the most pronounced modulator-induced receptor rearrangements occur at the transmembrane helix 7 (TM7). In particular, TM7 bending at the conserved P7.50 and G7.42 was observed. The latter resides next to the Y7.43, which is a significant part of the orthosteric binding site. Moreover, the W7.40 conformation seems to be affected by the presence of the positive allosteric modulator.

Conclusions

Our work reveals that allosteric modulation of the D_2L receptor can be affected by subtle ligand modifications. A change in configuration of a chiral carbon and/or minor structural modulator modifications are solely responsible for the functional outcome of the allosteric modulator.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s43440-021-00352-x.

Delineating the activation mechanism and conformational landscape of a class B G protein-coupled receptor glucagon receptor

Class B G protein-coupled receptors (GPCRs) are important targets in the treatment of metabolic syndrome and diabetes. Although multiple structures of class B GPCRs-G protein complexes have been elucidated, the detailed activation mechanism of the receptors remains unclear. Here, we combine Gaussian accelerated molecular dynamics simulations and Markov state models (MSM) to investigate the activation mechanism of a canonical class B GPCR, human glucagon receptor-GCGR, including the negative allosteric modulator-bound inactive state, the agonist glucagon-bound active state, and both glucagon- and Gs-bound fully active state. The free-energy landscapes of GCGR show the conformational ensemble consisting of three activation-associated states: inactive, active, and fully active. The structural analysis indicates the high dynamics of GCGR upon glucagon binding with both active and inactive conformations in the ensemble. Significantly, the H8 and TM6 exhibits distinct features from the inactive to the active states. The additional simulations demonstrate the role of H8 in the recruitment of Gs. Gs binding presents a crucial function of stabilizing the glucagon binding site and MSM highlights the absolute requirement of Gs to help the GCGR reach the fully active state. Together, our results reveal the detailed activation mechanism of GCGR from the view of conformational dynamics.

New insights into the molecular mechanism of rhodopsin retinitis pigmentosa from the biochemical and functional characterization of G90V, Y102H and I307N mutations

Mutations in the photoreceptor protein rhodopsin are known as one of the leading causes of retinal degeneration in humans. Two rhodopsin mutations, Y102H and I307N, obtained in chemically mutagenized mice, are currently the subject of increased interest as relevant models for studying the process of retinal degeneration in humans. Here, we report on the biochemical and functional characterization of the structural and functional alterations of these two rhodopsin mutants and we compare them with the G90V mutant previously analyzed, as a basis for a better understanding of in vivo studies. This mechanistic knowledge is fundamental to use it for developing novel therapeutic approaches for the treatment of inherited retinal degeneration in retinitis pigmentosa. We find that Y102H and I307N mutations affect the inactive–active equilibrium of the receptor. In this regard, the mutations reduce the stability of the inactive conformation but increase the stability of the active conformation. Furthermore, the initial rate of the functional activation of transducin, by the I307N mutant is reduced, but its kinetic profile shows an unusual increase with time suggesting a profound effect on the signal transduction process. This latter effect can be associated with a change in the flexibility of helix 7 and an indirect effect of the mutation on helix 8 and the C-terminal tail of rhodopsin, whose potential role in the functional activation of the receptor has been usually underestimated. In the case of the Y102H mutant, the observed changes can be associated with conformational alterations affecting the folding of the rhodopsin intradiscal domain, and its presumed involvement in the retinal binding process by the receptor.

Recurrent high-impact mutations at cognate structural positions in class A G protein-coupled receptors expressed in tumors

GPCRs and GPCR pathways are increasingly being implicated in human malignancies, placing them among the most promising cancer drug candidates. Our results reveal enrichment of highly impactful, recurrent GPCR mutations within cancers. We found that cognate mutations in selected class A GPCRs have deleterious effects on signaling function. The results also suggest that olfactory receptors, often considered inconsequential, display a nonrandom mutation pattern in tumors in which they are expressed. These findings support the idea that protein paralogs can act in parallel as members of an onco-group.

G protein-coupled receptors (GPCRs) are the largest family of human proteins. They have a common structure and, signaling through a much smaller set of G proteins, arrestins, and effectors, activate downstream pathways that often modulate hallmark mechanisms of cancer. Because there are many more GPCRs than effectors, mutations in different receptors could perturb signaling similarly so as to favor a tumor. We hypothesized that somatic mutations in tumor samples may not be enriched within a single gene but rather that cognate mutations with similar effects on GPCR function are distributed across many receptors. To test this possibility, we systematically aggregated somatic cancer mutations across class A GPCRs and found a nonrandom distribution of positions with variant amino acid residues. Individual cancer types were enriched for highly impactful, recurrent mutations at selected cognate positions of known functional motifs. We also discovered that no single receptor drives this pattern, but rather multiple receptors contain amino acid substitutions at a few cognate positions. Phenotypic characterization suggests these mutations induce perturbation of G protein activation and/or β-arrestin recruitment. These data suggest that recurrent impactful oncogenic mutations perturb different GPCRs to subvert signaling and promote tumor growth or survival. The possibility that multiple different GPCRs could moonlight as drivers or enablers of a given cancer through mutations located at cognate positions across GPCR paralogs opens a window into cancer mechanisms and potential approaches to therapeutics.

GPCR activation mechanisms across classes and macro/microscales

Two-thirds of human hormones and one-third of clinical drugs activate ~350 G-protein-coupled receptors (GPCR) belonging to four classes: A, B1, C and F. Whereas a model of activation has been described for class A, very little is known about the activation of the other classes, which differ by being activated by endogenous ligands bound mainly or entirely extracellularly. Here we show that, although they use the same structural scaffold and share several 'helix macroswitches', the GPCR classes differ in their 'residue microswitch' positions and contacts. We present molecular mechanistic maps of activation for each GPCR class and methods for contact analysis applicable for any functional determinants. This provides a superfamily residue-level rationale for conformational selection and allosteric communication by ligands and G proteins, laying the foundation for receptor-function studies and drugs with the desired modality.

Deconstructing the Pharmacological Contribution of Sphingosine-1 Phosphate Receptors to Mouse Models of Multiple Sclerosis Using the Species Selectivity of Ozanimod, a Dual Modulator of Human Sphingosine-1 Phosphate Receptor Subtypes 1 and 5

Ozanimod, a sphingosine-1 phosphate (S1P) receptor modulator that binds with high affinity selectively to S1P receptor subtypes 1 (S1P₁) and 5 (S1P₅), is approved for the treatment of relapsing multiple sclerosis (MS) in multiple countries. Ozanimod profiling revealed a species difference in its potency for S1P₅ in mouse, rat, and canine compared with that for human and monkey. Site-directed mutagenesis identified amino acid alanine at position 120 to be responsible for loss of activity for mouse, rat, and canine S1P₅ and mutation back to threonine as in human/monkey S1P₅ restored activity. Radioligand binding analysis performed with mouse S1P₅ confirmed the potency loss is a consequence of a loss of affinity of ozanimod for mouse S1P₅ and was restored with mutation of alanine 120 to threonine. Study of ozanimod in preclinical mouse models of MS can now determine the S1P receptor(s) responsible for observed efficacies with receptor engagement as measured using pharmacokinetic exposures of free drug. Hence, in the experimental autoimmune encephalomyelitis model, ozanimod exposures sufficient to engage S1P₁, but not S1P₅, resulted in reduced circulating lymphocytes, disease scores, and body weight loss; reduced inflammation, demyelination, and apoptotic cell counts in the spinal cord; and reduced circulating levels of the neuronal degeneration marker, neurofilament light. In the demyelinating cuprizone model, ozanimod prevented axonal degradation and myelin loss during toxin challenge but did not facilitate enhanced remyelination post-intoxication. Since free drug levels in this model only engaged S1P_1, we concluded that S1P₁ activation is neuroprotective but does not appear to affect remyelination. Significance Statement Ozanimod, a selective human S1P_1/5 modulator, displays reduced potency for rodent and dog S1P₅ compared with human, which results from mutation of threonine to alanine at position 120. Ozanimod can thus be used as a selective S1P₁ agonist in mouse models of multiple sclerosis to define efficacies driven by S1P₁ but not S1P₅ Based on readouts for experimental autoimmune encephalomyelitis and cuprizone intoxication, S1P₁ modulation is neuroprotective but S1P₅ activity may be required for remyelination.

Activation pathway of a G protein-coupled receptor uncovers conformational intermediates as targets for allosteric drug design

G protein-coupled receptors (GPCRs) are the most common proteins targeted by approved drugs. A complete mechanistic elucidation of large-scale conformational transitions underlying the activation mechanisms of GPCRs is of critical importance for therapeutic drug development. Here, we apply a combined computational and experimental framework integrating extensive molecular dynamics simulations, Markov state models, site-directed mutagenesis, and conformational biosensors to investigate the conformational landscape of the angiotensin II (AngII) type 1 receptor (AT1 receptor) - a prototypical class A GPCR-activation. Our findings suggest a synergistic transition mechanism for AT1 receptor activation. A key intermediate state is identified in the activation pathway, which possesses a cryptic binding site within the intracellular region of the receptor. Mutation of this cryptic site prevents activation of the downstream G protein signaling and β-arrestin-mediated pathways by the endogenous AngII octapeptide agonist, suggesting an allosteric regulatory mechanism. Together, these findings provide a deeper understanding of AT1 receptor activation at an atomic level and suggest avenues for the design of allosteric AT1 receptor modulators with a broad range of applications in GPCR biology, biophysics, and medicinal chemistry.

Bitter Taste and Olfactory Receptors: Beyond Chemical Sensing in the Tongue and the Nose

Abstract

The Up-and-Coming-Scientist section of the current issue of the Journal of Membrane Biology features the invited essay by Dr. Mercedes Alfonso-Prieto, Assistant Professor at the Forschungszentrum Jülich (FZJ), Germany, and the Heinrich-Heine University Düsseldorf, Vogt Institute for Brain Research. Dr. Alfonso-Prieto completed her doctoral degree in chemistry at the Barcelona Science Park, Spain, in 2009, pursued post-doctoral research in computational molecular sciences at Temple University, USA, and then, as a Marie Curie post-doctoral fellow at the University of Barcelona, worked on computations of enzyme reactions and modeling of photoswitchable ligands targeting neuronal receptors. In 2016, she joined the Institute for Advanced Science and the Institute for Computational Biomedicine at the FZJ, where she pursues research on modeling and simulation of chemical senses. The invited essay by Dr. Alfonso-Prieto discusses state-of-the-art modeling of molecular receptors involved in chemical sensing – the senses of taste and smell. These receptors, and computational methods to study them, are the focus of Dr. Alfonso-Prieto’s research. Recently, Dr. Alfonso-Prieto and colleagues have presented a new methodology to predict ligand binding poses for GPCRs, and extensive computations that deciphered the ligand selectivity determinants of bitter taste receptors. These developments inform our current understanding of how taste occurs at the molecular level.

Graphic Abstract

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.

Structures of active-state orexin receptor 2 rationalize peptide and small-molecule agonist recognition and receptor activation

Narcolepsy type 1 (NT1) is a chronic neurological disorder that impairs the brain’s ability to control sleep-wake cycles. Current therapies are limited to the management of symptoms with modest effectiveness and substantial adverse effects. Agonists of the orexin receptor 2 (OX₂R) have shown promise as novel therapeutics that directly target the pathophysiology of the disease. However, identification of drug-like OX₂R agonists has proven difficult. Here we report cryo-electron microscopy structures of active-state OX₂R bound to an endogenous peptide agonist and a small-molecule agonist. The extended carboxy-terminal segment of the peptide reaches into the core of OX₂R to stabilize an active conformation, while the small-molecule agonist binds deep inside the orthosteric pocket, making similar key interactions. Comparison with antagonist-bound OX₂R suggests a molecular mechanism that rationalizes both receptor activation and inhibition. Our results enable structure-based discovery of therapeutic orexin agonists for the treatment of NT1 and other hypersomnia disorders.

Agonists of the orexin receptor 2 (OX₂R) show promise in the treatment of narcolepsy. Cryo-EM structures of active-state OX₂R bound to an endogenous peptide agonist and a small-molecule agonist suggest a molecular mechanism that rationalizes both receptor activation and inhibition.

Structural Assessment of Agonist Efficacy in the μ-Opioid Receptor: Morphine and Fentanyl Elicit Different Activation Patterns

Over the past two decades, the opioid epidemic in the United States and Canada has evidenced the need for a better understanding of the molecular mechanisms of medications used to fight pain. Morphine and fentanyl are widely used in opiate-mediated analgesia for the treatment of chronic pain. These compounds target the μ-opioid receptor (MOR), a class A G protein-coupled receptor (GPCR). In light of described higher efficacy of fentanyl with respect to morphine, we have performed independent μs-length unbiased molecular dynamics (MD) simulations of MOR complexes with each of these ligands, including the MOR antagonist naltrexone as a negative control. Consequently, MD simulations totaling 58 μs have been conducted to elucidate at the atomic level ligand-specific receptor activity and signal transmission in the MOR. In particular, we have identified stable binding poses of morphine and fentanyl, which interact differently with the MOR. Different ligand-receptor interaction landscapes directly induce sidechain conformational changes of orthosteric pocket residues: Asp149^3.32, Tyr150^3.33, Gln126^2.60, and Lys235^5.39. The induced conformations determine Asp149^3.32-Tyr328^7.43 sidechain-sidechain interactions and Trp295^6.48-Ala242^5.46 sidechain-backbone H-bond formations, as well as Met153^3.36 conformational changes. In addition to differences in ligand binding, different intracellular receptor conformational changes are observed as morphine preferentially activates transmembrane (TM) helices: TM3 and TM5, while fentanyl preferentially activates TM6 and TM7. As conformational changes in TM6 and TM7 are widely described as being the most crucial aspect in GPCR activation, this may contribute to the greater efficacy of fentanyl over morphine. These computationally observed functional differences between fentanyl and morphine may provide new avenues for the design of safer but not weaker opioid drugs because it is desirable to increase the safety of medicines without sacrificing their efficacy.

Probing biased activation of mu-opioid receptor by the biased agonist PZM21 using all atom molecular dynamics simulation

Morphine is a commonly used opioid drug to treat acute pain by binding to the mu-opioid receptor (MOR), but its effective analgesic efficacy via triggering of the heterotrimeric G_i protein pathway is accompanied by a series of adverse side effects via triggering of the β-arrestin pathway. Recently, PZM21, a recently developed MOR biased agonist, shows preferentially activating the G protein pathway over β-arrestin pathway. However, there is no high-resolution receptor structure in complex with PZM21 and its action mechanism remains elusive. In this study, PZM21 and Morphine were docked to the active human MOR-1 homology structure and then subjected to the molecular dynamics (MD) simulations in two different situations (i.e., one situation includes the crystal waters but another does not). Detailed comparisons between the two systems were made to characterize the differences in protein-ligand interactions, protein secondary and tertiary structures and dynamics networks. PZM21 could strongly interact with Y328^7.43 of TM7, besides the residues (Asp149^3.32 and Tyr150^3.33) of TM3. The two systems' network paths to the intracellular end of TM6 were roughly similar but the paths to the end of TM7 were different. The PZM21-bound MOR's intracellular ends of TM5-7 bent outward more along with the distance changes of the three key molecular switches (ionic lock, transmission and Tyr toggle) and the distance increase of some conserved inter-helical residue pairs. The larger intracellular opening of the receptor could potentially facilitate G protein binding.

Probing Allosteric Mechanism with Long-Range Rigidity Transmission Across Protein Networks

Allosteric transmission refers to regulation of protein function at a distance. "Allostery" involves regulation and/or signal transduction induced by a perturbation event. Allostery, which has been coined the "second secret of life," is a fundamental property of most dynamics proteins. Most of critical questions surrounding allostery are largely unresolved. One of the key puzzles is to describe the physical mechanism of distant coupled conformational change. Another hot research area surrounding allostery is detection of allosteric pathways or regions (residues) in the protein that are the most critical for transmission of allosteric information. Using techniques inspired by mathematical rigidity theory and mechanical linkages, we have previously proposed a mechanistic model and description of allosteric transmission and an accompanying computational method, the Rigidity Transmission Allostery (RTA) algorithm. The RTA algorithm and method are designed to predict if mechanical perturbation of rigidity, for example, due to ligand binding, at one site of the protein can transmit and propagate across a protein structure and in turn cause a change in available conformational degrees of freedom and a change in conformation at a second distant site, equivalently resulting in allosteric transmission. The RTA algorithm is computationally very fast and can rapidly scan many unknown sites for allosteric transmission, identifying potential novel allosteric sites and quantify their allosteric effect. In this chapter we will discuss the rigidity-based mechanistic model of allosteric communication. As a case illustrative study, we will demonstrate RTA analysis on a G protein coupled receptor (GPCR) human adenosine A_2A receptor. Our method gives important implications and a novel prospective for general mechanistic description of allosteric communication.

Ligand based conformational space studies of the μ-opioid receptor

BACKGROUND

G protein-coupled receptors (GPCRs) comprise a family of membrane proteins that can be activated by a variety of external factors. The μ-opioid receptor (MOR), a class A GPCR, is the main target of morphine. Recently, enhanced sampling molecular dynamics simulations of a constitutively active mutant of MOR in its apo form allowed us to capture the novel intermediate states of activation, as well as the active state. This prompted us to apply the same techniques to wild type MOR in complex with ligands, in order to explore their contributions to the receptor conformational changes in the activation process.

METHODS

MOR was modeled in complex with agonists (morphine, BU72), a partial agonist (naloxone benzoylhydrazone) and an antagonist (naloxone). Replica exchange with solute tempering (REST2) molecular dynamics simulations were carried out for all systems. Trajectory frames were clustered, and the activation state of each cluster was assessed by two different methods.

RESULTS

Cluster sizes and activation indices show that while agonists stabilized structures in a higher activation state, the antagonist behaved oppositely. Morphine tends to drive the receptor towards increasing R165-T279 distances, while naloxone tends to increase the NPxxYA motif conformational change.

CONCLUSIONS

Despite not observing a full transition between inactive and active states, an important conformational change of transmembrane helix 5 was observed and associated with a ligand-driven step of the process.

GENERAL SIGNIFICANCE

The activation process of GPCRs is widely studied but still not fully understood. Here we carried out a step forward in the direction of gaining more details of this process.

G protein-coupled receptor-G protein interactions: a single-molecule perspective

G protein-coupled receptors (GPCRs) regulate many cellular and physiological processes, responding to a diverse range of extracellular stimuli including hormones, neurotransmitters, odorants, and light. Decades of biochemical and pharmacological studies have provided fundamental insights into the mechanisms of GPCR signaling. Thanks to recent advances in structural biology, we now possess an atomistic understanding of receptor activation and G protein coupling. However, how GPCRs and G proteins interact in living cells to confer signaling efficiency and specificity remains insufficiently understood. The development of advanced optical methods, including single-molecule microscopy, has provided the means to study receptors and G proteins in living cells with unprecedented spatio-temporal resolution. The results of these studies reveal an unexpected level of complexity, whereby GPCRs undergo transient interactions among themselves as well as with G proteins and structural elements of the plasma membrane to form short-lived signaling nanodomains that likely confer both rapidity and specificity to GPCR signaling. These findings may provide new strategies to pharmaceutically modulate GPCR function, which might eventually pave the way to innovative drugs for common diseases such as diabetes or heart failure.

Activation Microswitches in Adenosine Receptor A2A Function as Rheostats in the Cell Membrane

Although multiple components of the cell membrane modulate the stability and activation of G protein-coupled receptors (GPCRs), insights into the dynamics of GPCR structures come from biophysical studies conducted in detergents. This is because of the challenges of studying activation in a multicomponent lipid bilayer. To understand the role of cellular membrane lipids and cations in GPCR activation, we performed multiscale molecular dynamics simulations (56 μs) on three different conformational states of adenosine receptor A2AR, in both the cell membrane-like lipid bilayer and in detergent micelles. Molecular dynamics (MD) simulations show that the phosphatidylinositol bisphosphate (PIP2) interacts with the basic residues in the intracellular regions of A2AR, thereby reducing the flexibility of the receptor in the inactive state and limiting the transition to the active-intermediate state. In the G protein-coupled fully active state, PIP2 stabilizes the GPCR:G protein complex. Such stiffening effects are absent in non-ionic detergent micelles, and therefore, more transitions have been observed in detergents. The inter-residue distances that change significantly upon GPCR activation are known as activation microswitches. The activation microswitches show different levels of activation in the cell membrane, in the pure POPC bilayer, and in detergents. Thus, the temporal heat map of different activation microswitches calculated from the MD simulations suggests a rheostat model of GPCR activation microswitches rather than the binary switch model. These simulation results connect the chemistry of cell membrane lipids to receptor activity, which is useful for the design of detergents mimicking the cell membrane.

HomolWat: a web server tool to incorporate 'homologous' water molecules into GPCR structures

Internal water molecules play an essential role in the structure and function of membrane proteins including G protein-coupled receptors (GPCRs). However, technical limitations severely influence the number and certainty of observed water molecules in 3D structures. This may compromise the accuracy of further structural studies such as docking calculations or molecular dynamics simulations. Here we present HomolWat, a web application for incorporating water molecules into GPCR structures by using template-based modelling of homologous water molecules obtained from high-resolution structures. While there are various tools available to predict the positions of internal waters using energy-based methods, the approach of borrowing lacking water molecules from homologous GPCR structures makes HomolWat unique. The tool can incorporate water molecules into a protein structure in about a minute with around 85% of water recovery. The web server is freely available at http://lmc.uab.es/homolwat.

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.

Coevolutionary data-based interaction networks approach highlighting key residues across protein families: The case of the G-protein coupled receptors

We present an approach that, by integrating structural data with Direct Coupling Analysis, is able to pinpoint most of the interaction hotspots (i.e. key residues for the biological activity) across very sparse protein families in a single run. An application to the Class A G-protein coupled receptors (GPCRs), both in their active and inactive states, demonstrates the predictive power of our approach. The latter can be easily extended to any other kind of protein family, where it is expected to highlight most key sites involved in their functional activity.

Truncated (N)-Methanocarba Nucleosides as Partial Agonists at Mouse and Human A3 Adenosine Receptors: Affinity Enhancement by N6-(2-Phenylethyl) Substitution

Dopamine-derived N6-substituents, compared to N6-(2-phenylethyl), in truncated (N)-methanocarba (bicyclo[3.1.0]hexyl) adenosines favored high A3 adenosine receptor (AR) affinity/selectivity, e.g., C2-phenylethynyl analogue 15 (MRS7591, Ki = 10.9/17.8 nM, at human/mouse A3AR). 15 was a partial agonist in vitro (hA3AR, cAMP inhibition, 31% Emax; mA3AR, [35S]GTP-γ-S binding, 16% Emax) and in vivo and also antagonized hA3AR in vitro. Distal H-bonding substitutions of the N6-(2-phenylethyl) moiety particularly enhanced mA3AR affinity by polar interactions with the extracellular loops, predicted using docking and molecular dynamics simulation with newly constructed mA3AR and hA3AR homology models. These hybrid models were based on an inactive antagonist-bound hA1AR structure for the upper part of TM2 and an agonist-bound hA2AAR structure for the remaining TM portions. These species-independent A3AR-selective nucleosides are low efficacy partial agonists and novel, nuanced modulators of the A3AR, a drug target of growing interest.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Advances in therapeutic peptides targeting G protein-coupled receptors

Dysregulation of peptide-activated pathways causes a range of diseases, fostering the discovery and clinical development of peptide drugs. Many endogenous peptides activate G protein-coupled receptors (GPCRs) - nearly 50 GPCR peptide drugs have been approved to date, most of them for metabolic disease or oncology, and more than 10 potentially first-in-class peptide therapeutics are in the pipeline. The majority of existing peptide therapeutics are agonists, which reflects the currently dominant strategy of modifying the endogenous peptide sequence of ligands for peptide-binding GPCRs. Increasingly, novel strategies are being employed to develop both agonists and antagonists, to both introduce chemical novelty and improve drug-like properties. Pharmacodynamic improvements are evolving to allow biasing ligands to activate specific downstream signalling pathways, in order to optimize efficacy and reduce side effects. In pharmacokinetics, modifications that increase plasma half-life have been revolutionary. Here, we discuss the current status of the peptide drugs targeting GPCRs, with a focus on evolving strategies to improve pharmacokinetic and pharmacodynamic properties.

Exploring the Activation Mechanism of the mGlu5 Transmembrane Domain

As a class C GPCR and regulator of synaptic activity, mGlu5 is an attractive drug target, potentially offering treatment for several neurologic and psychiatric disorders. As little is known about the activation mechanism of mGlu5 at a structural level, potential of mean force calculations linked to molecular dynamics simulations were performed on the mGlu5 transmembrane domain crystal structure to explore various internal mechanisms responsible for its activation. Our results suggest that the hydrophilic interactions between intracellular loop 1 and the intracellular side of TM6 have to be disrupted to reach a theoretically active-like conformation. In addition, interactions between residues that are key for mGlu5 activation (Tyr659^3.44 and Ile751^5.51) and mGlu5 inactivation (Tyr659^3.44 and Ser809^7.39) have been identified. Inasmuch as mGlu5 receptor signaling is poorly understood, potentially showing a complex network of micro-switches and subtle structure-activity relationships, the present study represents a step forward in the understanding of mGlu5 transmembrane domain activation.

Relevance of Hydrogen Bonds for the Histamine H2 Receptor-Ligand Interactions: A Lesson from Deuteration

We used a combination of density functional theory (DFT) calculations and the implicit quantization of the acidic N–H and O–H bonds to assess the effect of deuteration on the binding of agonists (2-methylhistamine and 4-methylhistamine) and antagonists (cimetidine and famotidine) to the histamine H2 receptor. The results show that deuteration significantly increases the affinity for 4-methylhistamine and reduces it for 2-methylhistamine, while leaving it unchanged for both antagonists, which is found in excellent agreement with experiments. The revealed trends are interpreted in the light of the altered strength of the hydrogen bonding upon deuteration, known as the Ubbelohde effect, which affects ligand interactions with both active sites residues and solvent molecules preceding the binding, thus providing strong evidence for the relevance of hydrogen bonding for this process. In addition, computations further underline an important role of the Tyr250 residue for the binding. The obtained insight is relevant for the therapy in the context of (per)deuterated drugs that are expected to enter therapeutic practice in the near future, while this approach may contribute towards understanding receptor activation and its discrimination between agonists and antagonists.

Exploring the Activation Mechanism of a Metabotropic Glutamate Receptor Homodimer via Molecular Dynamics Simulation

Metabotropic glutamate receptors of class C GPCRs exist as constitutive dimers, which play important roles in activating excitatory synapses of the central nervous system. However, the activation mechanism induced by agonists has not been clarified in experiments. To address the problem, we used microsecond all-atom molecular dynamics (MD) simulation couple with protein structure network (PSN) to explore the glutamate-induced activation for the mGluR1 homodimer. The results indicate that glutamate binding stabilizes not only the closure of Venus flytrap domains but also the polar interaction of LB2-LB2, in turn keeping the extracelluar domain in the active state. The activation of the extracelluar domain drives transmembrane domains (TMDs) of the two protomers closer and induces asymmetric activation for the TMD domains of the two protomers. One protomer with lower binding affinity to the agonist is activated, while the other protomer with higher binding energy is still in the inactive state. The PSN analysis identifies the allosteric regulation pathway from the ligand-binding pocket in the extracellular domain to the G-protein binding site in the intracellular TMD region and further reveals that the asymmetric activation is attributed to a combination of trans-pathway and cis-pathway regulations from two glumatates, rather than a single activation pathway. These observations could provide valuable molecular information for understanding of the structure and the implications in drug efficacy for the class C GPCR dimers.

Preferential Coupling of Dopamine D2S and D2L Receptor Isoforms with Gi1 and Gi2 Proteins-In Silico Study

The dopamine D₂ receptor belongs to rhodopsin-like G protein-coupled receptors (GPCRs) and it is an important molecular target for the treatment of many disorders, including schizophrenia and Parkinson's disease. Here, computational methods were used to construct the full models of the dopamine D₂ receptor short (D_2S) and long (D_2L) isoforms (differing with 29 amino acids insertion in the third intracellular loop, ICL3) and to study their coupling with G_i1 and G_i2 proteins. It was found that the D_2L isoform preferentially couples with the G_i2 protein and D_2S isoform with the G_i1 protein, which is in accordance with experimental data. Our findings give mechanistic insight into the interplay between isoforms of dopamine D₂ receptors and G_i proteins subtypes, which is important to understand signaling by these receptors and their mediation by pharmaceuticals, in particular psychotic and antipsychotic agents.