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.

Luminescence-based assay methods in drug discovery

The drug discovery process is facing new challenges in the evaluation process of the lead compounds as the number of new compounds synthesized is increasing. The potentiality of test compounds is most frequently assayed through the binding of the test compound to the target molecule or receptor, or measuring functional secondary effects caused by the test compound in the target model cells, tissues or organism. Modern homogeneous high-throughput-screening (HTS) assays for purified estrogen receptors (ER) utilize various luminescence based detection methods. Fluorescence polarization (FP) is a standard method for ER ligand binding assay. It was used to demonstrate the performance of two-photon excitation of fluorescence (TPFE) vs. the conventional one-photon excitation method. As result, the TPFE method showed improved dynamics and was found to be comparable with the conventional method. It also held potential for efficient miniaturization. Other luminescence based ER assays utilize energy transfer from a long-lifetime luminescent label e.g. lanthanide chelates (Eu, Tb) to a prompt luminescent label, the signal being read in a time-resolved mode. As an alternative to this method, a new single-label (Eu) time-resolved detection method was developed, based on the quenching of the label by a soluble quencher molecule when displaced from the receptor to the solution phase by an unlabeled competing ligand. The new method was paralleled with the standard FP method. It was shown to yield comparable results with the FP method and found to hold a significantly higher signal-tobackground ratio than FP. Cell-based functional assays for determining the extent of cell surface adhesion molecule (CAM) expression combined with microscopy analysis of the target molecules would provide improved information content, compared to an expression level assay alone. In this work, immune response was simulated by exposing endothelial cells to cytokine stimulation and the resulting increase in the level of adhesion molecule expression was analyzed on fixed cells by means of immunocytochemistry utilizing specific long-lifetime luminophore labeled antibodies against chosen adhesion molecules. Results showed that the method was capable of use in amulti-parametric assay for protein expression levels of several CAMs simultaneously, combined with analysis of the cellular localization of the chosen adhesion molecules through time-resolved luminescence microscopy inspection.

Homogeneous assay technologies in drug screening: Quenching Resonance Energy Transfer (QRET) technique

Information gained from the human genome project and improvements in compound synthesizing have increased the number of both therapeutic targets and potential lead compounds. This has evolved a need for better screening techniques to have a capacity to screen number of compound libraries against increasing amount of targets. Radioactivity based assays have been traditionally used in drug screening but the fluorescence based assays have become more popular in high throughput screening (HTS) as they avoid safety and waste problems confronted with radioactivity. In comparison to conventional fluorescence more sensitive detection is obtained with time-resolved luminescence which has increased the popularity of time-resolved fluorescence resonance energy transfer (TR-FRET) based assays. To simplify the current TR-FRET based assay concept the luminometric homogeneous single-label utilizing assay technique, Quenching Resonance Energy Transfer (QRET), was developed. The technique utilizes soluble quencher to quench non-specifically the signal of unbound fraction of lanthanide labeled ligand. One labeling procedure and fewer manipulation steps in the assay concept are saving resources. The QRET technique is suitable for both biochemical and cell-based assays as indicated in four studies:1) ligand screening study of β2 -adrenergic receptor (cell-based), 2) activation study of Gs-/Gi-protein coupled receptors by measuring intracellular concentration of cyclic adenosine monophosphate (cell-based), 3) activation study of G-protein coupled receptors by observing the binding of guanosine-5’-triphosphate (cell membranes), and 4) activation study of small GTP binding protein Ras (biochemical). Signal-to-background ratios were between 2.4 to 10 and coefficient of variation varied from 0.5 to 17% indicating their suitability to HTS use.

Metal-Ion-Binding Oligonucleotides: High-Affinity Probes for Nucleic Acid Sequences

Small non-coding RNAs have numerous biological functions in cell and are divided into different classes such as: microRNA, snoRNA, snRNA and siRNA. MicroRNA (miRNA) is the most studied non-coding RNA to date and is found in plants, animals and some viruses. miRNA with short sequences is involved in suppressing translation of target genes by binding to their mRNA post-transcriptionally and silencing it. Their function besides silencing of the viral gene, can be oncogenic and therefore the cause of cancer. Hence, their roles are highlighted in human diseases, which increases the interest in using them as biomarkers and drug targets. One of the major problems to overcome is recognition of miRNA. Owing to a stable hairpin structure, chain invasion by conventional Watson-Crick base-pairing is difficult. One way to enhance the hybridization is exploitation of metal-ion mediated base-pairing, i. e. oligonucleotide probes that tightly bind a metal ions and are able to form a coordinative bonds between modified and natural nucleobases. This kind of metallo basepairs containing short modified oligonucleotides can also be useful for recognition of other RNA sequences containing hairpin-like structural motives, such as the TAR sequence of HIV. In addition, metal-ion-binding oligonucleotides will undoubtedly find applications in DNA-based nanotechnology. In this study, the 3,5-dimethylpyrazol-1-yl substituted purine derivatives were successfully incorporated within oligonucleotides, into either a terminal or non-terminal position. Among all of the modified oligonucleotides studied, a 2-(3,5-dimethylpyrazol-1-yl)-6-oxopurine base containing oligonucleotide was observed to bind most efficiently to their unmodified complementary sequences in the presence of both Cu2+ or Zn2+. The oligonucleotide incorporating 2,6-bis(3,5-dimethylpyrazol-1-yl)purine base also markedly increased the stability of duplexes in the presence of Cu2+ without losing the selectivity.

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.