Study of ligand-based virtual screening tools in computer-aided drug design

Virtual screening is a central technique in drug discovery today. Millions of molecules can be tested in silico with the aim to only select the most promising and test them experimentally. The topic of this thesis is ligand-based virtual screening tools which take existing active molecules as starting point for finding new drug candidates. One goal of this thesis was to build a model that gives the probability that two molecules are biologically similar as function of one or more chemical similarity scores. Another important goal was to evaluate how well different ligand-based virtual screening tools are able to distinguish active molecules from inactives. One more criterion set for the virtual screening tools was their applicability in scaffold-hopping, i.e. finding new active chemotypes. In the first part of the work, a link was defined between the abstract chemical similarity score given by a screening tool and the probability that the two molecules are biologically similar. These results help to decide objectively which virtual screening hits to test experimentally. The work also resulted in a new type of data fusion method when using two or more tools. In the second part, five ligand-based virtual screening tools were evaluated and their performance was found to be generally poor. Three reasons for this were proposed: false negatives in the benchmark sets, active molecules that do not share the binding mode, and activity cliffs. In the third part of the study, a novel visualization and quantification method is presented for evaluation of the scaffold-hopping ability of virtual screening tools.

PPARγ induces growth inhibition and apoptosis through upregulation of insulin-like growth factor-binding protein-3 in gastric cancer cells

Peroxisome proliferator activator receptor-gamma (PPARγ) is a ligand-activated transcriptional factor involved in the carcinogenesis of various cancers. Insulin-like growth factor-binding protein-3 (IGFBP-3) is a tumor suppressor gene that has anti-apoptotic activity. The purpose of this study was to investigate the anticancer mechanism of PPARγ with respect to IGFBP-3. PPARγ was overexpressed in SNU-668 gastric cancer cells using an adenovirus gene transfer system. The cells in which PPARγ was overexpressed exhibited growth inhibition, induction of apoptosis, and a significant increase in IGFBP-3 expression. We investigated the underlying molecular mechanisms of PPARγ in SNU-668 cells using an IGFBP-3 promoter/luciferase reporter system. Luciferase activity was increased up to 15-fold in PPARγ transfected cells, suggesting that PPARγ may directly interact with IGFBP-3 promoter to induce its expression. Deletion analysis of the IGFBP-3 promoter showed that luciferase activity was markedly reduced in cells without putative p53-binding sites (-Δ1755, -Δ1795). This suggests that the critical PPARγ-response region is located within the p53-binding region of the IGFBP-3 promoter. We further demonstrated an increase in PPARγ-induced luciferase activity even in cells treated with siRNA to silence p53 expression. Taken together, these data suggest that PPARγ exhibits its anticancer effect by increasing IGFBP-3 expression, and that IGFBP-3 is a significant tumor suppressor.

VAP-1 as drug target in inflammation : structural features determining ligand interactions

