Gene Regulation in The Human Immune System. Gene Expression Signatures of Th2 Cell Differentiation, Type 1 Diabetes, and Intrauterine Immune Adaptation

The human immune system is constantly interacting with the surrounding stimuli and microorganisms. However, when directed against self or harmless antigens, these vital defense mechanisms can cause great damage. In addition, the understanding the underlying mechanism of several human diseases caused by aberrant immune cell functions, for instance type 1 diabetes and allergies, remains far from being complete. In this Ph.D. study these questions were addressed using genome-wide transcriptomic analyses. Asthma and allergies are characterized by a hyperactive response of the T helper 2 (Th2) immune cells. In this study, the target genes of the STAT6 transcription factor in naïve human T cells were identified with RNAi for the first time. STAT6 was shown to act as a central activator of the genes expression upon IL-4 signaling, with both direct and indirect effects on Th2 cell transcriptome. The core transcription factor network induced by IL-4 was identified from a kinetic analysis of the transcriptome. Type 1 diabetes is an autoimmune disease influenced by both the genetic susceptibility of an individual and the disease-triggering environmental factors. To improve understanding of the autoimmune processes driving pathogenesis in the prediabetic phase in humans, a unique series of prospective whole-blood RNA samples collected from HLA-susceptible children in the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) study was studied. Changes in different timewindows of the pathogenesis process were identified, and especially the type 1 interferon response was activated early and throughout the preclinical T1D. The hygiene hypothesis states that allergic diseases, and lately also autoimmune diseases, could be prevented by infections and other microbial contacts acquired in early childhood, or even prenatally. To study the effects of the standard of hygiene on the development of neonatal immune system, cord blood samples from children born in Finland (high standard of living), Estonia (rapid economic growth) and Russian Karelia (low standard of living) were compared. Children born in Russian Karelia deviated from Finnish and Estonian children in many aspects of the neonatal immune system, which was developmentally more mature in Karelia, resembling that of older infants. The results of this thesis offer significant new information on the regulatory networks associated with immune-mediated diseases in human. The results will facilitate understanding and further research on the role of the identified target genes and mechanisms driving the allergic inflammation and type 1 diabetes, hopefully leading to a new era of drug development.

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.

The Structure and Function of αI Domains in Collagen Receptor and Leukocyte Integrins

Integrins are heterodimeric, signaling transmembrane adhesion receptors that connect the intracellular actin microfilaments to the extracellular matrix composed of collagens and other matrix molecules. Bidirectional signaling is mediated via drastic conformational changes in integrins. These changes also occur in the integrin αI domains, which are responsible for ligand binding by collagen receptor and leukocyte specific integrins. Like intact integrins, soluble αI domains exist in the closed, low affinity form and in the open, high affinity form, and so it is possible to use isolated αI domains to study the factors and mechanisms involved in integrin activation/deactivation. Integrins are found in all mammalian tissues and cells, where they play crucial roles in growth, migration, defense mechanisms and apoptosis. Integrins are involved in many human diseases, such as inflammatory, cardiovascular and metastatic diseases, and so plenty of effort has been invested into developing integrin specific drugs. Humans have 24 different integrins, four of which are collagen receptor (α1β1, α2β1, α10β1, α11β1) and five leukocyte specific integrins (αLβ2, αMβ2, αXβ2, αDβ2, αEβ7). These two integrin groups are quite unselective having both primary and secondary ligands. This work presents the first systematic studies performed on these integrin groups to find out how integrin activation affects ligand binding and selectivity. These kinds of studies are important not only for understanding the partially overlapping functions of integrins, but also for drug development. In general, our results indicated that selectivity in ligand recognition is greatly reduced upon integrin activation. Interestingly, in some cases the ligand binding properties of integrins have been shown to be cell type specific. The reason for this is not known, but our observations suggest that cell types with a higher integrin activation state have lower ligand selectivity, and vice versa. Furthermore, we solved the three-dimensional structure for the activated form of the collagen receptor α1I domain. This structure revealed a novel intermediate conformation not previously seen with any other integrin αI domain. This is the first 3D structure for an activated collagen receptor αI domain without ligand. Based on the differences between the open and closed conformation of the αI domain we set structural criteria for a search for effective collagen receptor drugs. By docking a large number of molecules into the closed conformation of the α2I domain we discovered two polyketides, which best fulfilled the set structural criteria, and by cell adhesion studies we showed them to be specific inhibitors of the collagen receptor integrins.

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.

