[19] Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors

Positive allosteric modulation of the cannabinoid CB1 receptor potentiates endocannabinoid signalling and changes ERK1/2 phosphorylation kinetics

BACKGROUND AND PURPOSE

Activation of CB1 by exogenous agonists causes adverse effects in vivo. Positive allosteric modulation may offer improved therapeutic potential and a reduced on-target adverse effect profile compared with orthosteric agonists, due to reduced desensitisation/tolerance, but this has not been directly tested. This study investigated the ability of PAMs/ago-PAMs to induce receptor regulation pathways, including desensitisation and receptor internalisation.

EXPERIMENTAL APPROACH

Bioluminescence resonance energy transfer (BRET) assays in HEK293 cells were performed to investigate G protein dissociation, ERK1/2 phosphorylation and β-arrestin 2 translocation, while immunocytochemistry was performed to measure internalisation of CB1 in response to the PAMs ZCZ011, GAT229 and ABD1236 alone and in combination with the orthosteric agonists AEA, 2-AG, and AMB-FUBINACA.

KEY RESULTS

ZCZ011, GAT229 and ABD1236 were allosteric agonists in all pathways tested. The ago-PAM ZCZ011 induced a biphasic ERK1/2 phosphorylation time course compared to transient activation by orthosteric agonists. In combination with 2-AG but not AEA or AMB-FUBINACA, ZCZ011 and ABD1236 caused the transient peak of ERK1/2 phosphorylation to become sustained. All PAMs increased the potency and efficacy of AEA-induced signalling in all pathways tested; however, no notable potentiation of 2-AG or AMB-FUBINACA was observed.

CONCLUSION AND IMPLICATIONS

Ago-PAMs can potentiate endocannabinoid CB1 agonism by AEA to a larger extent compared with 2-AG. However, all compounds were found to be allosteric agonists and induce activation of CB1 in the absence of endocannabinoid, including β-arrestin 2 recruitment and internalisation. Thus, the spatiotemporal signalling of endogenous cannabinoids will not be retained in vivo.

Scutellarein derivatives with histamine H3 receptor antagonism and cholinesterase inhibitory potency as multi target-directed ligands for possible Alzheimer's disease therapy

A series of scutellarein 7-l-amino acid carbamate-4'-cycloalkylamine propyl ether conjugates were designed and synthesized for the first time as multifunctional agents for Alzheimer's disease (AD) therapy. The designed compounds exhibited more balanced and effective multi-target potency. Among them, compound 11l, l-Valine carbamate derivative of scutellarein cycloheptylamine ether, exhibited the most potent inhibition of electric eel AChE enzymes and human AChE enzymes, with an IC50 values of 7.04 μM and 9.73 μM, respectively. Moreover, 11l exhibited more potent H3R antagonistic activities than clobenpropit, with an IC50 value of 1.09 nM. Compound 11l not only displayed excellent inhibition of self- and Cu2+-induced Aβ1-42 aggregation (95.48 % and 88.63 % inhibition, respectively) but also induced the disassembly of self- and Cu2+-induced Aβ fibrils (80.16 % and 89.30 % disaggregation, respectively). Moreover, 11l significantly reduced tau protein hyperphosphorylation induced by Aβ25-35. It exhibited effective antioxidant activity and neuroprotective potency, and inhibited RSL3-induced PC12 cell ferroptosis. Assays of hCMEC/D3 and hPepT1-MDCK cell line permeability indicated that 11l would have optimal blood-brain barrier permeability and intestinal absorption characteristics. In addition, in vivo studies revealed that compound 11l significantly attenuated learning and memory impairment in an AD mouse model. Finally, a pharmacokinetic characterization of 11l indicated favorable druggability and pharmacokinetic properties. Taken together, our results suggest that 11l is a potential candidate for AD treatment and merits further investigation.

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.

Pharmacological characterization of allosteric modulators: A case for chemokine receptors

Chemokine receptors are relevant targets for a multitude of immunological diseases, but drug attrition for these receptors is remarkably high. While many drug discovery programs have been pursued, most prospective drugs failed in the follow-up studies due to clinical inefficacy, and hence there is a clear need for alternative approaches. Allosteric modulators of receptor function represent an excellent opportunity for novel drugs, as they modulate receptor activation in a controlled manner and display increased selectivity, and their pharmacological profile can be insurmountable. Here, we discuss allosteric ligands and their pharmacological characterization for modulation of chemokine receptors. Ligands are included if (1) they show clear signs of allosteric modulation in vitro and (2) display evidence of binding in a topologically distinct manner compared to endogenous chemokines. We discuss how allosteric ligands affect binding of orthosteric (endogenous) ligands in terms of affinity as well as binding kinetics in radioligand binding assays. Moreover, their effects on signaling events in functional assays and how their binding site can be elucidated are specified. We substantiate this with examples of published allosteric ligands targeting chemokine receptors and hypothetical graphs of pharmacological behavior. This review should serve as an effective starting point for setting up assays for characterizing allosteric ligands to develop safer and more efficacious drugs for chemokine receptors and, ultimately, other G protein-coupled receptors.

Molecular recognition of an odorant by the murine trace amine-associated receptor TAAR7f

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin, dopamine and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein Gs and bound to the odorant N,N-dimethylcyclohexylamine (DMCHA) to an overall resolution of 2.9 Å. DMCHA is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the β-adrenoceptors, with the significant difference that DMCHA interacts directly with one of the main activation microswitch residues, Trp6.48.

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.

Structural Basis for the Recognition of GPRC5D by Talquetamab, a Bispecific Antibody for Multiple Myeloma

