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.

Epidemiology of Cytochrome P450-mediated Drug-Drug Interactions

Drug-drug interactions (DDIs) comprise an important cause of adverse drug reactions leading to excess hospitalizations. Drug metabolism is catalyzed by 75% by cytochrome P450 (CYP) enzymes and thus they are often involved in pharmacokinetic DDIs. In general, DDIs are studied in randomized controlled clinical trials in selected study populations. The overall aim of the present studies was to perform observational pharmacoepidemiological surveys on CYP-mediated DDIs in diseases important at the population level. The prevalence of co-administrations of four prodrugs (losartan, codeine, tramadol, and clopidogrel), three sulphonylureas (glibenclamide, glimepiride, and glipizide), or two statins (lovastatin and simvastatin) with well established agents altering CYP activity, as well as of statins with fibrates, was studied in Finland utilizing data from a university hospital medication database (inpatients) and the National Prescription Register of the Social Insurance Institution of Finland, Kela (outpatients). Clinical consequences of potential DDIs were estimated by reviewing laboratory data, and information from hospital care and cause-of-death registers. Concomitant use of study substrates with interacting medication was detected in up to one fifth of patients in both hospital and community settings. Potential CYP3A4 interactions in statin users did not manifest in clear adverse laboratory values but pharmacodynamic DDIs between statins and fibrates predisposed patients to muscular toxicity. Sulphonylurea DDIs with CYP2C9 inhibitors increased the risk of hypoglycaemia. CYP3A4 inhibitor use with clopidogrel was not associated with significant changes in mortality but non-fatal thrombosis and haemorrhage complications were seen less often in this group. Concomitant administration of atorvastatin with clopidogrel moderately attenuated the antithrombotic effect by clopidogrel. The overall mortality was increased in CYP3A4 inducer and clopidogrel co-users. Atorvastatin used concomitantly with prodrug clopidogrel seems to be beneficial in terms of total and LDL cholesterol concentrations, and overall mortality compared with clopidogrel use without interacting medication. In conclusion, CYP-mediated DDIs are a common and often unrecognized consequence of irrational drug prescribing.

The spindle assembly checkpoint as a drug target - Novel small-molecule inhibitors of Aurora kinases

Cell division (mitosis) is a fundamental process in the life cycle of a cell. Equal distribution of chromosomes between the daughter cells is essential for the viability and well-being of an organism: loss of fidelity of cell division is a contributing factor in human cancer and also gives rise to miscarriages and genetic birth defects. For maintaining the proper chromosome number, a cell must carefully monitor cell division in order to detect and correct mistakes before they are translated into chromosomal imbalance. For this purpose an evolutionarily conserved mechanism termed the spindle assembly checkpoint (SAC) has evolved. The SAC comprises a complex network of proteins that relay and amplify mitosis-regulating signals created by assemblages called kinetochores (KTs). Importantly, minor defects in SAC signaling can cause loss or gain of individual chromosomes (aneuploidy) which promotes tumorigenesis while complete failure of SAC results in cell death. The latter event has raised interest in discovery of low molecular weight (LMW) compounds targeting the SAC that could be developed into new anti-cancer therapeutics. In this study, we performed a cell-based, phenotypic high-throughput screen (HTS) to identify novel LMW compounds that inhibit SAC function and result in loss of cancer cell viability. Altogether, we screened 65 000 compounds and identified eight that forced the cells prematurely out of mitosis. The flavonoids fisetin and eupatorin, as well as the synthetic compounds termed SACi2 and SACi4, were characterized in more detail utilizing versatile cell-based and biochemical assays. To identify the molecular targets of these SAC-suppressing compounds, we investigated the conditions in which SAC activity became abrogated. Eupatorin, SACi2 and SACi4 preferentially abolished the tensionsensitive arm of the SAC, whereas fisetin lowered also the SAC activity evoked by lack of attachments between microtubules (MTs) and KTs. Consistent with the abrogation of SAC in response to low tension, our data indicate that all four compounds inhibited the activity of Aurora B kinase. This essential mitotic protein is required for correction of erratic MT-KT attachments, normal SAC signaling and execution of cytokinesis. Furthermore, eupatorin, SACi2 and SACi4 also inhibited Aurora A kinase that controls the centrosome maturation and separation and formation of the mitotic spindle apparatus. In line with the established profound mitotic roles of Aurora kinases, these small compounds perturbed SAC function, caused spindle abnormalities, such as multi- and monopolarity and fragmentation of centrosomes, and resulted in polyploidy due to defects in cytokinesis. Moreover, the compounds dramatically reduced viability of cancer cells. Taken together, using a cell-based HTS we were able to identify new LMW compounds targeting the SAC. We demonstrated for the first time a novel function for flavonoids as cellular inhibitors of Aurora kinases. Collectively, our data support the concept that loss of mitotic fidelity due to a non-functional SAC can reduce the viability of cancer cells, a phenomenon that may possess therapeutic value and fuel development of new anti-cancer drugs.

Cell Model Systems in Characterizing Alpha2-Adrenoceptors as Drug Targets.

