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.

Structural basis for the ligand recognition and signaling of free fatty acid receptors.

Free fatty acid receptors 1 to 4 (FFA1 to FFA4) are class A G protein-coupled receptors (GPCRs). FFA1 to FFA3 share substantial sequence similarity, whereas FFA4 is unrelated. However, FFA1 and FFA4 are activated by long-chain fatty acids, while FFA2 and FFA3 respond to short-chain fatty acids generated by intestinal microbiota. FFA1, FFA2, and FFA4 are potential drug targets for metabolic and inflammatory conditions. Here, we determined the active structures of FFA1 and FFA4 bound to docosahexaenoic acid, FFA4 bound to the synthetic agonist TUG-891, and butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands to their respective GPCRs. Our findings unveiled distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research on this group of GPCRs.

ß-Arrestin-dependent and -independent endosomal G protein activation by the vasopressin type 2 receptor.

The vasopressin type 2 receptor (V2R) is an essential G protein-coupled receptor (GPCR) in renal regulation of water homeostasis. Upon stimulation, the V2R activates Gαs and Gαq/11, which is followed by robust recruitment of ß-arrestins and receptor internalization into endosomes. Unlike canonical GPCR signaling, the ß-arrestin association with the V2R does not terminate Gαs activation, and thus, Gαs-mediated signaling is sustained while the receptor is internalized. Here, we demonstrate that this V2R ability to co-interact with G protein/ß-arrestin and promote endosomal G protein signaling is not restricted to Gαs, but also involves Gαq/11. Furthermore, our data imply that ß-arrestins potentiate Gαs/Gαq/11 activation at endosomes rather than terminating their signaling. Surprisingly, we found that the V2R internalizes and promote endosomal G protein activation independent of ß-arrestins to a minor degree. These new observations challenge the current model of endosomal GPCR signaling and suggest that this event can occur in both ß-arrestin-dependent and -independent manners.

Structural basis for the ligand recognition and signaling of free fatty acid receptors.

Free fatty acid receptors 1-4 (FFA1-4) are class A G protein-coupled receptors (GPCRs). FFA1-3 share substantial sequence similarity whereas FFA4 is unrelated. Despite this FFA1 and FFA4 are activated by the same range of long chain fatty acids (LCFAs) whilst FFA2 and FFA3 are instead activated by short chain fatty acids (SCFAs) generated by the intestinal microbiota. Each of FFA1, 2 and 4 are promising targets for novel drug development in metabolic and inflammatory conditions. To gain insights into the basis of ligand interactions with, and molecular mechanisms underlying activation of, FFAs by LCFAs and SCFAs, we determined the active structures of FFA1 and FFA4 bound to the polyunsaturated LCFA docosahexaenoic acid (DHA), FFA4 bound to the synthetic agonist TUG-891, as well as SCFA butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands with varying chain lengths to their respective GPCRs. Our findings unveil distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research to understand allosteric modulation and biased signaling of this group of GPCRs.

ß-arrestin-dependent and -independent endosomal G protein activation by the vasopressin type 2 receptor.

The vasopressin type 2 receptor (V2R) is an essential GPCR in renal regulation of water homeostasis. Upon stimulation, the V2R activates Gαs and Gαq/11, which is followed by robust recruitment of ß-arrestins and receptor internalization into endosomes. Unlike canonical GPCR signaling, the ß-arrestin association with the V2R does not terminate Gαs activation, and thus, Gαs-mediated signaling is sustained while the receptor is internalized. Here, we demonstrate that this V2R ability to co-interact with G protein/ß-arrestin and promote endosomal G protein signaling is not restricted to Gαs, but also involves Gαq/11. Furthermore, our data implies that ß-arrestins potentiate Gαs/Gαq/11 activation at endosomes rather than terminating their signaling. Surprisingly, we found that the V2R internalizes and promote endosomal G protein activation independent of ß-arrestins to a minor degree. These new observations challenge the current model of endosomal GPCR signaling and suggest that this event can occur in both ß-arrestin-dependent and -independent manners.

Characterization of the mechanism of bile salt hydrolase substrate specificity by experimental and computational analyses.

Bile salt hydrolases (BSHs) are currently being investigated as target enzymes for metabolic regulators in humans and as growth promoters in farm animals. Understanding structural features underlying substrate specificity is necessary for inhibitor design. Here, we used a multidisciplinary workflow including mass spectrometry, mutagenesis, molecular dynamic simulations, machine learning, and crystallography to demonstrate substrate specificity in Lactobacillus salivarius BSH, the most abundant enzyme in human and farm animal intestines. We show the preference of substrates with a taurine head and a dehydroxylated sterol ring for hydrolysis. A regression model that correlates the relative rates of hydrolysis of various substrates in various enzyme mutants with the residue-substrate interaction energies guided the identification of structural determinants of substrate binding and specificity. In addition, we found T208 from another BSH protomer regulating the hydrolysis. The designed workflow can be used for fast and comprehensive characterization of enzymes with a broad range of substrates.