Multiple myeloma (MM) is a complex hematological malignancy characterized by abnormal antibody production from plasma cells. Despite advances in the treatment, many patients experience disease relapse or become refractory to treatment. G-protein-coupled receptor class C group 5 member D (GPRC5D), an orphan GPCR predominantly expressed in MM cells, is emerging as a promising target for MM immunotherapy. Talquetamab, a Food and Drug Administration-approved T-cell-directing bispecific antibody developed for treatment of MM, targets GPRC5D. Here, we elucidate the structure of GPRC5D complexed with the Fab fragment of talquetamab, using cryo-electron microscopy, providing the basis for recognition of GPRC5D by the bispecific antibody. GPRC5D forms a symmetric homodimer with the interface between transmembrane helix (TM) 4 of one protomer and TM4/5 of the other protomer. A single talquetamab Fab interacts with the GPRC5D dimer with its orientation toward the dimer interface. All six complementarity-determining regions of talquetamab engage with extracellular loops and TM3/5/7. In particular, the side-chain of an arginine residue from the antibody penetrates into a shallow pocket on the extracellular surface of GPRC5D. The structure offers insights for optimizing antibody design against GPRC5D for relapsed or refractory MM therapy.

Biased receptor signalling and intracellular trafficking profiles of structurally distinct formylpeptide receptor 2 agonists

BACKGROUND

There is increasing interest in developing FPR2 agonists (compound 43, ACT-389949 and BMS-986235) as potential pro-resolving therapeutics, with ACT-389949 and BMS-986235 having entered phase I clinical development. FPR2 activation leads to diverse downstream outputs. ACT-389949 was observed to cause rapid tachyphylaxis, while BMS-986235 and compound 43 induced cardioprotective effects in preclinical models. We aim to characterise the differences in ligand-receptor engagement and downstream signalling and trafficking bias profile.

EXPERIMENTAL APPROACH

Concentration-response curves to G protein dissociation, β-arrestin recruitment, receptor trafficking and second messenger signalling were generated using FPR2 ligands (BMS-986235, ACT-389949, compound 43 and WKYMVm), in HEK293A cells. Log(τ/KA) was obtained from the operational model for bias analysis using WKYMVm as a reference ligand. Docking of FPR2 ligands into the active FPR2 cryoEM structure (PDBID: 7T6S) was performed using ICM pro software.

KEY RESULTS

Bias analysis revealed that WKYMVm and ACT-389949 shared a very similar bias profile. In comparison, BMS-986235 and compound 43 displayed approximately 5- to 50-fold bias away from β-arrestin recruitment and trafficking pathways, while being 35- to 60-fold biased towards cAMP inhibition and pERK1/2. Molecular docking predicted key amino acid interactions at the FPR2 shared between WKYMVm and ACT-389949, but not with BMS-986235 and compound 43.

CONCLUSION AND IMPLICATIONS

In vitro characterisation demonstrated that WKYMVm and ACT-389949 differ from BMS-986235 and compound 43 in their signalling and protein coupling profile. This observation may be explained by differences in the ligand-receptor interactions. In vitro characterisation provided significant insights into identifying the desired bias profile for FPR2-based pharmacotherapy.

Monitoring GPCR conformation with GFP-inspired dyes

Solvatochromic compounds have emerged as valuable environment-sensitive probes for biological research. Here we used thiol-reactive solvatochromic analogs of the green fluorescent protein (GFP) chromophore to track conformational changes in two proteins, recoverin and the A2A adenosine receptor (A2AAR). Two dyes showed Ca2+-induced fluorescence changes when attached to recoverin. Our best-performing dye, DyeC, exhibited agonist-induced changes in both intensity and shape of its fluorescence spectrum when attached to A2AAR; none of these effects were observed with other common environment-sensitive dyes. Molecular dynamics simulations showed that activation of the A2AAR led to a more confined and hydrophilic environment for DyeC. Additionally, an allosteric modulator of A2AAR induced distinct fluorescence changes in the DyeC spectrum, indicating a unique receptor conformation. Our study demonstrated that GFP-inspired dyes are effective for detecting structural changes in G protein-coupled receptors (GPCRs), offering advantages such as intensity-based and ratiometric tracking, redshifted fluorescence spectra, and sensitivity to allosteric modulation.

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).

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.

A versatile residue numbering scheme for Nav and Cav channels

Voltage-gated sodium (Nav) and calcium (Cav) channels are responsible for the initiation of electrical signals. They have long been targeted for the treatment of various diseases. The mounting number of cryoelectron microscopy (cryo-EM) structures for diverse subtypes of Nav and Cav channels from multiple organisms necessitates a generic residue numbering system to establish the structure-function relationship and to aid rational drug design or optimization. Here we suggest a structure-based residue numbering scheme, centering around the most conserved residues on each of the functional segments. We elaborate the generic numbers through illustrative examples, focusing on representative drug-binding sites of eukaryotic Nav and Cav channels. We also extend the numbering scheme to compare common disease mutations among different Nav subtypes. Application of the generic residue numbering scheme affords immediate insights into hotspots for pathogenic mutations and critical loci for drug binding and will facilitate drug discovery targeting Nav and Cav channels.

Exploring an Intracellular Allosteric Site of CC-Chemokine Receptor 4 from 3D Models, Probe Simulations, and Mutagenesis

We applied our previously developed probe confined dynamic mapping protocol, which combines enhanced sampling molecular dynamics (MD) simulations and fragment-based approaches, to identify the binding site of GSK2239633A (N-[[3-[[3-[(5-chlorothiophen-2-yl)sulfonylamino]-4-methoxyindazol-1-yl]methyl]phenyl]methyl]-2-hydroxy-2-methylpropanamide), a selective CC-chemokine receptor type 4 (CCR4) negative allosteric modulator, using CCR4 homology and AlphaFold models. By comparing the performance across five computational models, we identified conserved (K3108.49 and Y3047.53) and non-conserved (M2436.36) residue hotspots for GSK2239633A binding, which were validated by mutagenesis and bioluminescence resonance energy transfer assay. Further analysis of 3D models and MD simulations highlighted the pair of residues 6.36 and 7.56 that might account for antagonist selectivity among chemokine receptors. Our in silico protocol provides a promising approach for characterizing ligand binding sites in membrane proteins, considering receptor dynamics and adaptability and guiding protein template selection for ligand design.

AlphaFold accelerated discovery of psychotropic agonists targeting the trace amine-associated receptor 1