Vascular adhesion protein-1 (VAP-1), which belongs to the copper amine oxidases (CAOs), is a validated drug target in inflammatory diseases. Inhibition of VAP-1 blocks the leukocyte trafficking to sites of inflammation and alleviates inflammatory reactions. In this study, a novel set of potent pyridazinone inhibitors is presented together with their X-ray structure complexes with VAP-1. The crystal structure of serum VAP-1 (sVAP-1) revealed an imidazole binding site in the active site channel and, analogously, the pyridazinone inhibitors were designed to bind into the channel. This is the first time human VAP-1 has been crystallized with a reversible inhibitor and the structures reveal detailed information of the binding mode on the atomic level. Similarly to some earlier studied inhibitors of human VAP-1, the designed pyridazinone inhibitors bind rodent VAP-1 with a lower affinity than human VAP-1. Therefore, we made homology models of rodent VAP-1 and compared human and rodent enzymes to determine differences that might affect the inhibitor binding. The comparison of the crystal structures of the human VAP-1 and the mouse VAP-1 homology model revealed key differences important for the species specific binding properties. In general, the channel in mouse VAP-1 is more narrow and polar than the channel in human VAP-1, which is wider and more hydrophobic. The differences are located in the channel leading to the active site, as well as, in the entrance to the active site channel. The information obtained from these studies is of great importance for the development and design of drugs blocking the activity of human VAP-1, as rodents are often used for in vivo testing of candidate drugs. In order to gain more insight into the selective binding properties of the different CAOs in one species a comprehensive evolutionary study of mammalian CAOs was performed. We found that CAOs can be classified into sub-families according to the residues X1 and X2 of the Thr/Ser-X1-X2-Asn-Tyr-Asp active site motif. In the phylogenetic tree, CAOs group into diamine oxidase, retina specific amine oxidase and VAP-1/serum amine oxidase clades based on the residue in the position X2. We also found that VAP-1 and SAO can be further differentiated based on the residue in the position X1. This is the first large-scale comparison of CAO sequences, which explains some of the reasons for the unique substrate specificities within the CAO family.

Unconventional integrin-ligand interactions

Integrins are cell surface adhesion and signaling receptors. Cells use integrins to attach to the extracellular matrix and to other cells, as well as for sensing their environment. In addition to adhesion and migration, integrins have been shown to be important for many biological processes including apoptosis, cell proliferation, and differentiation into specific tissues. Many important next generation biological drugs inhibit integrin functions. Thus, research into interactions between integrins and their ligands under different physiological and pathological conditions is not only of academic interest, but is also important for the field of drug discovery. In this Ph.D. project, the functions of integrin-ligand interactions were studied under different physiologically interesting conditions including 1) human echovirus 1 binding to integrin α2β1, 2) integrin α2β1 binding to collagen under flow conditions, 3) integrin α2β1 binding to a ligand in the presence of the angiogenesis inhibitor histidine rich glycoprotein (HRG) and 4) integrin binding to posttranslationally citrullinated ligands. As a result of the project, we could show that for each condition the integrin-ligand interaction is somewhat unconventional. 1) Echovirus 1 binds only to non-activated conformations of integrin α2β1. 2) Surprisingly, the non-activated conformation is also the primary conformation of integrin α2β1 when it binds to collagen under flow conditions, like when platelets adhere to subendothelial collagen in vascular injuries. In addition, the pre-activation of integrin α2β1 does not increase adhesion under flow. 3) HRG binds to integrin α2β1 through a low-affinity interaction that inhibits integrin binding to collagen. This shows that low affinity interactions could be biologically relevant and possibly regulate angiogenesis. 4) The citrullination of collagen, a posttranslational modification reported to occur in rheumatoid arthritis, specifically inhibits the binding of integrin α10β1 and α11β1, but does not affect the binding of α1β1 ja α2β1. On the other hand, the citrullination of isoDGR in fibronectin and RGD in pro-TGF- β:n inhibit integrin binding completely. Citrullination seems to be an inflammation related process and integrin ligands become citrullinated frequently in vivo. This Ph.D. thesis suggests that unconventional interaction mechanisms between integrins and their ligands, such as posttranslational modifications, low affinity interactions, and non-activated integrin conformations, can have an important role in pathological processes. The study of these kinds of integrin-ligand interactions is important for understanding biological phenomena more deeply. The research might also be beneficial for the development of integrin based therapies for treating diseases.

The effect of lysobisphosphatidic acid (LBPA) on α-tocopherol transfer protein (α-TTP) binding to lipid membranes

