Computational models for and from biology : simple gene assembly and reaction systems

The advancement of science and technology makes it clear that no single perspective is any longer sufficient to describe the true nature of any phenomenon. That is why the interdisciplinary research is gaining more attention overtime. An excellent example of this type of research is natural computing which stands on the borderline between biology and computer science. The contribution of research done in natural computing is twofold: on one hand, it sheds light into how nature works and how it processes information and, on the other hand, it provides some guidelines on how to design bio-inspired technologies. The first direction in this thesis focuses on a nature-inspired process called gene assembly in ciliates. The second one studies reaction systems, as a modeling framework with its rationale built upon the biochemical interactions happening within a cell. The process of gene assembly in ciliates has attracted a lot of attention as a research topic in the past 15 years. Two main modelling frameworks have been initially proposed in the end of 1990s to capture ciliates’ gene assembly process, namely the intermolecular model and the intramolecular model. They were followed by other model proposals such as templatebased assembly and DNA rearrangement pathways recombination models. In this thesis we are interested in a variation of the intramolecular model called simple gene assembly model, which focuses on the simplest possible folds in the assembly process. We propose a new framework called directed overlap-inclusion (DOI) graphs to overcome the limitations that previously introduced models faced in capturing all the combinatorial details of the simple gene assembly process. We investigate a number of combinatorial properties of these graphs, including a necessary property in terms of forbidden induced subgraphs. We also introduce DOI graph-based rewriting rules that capture all the operations of the simple gene assembly model and prove that they are equivalent to the string-based formalization of the model. Reaction systems (RS) is another nature-inspired modeling framework that is studied in this thesis. Reaction systems’ rationale is based upon two main regulation mechanisms, facilitation and inhibition, which control the interactions between biochemical reactions. Reaction systems is a complementary modeling framework to traditional quantitative frameworks, focusing on explicit cause-effect relationships between reactions. The explicit formulation of facilitation and inhibition mechanisms behind reactions, as well as the focus on interactions between reactions (rather than dynamics of concentrations) makes their applicability potentially wide and useful beyond biological case studies. In this thesis, we construct a reaction system model corresponding to the heat shock response mechanism based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We also introduce for RS various concepts inspired by biology, e.g., mass conservation, steady state, periodicity, etc., to do model checking of the reaction systems based models. We prove that the complexity of the decision problems related to these properties varies from P to NP- and coNP-complete to PSPACE-complete. We further focus on the mass conservation relation in an RS and introduce the conservation dependency graph to capture the relation between the species and also propose an algorithm to list the conserved sets of a given reaction system.

The Dual Role of HIF Hydroxylase PHD3 in Cancer Cell Survival

Most advanced tumours face periods of reduced oxygen availability i.e. hypoxia. During these periods tumour cells undergo adaptive changes enabling their survival under adverse conditions. In cancer hypoxia-induced cellular changes cause tumour progression, hinder cancer treatment and are indicative of poor prognosis. Within cells the main regulator of hypoxic responses is the hypoxia-inducible factor (HIF). HIF governs the expression of over a hundred hypoxia-inducible genes that regulate a number of cellular functions such as angiogenesis, glucose metabolism and cell migration. Therefore the activity of HIF must be tightly governed. HIF is regulated by a family of prolyl hydroxylase enzymes, PHDs, which mark HIF for destruction in normoxia. Under hypoxic conditions PHDs lose much of their enzymatic activity as they need molecular oxygen as a cofactor. Out of the three PHDs (PHD1, 2 and 3) PHD2 has been considered to be the main HIF-1 regulator in normoxic conditions. PHD3 on the other hand shows the most robust induction in response to oxygen deprivation and it has been implied as the main HIF-1 regulator under prolonged hypoxia. SQSTM1/p62 (p62) is an adaptor protein that functions through its binding motifs to bring together proteins in order to regulate signal transduction. In non-stressed situations p62 levels are kept low but its expression has been reported to be upregulated in many cancers. It has a definitive role as an autophagy receptor and as such it serves a key function in cancer cell survival decisions. In my thesis work I evaluated the significance of PHD3 in cancer cell and tumour biology. My results revealed that PHD3 has a dual role in cancer cell fate. First, I demonstrated that PHD3 forms subcellular protein aggregates in oxygenated carcinoma cells and that this aggregation promotes apoptosis induction in a subset of cancer cells. In these aggregates an adaptor protein SQSTM1/p62 interacts with PHD3 and in so doing regulates PHD3 expression. SQSTM1/p62 expression is needed to keep PHD3 levels low in normoxic conditions. Its levels rapidly decrease in response to hypoxia allowing PHD3 protein levels to be upregulated and the protein to be diffusely expressed throughout the cell. The interaction between PHD3 and SQSTM1/p62 limits the ability of PHD3 to function on its hydroxylation target protein HIF-1alpha. Second, the results indicate that when PHD3 is upregulated under hypoxia it protects cancer cells by allowing cell cycle to proceed from G1 to S-phase. My data demonstrates that PHD3 may either cause cell death or protect the cells depending on its expression pattern and the oxygen availability of tumours.