Artificial intelligence is revolutionizing protein structure prediction, providing unprecedented opportunities for drug design. To assess the potential impact on ligand discovery, we compared virtual screens using protein structures generated by the AlphaFold machine learning method and traditional homology modeling. More than 16 million compounds were docked to models of the trace amine-associated receptor 1 (TAAR1), a G protein-coupled receptor of unknown structure and target for treating neuropsychiatric disorders. Sets of 30 and 32 highly ranked compounds from the AlphaFold and homology model screens, respectively, were experimentally evaluated. Of these, 25 were TAAR1 agonists with potencies ranging from 12 to 0.03 μM. The AlphaFold screen yielded a more than twofold higher hit rate (60%) than the homology model and discovered the most potent agonists. A TAAR1 agonist with a promising selectivity profile and drug-like properties showed physiological and antipsychotic-like effects in wild-type but not in TAAR1 knockout mice. These results demonstrate that AlphaFold structures can accelerate drug discovery.

Hit-to-Lead Optimization of Heterocyclic Carbonyloxycarboximidamides as Selective Antagonists at Human Adenosine A3 Receptor

Antagonism of the human adenosine A3 receptor (hA3R) has potential therapeutic application. Alchemical relative binding free energy calculations of K18 and K32 suggested that the combination of a 3-(2,6-dichlorophenyl)-isoxazolyl group with 2-pyridinyl at the ends of a carbonyloxycarboximidamide group should improve hA3R affinity. Of the 25 new analogues synthesized, 37 and 74 showed improved hA3R affinity compared to K18 (and K32). This was further improved through the addition of a bromine group to the 2-pyridinyl at the 5-position, generating compound 39. Alchemical relative binding free energy calculations, mutagenesis studies and MD simulations supported the compounds' binding pattern while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, so highlighting the importance of rigidification of the carbonyloxycarboximidamide moiety. MD simulations highlighted the importance of rigidification of the carbonyloxycarboximidamide, while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, which was supported through mutagenesis studies 39 also selectively antagonized endogenously expressed hA3R in nonsmall cell lung carcinoma cells, while pharmacokinetic studies indicated low toxicity enabling in vivo evaluation. We therefore suggest that 39 has potential for further development as a high-affinity hA3R antagonist.

G protein-coupled receptor (GPCR) gene variants and human genetic disease

Genetic variations in the genes encoding G protein-coupled receptors (GPCRs) can disrupt receptor structure and function, which can result in human genetic diseases. Disease-causing mutations have been reported in at least 55 GPCRs for more than 66 monogenic diseases in humans. The spectrum of pathogenic and likely pathogenic variants includes loss of function variants that decrease receptor signaling on one extreme and gain of function that may result in biased signaling or constitutive activity, originally modeled on prototypical rhodopsin GPCR variants identified in retinitis pigmentosa, on the other. GPCR variants disrupt ligand binding, G protein coupling, accessory protein function, receptor desensitization and receptor recycling. Next generation sequencing has made it possible to identify variants of uncertain significance (VUS). We discuss variants in receptors known to result in disease and in silico strategies for disambiguation of VUS such as sorting intolerant from tolerant and polymorphism phenotyping. Modeling of variants has contributed to drug development and precision medicine, including drugs that target the melanocortin receptor in obesity and interventions that reverse loss of gonadotropin-releasing hormone receptor from the cell surface in idiopathic hypogonadotropic hypogonadism. Activating and inactivating variants of the calcium sensing receptor (CaSR) gene that are pathogenic in familial hypocalciuric hypercalcemia and autosomal dominant hypocalcemia have enabled the development of calcimimetics and calcilytics. Next generation sequencing has continued to identify variants in GPCR genes, including orphan receptors, that contribute to human phenotypes and may have therapeutic potential. Variants of the CaSR gene, some encoding an arginine-rich region that promotes receptor phosphorylation and intracellular retention, have been linked to an idiopathic epilepsy syndrome. Agnostic strategies have identified variants of the pyroglutamylated RF amide peptide receptor gene in intellectual disability and G protein-coupled receptor 39 identified in psoriatic arthropathy. Coding variants of the G protein-coupled receptor L1 (GPR37L1) orphan receptor gene have been identified in a rare familial progressive myoclonus epilepsy. The study of the role of GPCR variants in monogenic, Mendelian phenotypes has provided the basis of modeling the significance of more common variants of pharmacogenetic significance.

Structure-guided discovery of orexin receptor-binding PET ligands

Molecular imaging using positron emission tomography (PET) can serve as a promising tool for visualizing biological targets in the brain. Insights into the expression pattern and the in vivo imaging of the G protein-coupled orexin receptors OX1R and OX2R will further our understanding of the orexin system and its role in various physiological and pathophysiological processes. Guided by crystal structures of our lead compound JH112 and the approved hypnotic drug suvorexant bound to OX1R and OX2R, respectively, we herein describe the design and synthesis of two novel radioligands, [18F]KD23 and [18F]KD10. Key to the success of our structural modifications was a bioisosteric replacement of the triazole moiety with a fluorophenyl group. The 19F-substituted analog KD23 showed high affinity for the OX1R and selectivity over OX2R, while the high affinity ligand KD10 displayed similar Ki values for both subtypes. Radiolabeling starting from the respective pinacol ester precursors resulted in excellent radiochemical yields of 93% and 88% for [18F]KD23 and [18F]KD10, respectively, within 20 min. The new compounds will be useful in PET studies aimed at subtype-selective imaging of orexin receptors in brain tissue.

The Inferential Binding Sites of GCGR for Small Molecules Using Protein Dynamic Conformations and Crystal Structures