The a-tocopherol transfer protein (a-TTP) is responsible for the retention of the atocopherol form of vitamin E in living organisms. The detailed ligand transfer mechanism by a-TTP is still yet to be fully elucidated. To date, studies show that a-TTP transfers a-tocopherol from late endosomes in liver cells to the plasma membrane where it is repackaged into very low density lipoprotein (VLDL) and released into the circulation. Late endosomes have been shown to contain a lipid known as lysobisphosphatidic acid (LBP A) that is unique to this cellular compartment. LBPA plays a role in intracellular trafficking and controlling membrane curvature. Taking these observations into account plus the fact that certain proteins are recruited to membranes based on membrane curvature, the specific aim of this project was to examine the effect of LBP A on a-TTP binding to lipid membranes. To achieve this objective, dual polarization interferometry (DPI) and a vesicle binding assay were employed. Whilst DPI allows protein binding affinity to be measured on a flat lipid surface, the vesicle binding assay determines protein binding affinity to lipid vesicles mimicking curved membranes. DPI analysis revealed that the amount of a-TTP bound to lipid membranes is higher when LBPA is present. Using the vesicle binding assay, a similar result was seen where a greater amount of protein is bound to large unilamellar vesicles (LUV s) containing LBP A. However, the effect of LBP A was attenuated when small unilamellar vesicles (SUVs) were replaced with LUVs. The outcome of this project suggests that aTTP binding to membranes is influenced by membrane curvature, which in turn is induced by the presence of LBP A.

Directed evolution of human dihydrofolate reductase: towards a better understanding of binding at the active site

La dihydrofolate réductase humaine (DHFRh) est une enzyme essentielle à la prolifération cellulaire, ce qui en fait une cible de choix pour le traitement de différents cancers. À cet effet, plusieurs inhibiteurs spécifiques de la DHFRh, les antifolates, ont été mis au point : le méthotrexate (MTX) et le pemetrexed (PMTX) en sont de bons exemples. Malgré l’efficacité clinique certaine de ces antifolates, le développement de nouveaux traitements s’avère nécessaire afin de réduire les effets secondaires liés à leur utilisation. Enfin, dans l’optique d’orienter la synthèse de nouveaux composés inhibiteurs des DHFRh, une meilleure connaissance des interactions entre les antifolates et leur enzyme cible est primordiale. À l’aide de l’évolution dirigée, il a été possible d’identifier des mutants de la DHFRh pour lesquels l’affinité envers des antifolates cliniquement actifs se voyait modifiée. La mutagenèse dite ¬¬de saturation a été utilisée afin de générer des banques de mutants présentant une diversité génétique au niveau des résidus du site actif de l’enzyme d’intérêt. De plus, une nouvelle méthode de criblage a été mise au point, laquelle s’est avérée efficace pour départager les mutations ayant entrainé une résistance aux antifolates et/ou un maintient de l’activité enzymatique envers son substrat natif, soient les phénotypes d’activité. La méthode de criblage consiste dans un premier temps en une sélection bactérienne à haut débit, puis dans un second temps en un criblage sur plaques permettant d’identifier les meilleurs candidats. Plusieurs mutants actifs de la DHFRh, résistants aux antifolates, ont ainsi pu être identifiés et caractérisés lors d’études de cinétique enzymatique (kcat et IC50). Sur la base de ces résultats cinétiques, de la modélisation moléculaire et des données structurales de la littérature, une étude structure-activité a été effectuée. En regardant quelles mutations ont les effets les plus significatif sur la liaison, nous avons commencé à construire un carte moléculaire des contacts impliqués dans la liaison des ligands. Enfin, des connaissances supplémentaires sur les propriétés spécifiques de liaison ont put être acquises en variant l’inhibiteur testé, permettant ainsi une meilleure compréhension du phénomène de discrimination du ligand.

Combinatorial peptide ligand libraries and plant proteomics: a winning strategy at a price