Oxidation modulates LINGO2-induced inactivation of large conductance, Ca2+-activated potassium channels.

Ca2+ and voltage-activated K+ (BK) channels are ubiquitous ion channels that can be modulated by accessory proteins, including ß, γ, and LINGO1 BK subunits. In this study, we utilized a combination of site-directed mutagenesis, patch clamp electrophysiology, and molecular modeling to investigate if the biophysical properties of BK currents were affected by coexpression of LINGO2 and to examine how they are regulated by oxidation. We demonstrate that LINGO2 is a regulator of BK channels, since its coexpression with BK channels yields rapid inactivating currents, the activation of which is shifted ∼-30 mV compared to that of BKα currents. Furthermore, we show the oxidation of BK:LINGO2 currents (by exposure to epifluorescence illumination or chloramine-T) abolished inactivation. The effect of illumination depended on the presence of GFP, suggesting that it released free radicals which oxidized cysteine or methionine residues. In addition, the oxidation effects were resistant to treatment with the cysteine-specific reducing agent DTT, suggesting that methionine rather than cysteine residues may be involved. Our data with synthetic LINGO2 tail peptides further demonstrate that the rate of inactivation was slowed when residues M603 or M605 were oxidized, and practically abolished when both were oxidized. Taken together, these data demonstrate that both methionine residues in the LINGO2 tail mediate the effect of oxidation on BK:LINGO2 channels. Our molecular modeling suggests that methionine oxidation reduces the lipophilicity of the tail, thus preventing it from occluding the pore of the BK channel.

Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations.

Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.

Investigating the Structure-Activity Relationship of 1,2,4-Triazine G-Protein-Coupled Receptor 84 (GPR84) Antagonists.

G-protein-coupled receptor 84 (GPR84) is a proinflammatory orphan G-protein-coupled receptor implicated in several inflammatory and fibrotic diseases. Several agonist and antagonist ligands have been developed that target GPR84; however, a noncompetitive receptor blocker that was progressed to phase II clinical trials failed to demonstrate efficacy. New high-quality antagonists are required to investigate the pathophysiological role of GPR84 and to validate GPR84 as a therapeutic target. We previously reported the discovery of a novel triazine GPR84 competitive antagonist 1. Here, we describe an extensive structure-activity relationship (SAR) of antagonist 1 and also present in silico docking with supporting mutagenesis studies that reveals a potential binding pose for this type of orthosteric antagonist. Lead compound 42 is a potent GPR84 antagonist with a favorable pharmacokinetic (PK) profile suitable for further drug development.

Luminogenic HiBiT Peptide-Based NanoBRET Ligand Binding Assays for Melatonin Receptors.

The two human melatonin receptors MT1 and MT2, which belong to the G protein-coupled receptor (GPCR) family, are important drug targets with approved indications for circadian rhythm- and sleep-related disorders and major depression. Currently, most of the pharmacological studies were performed using [3H]melatonin and 2-[125I]iodomelatonin (2-[125I]-MLT) radioligands. Recently, NanoLuc-based bioluminescence resonance energy transfer (NanoBRET) monitoring competitive binding between fluorescent tracers and unmodified test compounds has emerged as a sensitive, nonradioactive alternative for quantifying GPCR ligand engagement on the surface of living cells in equilibrium and real time. However, developing such assays for the two melatonin receptors depends on the availability of fluorescent tracers, which has been challenging predominantly owing to their narrow ligand entry channel and small ligand binding pocket. Here, we generated a set of melatonergic fluorescent tracers and used NanoBRET to evaluate their engagement with MT1 and MT2 receptors that are genetically fused to an N-terminal luminogenic HiBiT-peptide. We identified several nonselective and subtype-selective tracers. Among the selective tracers, PBI-8238 exhibited high nanomolar affinity to MT1, and PBI-8192 exhibited low nanomolar affinity to MT2. The pharmacological profiles of both tracers were in good agreement with those obtained with the current standard 2-[125I]-MLT radioligand. Molecular docking and mutagenesis studies suggested the binding mode of PBI-8192 in MT2 and its selectivity over MT1. In conclusion, we describe the development of the first nonradioactive, real-time binding assays for melatonin receptors expressed at the cell surface of living cells that are likely to accelerate drug discovery for melatonin receptors.