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.

Conformational Regulation of alpha2beta1 Integrin in Ligand Binding

Integrins are heterodimeric cell adhesion receptors involved in cell-cell and cell-extracellular matrix (ECM) interactions. They transmit bidirectional signals across the cell membrane. This results in a wide range of biological events from cell differentiation to apoptosis. alpha2beta1 integrin is an abundant collagen receptor expressed on the surface of several cell types. In addition to ECM ligands, alpha2beta1 integrins are bound by echovirus 1 (EV1) which uses alpha2beta1 as a receptor to initiate its life cycle in the infected cell. The aim of this thesis project was to provide further insight into the mechanisms of alpha2beta1 integrin ligand recognition and receptor activation. Collagen fibrils are the principal tensile elements of the ECM. Yet, the interaction of alpha2beta1 integrin with the fibrillar form of collagen I has received relatively little attention. This research focused on the ability of alpha2beta1 integrin to act as a receptor for type I collagen fibrils. Also the molecular requirements of the EV1 interaction with alpha2beta1 were studied. Conventionally, ligand binding has been suggested to require integrin activation and the binding may further trigger integrin signalling. Another main objective of this study was to elucidate both the inside-out and outside-in signalling mechanisms of alpha2beta1 integrin in adherent cells. The results indicated that alpha2beta1 integrin is the principal integrin-type collagen receptor for type I collagen fibrils, and alpha2beta1 may participate in the regulation of pericellular collagen fibrillogenesis. Furthermore, alpha2beta1 integrin inside-out activation appeared to be synergistically regulated by integrin clustering and conformational activation. The triggering of alpha2beta1 integrin outside-in signalling, however, was shown to require both conformational changes and clustering. In contrast to ECM ligands, EV1 appeared to take advantage of the bent, inactive form of alpha2beta1 integrin in initiating its life cycle in the cell. This research together with other recent studies, has shed light on the molecular mechanisms of integrin activation. It is becoming evident that large ligands are able to bind to the bent form of integrin, which has been previously considered to be physiologically inactive. Consequently, our understanding of the conformational modulation of integrins upon activation is changing.

Alpha2-adrenoceptors: Structure and ligand binding properties at the molecular level

Alpha₂-Adrenoceptors: structure and ligand binding properties at the molecular level The mouse is the most frequently used animal model in biomedical research, but the use of zebrafish as a model organism to mimic human diseases is on the increase. Therefore it is considered important to understand their pharmacological differences from humans also at the molecular level. The zebrafish Alpha_2-adrenoceptors were expressed in mammalian cells and the binding affinities of 20 diverse ligands were determined and compared to the corresponding human receptors. The pharmacological properties of the human and zebrafish Alpha_2--adrenoceptors were found to be quite well conserved. Receptor models based on the crystal structures of bovine rhodopsin and the human Beta₂-adrenoceptor revealed that most structural differences between the paralogous and orthologous Alpha_2--adrenoceptors were located within the second extracellular loop (XL2). Reciprocal mutations were generated in the mouse and human Alpha_2--adrenoceptors. Ligand binding experiments revealed that substitutions in XL2 reversed the binding profiles of the human and mouse Alpha_2--adrenoceptors for yohimbine, rauwolscine and RS-79948-197, evidence for a role for XL2 in the determination of species-specific ligand binding. Previous mutagenesis studies had not been able to explain the subtype preference of several large Alpha_2--adrenoceptor antagonists. We prepared chimaeric Alpha_2--adrenoceptors where the first transmembrane (TM1) domain was exchanged between the three human Alpha_2--adrenoceptor subtypes. The binding affinities of spiperone, spiroxatrine and chlorpromazine were observed to be significantly improved by TM1 substitutions of the Alpha_2a--adrenoceptor. Docking simulations indicated that indirect effects, such as allosteric modulation, are more likely to be involved in this phenomenon rather than specific side-chain interactions between ligands and receptors.

Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.

Probabilistic contact preferences in protein-ligand and protein-protein complexes