Pluripotency and genetic stability of human pluripotent stem cells

Pluripotent cells have the potential to differentiate into all somatic cell types. As the adult human body is unable to regenerate various tissues, pluripotent cells provide an attractive source for regenerative medicine. Human embryonic stem cells (hESCs) can be isolated from blastocyst stage embryos and cultured in the laboratory environment. However, their use in regenerative medicine is restricted due to problems with immunosuppression by the host and ethical legislation. Recently, a new source of pluripotent cells was established via the direct reprogramming of somatic cells. These human induced pluripotent stem cells (hiPSCs) enable the production of patient specific cell types. However, numerous challenges, such as efficient reprogramming, optimal culture, directed differentiation, genetic stability and tumor risk need to be solved before the launch of therapeutic applications. The main objective of this thesis was to understand the unique properties of human pluripotent stem cells. The specific aims were to identify novel factors involved in maintaining pluripotency, characterize the effects of low oxygen culture on hESCs, and determine the high resolution changes in hESCs and hiPSCs during culture and reprogramming. As a result, the previously uncharacterized protein L1TD1 was determined to be specific for pluripotent cells and essential for the maintenance of pluripotency. The low oxygen culture supported undifferentiated growth and affected expression of stem cell associated transcripts. High resolution screening of hESCs identified a number of culture induced copy number variations and loss of heterozygosity changes. Further, screening of hiPSCs revealed that reprogramming induces high resolution alterations. The results obtained in this thesis have important implications for stem cell and cancer biology and the therapeutic potential of pluripotent cells.

International adoption: Symptoms of attachment disorders and their associations with the child’s background and developmental outcome

Many internationally adopted children have lived their first years of life in an environment with limited opportunities for primary caregiving. The lack of consistent care increases the prevalence of attachment disorders among them. Less is known about the influences of attachment disorders on a child’s later course of life. This study is part of the Finnish Adoption Study. Parents of all Finnish children who had been internationally adopted by legal adoption organisations between 1985 and 2007 were sent questionnaires (N=1450). Parental evaluations of the children’s symptoms of reactive attachment disorder (RAD) at the time of adoption, their later learning or language problems using a screening scale, and children’s self-reported school bullying experiences were evaluated. Each child’s attachment-related behavioural problems were requested in a follow-up survey 1.9 and 3.8 years after adoption and compared with a Finnish reference group. This study indicated that Finnish internationally adopted children have at least three-fold prevalence of learning and language problems compared with their age-mates. A child’s symptoms of attachment disorders were associated with learning or language problems at school age as well as with his/her school bullying experiences. The adopted children had more attachment-related behavioural problems two years after adoption than their age-mates, but the difference was no longer evident four years after adoption. In conclusion, this study showed that the symptoms of attachment disorder indicate a risk for an adopted child’s later developmental outcome. The findings demonstrate the need for comprehensive clinical examinations and planning of treatment strategies for children with symptoms of RAD.

Modeling and Control of the Power Conversion Unit in a Solid Oxide Fuel Cell Environment