Models to study intermediate filament dynamics and functions

Intermediate filaments are part of the cytoskeleton and nucleoskeleton; they provide cells with structure and have important roles in cell signalling. The IFs are a large protein family with more than 70 members; each tightly regulated and expressed in a cell type-specific manner. Although the IFs have been known and studied for decades, our knowledge about their specific functions is still limited, despite the fact that mutations in IF genes cause numerous severe human diseases. In this work, three IF proteins are examined more closely; the nuclear lamin A/C and the cytoplasmic nestin and vimentin. In particular the regulation of lamin A/C dynamics, the role of nestin in muscle and body homeostasis as well as the functions and evolutionary aspects of vimentin are investigated. Together this data highlights some less well understood functions of these IFs. We used mass-spectrometry to identify inter-phase specific phosphorylation sites on lamin A. With the use of genetically engineered lamin A protein in combination with high resolution microscopy and biochemical methods we discovered novel roles for this phosphorylation in regulation of lamin dynamics. More specifically, our data suggests that the phosphorylation of certain amino acids in lamin A determines the localization and dynamics of the protein. In addition, we present results demonstrating that lamin A regulates Cdk5-activity. In the second study we use mice lacking nestin to gain more knowledge of this seldom studied protein. Our results show that nestin is essential for muscle regeneration; mice lacking nestin recover more slowly from muscle injury and show signs of spontaneous muscle regeneration, indicating that their muscles are more sensitive to stresses and injury. The absence of nestin also leads to decreased over-all muscle mass and slower body growth. Furthermore, nestin has a role in controlling testicle homeostasis as nestin-/- male mice show a greater variation in testicle size. The common fruit fly Drosophila melanogaster lacks cytoplasmic IFs as most insects do. By creating a fly that expresses human vimentin we establish a new research platform for vimentin studies, as well as provide a new tool for the studies of IF evolution.

The role of keratin intermediate filaments in the colon epithelial cells

Keratins (K) are cytoskeletal proteins mainly expressed in the epithelium and constitute the largest subgroup of intermediate filaments (IFs). Simple epithelial keratins (SEKs) K7-K8 and K18-K20 are the major IF elements in the colon. SEK mutations are known to cause around 30 human diseases, mainly affecting liver and skin. However, so far no strong associations between K8 mutations and the development of human colitis have been found. The keratin contribution to colonic health comes from the K8 knock-out (K8-/-) mouse model, which develops an early chronic inflammation and hyperproliferation in the colon. The aim of this thesis was to investigate how keratins contribute to intestinal health and disease mainly by the experimental analysis using the K8-/- mouse colon and cell culture models. The work described here is divided into three studies. The first study revealed involvement of keratins in Notch1 signaling, which is the master regulator of cell fate in the colon. Immunoprecipitation and immunostaining, both in vitro and in vivo showed that K8 binds and co-localizes with Notch1. Interestingly, overexpression of keratins enhanced Notch1 levels and stabilized Notch intracellular domain (NICD), leading to higher activity of Notch signaling. The dramatic decrease in Notch activity in the K8-/- colon resulted in a differentiation shift towards goblet and enteroendocrine cells. The second study focused on the involvement of keratins in colitis-associated cancer (CAC). Although, the K8-/- inflamed colon did not develop colorectal cancer (CRC) spontaneously, it was dramatically more susceptible to induced CRC in two CRC models: azoxymethane (AOM) and multiple intestinal neoplasia (ApcMin/+). To understand how the loss of K8 contributes to CAC, the epithelial inflammasome signaling pathway was analyzed. The released component of active inflammasome, cleaved caspase-1 and its downstream protein, interleukin (IL)-18, were significantly increased in K8-/- and K8-/-ApcMin/+ colons. The inflammasome pathway has recently been suggested to control the levels of IL-22 binding protein (IL-22BP), which is a negative regulator of IL-22 activity. Interestingly, the activated inflammasome correlated with an upregulation of IL-22 and a complete loss of IL-22BP in the K8-null colons. The activation of IL-22 was confirmed by increased levels of downstream signaling, which is phosphorylated signal transducer and activator of transcription 3 (P-STAT3), a transcription factor promoting proliferation and tissue regeneration in the colon. The objective of the third study, was to examine the role of keratins in colon energy metabolism. A proteomic analysis identified mitochondrial 3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2) as the major ownregulated protein in the K8-/- colonocytes. HMGCS2 is the rate-limiting enzyme in ketogenesis, where energy from bacterially produced short chain fatty acids (SCFAs), mainly butyrate, is converted into ketone bodies in colonic epithelium. Lower levels and activity of HMGCS2 in the K8-/- colon resulted in a blunted ketogenesis. The studies upstream from HMGCS2, identified decreased levels of the SCFA-transporter monocarboxylate transporter 1 (MCT1), which led to increased SCFA content in the stool suggesting impaired butyrate transport through the colonic epithelium. Taken together, the results of the herein thesis indicate that keratins are essential regulators of colon homeostasis, in particular epithelial differentiation, tumorigenesis and energy metabolism.