Glucagon receptor (GCGR) is a class B1 G-protein-coupled receptor that plays a crucial role in maintaining human blood glucose homeostasis and is a significant target for the treatment of type 2 diabetes mellitus (T2DM). Currently, six small molecules (Bay 27-9955, MK-0893, MK-3577, LY2409021, PF-06291874, and LGD-6972) have been tested or are undergoing clinical trials, but only the binding site of MK-0893 has been resolved. To predict binding sites for other small molecules, we utilized both the crystal structure of the GCGR and MK-0893 complex and dynamic conformations. We docked five small molecules and selected the best conformation based on binding mode, docking score, and binding free energy. We performed MD simulations to verify the binding mode of the selected small molecules. Moreover, when selecting conformations, results of competitive binding were referred to. MD simulation indicated that Bay 27-9955 exhibits moderate binding stability in Pocket 3. MK-3577, LY2409021, and PF-06291874 exhibited highly stable binding to Pocket 2, consistent with experimental results. However, LY2409021 may also bind to Pocket 5. Additionally, LGD-6972 exhibited relatively stable binding in Pocket 5. We also conducted structural modifications of LGD-6972 based on the results of MD simulations and predicted its analogues' bioavailability, providing a reference for the study of GCGR small molecules.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

Lipid Trolling to Optimize A3 Adenosine Receptor-Positive Allosteric Modulators (PAMs)

A3 adenosine receptor (A3AR) positive allosteric modulators (PAMs) (2,4-disubstituted-1H-imidazo[4,5-c]quinolin-4-amines) allosterically increase the Emax of A3AR agonists, but not potency, due to concurrent orthosteric antagonism. Following mutagenesis/homology modeling of the proposed lipid-exposed allosteric binding site on the cytosolic side, we functionalized the scaffold, including heteroatom substitutions and exocyclic phenylamine extensions, to increase allosteric binding. Strategically appended linear alkyl-alkynyl chains with terminal amino/guanidino groups improved allosteric effects at both human and mouse A3ARs. The chain length, functionality, and attachment position were varied to modulate A3AR PAM activity. For example, 26 (MRS8247, p-alkyne-linked 8 methylenes) and homologues increased agonist Cl-IB-MECA's Emax and potency ([35S]GTPγS binding). The putative mechanism involves a flexible, terminally cationic chain penetrating the lipid environment for stable electrostatic anchoring to cytosolic phospholipid head groups, suggesting "lipid trolling", supported by molecular dynamic simulation of the active-state model. Thus, we have improved A3AR PAM activity through rational design based on an extrahelical, lipidic binding site.

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.

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.

Structural insights into ligand recognition and activation of the succinate receptor SUCNR1

Succinate, a citric acid cycle intermediate, serves important functions in energy homeostasis and metabolic regulation. Extracellular succinate acts as a stress signal through succinate receptor (SUCNR1), a class A G protein-coupled receptor. Research on succinate signaling is hampered by the lack of high-resolution structures of the agonist-bound receptor. We present cryoelectron microscopy (cryo-EM) structures of SUCNR1-Gi complexes bound to succinate and its non-metabolite derivative cis-epoxysuccinate. Key determinants for the recognition of succinate in cis conformation include R2817.39 and Y832.64, while Y301.39 and R993.29 participate in the binding of both succinate and cis-epoxysuccinate. Extracellular loop 2, through F175ECL2 in its β-hairpin, forms a hydrogen bond with succinate and caps the binding pocket. At the receptor-Gi interface, agonist binding induces the rearrangement of a hydrophobic network on transmembrane (TM)5 and TM6, leading to TM signaling through TM3 and TM7. These findings extend our understanding of succinate recognition by SUCNR1, aiding the development of therapeutics for the succinate receptor.

Conformational dynamics underlying atypical chemokine receptor 3 activation

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

AR71, Histamine H3 Receptor Ligand-In Vitro and In Vivo Evaluation (Anti-Inflammatory Activity, Metabolic Stability, Toxicity, and Analgesic Action)

The future of therapy for neurodegenerative diseases (NDs) relies on new strategies targeting multiple pharmacological pathways. Our research led to obtaining the compound AR71 [(E)-3-(3,4,5-trimethoxyphenyl)-1-(4-(3-(piperidin-1-yl)propoxy)phenyl)prop-2-en-1-one], which has high affinity for human H3R (Ki = 24 nM) and selectivity towards histamine H1 and H4 receptors (Ki > 2500 nM), and showed anti-inflammatory activity in a model of lipopolysaccharide-induced inflammation in BV-2 cells. The presented tests confirmed its antagonist/inverse agonist activity profile and good metabolic stability while docking studies showed the binding mode to histamine H1, H3, and H4 receptors. In in vitro tests, cytotoxicity was evaluated at three cell lines (neuroblastoma, astrocytes, and human peripheral blood mononuclear cells), and a neuroprotective effect was observed in rotenone-induced toxicity. In vivo experiments in a mouse neuropathic pain model demonstrated the highest analgesic effects of AR71 at the dose of 20 mg/kg body weight. Additionally, AR71 showed antiproliferative activity in higher concentrations. These findings suggest the need for further evaluation of AR71's therapeutic potential in treating ND and CNS cancer using animal experimental models.

Functional consequences of spatial, temporal and ligand bias of G protein-coupled receptors

G protein-coupled receptors (GPCRs) regulate every aspect of kidney function by mediating the effects of various endogenous and exogenous substances. A key concept in GPCR function is biased signalling, whereby certain ligands may selectively activate specific pathways within the receptor's signalling repertoire. For example, different agonists may induce biased signalling by stabilizing distinct active receptor conformations - a concept that is supported by advances in structural biology. However, the processes underlying functional selectivity in receptor signalling are extremely complex, involving differences in subcellular compartmentalization and signalling dynamics. Importantly, the molecular mechanisms of spatiotemporal bias, particularly its connection to ligand binding kinetics, have been detailed for GPCRs critical to kidney function, such as the AT1 angiotensin receptor (AT1R), V2 vasopressin receptor (V2R) and the parathyroid hormone 1 receptor (PTH1R). This expanding insight into the multifaceted nature of biased signalling paves the way for innovative strategies for targeting GPCR functions; the development of novel biased agonists may represent advanced pharmacotherapeutic approaches to the treatment of kidney diseases and related systemic conditions, such as hypertension, diabetes and heart failure.

Geosmin, a Food- and Water-Deteriorating Sesquiterpenoid and Ambivalent Semiochemical, Activates Evolutionary Conserved Receptor OR11A1

