Irina G. Tikhonova

Discovery of novel Agonists and antagonists of the free fatty acid receptor 1 (FFAR1) using virtual screening

The G-protein-coupled receptor free fatty acid receptor 1 (FFAR1), previously named GPR40, is a possible novel target for the treatment of type 2 diabetes. In an attempt to identify new ligands for this receptor, we performed virtual screening (VS) based on two-dimensional (2D) similarity, three-dimensional (3D) pharmacophore searches, and docking studies by using the structure of known agonists and our model of the ligand binding site, which was validated by mutagenesis. VS of a database of 2.6 million compounds followed by extraction of structural neighbors of functionally confirmed hits resulted in identification of 15 compounds active at FFAR1 either as full agonists, partial agonists, or pure antagonists. Site-directed mutagenesis and docking studies revealed different patterns of ligand-receptor interactions and provided important information on the role of specific amino acids in binding and activation of FFAR1.

Identification of residues important for agonist recognition and activation in GPR40

GPR40 was formerly an orphan G protein-coupled receptor whose endogenous ligands have recently been identified as free fatty acids (FFAs). The receptor, now named FFA receptor 1, has been implicated in the pathophysiology of type 2 diabetes and is a drug target because of its role in FFA-mediated enhancement of glucose-stimulated insulin release. Guided by molecular modeling, we investigated the molecular determinants contributing to binding of linoleic acid, a C18 polyunsaturated FFA, and GW9508, a synthetic small molecule agonist. Twelve residues within the putative GPR40-binding pocket including hydrophilic/positively charged, aromatic, and hydrophobic residues were identified and were subjected to site-directed mutagenesis. Our results suggest that linoleic acid and GW9508 are anchored on their carboxylate groups by Arg183, Asn244, and Arg258. Moreover, His86, Tyr91, and His137 may contribute to aromatic and/or hydrophobic interactions with GW9508 that are not present, or relatively weak, with linoleic acid. The anchor residues, as well as the residues Tyr12, Tyr91, His137, and Leu186, appear to be important for receptor activation also. Interestingly, His137 and particularly His86 may interact with GW9508 in a manner dependent on its protonation status. The greater number of putative interactions between GPR40 and GW9508 compared with linoleic acid may explain the higher potency of GW9508.

Selective Orthosteric Free Fatty Acid Receptor 2 (FFA2) Agonists IDENTIFICATION OF THE STRUCTURAL AND CHEMICAL REQUIREMENTS FOR SELECTIVE ACTIVATION OF FFA2 VERSUS FFA3

Free fatty acid receptor 2 (FFA2; GPR43) is a G protein-coupled seven-transmembrane receptor for short-chain fatty acids (SCFAs) that is implicated in inflammatory and metabolic disorders. The SCFA propionate has close to optimal ligand efficiency for FFA2 and can hence be considered as highly potent given its size. Propionate, however, does not discriminate between FFA2 and the closely related receptor FFA3 (GPR41). To identify FFA2-selective ligands and understand the molecular basis for FFA2 selectivity, a targeted library of small carboxylic acids was examined using holistic, label-free dynamic mass redistribution technology for primary screening and the receptor-proximal G protein [35S]guanosine 5′-(3-O-thio)triphosphate activation, inositol phosphate, and cAMP accumulation assays for hit confirmation. Structure-activity relationship analysis allowed formulation of a general rule to predict selectivity for small carboxylic acids at the orthosteric binding site where ligands with substituted sp3-hybridized α-carbons preferentially activate FFA3, whereas ligands with sp2- or sp-hybridized α-carbons prefer FFA2. The orthosteric binding mode was verified by site-directed mutagenesis: replacement of orthosteric site arginine residues by alanine in FFA2 prevented ligand binding, and molecular modeling predicted the detailed mode of binding. Based on this, selective mutation of three residues to their non-conserved counterparts in FFA3 was sufficient to transfer FFA3 selectivity to FFA2. Thus, selective activation of FFA2 via the orthosteric site is achievable with rather small ligands, a finding with significant implications for the rational design of therapeutic compounds selectively targeting the SCFA receptors.

Extracellular Ionic Locks Determine Variation in Constitutive Activity and Ligand Potency between Species Orthologs of the Free Fatty Acid Receptors FFA2 and FFA3