The three alpha2-adrenoceptor (alpha2-AR) subtypes belong to the G protein-coupled receptor superfamily and represent potential drug targets. These receptors have many vital physiological functions, but their actions are complex and often oppose each other. Current research is therefore driven towards discovering drugs that selectively interact with a specific subtype. Cell model systems can be used to evaluate a chemical compound's activity in complex biological systems. The aim of this thesis was to optimize and validate cell-based model systems and assays to investigate alpha2-ARs as drug targets. The use of immortalized cell lines as model systems is firmly established but poses several problems, since the protein of interest is expressed in a foreign environment, and thus essential components of receptor regulation or signaling cascades might be missing. Careful cell model validation is thus required; this was exemplified by three different approaches. In cells heterologously expressing alpha2A-ARs, it was noted that the transfection technique affected the test outcome; false negative adenylyl cyclase test results were produced unless a cell population expressing receptors in a homogenous fashion was used. Recombinant alpha2C-ARs in non-neuronal cells were retained inside the cells, and not expressed in the cell membrane, complicating investigation of this receptor subtype. Receptor expression enhancing proteins (REEPs) were found to be neuronalspecific adapter proteins that regulate the processing of the alpha2C-AR, resulting in an increased level of total receptor expression. Current trends call for the use of primary cells endogenously expressing the receptor of interest; therefore, primary human vascular smooth muscle cells (SMC) expressing alpha2-ARs were tested in a functional assay monitoring contractility with a myosin light chain phosphorylation assay. However, these cells were not compatible with this assay due to the loss of differentiation. A rat aortic SMC cell line transfected to express the human alpha2B-AR was adapted for the assay, and it was found that the alpha2-AR agonist, dexmedetomidine, evoked myosin light chain phosphorylation in this model.

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.

The alpha2C-adrenoceptor as a neuropsychiatric drug tar-get - PET studies in human subjects

Positron emission tomography imaging has both academic and applied uses in revealing the distribution and density of different molecular targets in the central nervous system. Following the significant progress made with the dopamine D2 receptor, advances have been made in developing PET tracers to allow analysis of receptor occupancy of many other receptor types as well as evaluating changes in endogenous synaptic transmitter concentrations of transmitters e.g. serotonin and noradrenaline. Noradrenergic receptors are divided into α1-, α2- and β-adrenoceptor subfamilies, in humans each of which is composed of three receptor subtypes. The α2-adrenoceptors have an important presynaptic auto-inhibitory function on noradrenaline release but they also have postsynaptic roles in modulating the release of other neurotransmitters, such as serotonin and dopamine. One of the subtypes, the α2C-adrenoceptor, has been detected at distinct locations in the central nervous system, most notably the dorsal striatum. Several serious neurological conditions causing dementia, Alzheimer’s disease and Parkinson’s disease have been linked to disturbed noradrenergic signaling. Furthermore, altered noradrenergic signaling has also been implicated in conditions like ADHD, depression, anxiety and schizophrenia. In order to benefit future research into these central nervous system disorders as well as being useful in the clinical development of drugs affecting brain noradrenergic neurotransmission, validation work of a novel tracer for positron emission tomography studies in humans was performed. Altogether 85 PET imaging experiments were performed during four separate clinical trials. The repeatability of [11C]ORM-13070 binding was tested in healthy individuals, followed by a study to evaluate the dose-dependent displacement of [11C]ORM-13070 from α2C-adrenoceptors by a competing ligand, and the final two studies examined the sensitivity of [11C]ORM-13070 binding to reflect changes in endogenous noradrenaline levels. The repeatability of [11C]ORM-13070 binding was very high. The binding properties of the tracer allowed for a reliable estimation of α2C-AR occupancy by using the reference tissue ratio method with low test-retest variability. [11C]ORM-13070 was dose-dependently displaced from its specific binding sites by the subtype-nonselective α2-adrenoceptor antagonist atipamezole, and thus it proved suitable for use in clinical drug development of novel α2C-adrenoceptor ligands e.g. to determine the best doses and dosing intervals for clinical trials. Convincing experimental evidence was gained to support the suitability of [11C]ORM-13070 for detecting an increase in endogenous synaptic noradrenaline in the human brain. Tracer binding in the thalamus tended to increase in accordance with reduced activity of noradrenergic projections from the locus coeruleus, although statistical significance was not reached. Thus, the investigation was unable to fully validate [11C]ORM-13070 for the detection of pharmacologically evoked reductions in noradrenaline levels.

Probing MST1 Kinase Sequence and Structure for Rational Drug Discovery (Targeting diabetes by blocking protein-protein interactions)

Apoptotic beta cell death is an underlying cause majorly for type I and to a lesser extent for type II diabetes. Recently, MST1 kinase was identified as a key apoptotic agent in diabetic condition. In this study, I have examined MST1 and closely related kinases namely, MST2, MST3 and MST4, aiming to tackle diabetes by exploring ways to selectively block MST1 kinase activity. The first investigation was directed towards evaluating possibilities of selectively blocking the ATP binding site of MST1 kinase that is essential for the activity of the enzymes. Structure and sequence analyses of this site however revealed a near absolute conservation between the MSTs and very few changes with other kinases. The observed residue variations also displayed similar physicochemical properties making it hard for selective inhibition of the enzyme. Second, possibilities for allosteric inhibition of the enzyme were evaluated. Analysis of the recognized allosteric site also posed the same problem as the MSTs shared almost all of the same residues. The third analysis was made on the SARAH domain, which is required for the dimerization and activation of MST1 and MST2 kinases. MST3 and MST4 lack this domain, hence selectivity against these two kinases can be achieved. Other proteins with SARAH domains such as the RASSF proteins were also examined. Their interaction with the MST1 SARAH domain were evaluated to mimic their binding pattern and design a peptide inhibitor that interferes with MST1 SARAH dimerization. In molecular simulations the RASSF5 SARAH domain was shown to strongly interact with the MST1 SARAH domain and possibly preventing MST1 SARAH dimerization. Based on this, the peptidic inhibitor was suggested to be based on the sequence of RASSF5 SARAH domain. Since the MST2 kinase also interacts with RASSF5 SARAH domain, absolute selectivity might not be achieved.