Common activation mechanism of class A GPCRs

Insight into structural properties of viral G protein-coupled receptors and their role in the viral infection: IUPHAR Review 41

G protein-coupled receptors (GPCRs) are pivotal in cellular signalling and drug targeting. Herpesviruses encode GPCRs (vGPCRs) to manipulate cellular signalling, thereby regulating various aspects of the virus life cycle, such as viral spreading and immune evasion. vGPCRs mimic host chemokine receptors, often with broader signalling and high constitutive activity. This review focuses on the recent advancements in structural knowledge about vGPCRs, with an emphasis on molecular mechanisms of action and ligand binding. The structures of US27 and US28 from human cytomegalovirus (HCMV) are compared to their closest human homologue, CX3CR1. Contrasting US27 and US28, the homotrimeric UL78 structure (HCMV) reveals more distance to chemokine receptors. Open reading frame 74 (ORF74; Kaposi's sarcoma-associated herpesvirus) is compared to CXCRs, whereas BILF1 (Epstein-Barr virus) is discussed as a putative lipid receptor. Furthermore, the roles of vGPCRs in latency and lytic replication, reactivation, dissemination and immune evasion are reviewed, together with their potential as drug targets for virus infections and virus-related diseases.

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.

The effect of cancer-associated mutations on ligand binding and receptor function - A case for the 5-HT2C receptor

The serotonin 5-HT2C receptor is a G protein-coupled receptor (GPCR) mainly expressed in the central nervous system. Besides regulating mood, appetite, and reproductive behavior, it has been identified as a potential target for cancer treatment. In this study, we aimed to investigate the effects of cancer patient-derived 5-HT2C receptor mutations on ligand binding and receptor functionality. By filtering the sequencing data from the Genomic Data Commons data portal (GDC), we selected 12 mutations from multiple cancer types. We found that the affinity of the endogenous agonist serotonin (5-HT) and inverse agonist mesulergine were both drastically decreased by mutations L209HECL2 and F328S6.52, which are located in the orthosteric binding pocket. In the calcium-flux assay, the potency of 5-HT was decreased at F328S6.52, while a trend of increased efficacy was observed. In contrast, 5-HT displayed higher affinity at E306K6.30 and E306A6.30, while a trend of decreased efficacy was observed. These two mutations may disrupt the conserved ionic interaction between E6.30 and R3.50, and thus increase the constitutive activity of the receptor. The inhibitory potency of mesulergine was increased at E306A6.30 but not E306K6.30. Lastly, P365H7.50 decreased the expression level of the receptor by more than ten-fold, which prevented further functional analyses. This study shows that cancer-associated mutations of 5-HT2C receptor have diverse effects on ligand binding and function. Such mutations may affect serotonin-mediated signaling in tumor cells as well as treatment strategies targeting this receptor.

Targeting adhesion G protein-coupled receptors. Current status and future perspectives

G protein-coupled receptors (GPCRs) orchestrate many physiological functions and are a crucial target in drug discovery. Adhesion GPCRs (aGPCRs), the second largest family within this superfamily, are promising yet underexplored targets for treating various diseases, including obesity, psychiatric disorders, and cancer. However, the receptors' unique and complex structure and miscellaneous interactions complicate comprehensive pharmacological studies. Despite recent progress in determining structures and elucidation of the activation mechanism, the function of many receptors remains to be determined. This review consolidates current knowledge on aGPCR ligands, focusing on small molecule orthosteric ligands and allosteric modulators identified for the ADGRGs subfamily (subfamily VIII), (GPR56/ADGRG1, GPR64/ADGRG2, GPR97/ADGRG3, GPR114/ADGRG5, GPR126/ADGRG6, and GPR128/ADGRG7). We discuss challenges in hit identification, target validation, and drug discovery, highlighting molecular compositions and recent structural breakthroughs. ADGRG ligands can offer new insights into aGPCR modulation and have significant potential for novel therapeutic interventions targeting various diseases.

Divergent roles of DRY and NPxxY motifs in selective activation of downstream signalling by the apelin receptor

G protein-coupled receptors (GPCRs) serve as critical communication hubs, translating a wide range of extracellular signals into intracellular responses that govern numerous physiological processes. In class-A GPCRs, conserved motifs mediate conformational changes of the active states of the receptor, and signal transduction is achieved by selectively binding to Gα proteins and/or adapter protein, arrestin. Apelin receptor (APJR) is a class-A GPCR that regulates a wide range of intracellular signalling cascades in response to apelin and elabela peptide ligands. Understanding how conserved motifs within APJR mediate activation and signal specificity remains unexplored. This study focuses on the functional roles of the DRY and NPxxY motifs within APJR by analyzing their impact on downstream signaling pathways across the receptor's conformational ensembles. Our findings provide compelling evidence that mutations within the conserved DRY and NPxxY motifs of APJR significantly alter its conformational preferences where modification of DRY motif leads to abrogation of G-protein coupling and mutation of NPxxY motif causing abolition of β-arrestin-2 recruitment. These observations shed light on the importance of these motifs in APJR activation and its potential for functional selectivity, highlighting the role of DRY/NPxxY as conformational switches of APJR signalling.

Exploring TAS2R46 biomechanics through molecular dynamics and network analysis

Understanding the intricate interplay between structural features and signal-processing events is crucial for unravelling the mechanisms of biomolecular systems. G protein-coupled receptors (GPCRs), a pervasive protein family in humans, serve a wide spectrum of vital functions. TAS2Rs, a subfamily of GPCRs, play a primary role in recognizing bitter molecules and triggering events leading to the perception of bitterness, a crucial defence mechanism against spoiled or poisonous food. Beyond taste, TAS2Rs function is associated with many diseases as they are expressed in several extra-oral tissues. Given that the precise functioning mechanisms of TAS2R remain poorly understood, this study employed molecular dynamics simulations combined with network-based analysis to investigate local conformational changes and global structural correlations in different states of the receptor. The focus was on the human TAS2R46 bitter taste receptor, recently resolved experimentally, both in the presence and absence of strychnine, a known bitter agonist. The results showed that the ligand-bound state of the receptor exhibited more correlated dynamics compared to the apo state, and the presence of the agonist mediated the allosteric network between two helices (TM3 and TM6) which mainly convey the signal transferring from the extracellular to the intracellular region. By elucidating the hallmarks of the conformational changes and allosteric network of TAS2R46 under varying conditions, this study has enabled the identification of the unique structural and dynamics features of this receptor, thereby establishing a foundation for a more profound characterisation of this intriguing class of receptors.

Exploring key features of selectivity in somatostatin receptors through molecular dynamics simulations

Somatostatin receptors (SSTRs) are widely distributed throughout the human body and play crucial roles in various physiological processes. They are recognized as key targets for both radiotherapy and radiodiagnosis due to their overexpression in several cancer types. However, the discovery and design of selective drugs for each of the five isoforms have been significantly hindered by the lack of complete structural information. In this study, we conducted a systematic computational analysis of all five SSTRs in complex with the endogenous ligand somatostatin to elucidate their structural and dynamic features. We thoroughly characterized each isoform using available experimental structures for SSTR2 and SSTR4, as well as AlphaFold2 models for SSTR1, SSTR3, and SSTR5. By performing multi-copy μs-long molecular dynamics simulations, we examined the differences and similarities in dynamical behavior and somatostatin binding among all SSTRs. Our analysis focused on understanding the opening and closing movements of the extracellular loop 2, which are crucial for ligand binding and recognition. Interestingly, we observed a unique conformation of somatostatin within the binding pocket of SSTR5 in which the loop can partially close, as compared to the other isoforms. Fingerprint analyses provided distinct interaction patterns of somatostatin with all receptors, thus enabling precise guidelines for the discovery and development of more selective somatostatin-based pharmaceuticals tailored for precision medicine therapies.

Evolutionary diversity of CXCL16-CXCR6: Convergent substitutions and recurrent gene loss in sauropsids

The CXCL16-CXCR6 axis is crucial for regulating the persistence of CD8 tissue-resident memory T cells (TRM). CXCR6 deficiency lowers TRM cell numbers in the lungs and depletes ILC3s in the lamina propria, impairing mucosal defence. This axis is linked to diseases like HIV/SIV, cancer, and COVID-19. Together, these highlight that the CXCL16-CXCR6 axis is pivotal in host immunity. Previous studies of the CXCL16-CXCR6 axis found genetic variation among species but were limited to primates and rodents. To understand the evolution and diversity of CXCL16-CXCR6 across vertebrates, we compared approximately 400 1-to-1 CXCR6 orthologs spanning diverse vertebrates. The unique DRF motif of CXCR6 facilitates leukocyte adhesion by interacting with cell surface-expressed CXCL16 and plays a key role in G-protein selectivity during receptor signalling; however, our findings show that this motif is not universal. The DRF motif is restricted to mammals, turtles, and frogs, while the DRY motif, typical in other CKRs, is found in snakes and lizards. Most birds exhibit the DRL motif. These substitutions at the DRF motif affect the receptor-Gi/o protein interaction. We establish recurrent CXCR6 gene loss in 10 out of 36 bird orders, including Galliformes and Passeriformes, Crocodilia, and Elapidae, attributed to segmental deletions and/or frame-disrupting changes. Notably, single-cell RNA sequencing of the lung shows a drop in TRM cells in species with CXCR6 loss, suggesting a possible link. The concurrent loss of ITGAE, CXCL16, and CXCR6 in chickens may have altered CD8 TRM cell abundance, with implications for immunity against viral diseases and vaccines inducing CD8 TRM cells.

Biased activation of the vasopressin V2 receptor probed by molecular dynamics simulations, NMR and pharmacological studies

G protein-coupled receptors (GPCRs) control critical cell signaling. Their response to extracellular stimuli involves conformational changes to convey signals to intracellular effectors, among which the most important are G proteins and β-arrestins (βArrs). Biased activation of one pathway is a field of intense research in GPCR pharmacology. Combining NMR, site-directed mutagenesis, molecular pharmacology, and molecular dynamics (MD) simulations, we studied the conformational diversity of the vasopressin V2 receptor (V2R) bound to different types of ligands: the antagonist Tolvaptan, the endogenous unbiased agonist arginine-vasopressin, and MCF14, a partial Gs protein-biased agonist. A double-labeling NMR scheme was developed to study the receptor conformational changes and ligand binding: V2R was subjected to lysine 13CH3 methylation for complementary NMR studies, whereas the agonists were tagged with a paramagnetic probe. Paramagnetic relaxation enhancements and site-directed mutagenesis validated the ligand binding modes in the MD simulations. We found that the bias for the Gs protein over the βArr pathway involves interactions between the conserved NPxxY motif in the transmembrane helix 7 (TM7) and TM3, compacting helix 8 (H8) toward TM1 and likely inhibiting βArr signaling. A similar mechanism was elicited for the pathogenic mutation I130N, which constitutively activates the Gs proteins without concomitant βArr recruitment. The findings suggest common patterns of biased signaling in class A GPCRs, as well as a rationale for the design of G protein-biased V2R agonists.

Early Events in β2AR Dimer Dynamics Mediated by Activation-Related Microswitches

G-Protein-Coupled Receptors (GPCRs) make up around 3-4% of the human genome and are the targets of one-third of FDA-approved drugs. GPCRs typically exist as monomers but also aggregate to form higher-order oligomers, including dimers. β2AR, a pharmacologically relevant GPCR, is known to be targeted for the treatment of asthma and cardiovascular diseases. The activation of β2AR at the dimer level remains under-explored. In the current study, molecular dynamics (MD) simulations have been performed to understand activation-related structural changes in β2AR at the dimer level. The transition from inactive to active and vice versa has been studied by starting the simulations in the apo, agonist-bound, and inverse agonist-bound β2AR dimers for PDB ID: 2RH1 and PDB ID: 3P0G, respectively. A cumulative total of around 21-μs simulations were performed. Residue-based distances, RMSD, and PCA calculations suggested that either of the one monomer attained activation-related features for the apo and agonist-bound β2AR dimers. The TM5 and TM6 helices within the two monomers were observed to be in significant variation in all the simulations. TM5 bulge and proximity of TM2 and TM7 helices may be contributing to one of the early events in activation. The dimeric interface between TM1 and helix 8 were observed to be well maintained in the apo and agonist-bound simulations. The presence of inverse agonists favored inactive features in both the monomers. These key features of activation known for monomers were observed to have an impact on β2AR dimers, thereby providing an insight into the oligomerization mechanism of GPCRs.

Structural basis for ligand recognition of the human hydroxycarboxylic acid receptor HCAR3

Hydroxycarboxylic acid receptor 3 (HCAR3), a class A G-protein-coupled receptor, is an important cellular energy metabolism sensor with a key role in the regulation of lipolysis in humans. HCAR3 is deeply involved in many physiological processes and serves as a valuable target for the treatment of metabolic diseases, tumors, and immune diseases. Here, we report four cryoelectron microscopy (cryo-EM) structures of human HCAR3-Gi1 complexes with or without agonists: the endogenous ligand 3-hydroxyoctanoic acid, the drug niacin, the highly subtype-specific agonist compound 5c (4-(n-propyl)amino-3-nitrobenzoic acid), and the apo form. Together with mutagenesis and functional analyses, we revealed the recognition mechanisms of HCAR3 for different agonists. In addition, the key residues that determine the ligand selectivity between HCAR2 and HCAR3 were also illuminated. Overall, these findings provide a structural basis for the ligand recognition, activation, and selectivity and G-protein coupling mechanisms of HCAR3, which contribute to the design of HCAR3-targeting drugs with high efficacy and selectivity.

Molecular basis for chemokine recognition and activation of XCR1

The X-C motif chemokine receptor XCR1, which selectively binds to the chemokine XCL1, is highly expressed in conventional dendritic cells subtype 1 (cDC1s) and crucial for their activation. Modulating XCR1 signaling in cDC1s could offer novel opportunities in cancer immunotherapy and vaccine development by enhancing the antigen presentation function of cDC1s. To investigate the molecular mechanism of XCL-induced XCR1 signaling, we determined a high-resolution structure of the human XCR1 and Gi complex with an engineered form of XCL1, XCL1 CC3, by cryoelectron microscopy. Through mutagenesis and structural analysis, we elucidated the molecular details for the binding of the N-terminal segment of XCL1 CC3, which is vital for activating XCR1. The unique arrangement within the XCL1 CC3 binding site confers specificity for XCL1 in XCR1. We propose an activation mechanism for XCR1 involving structural alterations of key residues at the bottom of the XCL1 binding pocket. These detailed insights into XCL1 CC3-XCR1 interaction and XCR1 activation pave the way for developing novel XCR1-targeted therapeutics.

A non-canonical mechanism of GPCR activation

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

Unraveling the Structural Basis of Biased Agonism in the β2-Adrenergic Receptor Through Molecular Dynamics Simulations

Biased agonism in G protein-coupled receptors is a phenomenon resulting in the selective activation of distinct intracellular signaling pathways by different agonists, which may exhibit bias toward either Gs, Gi, or arrestin-mediated pathways. This study investigates the structural basis of ligand-induced biased agonism within the context of the β2-adrenergic receptor (β2-AR). Atomistic molecular dynamics simulations were conducted for β2-AR complexes with two stereoisomers of methoxynaphtyl fenoterol (MNFen), that is, compounds eliciting qualitatively different cellular responses. The simulations reveal distinct interaction patterns within the binding cavity, dependent on the stereoisomer. These changes propagate to the intracellular parts of the receptor, triggering various structural responses: the dynamic structure of the intracellular regions of the (R,R)-MNFen complex more closely resembles the "Gs-compatible" and "β-arrestin-compatible" conformation of β2-AR, while both stereoisomers maintain structural responses equidistant from the inactive conformation. These findings are confirmed by independent coarse-grained simulations. In the context of deciphered molecular mechanisms, Trp313 plays a pivotal role, altering its orientation upon interactions with (R,R)-MNFen, along with the Lys305-Asp192 ionic bridge. This effect, accompanied by ligand interactions with residues on TM2, increases the strength of interactions within the extracellular region and the binding cavity, resulting in a slightly more open conformation and a minor (by ca. 0.2 nm) increase in the distance between the TM5-TM7, TM1-TM6, TM6-TM7, and TM1-TM5 pairs. On the other hand, an even slighter decrease in the distance between the TM1-TM4 and TM2-TM4 pairs is observed.

Investigating the Effect of GLU283 Protonation State on the Conformational Heterogeneity of CCR5 by Molecular Dynamics Simulations

CCR5 is a class A GPCR and serves as one of the coreceptors facilitating HIV-1 entry into host cells. This receptor has vital roles in the immune system and is involved in the pathogenesis of different diseases. Various studies were conducted to understand its activation mechanism, including structural studies in which inactive and active states of the receptor were determined in complex with various binding partners. These determined structures provided opportunities to perform molecular dynamics (MD) simulations and to analyze conformational changes observed in the protein structures. The atomic-level dynamic studies allow us to explore the effects of ionizable residues on the receptor. Here, our aim was to investigate the conformational changes in CCR5 when it forms a complex with either the inhibitor maraviroc (MRV), an approved anti-HIV drug, or HIV-1 envelope protein GP120, and compare these changes to the receptor's apo form. In our simulations, we considered both ionized and protonated states of ionizable binding site residue GLU2837.39 in CCR5 as the protonation state of this residue was considered ambiguously in previous studies. Our molecular simulations results suggested that in fact, the change in the protonation state of GLU2837.39 caused interaction profiles to be different between CCR5 and its binding partners, GP120 or MRV. We observed that when the protonated state of GLU2837.39 was considered in complex with the envelope protein GP120, there were substantial structural changes in CCR5, indicating that it adopts a more active-like conformation. On the other hand, CCR5 in complex with MRV always adopted an inactive conformation regardless of the protonation state. Hence, the CCR5 coreceptor displays conformational heterogeneity not only depending on its binding partner but also influenced by the protonation state of the binding site binding site residue GLU2837.39. This outcome is also in accordance with some studies showing that GP120 binding could activate signaling pathways. This outcome could also have significant implications for discovering novel CCR5 inhibitors as anti-HIV drugs using in silico methods such as molecular docking, as it may be necessary to consider the protonated state of GLU2837.39.

Mechanism of negative μ-opioid receptor modulation by sodium ions

Negative allosteric modulation of G-protein coupled receptors (GPCRs) by Na+ ions was first described in the 1970s for opioid receptors (ORs) and has subsequently been detected for most class A GPCRs. In high-resolution structures of inactive-state class A GPCRs, a Na+ ion binds to a conserved pocket near residue D2.50, whereas active-state structures of GPCRs are incompatible with Na+ binding. Correspondingly, Na+ diminishes agonist affinity, stabilizes the receptors in the inactive state, and reduces basal signaling. We applied a mutual-information based analysis to μs-timescale biomolecular simulations of the μ-opioid receptor (μ-OR). Our results reveal that Na+ binding is coupled to a water wire linking the Na+ binding site with the agonist binding pocket and to rearrangements in polar networks propagating conformational changes to the agonist and G-protein binding sites. These findings provide a new mechanistic link between the presence of the ion, altered agonist affinity, receptor deactivation, and lowered basal signaling levels.

Conserved pattern-based classification of human odorant receptor multigene family

Conserved protein-coding sequences are critical for maintaining protein function across species. Odorant receptors (ORs), a large poorly understood multigene family responsible for odor detection, lack comprehensive classification methods that reflect their functional diversity. In this study, we propose a new approach called conserved motif-based classification (CMC) for classifying ORs based on amino acid sequence similarities within conserved motifs. Specifically, we focused on three well-conserved motifs: MAYDRYVAIC in TM3, KAFSTCASH in TM6, and PMLNPFIY in TM7. Using an unsupervised clustering technique, we classified human ORs (hORs) into two main clusters with six sub-clusters. CMC partly reflects previously identified subfamilies, revealing altered residue positions among the sub-clusters. These altered positions interacted with specific residues within or adjacent to the transmembrane domain, suggesting functional implications. Furthermore, we found that the CMC correlated with both ligand responses and ectopic expression patterns, highlighting its relevance to OR function. This conserved motif-based classification will help in understanding the functions and features that are not understood by classification based solely on entire amino acid sequence similarity.

BaNDyT: Bayesian Network modeling of molecular Dynamics Trajectories

Bayesian network modeling (BN modeling, or BNM) is an interpretable machine learning method for constructing probabilistic graphical models from the data. In recent years, it has been extensively applied to diverse types of biomedical datasets. Concurrently, our ability to perform long-timescale molecular dynamics (MD) simulations on proteins and other materials has increased exponentially. However, the analysis of MD simulation trajectories has not been data-driven but rather dependent on the user's prior knowledge of the systems, thus limiting the scope and utility of the MD simulations. Recently, we pioneered using BNM for analyzing the MD trajectories of protein complexes. The resulting BN models yield novel fully data-driven insights into the functional importance of the amino acid residues that modulate proteins' function. In this report, we describe the BaNDyT software package that implements the BNM specifically attuned to the MD simulation trajectories data. We believe that BaNDyT is the first software package to include specialized and advanced features for analyzing MD simulation trajectories using a probabilistic graphical network model. We describe here the software's uses, the methods associated with it, and a comprehensive Python interface to the underlying generalist BNM code. This provides a powerful and versatile mechanism for users to control the workflow. As an application example, we have utilized this methodology and associated software to study how membrane proteins, specifically the G protein-coupled receptors, selectively couple to G proteins. The software can be used for analyzing MD trajectories of any protein as well as polymeric materials.

Structural insights into the agonist selectivity of the adenosine A3 receptor

Adenosine receptors play pivotal roles in physiological processes. Adenosine A3 receptor (A3R), the most recently identified adenosine receptor, is expressed in various tissues, exhibiting important roles in neuron, heart, and immune cells, and is often overexpressed in tumors, highlighting the therapeutic potential of A3R-selective agents. Recently, we identified RNA-derived N6-methyladenosine (m6A) as an endogenous agonist for A3R, suggesting the relationship between RNA-derived modified adenosine and A3R. Despite extensive studies on the other adenosine receptors, the selectivity mechanism of A3R, especially for A3R-selective agonists such as m6A and namodenoson, remained elusive. Here, we identify tRNA-derived N6-isopentenyl adenosine (i6A) as an A3R-selective ligand via screening of modified nucleosides against the adenosine receptors. Like m6A, i6A is found in the human body and may be an endogenous A3R ligand. Our cryo-EM analyses elucidate the A3R-Gi complexes bound to adenosine, 5'-N-ethylcarboxamidoadenosine (NECA), m6A, i6A, and namodenoson at overall resolutions of 3.27 Å (adenosine), 2.86 Å (NECA), 3.19 Å (m6A), 3.28 Å (i6A), and 3.20 Å (namodenoson), suggesting the selectivity and activation mechanism of A3R. We further conduct structure-guided engineering of m6A-insensitive A3R, which may aid future research targeting m6A and A3R, providing a molecular basis for future drug discovery.

Constitutive activity of an atypical chemokine receptor revealed by inverse agonistic nanobodies

Chemokine stimulation of atypical chemokine receptor 3 (ACKR3) does not activate G proteins but recruits arrestins. It is a chemokine scavenger that indirectly influences responses by restricting the availability of CXCL12, an agonist shared with the canonical receptor CXCR4. ACKR3 is upregulated in numerous disorders. Due to limited insights in chemokine-activated ACKR3 signaling, it is unclear how ACKR3 contributes to pathological phenotypes. One explanation may be that high constitutive activity of ACKR3 drives non-canonical signaling through a basal receptor state. Here we characterize the constitutive action of ACKR3 using novel inverse agonistic nanobodies to suppress basal activity. These new tools promote an inactive receptor conformation which decreased arrestin engagement and inhibited constitutive internalization. Basal, non-chemotactic, breast cancer cell motility was also suppressed, suggesting a role for ACKR3 in this process. The basal receptor activity in pathophysiology may provide a new therapeutic approach for targeting ACKR3.

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

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

Non-canonical olfactory pathway activation induces cell fusion of cervical cancer cells

Multinucleation occurs in various types of advanced cancers and contributes to their malignant characteristics, including anticancer drug resistance. Therefore, inhibiting multinucleation can improve cancer prognosis; however, the molecular mechanisms underlying multinucleation remain elusive. Here, we introduced a genetic mutation in cervical cancer cells to induce cell fusion-mediated multinucleation. The olfactory receptor OR1N2 was heterozygously mutated in these fused cells; the same OR1N2 mutation was detected in multinucleated cells from clinical cervical cancer specimens. The mutation-induced structural change in the OR1N2 protein activated protein kinase A (PKA), which, in turn, mediated the non-canonical olfactory pathway. PKA phosphorylated and activated furin protease, resulting in the cleavage of the fusogenic protein syncytin-1. Because this cleaved form of syncytin-1, processed by furin, participates in cell fusion, furin inhibitors could suppress multinucleation and reduce surviving cell numbers after anticancer drug treatment. The improved anticancer drug efficacy indicates a promising therapeutic approach for advanced cervical cancers.

Evolutionary history of calcium-sensing receptors unveils hyper/hypocalcemia-causing mutations

Despite advancements in understanding the structure and functions of the Calcium Sensing Receptor (CaSR), gaps persist in our knowledge of the specific functions of its residues. In this study, we used phylogeny-based techniques to identify functionally equivalent orthologs of CaSR, predict residue significance, and compute specificity-determining position (SDP) scores to understand its evolutionary basis. The analysis revealed exceptional conservation of the CaSR subfamily, emphasizing the critical role of residues with high SDP scores in receptor activation and pathogenicity. To further enhance the findings, gradient-boosting trees were applied to differentiate between gain- and loss-of-function mutations responsible for hypocalcemia and hypercalcemia. Lastly, we investigated the importance of these mutations in the context of receptor activation dynamics. In summary, through comprehensive exploration of the evolutionary history of the CaSR subfamily, coupled with innovative phylogenetic methodologies, we identified activating and inactivating residues, providing valuable insights into the regulation of calcium homeostasis and its connections to associated disorders.

Dynamic Mechanism for Subtype Selectivity of Endocannabinoids

Endocannabinoids are naturally occurring lipid-like molecules that bind to cannabinoid receptors (CB1 and CB2) and regulate many of human bodily functions via the endocannabinoid system. There is a tremendous interest in developing selective drugs that target the CB receptors. However, the biophysical mechanisms responsible for the subtype selectivity for endocannbinoids have not been established. Recent experimental structures of CB receptors show that endocannbinoids potentially bind via membrane using the lipid access channel in the transmembrane region of the receptors. Furthermore, the N-terminus of the receptor could move in and out of the binding pocket thereby modulating both the pocket volume and its residue composition. On the basis of these observations, we propose two hypothesis to explain the selectivity of the endocannabinoid, anandamide for CB1 receptor. First, the selectivity arises from distinct enthalpic ligand-protein interactions along the ligand binding pathway formed due to the movement of N-terminus and subsequent shifts in the binding pocket composition. Second, selectivity arises from the volumetric differences in the binding pocket allowing for differences in ligand conformational entropy. To quantitatively test these hypothesis, we perform extensive molecular dynamics simulations (∼0.9 milliseconds) along with Markov state modeling and deep learning-based VAMP-nets to provide an interpretable characterization of the anandamide binding process to cannabinoid receptors and explain its selectivity for CB1. Our findings reveal that the distinct N-terminus positions along lipid access channels between TM1 and TM7 lead to different binding mechanisms and interactions between anandamide and the binding pocket residues. To validate the critical stabilizing interactions along the binding pathway, relative free energy calculations of anandamide analogs are used. Moreover, the larger CB2 pocket volume increases the entropic effects of ligand binding by allowing higher ligand fluctuations but reduced stable interactions. Therefore, the opposing enthalpy and entropy effects between the receptors shape the endocannabinoid selectivity. Overall, the CB1 selectivity of anandamide is explained by the dominant enthalpy contributions due to ligand-protein interactions in stable binding poses. This study shed lights on potential selectivity mechanisms for endocannabinoids that would aid in the discovery of CB selective drugs.

A Multi-Angle Approach to Predict Peptide-GPCR Complexes: The N/OFQ-NOP System as a Successful AlphaFold Application Case Study

With nearly 700 structures solved and a growing number of customized structure prediction algorithms being developed at a fast pace, G protein-coupled receptors (GPCRs) are an optimal test case for validating new approaches for the prediction of receptor active state and ligand bioactive conformation complexes. In this study, we leveraged the availability of hundreds of peptide GPCRs in the active state and both classical homology and artificial intelligence (AI) based protein modeling combined with docking and AI-based peptide structure prediction approaches to predict the nociceptin/orphanin FQ-NOP receptor active state complex (N/OFQ-NOPa). The In Silico generated hypotheses were validated via the design, synthesis, and pharmacological characterization of novel linear N/OFQ(1-13)-NH2 analogues, leading to the discovery of a novel antagonist (3B; pKB = 6.63) bearing a single ring-constrained residue in place of the Gly2-Gly3 motif of the N/OFQ message sequence (FGGF). While the experimental validation was ongoing, the availability of the Cryo-EM structure of the predicted complex enabled us to unambiguously validate the generated hypotheses. To the best of our knowledge, this is the first example of a peptide-GPCR complex predicted with atomistic accuracy (full complex Cα RMSD < 1.0 Å) and of the N/OFQ message moiety being successfully modified with a rigid scaffold.

Exploring Endocannabinoid System: Unveiling New Roles in Modulating ER Stress

The endoplasmic reticulum (ER) is the organelle mainly involved in maintaining cellular homeostasis and driving correct protein folding. ER-dependent defects or dysfunctions are associated with the genesis/progression of several pathological conditions, including cancer, inflammation, and neurodegenerative disorders, that are directly or indirectly correlated to a wide set of events collectively named under the term "ER stress". Despite the recent increase in interest concerning ER activity, further research studies are needed to highlight all the mechanisms responsible for ER failure. In this field, recent discoveries paved the way for the comprehension of the strong interaction between ER stress development and the endocannabinoid system. The activity of the endocannabinoid system is mediated by the activation of cannabinoid receptors (CB), G protein-coupled receptors that induce a decrease in cAMP levels, with downstream anti-inflammatory effects. CB activation drives, in most cases, the recovery of ER homeostasis through the regulation of ER stress hallmarks PERK, ATF6, and IRE1. In this review, we focus on the CB role in modulating ER stress, with particular attention to the cellular processes leading to UPR activation and oxidative stress response extinguishment, and to the mechanisms underlying natural cannabinoids' modulation of this complex cellular machine.

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

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

Ligand-induced conformational changes in the β1-adrenergic receptor revealed by hydrogen-deuterium exchange mass spectrometry

G Protein Coupled Receptors (GPCRs) constitute the largest family of signalling proteins responsible for translating extracellular stimuli into intracellular functions. They play crucial roles in numerous physiological processes and are major targets for drug discovery. Dysregulation of GPCRs is implicated in various diseases, making understanding their structural dynamics critical for therapeutic development. Here, we use Hydrogen Deuterium Exchange Mass Spectrometry (HDX-MS) to explore the structural dynamics of the turkey β1-adrenergic receptor (tβ1AR) bound with nine different ligands, including agonists, partial agonists, and antagonists. We find that these ligands induce distinct dynamic patterns across the receptor, which can be grouped by compound modality. Notably, full agonist binding destabilises the intracellular loop 1 (ICL1), while antagonist binding stabilises it, highlighting ICL1's role in G protein recruitment. Our findings indicate that the conserved L72 residue in ICL1 is crucial for maintaining receptor structural integrity and stabilising the GDP-bound state. Overall, our results provide a platform for determining drug modality and highlight how HDX-MS can be used to dissect receptor ligand interaction properties and GPCR mechanism.

Active state structures of a bistable visual opsin bound to G proteins

Opsins are G protein-coupled receptors (GPCRs) that have evolved to detect light stimuli and initiate intracellular signaling cascades. Their role as signal transducers is critical to light perception across the animal kingdom. Opsins covalently bind to the chromophore 11-cis retinal, which isomerizes to the all-trans isomer upon photon absorption, causing conformational changes that result in receptor activation. Monostable opsins, responsible for vision in vertebrates, release the chromophore after activation and must bind another retinal molecule to remain functional. In contrast, bistable opsins, responsible for non-visual light perception in vertebrates and for vision in invertebrates, absorb a second photon in the active state to return the chromophore and protein to the inactive state. Structures of bistable opsins in the activated state have proven elusive, limiting our understanding of how they function as bidirectional photoswitches. Here we present active state structures of a bistable opsin, jumping spider rhodopsin isoform-1 (JSR1), in complex with its downstream signaling partners, the Gi and Gq heterotrimers. These structures elucidate key differences in the activation mechanisms between monostable and bistable opsins, offering essential insights for the rational engineering of bistable opsins into diverse optogenetic tools to control G protein signaling pathways.

Structural basis of psychedelic LSD recognition at dopamine D1 receptor

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

Selective ligand recognition and activation of somatostatin receptors SSTR1 and SSTR3

Somatostatin receptors (SSTRs) exert critical biological functions such as negatively regulating hormone release and cell proliferation, making them popular targets for developing therapeutics to treat endocrine disorders, especially neuroendocrine tumors. Although several panagonists mimicking the endogenous ligand somatostatin are available, the development of more effective and safer somatostatinergic therapies is limited due to a lack of molecular understanding of the ligand recognition and regulation of divergent SSTR subtypes. Here, we report four cryoelectron microscopy structures of Gi-coupled SSTR1 and SSTR3 activated by distinct agonists, including the FDA-approved panagonist pasireotide as well as their selective small molecule agonists L-797591 and L-796778. Our structures reveal a conserved recognition pattern of pasireotide in SSTRs attributed to the binding with a conserved extended binding pocket, distinct from SST14, octreotide, and lanreotide. Together with mutagenesis analyses, our structures further reveal the dynamic feature of ligand binding pockets in SSTR1 and SSTR3 to accommodate divergent agonists, the key determinants of ligand selectivity lying across the orthosteric pocket of different SSTR subtypes, as well as the molecular mechanism underlying diversity and conservation of receptor activation. Our work provides a framework for rational design of subtype-selective SSTR ligands and may facilitate drug development efforts targeting SSTRs with improved therapeutic efficacy and reduced side effects.

Homodimerization of CB2 cannabinoid receptor triggered by a bivalent ligand enhances cellular signaling

G protein-coupled receptors (GPCRs) exist within a landscape of interconvertible conformational states and in dynamic equilibrium between monomers and higher-order oligomers, both influenced by ligand binding. Here, we show that a homobivalent ligand formed by equal chromenopyrazole moieties as pharmacophores, connected by 14 methylene units, can modulate the dynamics of the cannabinoid CB2 receptor (CB2R) homodimerization by simultaneously binding both protomers of the CB2R-CB2R homodimer. Computational and pharmacological experiments showed that one of the ligand pharmacophores binds to the orthosteric site of one protomer, and the other pharmacophore to a membrane-oriented pocket between transmembranes 1 and 7 of the partner protomer. This results in unique pharmacological properties, including increased potency in Gi-mediated signaling and enhanced recruitment of β-arrestin. Thus, by modulating dimerization dynamics, it may be possible to fine-tune CB2R activity, potentially leading to improved therapeutic outcomes.

Structural basis of tethered agonism and G protein coupling of protease-activated receptors

Protease-activated receptors (PARs) are a unique group within the G protein-coupled receptor superfamily, orchestrating cellular responses to extracellular proteases via enzymatic cleavage, which triggers intracellular signaling pathways. Protease-activated receptor 1 (PAR1) is a key member of this family and is recognized as a critical pharmacological target for managing thrombotic disorders. In this study, we present cryo-electron microscopy structures of PAR1 in its activated state, induced by its natural tethered agonist (TA), in complex with two distinct downstream proteins, the Gq and Gi heterotrimers, respectively. The TA peptide is positioned within a surface pocket, prompting PAR1 activation through notable conformational shifts. Contrary to the typical receptor activation that involves the outward movement of transmembrane helix 6 (TM6), PAR1 activation is characterized by the simultaneous downward shift of TM6 and TM7, coupled with the rotation of a group of aromatic residues. This results in the displacement of an intracellular anion, creating space for downstream G protein binding. Our findings delineate the TA recognition pattern and highlight a distinct role of the second extracellular loop in forming β-sheets with TA within the PAR family, a feature not observed in other TA-activated receptors. Moreover, the nuanced differences in the interactions between intracellular loops 2/3 and the Gα subunit of different G proteins are crucial for determining the specificity of G protein coupling. These insights contribute to our understanding of the ligand binding and activation mechanisms of PARs, illuminating the basis for PAR1's versatility in G protein coupling.

Distinct binding hotspots for natural and synthetic agonists of FFA4 from in silico approaches

FFA4 has gained interest in recent years since its deorphanization in 2005 and the characterization of the Free Fatty Acids receptors family for their therapeutic potential in metabolic disorders. The expression of FFA4 (also known as GPR120) in numerous organs throughout the human body makes this receptor a highly potent target, particularly in fat sensing and diet preference. This offers an attractive approach to tackle obesity and related metabolic diseases. Recent cryo-EM structures of the receptor have provided valuable information for a potential active state although the previous studies of FFA4 presented diverging information. We performed molecular docking and molecular dynamics simulations of four agonist ligands, TUG-891, Linoleic acid, α-Linolenic acid, and Oleic acid, based on a homology model. Our simulations, which accumulated a total of 2 μs of simulation, highlighted two binding hotspots at Arg992.64 and Lys293 (ECL3). The results indicate that the residues are located in separate areas of the binding pocket and interact with various types of ligands, implying different potential active states of FFA4 and a highly adaptable binding intra-receptor pocket. This article proposes additional structural characteristics and mechanisms for agonist binding that complement the experimental structures.

Cyclopamine modulates smoothened receptor activity in a binding position dependent manner

Cyclopamine, a natural alkaloid, can act as an agonist when it binds to the Cysteine-Rich Domain (CRD) of Smoothened receptor and as an antagonist when it binds to the Transmembrane Domain (TMD). To study the effect of cyclopamine binding to each site experimentally, mutations in the other site are required. Hence, simulations are critical for understanding the WT activity due to binding at different sites. Using multi-milliseconds long aggregate MD simulations combined with Markov state models and machine learning, we explore the dynamic behavior of cyclopamine's interactions with different domains of WT SMO. A higher population of the active state at equilibrium, a lower free energy barrier of ~2 kcal/mol, and expansion of hydrophobic tunnel to facilitate cholesterol transport agrees with cyclopamine's agonistic behavior when bound to CRD. A higher population of the inactive state at equilibrium, a higher free energy barrier of ~4 kcal/mol and restricted hydrophobic tunnel shows cyclopamine's antagonistic behavior when bound to TMD. With cyclopamine bound to both sites, there is a slightly larger inactive population at equilibrium and an increased free energy barrier (~3.5 kcal/mol) exhibiting an overall weak antagonistic effect. These findings show cyclopamine's domain-specific modulation of SMO regulates Hedgehog signaling and cholesterol transport.

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

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

The high-light-sensitivity mechanism and optogenetic properties of the bacteriorhodopsin-like channelrhodopsin GtCCR4

Channelrhodopsins are microbial light-gated ion channels that can control the firing of neurons in response to light. Among several cation channelrhodopsins identified in Guillardia theta (GtCCRs), GtCCR4 has higher light sensitivity than typical channelrhodopsins. Furthermore, GtCCR4 shows superior properties as an optogenetic tool, such as minimal desensitization. Our structural analyses of GtCCR2 and GtCCR4 revealed that GtCCR4 has an outwardly bent transmembrane helix, resembling the conformation of activated G-protein-coupled receptors. Spectroscopic and electrophysiological comparisons suggested that this helix bend in GtCCR4 omits channel recovery time and contributes to high light sensitivity. An electrophysiological comparison of GtCCR4 and the well-characterized optogenetic tool ChRmine demonstrated that GtCCR4 has superior current continuity and action-potential spike generation with less invasiveness in neurons. We also identified highly active mutants of GtCCR4. These results shed light on the diverse structures and dynamics of microbial rhodopsins and demonstrate the strong optogenetic potential of GtCCR4.

Integrative residue-intuitive machine learning and MD Approach to Unveil Allosteric Site and Mechanism for β2AR

Allosteric drugs offer a new avenue for modern drug design. However, the identification of cryptic allosteric sites presents a formidable challenge. Following the allostery nature of residue-driven conformation transition, we propose a state-of-the-art computational pipeline by developing a residue-intuitive hybrid machine learning (RHML) model coupled with molecular dynamics (MD) simulation, through which we can efficiently identify the allosteric site and allosteric modulator as well as reveal their regulation mechanism. For the clinical target β2-adrenoceptor (β2AR), we discover an additional allosteric site located around residues D792.50, F2826.44, N3187.45 and S3197.46 and one putative allosteric modulator ZINC5042. Using Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and protein structure network (PSN), the allosteric potency and regulation mechanism are probed to further improve identification accuracy. Benefiting from sufficient computational evidence, the experimental assays then validate our predicted allosteric site, negative allosteric potency and regulation pathway, showcasing the effectiveness of the identification pipeline in practice. We expect that it will be applicable to other target proteins.

Progress on the development of Class A GPCR-biased ligands

Class A G protein-coupled receptors (GPCRs) continue to garner interest for their essential roles in cell signalling and their importance as drug targets. Although numerous drugs in the clinic target these receptors, over 60% GPCRs remain unexploited. Moreover, the adverse effects triggered by the available unbiased GPCR modulators, limit their use and therapeutic value. In this context, the elucidation of biased signalling has opened up new pharmacological avenues holding promise for safer therapeutics. Functionally selective ligands favour receptor conformations facilitating the recruitment of specific effectors and the modulation of the associated pathways. This review surveys the current drug discovery landscape of GPCR-biased modulators with a focus on recent advances. Understanding the biological effects of this preferential coupling is at different stages depending on the Class A GPCR family. Therefore, with a focus on individual GPCR families, we present a compilation of the functionally selective modulators reported over the past few years. In doing so, we dissect their therapeutic relevance, molecular determinants and potential clinical applications.

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.

Allosteric modulation of the CXCR4:CXCL12 axis by targeting receptor nanoclustering via the TMV-TMVI domain

CXCR4 is a ubiquitously expressed chemokine receptor that regulates leukocyte trafficking and arrest in both homeostatic and pathological states. It also participates in organogenesis, HIV-1 infection, and tumor development. Despite the potential therapeutic benefit of CXCR4 antagonists, only one, plerixafor (AMD3100), which blocks the ligand-binding site, has reached the clinic. Recent advances in imaging and biophysical techniques have provided a richer understanding of the membrane organization and dynamics of this receptor. Activation of CXCR4 by CXCL12 reduces the number of CXCR4 monomers/dimers at the cell membrane and increases the formation of large nanoclusters, which are largely immobile and are required for correct cell orientation to chemoattractant gradients. Mechanistically, CXCR4 activation involves a structural motif defined by residues in TMV and TMVI. Using this structural motif as a template, we performed in silico molecular modeling followed by in vitro screening of a small compound library to identify negative allosteric modulators of CXCR4 that do not affect CXCL12 binding. We identified AGR1.137, a small molecule that abolishes CXCL12-mediated receptor nanoclustering and dynamics and blocks the ability of cells to sense CXCL12 gradients both in vitro and in vivo while preserving ligand binding and receptor internalization.

pKa Calculations of GPCRs: Understanding Protonation States in Receptor Activation

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

Utilizing Molecular Dynamics Simulations, Machine Learning, Cryo-EM, and NMR Spectroscopy to Predict and Validate Protein Dynamics

Protein dynamics play a crucial role in biological function, encompassing motions ranging from atomic vibrations to large-scale conformational changes. Recent advancements in experimental techniques, computational methods, and artificial intelligence have revolutionized our understanding of protein dynamics. Nuclear magnetic resonance spectroscopy provides atomic-resolution insights, while molecular dynamics simulations offer detailed trajectories of protein motions. Computational methods applied to X-ray crystallography and cryo-electron microscopy (cryo-EM) have enabled the exploration of protein dynamics, capturing conformational ensembles that were previously unattainable. The integration of machine learning, exemplified by AlphaFold2, has accelerated structure prediction and dynamics analysis. These approaches have revealed the importance of protein dynamics in allosteric regulation, enzyme catalysis, and intrinsically disordered proteins. The shift towards ensemble representations of protein structures and the application of single-molecule techniques have further enhanced our ability to capture the dynamic nature of proteins. Understanding protein dynamics is essential for elucidating biological mechanisms, designing drugs, and developing novel biocatalysts, marking a significant paradigm shift in structural biology and drug discovery.

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

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

A bitopic agonist bound to the dopamine 3 receptor reveals a selectivity site

Although aminergic GPCRs are the target for ~25% of approved drugs, developing subtype selective drugs is a major challenge due to the high sequence conservation at their orthosteric binding site. Bitopic ligands are covalently joined orthosteric and allosteric pharmacophores with the potential to boost receptor selectivity and improve current medications by reducing off-target side effects. However, the lack of structural information on their binding mode impedes rational design. Here we determine the cryo-EM structure of the hD3R:GαOβγ complex bound to the D3R selective bitopic agonist FOB02-04A. Structural, functional and computational analyses provide insights into its binding mode and point to a new TM2-ECL1-TM1 region, which requires the N-terminal ordering of TM1, as a major determinant of subtype selectivity in aminergic GPCRs. This region is underexploited in drug development, expands the established secondary binding pocket in aminergic GPCRs and could potentially be used to design novel and subtype selective drugs.

Distinct binding conformations of epinephrine with α- and β-adrenergic receptors

Agonists targeting α2-adrenergic receptors (ARs) are used to treat diverse conditions, including hypertension, attention-deficit/hyperactivity disorder, pain, panic disorders, opioid and alcohol withdrawal symptoms, and cigarette cravings. These receptors transduce signals through heterotrimeric Gi proteins. Here, we elucidated cryo-EM structures that depict α2A-AR in complex with Gi proteins, along with the endogenous agonist epinephrine or the synthetic agonist dexmedetomidine. Molecular dynamics simulations and functional studies reinforce the results of the structural revelations. Our investigation revealed that epinephrine exhibits different conformations when engaging with α-ARs and β-ARs. Furthermore, α2A-AR and β1-AR (primarily coupled to Gs, with secondary associations to Gi) were compared and found to exhibit different interactions with Gi proteins. Notably, the stability of the epinephrine-α2A-AR-Gi complex is greater than that of the dexmedetomidine-α2A-AR-Gi complex. These findings substantiate and improve our knowledge on the intricate signaling mechanisms orchestrated by ARs and concurrently shed light on the regulation of α-ARs and β-ARs by epinephrine.

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.

Delineating the stepwise millisecond allosteric activation mechanism of the class C GPCR dimer mGlu5

Two-thirds of signaling hormones and one-third of approved drugs exert their effects by binding and modulating the G protein-coupled receptors (GPCRs) activation. While the activation mechanism for monomeric GPCRs has been well-established, little is known about GPCRs in dimeric form. Here, by combining transition pathway generation, extensive atomistic simulation-based Markov state models, and experimental signaling assays, we reveal an asymmetric, stepwise millisecond allosteric activation mechanism for the metabotropic glutamate receptor subtype 5 receptor (mGlu5), an obligate dimeric class C GPCR. The dynamic picture is presented that agonist binding induces dimeric ectodomains compaction, amplified by the precise association of the cysteine-rich domains, ultimately loosely bringing the intracellular 7-transmembrane (7TM) domains into proximity and establishing an asymmetric TM6-TM6 interface. The active inter-domain interface enhances their intra-domain flexibility, triggering the activation of micro-switches crucial for downstream signal transduction. Furthermore, we show that the positive allosteric modulator stabilizes both the active inter-domain 7TM interface and an open, extended intra-domain ICL2 conformation. This stabilization leads to the formation of a pseudo-cavity composed of the ICL2, ICL3, TM3, and C-terminus, which facilitates G protein coordination. Our strategy may be generalizable for characterizing millisecond events in other allosteric systems.

Structural basis of frizzled 7 activation and allosteric regulation

Frizzleds (ten paralogs: FZD1-10) belong to the class F of G protein-coupled receptors (GPCRs), which remains poorly understood despite its crucial role in multiple key biological functions including embryonic development, stem cell regulation, and homeostasis in the adult. FZD7, one of the most studied members of the family, is more specifically involved in the migration of mesendoderm cells during the development and renewal of intestinal stem cells in adults. Moreover, FZD7 has been highlighted for its involvement in tumor development predominantly in the gastrointestinal tract. This study reports the structure of inactive FZD7, without any stabilizing mutations, determined by cryo-electron microscopy (cryo-EM) at 1.9 Å resolution. We characterize a fluctuating water pocket in the core of the receptor important for FZD7 dynamics. Molecular dynamics simulations are used to investigate the temporal distribution of those water molecules and their importance for potential conformational changes in FZD7. Moreover, we identify lipids interacting with the receptor core and a conserved cholesterol-binding site, which displays a key role in FZD7 association with a transducer protein, Disheveled (DVL), and initiation of downstream signaling and signalosome formation.

Structural insights into ligand recognition, selectivity, and activation of bombesin receptor subtype-3

Bombesin receptor subtype-3 (BRS3) is an important orphan G protein-coupled receptor that regulates energy homeostasis and insulin secretion. As a member of the bombesin receptor (BnR) family, the lack of known endogenous ligands and high-resolution structure has hindered the understanding of BRS3 signaling and function. We present two cryogenic electron microscopy (cryo-EM) structures of BRS3 in complex with the heterotrimeric Gq protein in its active states: one bound to the pan-BnR agonist BA1 and the other bound to the synthetic BRS3-specific agonist MK-5046. These structures reveal the architecture of the orthosteric ligand pocket underpinning molecular recognition and provide insights into the structural basis for BRS3's selectivity and low affinity for bombesin peptides. Examination of conserved micro-switches suggests a shared activation mechanism among BnRs. Our findings shed light on BRS3's ligand selectivity and signaling mechanisms, paving the way for exploring its therapeutic potential for diabetes, obesity, and related metabolic disorders.