Stanislav Engel

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.

A Virtual Screen for Diverse Ligands : Discovery of Selective G Protein-Coupled Receptor Antagonists

Virtual screening has become a major focus of bioactive small molecule lead identification, and reports of agonists and antagonists discovered via virtual methods are becoming more frequent. G protein-coupled receptors (GPCRs) are the one class of protein targets for which success with this approach has been limited. This is likely due to the paucity of detailed experimental information describing GPCR structure and the intrinsic function-associated structural flexibility of GPCRs which present major challenges in the application of receptor-based virtual screening. Here we describe an in silico methodology that diminishes the effects of structural uncertainty, allowing for more inclusive representation of a potential docking interaction with exogenous ligands. Using this approach, we screened one million compounds from a virtual database, and a diverse subgroup of 100 compounds was selected, leading to experimental identification of five structurally diverse antagonists of the thyrotropin-releasing hormone receptors (TRH-R1 and TRH-R2). The chirality of the most potent chemotype was demonstrated to be important in its binding affinity to TRH receptors; the most potent stereoisomer was noted to have a 13-fold selectivity for TRH-R1 over TRH-R2. A comprehensive mutational analysis of key amino acid residues that form the putative binding pocket of TRH receptors further verified the binding modality of these small molecule antagonists. The described virtual screening approach may prove applicable in the search for novel small molecule agonists and antagonists of other GPCRs.

Atomistic Insights into Rhodopsin Activation from a Dynamic Model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Low Affinity Analogs of Thyrotropin-releasing Hormone Are Super-agonists

We show that several analogs of thyrotropin-releasing hormone (TRH) are more efficacious agonists at TRH receptors R1 and R2 than TRH itself. The apparent efficacies of the analogs were inversely related to their potencies and were independent of the nature of the modifications in TRH structure. In studies in intact cells, we showed that the differences in apparent efficacies were not due to differences in G-protein coupling, receptor desensitization, or recycling. Moreover, the differences in efficacies persisted in experiments using accessory protein-free membranes. We conclude that the efficacy differences of TRH analogs originated from the enhanced ability of TRH-R complexed to the low affinity agonists to directly activate G-protein(s), and not by a modulation of the activity of accessory proteins, and propose possible mechanisms for this phenomenon.

Understanding the structural and functional differences between mouse thyrotropin-releasing hormone receptors 1 and 2.

Multiple computational methods have been employed in a comparative study of thyrotropin‐releasing hormone receptors 1 and 2 (TRH‐R1 and TRH‐R2) to explore the structural bases for the different functional properties of these G protein‐coupled receptors. Three‐dimensional models of both murine TRH receptors have been built and optimized by means of homology modeling based on the crystal structure of bovine rhodopsin, molecular dynamics simulations, and energy minimizations in a membrane‐aqueous environment. The comparison between the two models showed a correlation between the higher flexibility and higher basal activity of TRH‐R2 versus the lesser flexibility and lower basal activity of TRH‐R1 and supported the involvement of the highly conserved W6.48 in the signaling process. A correlation between the level of basal activity and conformational changes of TM5 was detected also. Comparison between models of the wild type receptors and their W6.48A mutants, which have reversed basal activities compared with their respective wild types, further supported these correlations. A flexible molecular docking procedure revealed that TRH establishes a direct interaction with W6.48 in TRH‐R2 but not in TRH‐R1. We designed and performed new mutagenesis experiments that strongly supported these observations. Proteins 2008; 71:783–794. Published 2007 Wiley‐Liss, Inc.

Dimerization of NKp46 Receptor Is Essential for NKp46-Mediated Lysis: Characterization of the Dimerization Site by Epitope Mapping

NKp46 is a primary activating receptor of NK cells that is involved in lysis of target cells by NK cells. Previous studies showed that the membrane-proximal domain of NKp46 (NKp46D2) retained the binding of NKp46 to its ligands and is involved in lysis. We studied NKp46D2 by using a peptide-based epitope mapping approach and identified an NKp46D2-derived linear epitope that inhibited NKp46-mediated lysis. The epitope, designated as pep4 (aa 136–155), interacted with NKp46, and lysis by NK cells was inhibited by the presence of pep4. Through modeling and mutagenesis, we showed that pep4 could be involved in NKp46 homodimerization. R145 and D147 contribute to the function of pep4, and R145Q mutation in recombinant NKp46 reduced its binding to target cells. At the cellular level, fluorescent resonance energy transfer analysis revealed that pep4 is indeed involved in dimerization of cell membrane-associated NKp46. We suggest that the NKp46-derived pep4 site is part of the dimerization surface of NKp46 and that NKp46 dimerization contributes to NKp46-mediated lysis by NK cells.

Purification of acetohydroxy acid synthase by separation in an aqueous two-phase system

Extraction in a polyethylene glycol (PEG)-phosphate aqueous two-phase system was considered as a primary step in purification of the acetohydroxy acid synthase III large catalytic subunit from an E. coli extract. Extraction optimization was achieved by varying the system parameters. Two systems with the following weight compositions were chosen for purification: PEG-2000 (16%)-phosphate (6%) and PEG-4000 (14%)-phosphate (5.5%)-KCl (8%), both at pH 7.0 and 1 mg total protein per 1 g system. Significant purification was achieved by a single extraction step with 70% recovery of the enzyme. After an additional ion-exchange chromatography step, pure enzyme was obtained in a 50% overall yield.

Determination of the dissociation constant of valine from acetohydroxy acid synthase by equilibrium partition in an aqueous two-phase system

An aqueous polyethylene glycol/salt two-phase system was used to estimate the dissociation constant, K(dis), of the Escherichia coli isoenzyme AHAS III regulatory subunit, ilvH protein, from the feedback inhibitor valine. The amounts of the bound and free radioactive valine in the system were determined. A Scatchard plot of the data revealed a 1:1 valine-protein binding ratio and K(dis) of 133+/-14 microM. The protein did not bind leucine, and the ilvH protein isolated from a valine resistant mutant showed no valine binding. This method is very simple, rapid and requires only a small amounts of protein compared to the presently used equilibrium dialysis method.

Surface Dynamics in Allosteric Regulation of Protein-Protein Interactions: Modulation of Calmodulin Functions by Ca2+

Knowledge of the structural basis of protein-protein interactions (PPI) is of fundamental importance for understanding the organization and functioning of biological networks and advancing the design of therapeutics which target PPI. Allosteric modulators play an important role in regulating such interactions by binding at site(s) orthogonal to the complex interface and altering the proteins propensity for complex formation. In this work, we apply an approach recently developed by us for analyzing protein surfaces based on steered molecular dynamics simulation (SMD) to the study of the dynamic properties of functionally distinct conformations of a model protein, calmodulin (CaM), whose ability to interact with target proteins is regulated by the presence of the allosteric modulator Ca2+. Calmodulin is a regulatory protein that acts as an intracellular Ca2+ sensor to control a wide variety of cellular processes. We demonstrate that SMD analysis is capable of pinpointing CaM surfaces implicated in the recognition of both the allosteric modulator Ca2+ and target proteins. Our analysis of changes in the dynamic properties of the CaM backbone elicited by Ca2+ binding yielded new insights into the molecular mechanism of allosteric regulation of CaM-target interactions.