In this doctoral thesis, a power conversion unit for a 10 kWsolid oxide fuel cell is modeled, and a suitable control system is designed. The need for research was identified based on an observation that there was no information available about the characteristics of the solid oxide fuel cell from the perspective of power electronics and the control system, and suitable control methods had not previously been studied in the literature. In addition, because of the digital implementation of the control system, the inherent characteristics of the digital system had to be taken into account in the characteristics of the solid oxide fuel cell (SOFC). The characteristics of the solid oxide fuel cell as well the methods for the modeling and control of the DC/DC converter and the grid converter are studied by a literature survey. Based on the survey, the characteristics of the SOFC as an electrical power source are identified, and a solution to the interfacing of the SOFC in distributed generation is proposed. A mathematical model of the power conversion unit is provided, and the control design for the DC/DC converter and the grid converter is made based on the proposed interfacing solution. The limit cycling phenomenon is identified as a source of low-frequency current ripple, which is found to be insignificant when connected to a grid-tied converter. A method to mitigate a second harmonic originating from the grid interface is proposed, and practical considerations of the operation with the solid oxide fuel cell plant are presented. At the theoretical level, the thesis discusses and summarizes the methods to successfully derive a model for a DC/DC converter, a grid converter, and a power conversion unit. The results of this doctoral thesis can also be used in other applications, and the models and methods can be adopted to similar applications such as photovoltaic systems. When comparing the results with the objectives of the doctoral thesis, we may conclude that the objectives set for the work are met. In this doctoral thesis, theoretical and practical guidelines are presented for the successful control design to connect a SOFC-based distributed generation plant to the utility grid.

Regulation of self-renewal and detection of karyotypic changes of pluripotent human embryonic stem cells

Human embryonic stem cells are pluripotent cells capable of renewing themselves and differentiating to specialized cell types. Because of their unique regenerative potential, pluripotent cells offer new opportunities for disease modeling, development of regenerative therapies, and treating diseases. Before pluripotent cells can be used in any therapeutic applications, there are numerous challenges to overcome. For instance, the key regulators of pluripotency need to be clarified. In addition, long term culture of pluripotent cells is associated with the accumulation of karyotypic abnormalities, which is a concern regarding the safe use of the cells for therapeutic purposes. The goal of the work presented in this thesis was to identify new factors involved in the maintenance of pluripotency, and to further characterize molecular mechanisms of selected candidate genes. Furthermore, we aimed to set up a new method for analyzing genomic integrity of pluripotent cells. The experimental design applied in this study involved a wide range of molecular biology, genome-wide, and computational techniques to study the pluripotency of stem cells and the functions of the target genes. In collaboration with instrument and reagent company Perkin Elmer, KaryoliteTM BoBsTM was implemented for detecting karyotypic changes of pluripotent cells. Novel genes were identified that are highly and specifically expressed in hES cells. Of these genes, L1TD1 and POLR3G were chosen for further investigation. The results revealed that both of these factors are vital for the maintenance of pluripotency and self-renewal of the hESCs. KaryoliteTM BoBsTM was validated as a novel method to detect karyotypic abnormalities in pluripotent stem cells. The results presented in this thesis offer significant new information on the regulatory networks associated with pluripotency. The results will facilitate in understanding developmental and cancer biology, as well as creating stem cell based applications. KaryoliteTM BoBsTM provides rapid, high-throughput, and cost-efficient tool for screening of human pluripotent cell cultures.

Choosing the target loci : heat shock factors HSF1 and HSF2 as regulators of cell stress and development

Heat shock factors (HSFs) are an evolutionarily well conserved family of transcription factors that coordinate stress-induced gene expression and direct versatile physiological processes in eukaryote organisms. The essentiality of HSFs for cellular homeostasis has been well demonstrated, mainly through HSF1-induced transcription of heat shock protein (HSP) genes. HSFs are important regulators of many fundamental processes such as gametogenesis, metabolic control and aging, and are involved in pathological conditions including cancer progression and neurodegenerative diseases. In each of the HSF-mediated processes, however, the detailed mechanisms of HSF family members and their complete set of target genes have remained unknown. Recently, rapid advances in chromatin studies have enabled genome-wide characterization of protein binding sites in a high resolution and in an unbiased manner. In this PhD thesis, these novel methods that base on chromatin immunoprecipitation (ChIP) are utilized and the genome-wide target loci for HSF1 and HSF2 are identified in cellular stress responses and in developmental processes. The thesis and its original publications characterize the individual and shared target genes of HSF1 and HSF2, describe HSF1 as a potent transactivator, and discover HSF2 as an epigenetic regulator that coordinates gene expression throughout the cell cycle progression. In male gametogenesis, novel physiological functions for HSF1 and HSF2 are revealed and HSFs are demonstrated to control the expression of X- and Y-chromosomal multicopy genes in a silenced chromatin environment. In stressed human cells, HSF1 and HSF2 are shown to coordinate the expression of a wide variety of genes including genes for chaperone machinery, ubiquitin, regulators of cell cycle progression and signaling. These results highlight the importance of cell type and cell cycle phase in transcriptional responses, reveal the myriad of processes that are adjusted in a stressed cell and describe novel mechanisms that maintain transcriptional memory in mitotic cell division.