Geosmin, a ubiquitous volatile sesquiterpenoid of microbiological origin, is causative for deteriorating the quality of many foods, beverages, and drinking water, by eliciting an undesirable "earthy/musty" off-flavor. Moreover, and across species from worm to human, geosmin is a volatile, chemosensory trigger of both avoidance and attraction behaviors, suggesting its role as semiochemical. Volatiles typically are detected by chemosensory receptors of the nose, which have evolved to best detect ecologically relevant food-related odorants and semiochemicals. An insect receptor for geosmin was recently identified in flies. A human geosmin-selective receptor, however, has been elusive. Here, we report on the identification and characterization of a human odorant receptor for geosmin, with its function being conserved in orthologs across six mammalian species. Notably, the receptor from the desert-dwelling kangaroo rat showed a more than 100-fold higher sensitivity compared to its human ortholog and detected geosmin at low nmol/L concentrations in extracts from geosmin-producing actinomycetes.

GprC of the nematode-trapping fungus Arthrobotrys flagrans activates mitochondria and reprograms fungal cells for nematode hunting

Initiation of development requires differential gene expression and metabolic adaptations. Here we show in the nematode-trapping fungus, Arthrobotrys flagrans, that both are achieved through a dual-function G-protein-coupled receptor (GPCR). A. flagrans develops adhesive traps and recognizes its prey, Caenorhabditis elegans, through nematode-specific pheromones (ascarosides). Gene-expression analyses revealed that ascarosides activate the fungal GPCR, GprC, at the plasma membrane and together with the G-protein alpha subunit GasA, reprograms the cell. However, GprC and GasA also reside in mitochondria and boost respiration. This dual localization of GprC in A. flagrans resembles the localization of the cannabinoid receptor CB1 in humans. The C. elegans ascaroside-sensing GPCR, SRBC66 and GPCRs of many fungi are also predicted for dual localization, suggesting broad evolutionary conservation. An SRBC64/66-GprC chimaeric protein was functional in A. flagrans, and C. elegans SRBC64/66 and DAF38 share ascaroside-binding sites with the fungal GprC receptor, suggesting 400-million-year convergent evolution.

Molecular Basis of MC1R Activation: Mutation-Induced Alterations in Structural Dynamics

The MC1R protein is a receptor found in melanocytes that plays a role in melanin synthesis. Mutations in this protein can impact hair color, skin tone, tanning ability, and increase the risk of skin cancer. The MC1R protein is activated by the alpha-melanocyte-stimulating hormone (α-MSH). Previous studies have shown that mutations affect the interaction between MC1R and α-MSH; however, the mechanism behind this process is poorly understood. Our study aims to shed light on this mechanism using molecular dynamics (MD) simulations to analyze the Asp84Glu and Asp294His variants. We simulated both the wild-type (WT) protein and the mutants with and without ligand. Our results reveal that mutations induce unique conformations during state transitions, hindering the switch between active and inactive states and decreasing cellular levels of cAMP. Interestingly, Asp294His showed increased ligand affinity but decreased protein activity, highlighting that tighter binding does not always lead to increased activation. Our study provides insights into the molecular mechanisms underlying the impact of MC1R mutations on protein activity.

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Structural basis of tolvaptan binding to the vasopressin V2 receptor

The vasopressin V2 receptor (V2R) is a validated therapeutic target for autosomal dominant polycystic kidney disease (ADPKD), with tolvaptan being the first FDA-approved antagonist. Herein, we used Gaussian accelerated molecular dynamics simulations to investigate the spontaneous binding of tolvaptan to both active and inactive V2R conformations at the atomic-level. Overall, the binding process consists of two stages. Tolvaptan binds initially to extracellular loops 2 and 3 (ECL2/3) before overcoming an energy barrier to enter the pocket. Our simulations result highlighted key residues (e.g., R181, Y205, F287, F178) involved in this process, which were experimentally confirmed by site-directed mutagenesis. This work provides structural insights into tolvaptan-V2R interactions, potentially aiding the design of novel antagonists for V2R and other G protein-coupled receptors.

Structure and dynamics of the pyroglutamylated RF-amide peptide QRFP receptor GPR103

Pyroglutamylated RF-amide peptide (QRFP) is a peptide hormone with a C-terminal RF-amide motif. QRFP selectively activates a class A G-protein-coupled receptor (GPCR) GPR103 to exert various physiological functions such as energy metabolism and appetite regulation. Here, we report the cryo-electron microscopy structure of the QRFP26-GPR103-Gq complex at 3.19 Å resolution. QRFP26 adopts an extended structure bearing no secondary structure, with its N-terminal and C-terminal sides recognized by extracellular and transmembrane domains of GPR103 respectively. This movement, reminiscent of class B1 GPCRs except for orientation and structure of the ligand, is critical for the high-affinity binding and receptor specificity of QRFP26. Mutagenesis experiments validate the functional importance of the binding mode of QRFP26 by GPR103. Structural comparisons with closely related receptors, including RY-amide peptide-recognizing GPCRs, revealed conserved and diversified peptide recognition mechanisms, providing profound insights into the biological significance of RF-amide peptides. Collectively, this study not only advances our understanding of GPCR-ligand interactions, but also paves the way for the development of novel therapeutics targeting metabolic and appetite disorders and emergency medical care.

A Study on the Robustness and Stability of Explainable Deep Learning in an Imbalanced Setting: The Exploration of the Conformational Space of G Protein-Coupled Receptors

G-protein coupled receptors (GPCRs) are transmembrane proteins that transmit signals from the extracellular environment to the inside of the cells. Their ability to adopt various conformational states, which influence their function, makes them crucial in pharmacoproteomic studies. While many drugs target specific GPCR states to exert their effects-thereby regulating the protein's activity-unraveling the activation pathway remains challenging due to the multitude of intermediate transformations occurring throughout this process, and intrinsically influencing the dynamics of the receptors. In this context, computational modeling, particularly molecular dynamics (MD) simulations, may offer valuable insights into the dynamics and energetics of GPCR transformations, especially when combined with machine learning (ML) methods and techniques for achieving model interpretability for knowledge generation. The current study builds upon previous work in which the layer relevance propagation (LRP) technique was employed to interpret the predictions in a multi-class classification problem concerning the conformational states of the β2-adrenergic (β2AR) receptor from MD simulations. Here, we address the challenges posed by class imbalance and extend previous analyses by evaluating the robustness and stability of deep learning (DL)-based predictions under different imbalance mitigation techniques. By meticulously evaluating explainability and imbalance strategies, we aim to produce reliable and robust insights.