Background: Differences in ligand potency and selectivity are observed between human and mouse FFA2 and FFA3 orthologs. Results: Potency differences result from differential constitutive activity between species. Conclusion: An “ionic lock” between extracellular loop 2 and the ligand binding pocket regulates constitutive activity. Significance: Understanding species differences in FFA2 and FFA3 function is critical to future studies with these receptors. Free fatty acid receptors 2 and 3 (FFA2 and FFA3) are G protein-coupled receptors for short chain free fatty acids (SCFAs). They respond to the same set of endogenous ligands but with distinct rank-order of potency such that acetate (C2) has been described as FFA2-selective, whereas propionate (C3) is non-selective. Although C2 was confirmed to be selective for human FFA2 over FFA3, this ligand was not selective between the mouse orthologs. Moreover, although C3 was indeed not selective between the human orthologs, it displayed clear selectivity for mouse FFA3 over mouse FFA2. This altered selectivity to C2 and C3 resulted from broad differences in SCFAs potency at the mouse orthologs. In studies to define the molecular basis for these observations, marked variation in ligand-independent constitutive activity was identified using a [35S]GTPγS assay. The orthologs with higher potency for the SCFAs, human FFA2 and mouse FFA3, displayed high constitutive activity in this assay, whereas the orthologs with lower potency for the agonist ligands, mouse FFA2 and human FFA3, did not. Sequence alignments of the second extracellular loop identified single negatively charged residues in FFA2 and FFA3 not conserved between species and predicted to form ionic lock interactions with arginine residues within the FFA2 or FFA3 agonist binding pocket to regulate constitutive activity and SCFA potency. Reciprocal mutation of these residues between species orthologs resulted in the induction (or repression) of constitutive activity and in most cases also yielded corresponding changes in SCFA potency.

Two Arginine-Glutamate Ionic Locks Near the Extracellular Surface of FFAR1 Gate Receptor Activation

Activation of a number of class A G protein-coupled receptors (GPCRs) is thought to involve two molecular switches, a rotamer toggle switch within the transmembrane domain and an ionic lock at the cytoplasmic surface of the receptor; however, the mechanism by which agonist binding changes these molecular interactions is not understood. Importantly, 80% of GPCRs including free fatty acid receptor 1 (FFAR1) lack the complement of amino acid residues implicated in either or both of these two switches; the mechanism of activation of these GPCRs is therefore less clear. By homology modeling, we identified two Glu residues (Glu-145 and Glu-172) in the second extracellular loop of FFAR1 that form putative interactions individually with two transmembrane Arg residues (Arg-183(5.39) and Arg-258(7.35)) to create two ionic locks. Molecular dynamics simulations showed that binding of agonists to FFAR1 leads to breakage of these Glu-Arg interactions. In mutagenesis experiments, breakage of these two putative interactions by substituting Ala for Glu-145 and Glu-172 caused constitutive receptor activation. Our results therefore reveal a molecular switch for receptor activation present on the extracellular surface of FFAR1 that is broken by agonist binding. Similar ionic locks between the transmembrane domains and the extracellular loops may constitute a mechanism common to other class A GPCRs also.

Extracellular Loop 2 Of The Free Fatty Acid Receptor 2 Mediates Allosterism Of A Phenylacetamide Ago-Allosteric Modulator

Allosteric agonists are powerful tools for exploring the pharmacology of closely related G protein-coupled receptors that have nonselective endogenous ligands, such as the short chain fatty acids at free fatty acid receptors 2 and 3 (FFA2/GPR43 and FFA3/GPR41, respectively). We explored the molecular mechanisms mediating the activity of 4-chloro-α-(1-methylethyl)-N-2-thiazolylbenzeneacetamide (4-CMTB), a recently described phenylacetamide allosteric agonist and allosteric modulator of endogenous ligand function at human FFA2, by combining our previous knowledge of the orthosteric binding site with targeted examination of 4-CMTB structure-activity relationships and mutagenesis and chimeric receptor generation. Here we show that 4-CMTB is a selective agonist for FFA2 that binds to a site distinct from the orthosteric site of the receptor. Ligand structure-activity relationship studies indicated that the N-thiazolyl amide is likely to provide hydrogen bond donor/acceptor interactions with the receptor. Substitution at Leu173 or the exchange of the entire extracellular loop 2 of FFA2 with that of FFA3 was sufficient to reduce or ablate, respectively, allosteric communication between the endogenous and allosteric agonists. Thus, we conclude that extracellular loop 2 of human FFA2 is required for transduction of cooperative signaling between the orthosteric and an as-yet-undefined allosteric binding site of the FFA2 receptor that is occupied by 4-CMTB.

Defining the molecular basis for the first potent and selective orthosteric agonists of the FFA2 free fatty acid receptor.

Background: Understanding the function of FFA2 has been slowed by a lack of selective orthosteric ligands. Results: Residues within FFA2 that dictate the recognition and function of potent and selective orthosteric agonists are described. Conclusion: Key aspects of ligand interaction with the orthosteric binding pocket of FFA2 are defined. Significance: This work will be invaluable in future drug development at the FFA2 receptor. FFA2 is a G protein-coupled receptor that responds to short chain fatty acids and has generated interest as a therapeutic target for metabolic and inflammatory conditions. However, definition of its functions has been slowed by a dearth of selective ligands that can distinguish it from the closely related FFA3. At present, the only selective ligands described for FFA2 suffer from poor potency, altered signaling due to allosteric modes of action, or a lack of function at non-human orthologs of the receptor. To address the need for novel selective ligands, we synthesized two compounds potentially having FFA2 activity and examined the molecular basis of their function. These compounds were confirmed to be potent and selective orthosteric FFA2 agonists. A combination of ligand structure-activity relationship, pharmacological analysis, homology modeling, species ortholog comparisons, and mutagenesis studies were then employed to define the molecular basis of selectivity and function of these ligands. From this, we identified key residues within both extracellular loop 2 and the transmembrane domain regions of FFA2 critical for ligand function. One of these ligands was active with reasonable potency at rodent orthologs of FFA2 and demonstrated the role of FFA2 in inhibition of lipolysis and glucagon-like peptide-1 secretion in murine-derived 3T3-L1 and STC-1 cell lines, respectively. Together, these findings describe the first potent and selective FFA2 orthosteric agonists and demonstrate key aspects of ligand interaction within the binding site of FFA2 that will be invaluable in future ligand development at this receptor.