Immunopathogenesis of Asthma and Atopic Diseases – The specific Role of a selected Panel of Genes in human T helper Cell Differentiation

T helper cell (Th) functions are crucial for proper immune defence against various intra- and extracellular pathogens. According to the specific immune responses, Th cells can be classified into subtypes, Th1 and Th2 cells being the most frequently characterized classes. Th1 and Th2 cells interact with other immune cells by regulating their functions with specific cytokine production. IFN, IL-2 and TNF- are the cytokines predominantly produced by Th1 cells whereas Th2 cells produce Th2-type cytokines, such as IL-4, IL-5 and IL-13. Upon TCR activation and in the presence of polarizing cytokines, Th cells differentiate into effector subtypes from a common precursor cell. IFN and IL-12 are the predominant Th1 polarizing cytokines whereas IL-4 directs Th2 polarization. The cytokines mediate their effects through specific receptor signalling. The differentiation process is complex, involving various signalling molecules and routes, as well as functions of the specific transcription factors. The functions of the Th1/Th2 cells are tightly regulated; however, knowledge on human Th cell differentiation is, as yet, fairly poor. The susceptibility for many immune-mediated disorders often originates from disturbed Th cell responses. Thus, research is needed for defining the molecular mechanisms involved in the differentiation and balanced functions of the Th cells. Importantly, the new information obtained will be crucial for a better understanding of the pathogenesis of immune-mediated disorders, such as asthma or autoimmune diseases. In the first subproject of this thesis, the role of genetic polymorphisms in the human STAT6, GATA3 and STAT4 genes were investigated for asthma or atopy susceptibility in Finnish asthma families by association analysis. These genes code for key transcription factors regulating Th cell differentiation. The study resulted in the identification of a GATA3 haplotype that associated with asthma and related traits (high serum IgE level). In the second subproject, an optimized method for human primary T cell transfection and enrichment was established. The method can be utilized for functional studies for the selected genes of interest. The method was also utilized in the third subproject, which aimed at the identification of novel genes involved in early human Th cell polarization (0-48h) using genome-wide oligonucleotide arrays. As a result, numerous genes and ESTs with known or unknown functions were identified in the study. Using an shRNA knockdown approach, a panel of novel IL-4/STAT6 regulated genes were identified in the functional studies of the genes. Moreover, one of the genes, NDFIP2, with a previously uncharacterized role in the human Th differentiation, was observed to promote IFN production of the differentiated Th1 cells. Taken together, the results obtained have revealed potential new relevant candidate genes serving as a basis for further studies characterizing the detailed networks involved in the human Th cell differentiation as well as in the genetic susceptibility of Th-mediated immune disorders.

Role of human hydroxysteroid (17beta) dehydrogenase type 1 (HSD17B1) in steroid-dependent diseases in females –novel indications for HSD17B1 inhibitors. Phenotypic analysis of transgenic mice overexpressing human HSD17B1.