Mechanistic insights into G-protein coupling with an agonist-bound G-protein-coupled receptor

G-protein-coupled receptors (GPCRs) activate heterotrimeric G proteins by promoting guanine nucleotide exchange. Here, we investigate the coupling of G proteins with GPCRs and describe the events that ultimately lead to the ejection of GDP from its binding pocket in the Gα subunit, the rate-limiting step during G-protein activation. Using molecular dynamics simulations, we investigate the temporal progression of structural rearrangements of GDP-bound Gs protein (Gs·GDP; hereafter GsGDP) upon coupling to the β2-adrenergic receptor (β2AR) in atomic detail. The binding of GsGDP to the β2AR is followed by long-range allosteric effects that significantly reduce the energy needed for GDP release: the opening of α1-αF helices, the displacement of the αG helix and the opening of the α-helical domain. Signal propagation to the Gs occurs through an extended receptor interface, including a lysine-rich motif at the intracellular end of a kinked transmembrane helix 6, which was confirmed by site-directed mutagenesis and functional assays. From this β2AR-GsGDP intermediate, Gs undergoes an in-plane rotation along the receptor axis to approach the β2AR-Gsempty state. The simulations shed light on how the structural elements at the receptor-G-protein interface may interact to transmit the signal over 30 Å to the nucleotide-binding site. Our analysis extends the current limited view of nucleotide-free snapshots to include additional states and structural features responsible for signaling and G-protein coupling specificity.

Trace amine associated receptor 1: predicted effects of single nucleotide variants on structure-function in geographically diverse populations

Trace Amine Associated Receptor 1 (TAAR1) is a novel pharmaceutical target under investigation for the treatment of several neuropsychiatric conditions. TAAR1 single nucleotide variants (SNV) have been found in patients with schizophrenia and metabolic disorders. However, the frequency of variants in geographically diverse populations and the functional effects of such variants are unknown. In this study, we aimed to characterise the distribution of TAAR1 SNVs in five different WHO regions using the Database of Genotypes and Phenotypes (dbGaP) and conducted a critical computational analysis using available TAAR1 structural data to identify SNVs affecting ligand binding and/or functional regions. Our analysis shows 19 orthosteric, 9 signalling and 16 micro-switch SNVs hypothesised to critically influence the agonist induced TAAR1 activation. These SNVs may non-proportionally influence populations from discrete regions and differentially influence the activity of TAAR1-targeting therapeutics in genetically and geographically diverse populations. Notably, our dataset presented with orthosteric SNVs D1033.32N (found only in the South-East Asian Region and Western Pacific Region) and T1945.42A (found only in South-East Asian Region), and 2 signalling SNVs (V1253.54A/T2526.36A, found in African Region and commonly, respectively), all of which have previously demonstrated to influence ligand induced functions of TAAR1. Furthermore, bioinformatics analysis using SIFT4G, MutationTaster 2, PROVEAN and MutationAssessor predicted all 16 micro-switch SNVs are damaging and may further influence the agonist activation of TAAR1, thereby possibly impacting upon clinical outcomes. Understanding the genetic basis of TAAR1 function and the impact of common mutations within clinical populations is important for the safe and effective utilisation of novel and existing pharmacotherapies.

Benchmarking AlphaFold-Generated Structures of Chemokine-Chemokine Receptor Complexes