Identification of determinants of glucose-dependent insulinotropic polypeptide receptor that interact with N-terminal biologically active region of the natural ligand.

Glucose-dependent insulinotropic polypeptide receptor (GIPR), a member of family B of the G-protein coupled receptors, is a potential therapeutic target for which discovery of nonpeptide ligands is highly desirable. Structure-activity relationship studies indicated that the N-terminal part of glucose-dependent insulinotropic polypeptide (GIP) is crucial for biological activity. Here, we aimed at identification of residues in the GIPR involved in functional interaction with N-terminal moiety of GIP. A homology model of the transmembrane core of GIPR was constructed, whereas a three-dimensional model of the complex formed between GIP and the N-terminal extracellular domain of GIPR was taken from the crystal structure. The latter complex was docked to the transmembrane domains of GIPR, allowing in silico identification of putative residues of the agonist binding/activation site. All mutants were expressed at the surface of human embryonic kidney 293 cells as indicated by flow cytometry and confocal microscopy analysis of fluorescent GIP binding. Mutation of residues Arg183, Arg190, Arg300, and Phe357 caused shifts of 76-, 71-, 42-, and 16-fold in the potency to induce cAMP formation, respectively. Further characterization of these mutants, including tests with alanine-substituted GIP analogs, were in agreement with interaction of Glu3 in GIP with Arg183 in GIPR. Furthermore, they strongly supported a binding mode of GIP to GIPR in which the N-terminal moiety of GIP was sited within transmembrane helices (TMH) 2, 3, 5, and 6 with biologically crucial Tyr1 interacting with Gln224 (TMH3), Arg300 (TMH5), and Phe357 (TMH6). These data represent an important step toward understanding activation of GIPR by GIP, which should facilitate the rational design of therapeutic agents.

Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs

When it is difficult to develop selective ligands within a family of related G‐protein‐coupled receptors (GPCRs), chemically engineered receptors activated solely by synthetic ligands (RASSLs) are useful alternatives for probing receptor function. In the present work, we explored whether a RASSL of the free fatty acid receptor 2 (FFA2) could be developed on the basis of pharmacological variation between species orthologs. For this, bovine FFA2 was characterized, revealing distinct ligand selectivity compared with human FFA2. Homology modeling and mutational analysis demonstrated a single mutation in human FFA2 of C4.57G resulted in a human FFA2 receptor with ligand selectivity similar to the bovine receptor. This was exploited to generate human FFA2‐RASSL by the addition of a second mutation at a known orthosteric ligand interaction site, H6.55Q. The resulting FFA2‐RASSL displayed a >100‐fold loss of activity to endogenous ligands, while responding to the distinct ligand sorbic acid with pEC50 values for inhibition of cAMP, 5.83 ± 0.11; Ca2+ mobilization, 4.63 ± 0.05; ERK phosphorylation, 5.61 ± 0.06; and dynamic mass redistribution, 5.35 ± 0.06. This FFA2‐RASSL will be useful in future studies on this receptor and demonstrates that exploitation of pharmacological variation between species orthologs is a powerful method to generate novel chemically engineered GPCRs.—Hudson, B. D., Christiansen, E., Tikhonova, I. G., Grundmann, M., Kostenis, E., Adams, D. R., Ulven, T., Milligan, G. Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs. FASEB J. 26, 4951–4965 (2012). www.fasebj.org

Probing the mid-gorge of cholinesterases with spacer-modified bivalent quinazolinimines leads to highly potent and selective butyrylcholinesterase inhibitors.

The spacer structure of homobivalent quinazolinimes acting as potent acetyl-(AChE)- and butyrylcholinesterase (BChE) inhibitors was chemically modified introducing tertiary amine and acyl-amide moieties, and the activities at both ChEs were evaluated. Molecular docking was applied to explain the data and probe the capacity of the mid-gorge site of both ChEs. The novel spacer structures considerably alter the biological profile of bivalent quinazolinimines with regard to both inhibitory activity and selectivity. Mutual interaction of binding to the various sites of the enzymes was further investigated by applying also different spacer lengths and ring sizes of the alicycle of the tricyclic quinazolinimines. In order to achieve selectivity toward BChE and to improve inhibitory activities, the spacer structure was optimized and identified a highly potent and selective BChE inhibitor.