Hormone-dependent diseases, e.g. cancers, rank high in mortality in the modern world, and thus, there is an urgent need for new drugs to treat these diseases. Although the diseases are clearly hormone-dependent, changes in circulating hormone concentrations do not explain all the pathological processes observed in the diseased tissues. A more inclusive explanation is provided by intracrinology – a regulation of hormone concentrations at the target tissue level. This is mediated by the expression of a pattern of steroid-activating and -inactivating enzymes in steroid target tissues, thus enabling a concentration gradient between the blood circulation and the tissue. Hydroxysteroid (17beta) dehydrogenases (HSD17Bs) form a family of enzymes that catalyze the conversion between low active 17-ketosteroids and highly active 17beta-hydroxysteroids. HSD17B1 converts low active estrogen (E1) to highly active estradiol (E2) with high catalytic efficiency, and altered HSD17B1 expression has been associated with several hormone-dependent diseases, including breast cancer, endometriosis, endometrial hyperplasia and cancer, and ovarian epithelial cancer. Because of its putative role in E2 biosynthesis in ovaries and peripheral target tissues, HSD17B1 is considered to be a promising drug target for estrogen-dependent diseases. A few studies have indicated that the enzyme also has androgenic activity, but they have been ignored. In the present study, transgenic mice overexpressing human HSD17B1 (HSD17B1TG mice) were used to study the effects of the enzyme in vivo. Firstly, the substrate specificity of human HSD17B1 was determined in vivo. The results indicated that human HSD17B1 has significant androgenic activity in female mice in vivo, which resulted in increased fetal testosterone concentration and female disorder of sexual development appearing as masculinized phenotype (increased anogenital distance, lack of nipples, lack of vaginal opening, combination of vagina with urethra, enlarged Wolffian duct remnants in the mesovarium and enlarged female prostate). Fetal androgen exposure has been linked to polycystic ovary syndrome (PCOS) and metabolic syndrome during adulthood in experimental animals and humans, but the genes involved in PCOS are largely unknown. A putative mechanism to accumulate androgens during fetal life by HSD17B1 overexpression was shown in the present study. Furthermore, as a result of prenatal androgen exposure locally in the ovaries, HSD17B1TG females developed ovarian benign serous cystadenomas in adulthood. These benign lesions are precursors of low-grade ovarian serous tumors. Ovarian cancer ranks fifth in mortality of all female cancers in Finland, and most of the ovarian cancers arise from the surface epithelium. The formation of the lesions was prevented by prenatal antiandrogen treatment and by transplanting wild type (WT) ovaries prepubertally into HSD17B1TG females. The results obtained in our non-clinical TG mouse model, together with a literature analysis, suggest that HSD17B1 has a role in ovarian epithelial carcinogenesis, and especially in the development of serous tumors. The role of androgens in ovarian carcinogenesis is considered controversial, but the present study provides further evidence for the androgen hypothesis. Moreover, it directly links HSD17B1-induced prenatal androgen exposure to ovarian epithelial carcinogenesis in mice. As expected, significant estrogenic activity was also detected for human HSD17B1. HSD17B1TG mice had enhanced peripheral conversion of E1 to E2 in a variety of target tissues, including the uterus. Furthermore, this activity was significantly decreased by treatments with specific HSD17B1 inhibitors. As a result, several estrogen-dependent disorders were found in HSD17B1TG females. Here we report that HSD17B1TG mice invariably developed endometrial hyperplasia and failed to ovulate in adulthood. As in humans, endometrial hyperplasia in HSD17B1TG females was reversible upon ovulation induction, triggering a rise in circulating progesterone levels, and in response to exogenous progestins. Remarkably, treatment with a HSD17B1 inhibitor failed to restore ovulation, yet completely reversed the hyperplastic morphology of epithelial cells in the glandular compartment. We also demonstrate that HSD17B1 is expressed in normal human endometrium, hyperplasia, and cancer. Collectively, our non-clinical data and literature analysis suggest that HSD17B1 inhibition could be one of several possible approaches to decrease endometrial estrogen production in endometrial hyperplasia and cancer. HSD17B1 expression has been found in bones of humans and rats. The non-clinical data in the present study suggest that human HSD17B1 is likely to have an important role in the regulation of bone formation, strength and length during reproductive years in female mice. Bone density in HSD17B1TG females was highly increased in femurs, but in lesser amounts also in tibias. Especially the tibia growth plate, but not other regions of bone, was susceptible to respond to HSD17B1 inhibition by increasing bone length, whereas the inhibitors did not affect bone density. Therefore, HSD17B1 inhibitors could be safer than aromatase inhibitors in regard to bone in the treatment of breast cancer and endometriosis. Furthermore, diseases related to improper growth, are a promising new indication for HSD17B1 inhibitors.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Selective inhibition of human tankyrases