Modeller för intermolekulär växelvärkan utnyttjas brett inom biologin. Analys av kontakter mellan proteiner och läkemedelsforskning representerar typiska tillämpningsområden för dylika modeller. En modell som beskriver sådana molekylära växelverkningar kan utformas med hjälp av biofysisk teori, vilket tenderar att resultera i ytterst tung beräkningsbörda även för enkla tillämpningar. Ett alternativt sätt att formulera modeller är att utnyttja stora databaser som innehåller strukturmätningar gjorda med hjälp av till exempel röntgendiffraktion. Då man använder sig av empiriska mätdata direkt, möjliggör en statistisk modell att osäkerheten och inexaktheten i datat tas till hänsyn på ett adekvat sätt, samtidigt som beräkningsbördan håller sig på en rimligare nivå jämfört med kvantmekaniska metoder som i princip borde ge de optimala resultaten. I avhandlingen utvecklades en 3D modell för numerisk undersökning av intermolekulär växelverkan baserad på Bayesiansk statistik. Modellens syfte är att åstadkomma prognoser för det hurdana eller vilka molekylstrukturer prefereras i en given kontext, d.v.s. är mer sannolika inom ramen för interaktion. Modellen testades i essentiella molekyläromgivningar - en liten molekyl vid sin bindningsplats hos ett protein och en gränsyta mellan proteinerna i ett komplex. De erhållna numeriska resultaten motsvarar väl experimentella resultat som tidigare rapporterats i litteraturen, exempelvis kvalitativa bindningsaffiniteter och kemisk kännedom av vissa aminosyrors rumsliga förmågor att utgöra bindningar. I avhandlingen gjordes ytterligare preliminära tester av den statistiska ansatsen för modellering av den centrala molekylära strukturella anpassningsbarheten. I praktiken är den utvecklade modellen ämnad som ett led i en mer omfattande analysmetod, så som en s.k. farmakofor modell. Molekyylivuorovaikutusten mallintamista hyödynnetään laajasti biologisten kysymysten tarkastelussa. Tyypillisiä esimerkkejä sovelluskohteista ovat proteiinien väliset kontaktit ja lääkesuunnittelu. Vuorovaikutuksia kuvaavan mallin lähtökohta voi olla molekyyleihin liittyvä teoria, jolloin soveltamiseen liittyvä laskenta saattaa olla erityisen raskasta, tai suuri havaintojoukko joka on saatu aikaan esimerkiksi mittaamalla rakenteita röntgendiffraktio menetelmällä. Tilastollinen malli mahdollistaa havaintoaineistossa olevan epätarkkuuden ja epävarmuuden huomioimisen, samalla pitäen laskennallisen kuorman pienempänä verrattuna periaatteessa parhaan tuloksen antavaan kvanttimekaaniseen mallinnukseen. Väitöstyössä kehitettiin bayesiläiseen tilastotieteeseen perustuva 3D malli molekyylien välisten vuorovaikutusten laskennalliseen tarkasteluun. Mallin tehtävä on tuottaa ennusteita sen suhteen, minkä tai millaisten molekyylirakenteiden väliset kompleksit ovat etusijalla, toisin sanoen todennäköisempiä, vuorovaikutustilanteessa. Työssä kehitetyn menetelmän toimivuutta testattiin käyttötarkoituksen suhteen olennaisissa molekyyliympäristöissä - pieni molekyyli sitoutumiskohdassaan proteiinissa sekä rajapinta kahden proteiinin välilllä proteiinikompleksissa. Saadut laskennalliset tulokset vastasivat hyvin vertailuun käytettyjä kirjallisuudesta saatuja kokeellisia tuloksia, kuten laadullisia sitoutumisaffiniteetteja, sekä kemiallista tietoa esimerkiksi tiettyjen aminohappojen avaruudellisesta sidoksenmuodostuksesta. Väitöstyössä myös alustavasti testattiin tilastollista lähestymistapaa tärkeän molekyylien rakenteellisen mukautuvuuden mallintamiseen. Käytännössä malli on tarkoitettu osaksi jotakin laajempaa analyysimenetelmää, kuten farmakoforimallia.

Synthesis of Short Oligonucleotides on a Soluble Support by the Phosphoramidite Method

In the last decades, the chemical synthesis of short oligonucleotides has become an important aspect of study due to the discovery of new functions for nucleic acids such as antisense oligonucleotides (ASOs), aptamers, DNAzymes, microRNA (miRNA) and small interfering RNA (siRNA). The applications in modern therapies and fundamental medicine on the treatment of different cancer diseases, viral infections and genetic disorders has established the necessity to develop scalable methods for their cheaper and easier industrial manufacture. While small scale solid-phase oligonucleotide synthesis is the method of choice in the field, various challenges still remain associated with the production of short DNA and RNA-oligomers in very large quantities. On the other hand, solution phase synthesis of oligonucleotides offers a more predictable scaling-up of the synthesis and is amenable to standard industrial manufacture techniques. In the present thesis, various protocols for the synthesis of short DNA and RNA oligomers have been studied on a peracetylated and methylated β-cyclodextrin, and also on a pentaerythritol-derived support. On using the peracetylated and methylated β-cyclodextrin soluble supports, the coupling cycle was simplified by replacement of the typical 5′-O-(4,4′-dimethoxytrityl) protecting group with an acid-labile acetal-protected 5′-O-(1-methoxy-1-methylethyl) group, which upon acid-catalyzed methanolysis released easily removable volatile products. For this reason monomeric building blocks 5′-O-(1-methoxy-1-methylethyl) 3′-(2-cyano-ethyl-N,N-diisopropylphosphoramidite) were synthesized. Alternatively, on using the precipitative pentaerythritol support, novel 2´-O-(2-cyanoethyl)-5´-O-(1-methoxy-1-methylethyl) protected phosphoramidite building blocks for RNA synthesis have been prepared and their applicability by the synthesis of a pentamer was demonstrated. Similarly, a method for the preparation of short RNAs from commercially available 5´-O-(4,4´-dimethoxytrityl)-2´-O-(tert-butyldimethyl-silyl)ribonucleoside 3´-(2-cyanoethyl-N,N-diisopropylphosphoramidite) building blocks has been developed

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.