The mechanism of action and properties of a solid-phase ligand library made of hexapeptides (combinatorial peptide ligand libraries or CPLL), for capturing the "hidden proteome", i.e. the low- and very low-abundance proteins constituting the vast majority of species in any proteome, as applied to plant tissues, are reviewed here. Plant tissues are notoriously recalcitrant to protein extraction and to proteome analysis. Firstly, rigid plant cell walls need to be mechanically disrupted to release the cell content and, in addition to their poor protein yield, plant tissues are rich in proteases and oxidative enzymes, contain phenolic compounds, starches, oils, pigments and secondary metabolites that massively contaminate protein extracts. In addition, complex matrices of polysaccharides, including large amount of anionic pectins, are present. All these species compete with the binding of proteins to the CPLL beads, impeding proper capture and identification / detection of low-abundance species. When properly pre-treated, plant tissue extracts are amenable to capture by the CPLL beads revealing thus many new species among them low-abundance proteins. Examples are given on the treatment of leaf proteins, of corn seed extracts and of exudate proteins (latex from Hevea brasiliensis). In all cases, the detection of unique gene products via CPLL capture is at least twice that of control, untreated sample.

Nongenomic signaling of the retinoid X receptor through binding and inhibiting Gq in human platelets

Retinoid X receptors (RXRs) are important transcriptional nuclear hormone receptors, acting as either homodimers or the binding partner for at least one fourth of all the known human nuclear receptors. Functional nongenomic effects of nuclear receptors are poorly understood; however, recently peroxisome proliferator-activated receptor (PPAR) gamma, PPAR beta, and the glucocorticoid receptor have all been found active in human platelets. Human platelets express RXR alpha, and RXR beta. RXR ligands inhibit platelet aggregation and TXA(2) release to ADP and the TXA(2) receptors, but only weakly to collagen. ADP and TXA(2) both signal via the G protein, Gq. RXR rapidly binds Gq but not Gi/z/o/t/gust in a ligand-dependent manner and inhibits Gq-induced Rac activation and intracellular calcium release. We propose that RXR ligands may have beneficial clinical actions through inhibition of platelet activation. Furthermore, our results demonstrate a novel nongenomic mode for nuclear receptor action and a functional cross-talk between G-protein and nuclear receptor signaling families. (C) 2007 by The American Society of Hematology.

Investigation of cooperativity in the binding of ligands to the D-2 dopamine receptor

The D 2 dopamine receptor exists as dimers or as higher-order oligomers, as determined from data from physical experiments. In this study, we sought evidence that this oligomerization leads to cooperativity by examining the binding of three radioligands ([H-3] nemonapride, [H-3] raclopride, and [H-3] spiperone) to D 2 dopamine receptors expressed in membranes of Sf9 cells. In saturation binding experiments, the three radioligands exhibited different B-max values, and the B-max values could be altered by the addition of sodium ions to assays. Despite labeling different numbers of sites, the different ligands were able to achieve full inhibition in competition experiments. Some ligand pairs also exhibited complex inhibition curves in these experiments. In radioligand dissociation experiments, the rate of dissociation of [H-3] nemonapride or [H-3] spiperone depended on the sodium ion concentration but was independent of the competing ligand. Although some of the data in this study are consistent with the behavior of a cooperative oligomeric receptor, not all of the data are in agreement with this model. It may, therefore, be necessary to consider more complex models for the behavior of this receptor.

Insights into the effects on metal binding of the systematic substitution of five key glutamate ligands in the ferritin of Escherichia coli

Ferritins are nearly ubiquitous iron storage proteins playing a fundamental role in iron metabolism. They are composed of 24 subunits forming a spherical protein shell encompassing a central iron storage cavity. The iron storage mechanism involves the initial binding and subsequent O-2-dependent oxidation of two Fe2+ ions located at sites A and B within the highly conserved dinuclear "ferroxidase center" in individual subunits. Unlike animal ferritins and the heme-containing bacterioferritins, the Escherichia coli ferritin possesses an additional iron-binding site (site C) located on the inner surface of the protein shell close to the ferroxidase center. We report the structures of five E. coli ferritin variants and their Fe3+ and Zn2+ (a redox-stable alternative for Fe2+) derivatives. Single carboxyl ligand replacements in sites A, B, and C gave unique effects on metal binding, which explain the observed changes in Fe2+ oxidation rates. Binding of Fe2+ at both A and B sites is clearly essential for rapid Fe2+ oxidation, and the linking of Fe-B(2+) to Fe-C(2+) enables the oxidation of three Fe2+ ions. The transient binding of Fe2+ at one of three newly observed Zn2+ sites may allow the oxidation of four Fe2+ by one dioxygen molecule.