Tankyrases belong to the Diphtheria toxin-like ADP-ribosyltransferase (ARTD) enzyme superfamily, also known as poly(ADP-ribose) polymerases (PARPs). They catalyze a covalent post-translational modification reaction where they transfer ADP-ribose units from NAD+ to target proteins. Tankyrases are involved in many cellular processes and their roles in telomere homeostasis, Wnt signaling and in several diseases including cancers have made them interesting drug targets. In this thesis project, selective inhibition of human tankyrases was studied. A homogeneous fluorescence-based assay was developed to screen the compound libraries. The assay is inexpensive, operationally easy, and performs well according to the statistical analysis. Assay suitability was confirmed by screening a natural product library. Flavone was identified as the most potent inhibitor in the library and this motivated us to screen a larger flavonoid library. Results showed that flavones were indeed the best inhibitor of tankyrases among flavonoids. To further study the structure-activity relationship, a small library of flavones containing single substitution was screened and potency measurements allowed us to generate structure-activity relationship. Compounds containing substitutions at 4´-position were more potent in comparison to other substitutions, and importantly, hydrophobic groups improved isoenzyme selectivity as well as the potency. A flavone derivative containing a hydrophobic isopropyl group (compound 22), displayed 6 nM potency against TNKS1, excellent isoenzyme selectivity and Wnt signaling inhibition. Protein interactions with compounds were studied by solving complex crystal structures of the compounds with TNKS2 catalytic domain. A novel tankyrase inhibitor (IWR-1) was also crystallized in complex with TNKS2 catalytic domain. The crystal structure of TNKS2 in complex with IWR-1 showed that the compound binds to adenosine site and it was the first known ARTD inhibitor of this kind. To date, there is no structural information available about the substrate binding with any of the ARTD family members; therefore NAD+ was soaked with TNKS2 catalytic domain crystals. However, analysis of crystal structure showed that NAD+ was hydrolyzed to nicotinamide. Also, a co-crystal structure of NAD+ mimic compound, EB-47, was solved which was used to deduce some insights about the substrate interactions with the enzyme. Like EB-47, other ARTD1 inhibitors were also shown to inhibit tankyrases. It indicated that selectivity of the ARTD1 inhibitors should be considered as some of the effects in cells could come from tankyrase inhibition. In conclusion, the study provides novel information on tankyrase inhibition and presents new insight into the selectivity and potency of compounds.

Scalable Feature Selection Applications for Genome-Wide Association Studies of Complex Diseases

Personalized medicine will revolutionize our capabilities to combat disease. Working toward this goal, a fundamental task is the deciphering of geneticvariants that are predictive of complex diseases. Modern studies, in the formof genome-wide association studies (GWAS) have afforded researchers with the opportunity to reveal new genotype-phenotype relationships through the extensive scanning of genetic variants. These studies typically contain over half a million genetic features for thousands of individuals. Examining this with methods other than univariate statistics is a challenging task requiring advanced algorithms that are scalable to the genome-wide level. In the future, next-generation sequencing studies (NGS) will contain an even larger number of common and rare variants. Machine learning-based feature selection algorithms have been shown to have the ability to effectively create predictive models for various genotype-phenotype relationships. This work explores the problem of selecting genetic variant subsets that are the most predictive of complex disease phenotypes through various feature selection methodologies, including filter, wrapper and embedded algorithms. The examined machine learning algorithms were demonstrated to not only be effective at predicting the disease phenotypes, but also doing so efficiently through the use of computational shortcuts. While much of the work was able to be run on high-end desktops, some work was further extended so that it could be implemented on parallel computers helping to assure that they will also scale to the NGS data sets. Further, these studies analyzed the relationships between various feature selection methods and demonstrated the need for careful testing when selecting an algorithm. It was shown that there is no universally optimal algorithm for variant selection in GWAS, but rather methodologies need to be selected based on the desired outcome, such as the number of features to be included in the prediction model. It was also demonstrated that without proper model validation, for example using nested cross-validation, the models can result in overly-optimistic prediction accuracies and decreased generalization ability. It is through the implementation and application of machine learning methods that one can extract predictive genotype–phenotype relationships and biological insights from genetic data sets.