AlphaFold and AlphaFold-Multimer have become two essential tools for the modeling of unknown structures of proteins and protein complexes. In this work, we extensively benchmarked the quality of chemokine-chemokine receptor structures generated by AlphaFold-Multimer against experimentally determined structures. Our analysis considered both the global quality of the model, as well as key structural features for chemokine recognition. To study the effects of template and multiple sequence alignment parameters on the results, a new prediction pipeline called LIT-AlphaFold (https://github.com/LIT-CCM-lab/LIT-AlphaFold) was developed, allowing extensive input customization. AlphaFold-Multimer correctly predicted differences in chemokine binding orientation and accurately reproduced the unique binding orientation of the CXCL12-ACKR3 complex. Further, the predictions of the full receptor N-terminus provided insights into a putative chemokine recognition site 0.5. The accuracy of chemokine N-terminus binding mode prediction varied between complexes, but the confidence score permitted the distinguishing of residues that were very likely well positioned. Finally, we generated a high-confidence model of the unsolved CXCL12-CXCR4 complex, which agreed with experimental mutagenesis and cross-linking data.

Arrestins: A Small Family of Multi-Functional Proteins

The first member of the arrestin family, visual arrestin-1, was discovered in the late 1970s. Later, the other three mammalian subtypes were identified and cloned. The first described function was regulation of G protein-coupled receptor (GPCR) signaling: arrestins bind active phosphorylated GPCRs, blocking their coupling to G proteins. It was later discovered that receptor-bound and free arrestins interact with numerous proteins, regulating GPCR trafficking and various signaling pathways, including those that determine cell fate. Arrestins have no enzymatic activity; they function by organizing multi-protein complexes and localizing their interaction partners to particular cellular compartments. Today we understand the molecular mechanism of arrestin interactions with GPCRs better than the mechanisms underlying other functions. However, even limited knowledge enabled the construction of signaling-biased arrestin mutants and extraction of biologically active monofunctional peptides from these multifunctional proteins. Manipulation of cellular signaling with arrestin-based tools has research and likely therapeutic potential: re-engineered proteins and their parts can produce effects that conventional small-molecule drugs cannot.

Engineering G protein-coupled receptors for stabilization

G protein-coupled receptors (GPCRs) are one of the most important families of targets for drug discovery. One of the limiting steps in the study of GPCRs has been their stability, with significant and time-consuming protein engineering often used to stabilize GPCRs for structural characterization and drug screening. Unfortunately, computational methods developed using globular soluble proteins have translated poorly to the rational engineering of GPCRs. To fill this gap, we propose GPCR-tm, a novel and personalized structurally driven web-based machine learning tool to study the impacts of mutations on GPCR stability. We show that GPCR-tm performs as well as or better than alternative methods, and that it can accurately rank the stability changes of a wide range of mutations occurring in various types of class A GPCRs. GPCR-tm achieved Pearson's correlation coefficients of 0.74 and 0.46 on 10-fold cross-validation and blind test sets, respectively. We observed that the (structural) graph-based signatures were the most important set of features for predicting destabilizing mutations, which points out that these signatures properly describe the changes in the environment where the mutations occur. More specifically, GPCR-tm was able to accurately rank mutations based on their effect on protein stability, guiding their rational stabilization. GPCR-tm is available through a user-friendly web server at https://biosig.lab.uq.edu.au/gpcr_tm/.

Molecular Determinant Underlying Selective Coupling of Primary G-Protein by Class A GPCRs

G-protein-coupled receptors (GPCRs) transmit downstream signals predominantly via G-protein pathways. However, the conformational basis of selective coupling of primary G-protein remains elusive. Histamine receptors H2R and H3R couple with Gs- or Gi-proteins respectively. Here, three cryo-EM structures of H2R-Gs and H3R-Gi complexes are presented at a global resolution of 2.6-2.7 Å. These structures reveal the unique binding pose for endogenous histamine in H3R, wherein the amino group interacts with E2065.46 of H3R instead of the conserved D1143.32 of other aminergic receptors. Furthermore, comparative analysis of the H2R-Gs and H3R-Gi complexes reveals that the structural geometry of TM5/TM6 determines the primary G-protein selectivity in histamine receptors. Machine learning (ML)-based structuromic profiling and functional analysis of class A GPCR-G-protein complexes illustrate that TM5 length, TM5 tilt, and TM6 outward movement are key determinants of the Gs and Gi/o selectivity among the whole Class A family. Collectively, the findings uncover the common structural geometry within class A GPCRs that determines the primary Gs- and Gi/o-coupling selectivity.

Bayesian network models identify cooperative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces, it is challenging to determine which amino acid pair interactions are cooperative. In this work, we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with Gα subunits. Our results reveal six GPCR:Gα contacts that are common to the different Gα subtypes and show strong cooperativity in the formation of interface. Both the C terminus helix5 and the core of the G protein are codependent entities and play an important role in GPCR coupling. We show that a promiscuous GPCR coupling to different Gα subtypes, makes all the GPCR:Gα contacts that are specific to each Gα subtype (Gαs, Gαi, and Gαq). This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Structural pharmacology and therapeutic potential of 5-methoxytryptamines

Psychedelic substances such as lysergic acid diethylamide (LSD) and psilocybin show potential for the treatment of various neuropsychiatric disorders1-3. These compounds are thought to mediate their hallucinogenic and therapeutic effects through the serotonin (5-hydroxytryptamine (5-HT)) receptor 5-HT2A (ref. 4). However, 5-HT1A also plays a part in the behavioural effects of tryptamine hallucinogens5, particularly 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT), a psychedelic found in the toxin of Colorado River toads6. Although 5-HT1A is a validated therapeutic target7,8, little is known about how psychedelics engage 5-HT1A and which effects are mediated by this receptor. Here we map the molecular underpinnings of 5-MeO-DMT pharmacology through five cryogenic electron microscopy (cryo-EM) structures of 5-HT1A, systematic medicinal chemistry, receptor mutagenesis and mouse behaviour. Structure-activity relationship analyses of 5-methoxytryptamines at both 5-HT1A and 5-HT2A enable the characterization of molecular determinants of 5-HT1A signalling potency, efficacy and selectivity. Moreover, we contrast the structural interactions and in vitro pharmacology of 5-MeO-DMT and analogues to the pan-serotonergic agonist LSD and clinically used 5-HT1A agonists. We show that a 5-HT1A-selective 5-MeO-DMT analogue is devoid of hallucinogenic-like effects while retaining anxiolytic-like and antidepressant-like activity in socially defeated animals. Our studies uncover molecular aspects of 5-HT1A-targeted psychedelics and therapeutics, which may facilitate the future development of new medications for neuropsychiatric disorders.

Antibodies expand the scope of angiotensin receptor pharmacology

G-protein-coupled receptors (GPCRs) are key regulators of human physiology and are the targets of many small-molecule research compounds and therapeutic drugs. While most of these ligands bind to their target GPCR with high affinity, selectivity is often limited at the receptor, tissue and cellular levels. Antibodies have the potential to address these limitations but their properties as GPCR ligands remain poorly characterized. Here, using protein engineering, pharmacological assays and structural studies, we develop maternally selective heavy-chain-only antibody ('nanobody') antagonists against the angiotensin II type I receptor and uncover the unusual molecular basis of their receptor antagonism. We further show that our nanobodies can simultaneously bind to angiotensin II type I receptor with specific small-molecule antagonists and demonstrate that ligand selectivity can be readily tuned. Our work illustrates that antibody fragments can exhibit rich and evolvable pharmacology, attesting to their potential as next-generation GPCR modulators.

Structural and dynamic insights into the activation of the μ-opioid receptor by an allosteric modulator

G-protein-coupled receptors (GPCRs) play pivotal roles in various physiological processes. These receptors are activated to different extents by diverse orthosteric ligands and allosteric modulators. However, the mechanisms underlying these variations in signaling activity by allosteric modulators remain largely elusive. Here, we determine the three-dimensional structure of the μ-opioid receptor (MOR), a class A GPCR, in complex with the Gi protein and an allosteric modulator, BMS-986122, using cryogenic electron microscopy. Our results reveal that BMS-986122 binding induces changes in the map densities corresponding to R1673.50 and Y2545.58, key residues in the structural motifs conserved among class A GPCRs. Nuclear magnetic resonance analyses of MOR in the absence of the Gi protein reveal that BMS-986122 binding enhances the formation of the interaction between R1673.50 and Y2545.58, thus stabilizing the fully-activated conformation, where the intracellular half of TM6 is outward-shifted to allow for interaction with the Gi protein. These findings illuminate that allosteric modulators like BMS-986122 can potentiate receptor activation through alterations in the conformational dynamics in the core region of GPCRs. Together, our results demonstrate the regulatory mechanisms of GPCRs, providing insights into the rational development of therapeutics targeting GPCRs.