1-(2 '-pyridylazo)-2-naphtholate complexes of ruthenium: synthesis, characterization, and DNA binding properties

Reaction of 1-(2'-pyridylazo)-2 -naphthol (Hpan) with [Ru(dmso)(4)Cl-2] (dmso=dimethylsulfoxide), [Ru(trpy)Cl-3] (trpy=2,2',2 ''-terpyridine), [Ru(bpy)Cl-3] (bpy=2,2'-bipyridine) and [Ru(PPh3)(3)Cl-2] in refluxing ethanol in the presence of a base (NEt3) affords, respectively, the [Ru(pan)(2)], [Ru(trpy)(pan)](+) (isolated as perchlorate salt), [Ru(bpy)(pan)Cl] and [Ru(PPh3)(2)(pan)Cl] complexes. Structures of these four complexes have been determined by X-ray crystallography. in each of these complexes, the pan ligand is coordinated to the metal center as a monoanionic tridentate N,N,O-donor. Reaction of the [Ru(bpy)(pan)Cl] complex with pyridine (py) and 4-picoline (pic) in the presence of silver ion has yielded the [Ru(bpy)(pan)(py)](+) and [Ru(bpy)(pan)(pic)](+) complexes (isolated as perchlorate salts), respectively. All the complexes are diamagnetic (low-spin d(6), S = 0) and show characteristic H-1 NMR signals and intense MLCT transitions in the visible region. Cyclic voltammetry on all the complexes shows a Ru(II)-Ru(III) oxidation on the positive side of SCE. Except in the [Ru(pan)(2)] complex, a second oxidative response has been observed in the other five complexes. Reductions of the coordinated ligands have also been observed on the negative side of SCE. The [Ru(trpy)(pan)]ClO4, [Ru(bpy)(pan)(py)]ClO4 and [Ru(bpy) (pan)(pic)]ClO4 complexes have been observed to bind to DNA, but they have not been able to cleave super-coiled DNA on UV irradiation. (c) 2008 Elsevier Ltd. All rights reserved.

Combinatorial peptide ligand libraries and plant proteomics: a winning strategy at a price

The mechanism of action and properties of a solid-phase ligand library made of hexapeptides (combinatorial peptide ligand libraries or CPLL, for capturing the "hidden proteome", i.e. the low- and very low-abundance proteins Constituting the vast majority of species in any proteome. as applied to plant tissues, are reviewed here. Plant tissues are notoriously recalcitrant to protein extraction and to proteome analysis, Firstly, rigid plant cell walls need to be mechanically disrupted to release the cell content and, in addition to their poor protein yield, plant tissues are rich in proteases and oxidative enzymes, contain phenolic Compounds, starches, oils, pigments and secondary metabolites that massively contaminate protein extracts. In addition, complex matrices of polysaccharides, including large amount of anionic pectins, are present. All these species compete with the binding of proteins to the CPLL beads, impeding proper capture and identification I detection of low-abundance species. When properly pre-treated, plant tissue extracts are amenable to capture by the CPLL beads revealing thus many new species among them low-abundance proteins. Examples are given on the treatment of leaf proteins, of corn seed extracts and of exudate proteins (latex from Hevea brasiliensis). In all cases, the detection of unique gene products via CPLL Capture is at least twice that of control, untreated sample. (c) 2008 Elsevier B.V. All rights reserved.

Ligand bridging of the DNA Holliday junction: molecular recognition of a stacked-X four-way junction by a small molecule

Planning a Holliday: A new mode of binding to a stacked-X, four-way Holliday junction is described in which a chromophore molecule binds across the center of the junction and two adenine residues are replaced by the acridine chromophores at either side of the crossover. This binding mode is specific for the Holliday junction and does not cause unwinding of the DNA helices.

Carboxylate anions binding and sensing by a novel tetraazamacrocycle containing ferrocene as receptor

A tetraazamacrocycle containing ferrocene moieties has been synthesized and characterized. The tetraprotonated form of this compound was evaluated as a receptor (R) for anion recognition of several substrates (S), Cl-, PF6-, HSO4-, H2PO4- and carboxylates, such as p-nitrobenzoate (p-nbz(-)), phthalate (ph(2-)), isophthalate (iph(2-)) and dipicolinate (dipic(2-)). H-1 NMR titrations in CD3OD indicated that this receptor is not suitable for recognizing HSO4- and H2PO4-, but weakly binds p-nbz(-), and strongly interacts with ph(2-), dipic(2-), and iph(2-) anions forming 1 : 2 assembled species. The largest beta(2) binding constant was determined for ph(2-), followed by dipic(2-) and finally iph(2-). The effect of the anionic substrates on the electron-transfer process of the ferrocene units of R was evaluated using cyclic voltammetry (CV) and square wave voltammetry (SWV) in methanol solution and 0.1 mol dm(-3) (CH3)(4)NCl as the supporting electrolyte. Titrations of the receptor were undertaken by addition of anion solutions in their tetrabutylammonium or tetramethylammonium forms. The protonated ligand exhibits a reversible voltammogram, which shifts cathodically in the presence of the substrates. The data revealed kinetic constraints in the formation of the receptor/substrate entity for dipic(2-), ph(2-) and iph(2-) anions, but not for p-nbz(-). In spite of the slow kinetics of assembled species formation with the ph(2-) substrate, this anion provides the largest redox-response when the supramolecular entity is formed, followed by dipic(2-), iph(2-) and finally p-nbz(-) anions. This trend is in agreement with the H-1 NMR results and the values of the binding constants. Single crystal X-ray structures of the receptor with PF6-, ph(2-), iph(2-) and p-nbz(-) were carried out and showed that supermolecules with a RS2 stoichiometry are formed with the first three anions, but RS4 with p-nbz(-). In all cases the binding occurs outside the macrocyclic cavity via N-H center dot center dot center dot O=C hydrogen bonds for carboxylate anions and N - H center dot center dot center dot F hydrogen bonds for the PF6- anion, which is in agreement with the solution results. The macrocyclic framework adopts different conformations in order to interact with each substrate having Fe center dot center dot center dot Fe intramolecular distances ranging from 10.125(14) to 12.783(15) angstrom.

A fluxional (CuN2O2)-N-I core: binding of a keto oxygen atom to Cu-I and Ag-I

[Cu(2-acetylpyridine)(2)]ClO4 (1), characterised here, has a novel Cu'N202 core in the solid state. Variable-temperature H-1 NMR studies show that the two chelate rings open up in solution at room temperature and the keto oxygen atoms dangle freely. As the temperature is lowered, the 0 atoms tend to bind to the metal atom. The corresponding silver(I) complex, [Ag(2-acetylpyridine)2]ClO4 (4), characterised by single-crystal X-ray crystallography, has an (AgN2)-N-I core in the solid state as well as in solution. Thus, while 1 is fluxional, 4 is not. In cyclic voltammetry, complex 1 displays a quasireversible Cu-II/I couple with a half-wave potential of 0.40 V vs. SCE. Complex I is easily oxidised by air and H2O2 in methanol to give rise to a dinuclear copper(II) complex where the ligand framework is not simple acetylpyridine. ((c) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2005).