Membrane mimetic-dependence of GPCR energy landscapes

We leveraged variable-temperature 19F-NMR spectroscopy to compare the conformational equilibria of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), across a range of temperatures ranging from lower temperatures typically employed in 19F-NMR experiments to physiological temperature. A2AAR complexes with partial agonists and full agonists showed large increases in the population of a fully active conformation with increasing temperature. NMR data measured at physiological temperature were more in line with functional data. This was pronounced for complexes with partial agonists, where the population of active A2AAR was nearly undetectable at lower temperature but became evident at physiological temperature. Temperature-dependent behavior of complexes with either full or partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. Cellular signaling experiments correlated with the temperature-dependent conformational equilibria of A2AAR in lipid nanodiscs but not in some detergents, underscoring the importance of the membrane environment in studies of GPCR function.

Time-resolved cryo-EM of G-protein activation by a GPCR

G-protein-coupled receptors (GPCRs) activate heterotrimeric G proteins by stimulating guanine nucleotide exchange in the Gα subunit1. To visualize this mechanism, we developed a time-resolved cryo-EM approach that examines the progression of ensembles of pre-steady-state intermediates of a GPCR-G-protein complex. By monitoring the transitions of the stimulatory Gs protein in complex with the β2-adrenergic receptor at short sequential time points after GTP addition, we identified the conformational trajectory underlying G-protein activation and functional dissociation from the receptor. Twenty structures generated from sequential overlapping particle subsets along this trajectory, compared to control structures, provide a high-resolution description of the order of main events driving G-protein activation in response to GTP binding. Structural changes propagate from the nucleotide-binding pocket and extend through the GTPase domain, enacting alterations to Gα switch regions and the α5 helix that weaken the G-protein-receptor interface. Molecular dynamics simulations with late structures in the cryo-EM trajectory support that enhanced ordering of GTP on closure of the α-helical domain against the nucleotide-bound Ras-homology domain correlates with α5 helix destabilization and eventual dissociation of the G protein from the GPCR. These findings also highlight the potential of time-resolved cryo-EM as a tool for mechanistic dissection of GPCR signalling events.

Flexible scaffold-based cheminformatics approach for polypharmacological drug design

Effective treatments for complex central nervous system (CNS) disorders require drugs with polypharmacology and multifunctionality, yet designing such drugs remains a challenge. Here, we present a flexible scaffold-based cheminformatics approach (FSCA) for the rational design of polypharmacological drugs. FSCA involves fitting a flexible scaffold to different receptors using different binding poses, as exemplified by IHCH-7179, which adopted a "bending-down" binding pose at 5-HT2AR to act as an antagonist and a "stretching-up" binding pose at 5-HT1AR to function as an agonist. IHCH-7179 demonstrated promising results in alleviating cognitive deficits and psychoactive symptoms in mice by blocking 5-HT2AR for psychoactive symptoms and activating 5-HT1AR to alleviate cognitive deficits. By analyzing aminergic receptor structures, we identified two featured motifs, the "agonist filter" and "conformation shaper," which determine ligand binding pose and predict activity at aminergic receptors. With these motifs, FSCA can be applied to the design of polypharmacological ligands at other receptors.

Calcium-Driven In Silico Inactivation of a Human Olfactory Receptor

Conformational changes as well as molecular determinants related to the activation and inactivation of olfactory receptors are still poorly understood due to the intrinsic difficulties in the structural determination of this GPCR family. Here, we perform, for the first time, the in silico inactivation of human olfactory receptor OR51E2, highlighting the possible role of calcium in this receptor state transition. Using molecular dynamics simulations, we show that a divalent ion in the ion binding site, coordinated by two acidic residues at positions 2.50 and 3.39 conserved across most ORs, stabilizes the receptor in its inactive state. In contrast, protonation of the same two acidic residues is not sufficient to drive inactivation within the microsecond timescale of our simulations. Our findings suggest a novel molecular mechanism for OR inactivation, potentially guiding experimental validation and offering insights into the possible broader role of divalent ions in GPCR signaling.

Comparison of the function of two novel human dopamine D2 receptor variants identifies a likely mechanism for their pathogenicity

Two recently discovered DRD2 mutations, c.634A > T, p.Ile212Phe and c.1121T > G, p.Met374Arg, cause hyperkinetic movement disorders that have overlapping features but apparently differ in severity. The two known carriers of the Met374Arg variant had early childhood disease onset and more severe motor, cognitive, and neuropsychiatric deficits than any known carriers of the Ile212Phe variant, whose symptoms were first apparent in adolescence. Here, we evaluated if differences in the function of the two variants in cultured cells could explain differing pathogenicity. Both variants were expressed less abundantly than the wild type receptor and exhibited loss of agonist-induced arrestin binding, but differences in expression and arrestin binding between the variants were minor. Basal and agonist-induced activation of heterotrimeric Gi/o/z proteins, however, showed clear differences; agonists were generally more potent at Met374Arg than at the Ile212Phe or wild type variants. Furthermore, all Gα subtypes tested were constitutively activated more by Met374Arg than by Ile212Phe. Met374Arg produced greater constitutive inhibition of cyclic AMP accumulation than Ile212Phe or the wild type D2 receptor. Met374Arg and Ile212Phe were more sensitive to thermal inactivation than the wild type D2 receptor, as reported for other constitutively active receptors, but Ile212Phe was affected more than Met374Arg. Additional pharmacological characterization suggested that the mutations differentially affect the shape of the agonist binding pocket and the potency of dopamine, norepinephrine, and tyramine. Molecular dynamics simulations provided a structural rationale for enhanced constitutive activation and agonist potency. Enhanced constitutive and agonist-induced G protein-mediated signaling likely contributes to the pathogenicity of these novel variants.