Methods for construction and analysis of computational models in systems biology : applications to the modelling of the heat shock response and the self-assembly of intermediate filaments

Systems biology is a new, emerging and rapidly developing, multidisciplinary research field that aims to study biochemical and biological systems from a holistic perspective, with the goal of providing a comprehensive, system- level understanding of cellular behaviour. In this way, it addresses one of the greatest challenges faced by contemporary biology, which is to compre- hend the function of complex biological systems. Systems biology combines various methods that originate from scientific disciplines such as molecu- lar biology, chemistry, engineering sciences, mathematics, computer science and systems theory. Systems biology, unlike “traditional” biology, focuses on high-level concepts such as: network, component, robustness, efficiency, control, regulation, hierarchical design, synchronization, concurrency, and many others. The very terminology of systems biology is “foreign” to “tra- ditional” biology, marks its drastic shift in the research paradigm and it indicates close linkage of systems biology to computer science. One of the basic tools utilized in systems biology is the mathematical modelling of life processes tightly linked to experimental practice. The stud- ies contained in this thesis revolve around a number of challenges commonly encountered in the computational modelling in systems biology. The re- search comprises of the development and application of a broad range of methods originating in the fields of computer science and mathematics for construction and analysis of computational models in systems biology. In particular, the performed research is setup in the context of two biolog- ical phenomena chosen as modelling case studies: 1) the eukaryotic heat shock response and 2) the in vitro self-assembly of intermediate filaments, one of the main constituents of the cytoskeleton. The range of presented approaches spans from heuristic, through numerical and statistical to ana- lytical methods applied in the effort to formally describe and analyse the two biological processes. We notice however, that although applied to cer- tain case studies, the presented methods are not limited to them and can be utilized in the analysis of other biological mechanisms as well as com- plex systems in general. The full range of developed and applied modelling techniques as well as model analysis methodologies constitutes a rich mod- elling framework. Moreover, the presentation of the developed methods, their application to the two case studies and the discussions concerning their potentials and limitations point to the difficulties and challenges one encounters in computational modelling of biological systems. The problems of model identifiability, model comparison, model refinement, model inte- gration and extension, choice of the proper modelling framework and level of abstraction, or the choice of the proper scope of the model run through this thesis.

Multifactor Performance Study on Complex Systems for the Benefit of Management

The primary objective is to identify the critical factors that have a natural impact on the performance measurement system. It is important to make correct decisions related to measurement systems, which are based on the complex business environment. The performance measurement system is combined with a very complex non-linear factor. The Six Sigma methodology is seen as one potential approach at every organisational level. It will be linked to the performance and financial measurement as well as to the analytical thinking on which the viewpoint of management depends. The complex systems are connected to the customer relationship study. As the primary throughput can be seen in a new well-defined performance measurement structure that will also be facilitated as will an analytical multifactor system. These critical factors should also be seen as a business innovation opportunity at the same time. This master's thesis has been divided into two different theoretical parts. The empirical part consists of both action-oriented and constructive research approaches with an empirical case study. The secondary objective is to seek a competitive advantage factor with a new analytical tool and the Six Sigma thinking. Process and product capabilities will be linked to the contribution of complex system. These critical barriers will be identified by the performance measuring system. The secondary throughput can be recognised as the product and the process cost efficiencies which throughputs are achieved with an advantage of management. The performance measurement potential is related to the different productivity analysis. Productivity can be seen as one essential part of the competitive advantage factor.

Autologous Stem Cell Transplantation In Multiple Myeloma

Background. Multiple myeloma (MM) is the second most common hematologic malignancy after lymphomas In Finland: the annual incidence of MM is approximately 200. For three decades the median survival remained at 3 to 4 years from diagnosis until high-dose melphalan treatment supported by autologous stem cell transplantation (ASCT) became the standard of care for newly diagnosed MM since the mid 1990’s and the median survival increased to 5 – 6 years. This study focuses on three important aspects of ASCT, namely 1) stem cell mobilization, 2) single vs. double ASCT as initial treatment, and 3) the role of minimal residual disease (MRD) for longterm outcome. Aim. The aim of this series of studies was to evaluate the outcomes of MM patients and the ASCT procedure at the Turku University Central Hospital, Finland. First, we tried to identify which factors predict unsuccessful mobilization of autologous stem cells. Second, we compared the use of short-acting granulocyte-colony stimulating factor (GCSF) with long-acting G-CSF as mobilization agents. Third, one and two successive ASCTs were compared in 100 patients with MM. Fourth, for patients in complete response (CR) after stem cell transplantation (SCT), patient-specific probes for quantitative allele-specific oligonucleotide polymerase-chain reaction (qASO-PCR) measurements were designed to evaluate MRD and its importance for long-term outcome. Results. The quantity of previous chemotherapy and previous interferon use were significant pre-mobilization factors that predicted mobilization failure, together with some factors related to mobilization therapy itself, such as duration and degree of cytopenias and occurrence of sepsis. Short-acting and long-acting G-CSF combined with chemotherapy were comparable as stem cells mobilizers. The progression free (PFS) and overall survival (OS) tended to be longer after double ASCT than after single ASCT with a median follow-up time of 4 years, but this difference disappeared as the follow-up time increased. qASO-PCR was a good and sensitive divider of the CR patients into two prognostic groups: MRD low/negative (≤ 0.01%) and MRD high (>0.01%) groups with a significant difference in PFS and suggestively also in OS. Conclusions. When the factors prediciting a poor outcome of stem cell mobilization prevail, it is possible to identify those patients who need specific efforts to maximize the mobilization efficacy. Long-acting pegfilgrastim is a practical and effective alternative to short-acting filgrastim for mobilization therapy. There is no need to perform double ASCT on all eligible patients. MRD assessment with qASO-PCR is a sensitive method for evaluation of the depth of the CR response and can be used to predict long-term outcome after ACST.

Fraud detection in the banking sector : a multi-agent approach

Fraud is an increasing phenomenon as shown in many surveys carried out by leading international consulting companies in the last years. Despite the evolution of electronic payments and hacking techniques there is still a strong human component in fraud schemes. Conflict of interest in particular is the main contributing factor to the success of internal fraud. In such cases anomaly detection tools are not always the best instruments, since the fraud schemes are based on faking documents in a context dominated by lack of controls, and the perpetrators are those ones who should control possible irregularities. In the banking sector audit team experts can count only on their experience, whistle blowing and the reports sent by their inspectors. The Fraud Interactive Decision Expert System (FIDES), which is the core of this research, is a multi-agent system built to support auditors in evaluating suspicious behaviours and to speed up the evaluation process in order to detect or prevent fraud schemes. The system combines Think-map, Delphi method and Attack trees and it has been built around audit team experts and their needs. The output of FIDES is an attack tree, a tree-based diagram to ”systematically categorize the different ways in which a system can be attacked”. Once the attack tree is built, auditors can choose the path they perceive as more suitable and decide whether or not to start the investigation. The system is meant for use in the future to retrieve old cases in order to match them with new ones and find similarities. The retrieving features of the system will be useful to simplify the risk management phase, since similar countermeasures adopted for past cases might be useful for present ones. Even though FIDES has been built with the banking sector in mind, it can be applied in all those organisations, like insurance companies or public organizations, where anti-fraud activity is based on a central anti-fraud unit and a reporting system.

Agent-augmented process automation system in a panel board factory

Kiristyvä kansainvälinen kilpailu pakottaa automaatiojärjestelmien valmistajat ottamaan käyttöön uusia menetelmiä, joiden avulla järjestelmien suorituskykyä ja joustavuutta saadaan parannettua. Agenttiteknologiaa on esitetty käytettäväksi olemassa olevien automaatiojärjestelmien kanssa vastaamaan automaatiolle asetettaviin uusiin haasteisiin. Agentit ovat itsenäisiä yhteisöllisiä toimijoita, jotka suorittavat niille ennalta määrättyjä tehtäviä. Ne tarjoavat yhtenäisen kehyksen kehittyneiden toimintojen toteutukselle. Agenttiteknologian avulla automaatiojärjestelmä saadaan toimimaan joustavasti ja vikasietoisesti. Tässä työssä selostetaan agenttiteknologian ajatuksia ja käsitteitä. Lisäksi selvitetään sen soveltuvuutta monimutkaisten ohjausjärjestelmien kehittämiseen ja etsitään käyttökohteita sen soveltamiselle levytehtaassa. Työssä käsitellään myös aatteita, jotka ovat johtaneet agenttiteknologian käyttöön automaatiojärjestelmissä, sekä selostetaan agenttiavusteisen esimerkkisovelluksen rakenne ja testitulokset. Tutkimuksen tuloksena löydettiin useita kohteita agenttiteknologian käytölle levytehtaassa. Esimerkkisovellus osoittaa sen sopivan hyvin kehittyneiden toimintojen toteutukseen automaatiojärjestelmissä.

Genetic regulatory networks in avian B cells

Biology is turning into an information science. The science of systems biology seeks to understand the genetic networks that govern organism development and functions. In this study the chicken was used as a model organism in the study of B cell regulatory factors. These studies open new avenues for plasma cell research by connecting the down regulation of the B cell gene expression program directly to the initiation of plasma cell differentiation. The unique advantages of the DT40 avian B cell model system, specifically its high homologous recombination rate, were utilized to study gene regulation in Pax5 knock out cell lines and to gain new insights into the B cell to plasma cell transitions that underlie the secretion of antibodies as part of the adaptive immune response. The Pax5 transcription factor is central to the commitment, development and maintenance of the B cell phenotype. Mice lacking the Pax5 gene have an arrest in development at the pro-B lymphocyte stage while DT40 cells have been derived from cells at a more mature stage of development. The DT40 Pax5-/- cells exhibited gene expression similarities with primary chicken plasma cells. The expression of the plasma cell transcription factors Blimp-1 and XBP-1 were significantly upregulated while the expression of the germinal centre factor BCL6 was diminished in Pax5-/- cells, and this alteration was normalized by Pax5 re-introduction. The Pax5-deficient cells further manifested substantially elevated secretion of IgM into the supernatant, another characteristic of plasma cells. These results for the first time indicated that the downregulation of the Pax5 gene in B cells promotes plasma cell differentiation. Cross-species meta-analysis of chicken and mouse Pax5 gene knockout studies uncovers genes and pathways whose regulatory relationship to Pax5 has remained unchanged for over 300 million years. Restriction of the hematopoietic stem cell fate to produce T, B and NK cell lineages is dependent on the Ikaros and its molecular partners, the closely related Helios and Aiolos. Ikaros family members are zinc finger proteins which act as transcriptional repressors while helping to activate lymphoid genes. Helios in mice is expressed from the hematopoietic stem cell level onwards, although later in development its expression seems to predominate in the T cell lineage. This study establishes the emergence and sequence of the chicken Ikaros family members. Helios expression in the bursa of Fabricius, germinal centres and B cell lines suggested a role for Helios in the avian B-cell lineage, too. Phylogenetic studies of the Ikaros family connect the expansion of the Ikaros family, and thus possibly the emergence of the adaptive immune system, with the second round of genome duplications originally proposed by Ohno. Paralogs that have arisen as a result of genome-wide duplications are sometimes termed ohnologs – Ikaros family proteins appear to fit that definition. This study highlighted the opportunities afforded by the genome sequencing efforts and somatic cell reverse genetics approaches using the DT40 cell line. The DT40 cell line and the avian model system promise to remain a fruitful model for mechanistic insight in the post-genomic era as well.

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.

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.

Quantitative refinement of reaction-based biomodels

In the field of molecular biology, scientists adopted for decades a reductionist perspective in their inquiries, being predominantly concerned with the intricate mechanistic details of subcellular regulatory systems. However, integrative thinking was still applied at a smaller scale in molecular biology to understand the underlying processes of cellular behaviour for at least half a century. It was not until the genomic revolution at the end of the previous century that we required model building to account for systemic properties of cellular activity. Our system-level understanding of cellular function is to this day hindered by drastic limitations in our capability of predicting cellular behaviour to reflect system dynamics and system structures. To this end, systems biology aims for a system-level understanding of functional intraand inter-cellular activity. Modern biology brings about a high volume of data, whose comprehension we cannot even aim for in the absence of computational support. Computational modelling, hence, bridges modern biology to computer science, enabling a number of assets, which prove to be invaluable in the analysis of complex biological systems, such as: a rigorous characterization of the system structure, simulation techniques, perturbations analysis, etc. Computational biomodels augmented in size considerably in the past years, major contributions being made towards the simulation and analysis of large-scale models, starting with signalling pathways and culminating with whole-cell models, tissue-level models, organ models and full-scale patient models. The simulation and analysis of models of such complexity very often requires, in fact, the integration of various sub-models, entwined at different levels of resolution and whose organization spans over several levels of hierarchy. This thesis revolves around the concept of quantitative model refinement in relation to the process of model building in computational systems biology. The thesis proposes a sound computational framework for the stepwise augmentation of a biomodel. One starts with an abstract, high-level representation of a biological phenomenon, which is materialised into an initial model that is validated against a set of existing data. Consequently, the model is refined to include more details regarding its species and/or reactions. The framework is employed in the development of two models, one for the heat shock response in eukaryotes and the second for the ErbB signalling pathway. The thesis spans over several formalisms used in computational systems biology, inherently quantitative: reaction-network models, rule-based models and Petri net models, as well as a recent formalism intrinsically qualitative: reaction systems. The choice of modelling formalism is, however, determined by the nature of the question the modeler aims to answer. Quantitative model refinement turns out to be not only essential in the model development cycle, but also beneficial for the compilation of large-scale models, whose development requires the integration of several sub-models across various levels of resolution and underlying formal representations.

Targeting oncogenic signaling proteins : new insights using a FRET-based chemical biology approach

Cancer affects more than 20 million people each year and this rate is increasing globally. The Ras/MAPK-pathway is one of the best-studied cancer signaling pathways. Ras proteins are mutated in almost 20% of all human cancers and despite numerous efforts, no effective therapy that specifically targets Ras is available to date. It is now well established that Ras proteins laterally segregate on the plasma membrane into transient nanoscale signaling complexes called nanoclusters. These Ras nanoclusters are essential for the high-fidelity signal transmission. Disruption of nanoclustering leads to reduction in Ras activity and signaling, therefore targeting nanoclusters opens up important new therapeutic possibilities in cancer. This work describes three different studies exploring the idea of membrane protein nanoclusters as novel anti-cancer drug targets. It is focused on the design and implementation of a simple, cell-based Förster Resonance Energy Transfer (FRET)-biosensor screening platform to identify compounds that affect Ras membrane organization and nanoclustering. Chemical libraries from different sources were tested and a number of potential hit molecules were validated on full-length oncogenic proteins using a combination of imaging, biochemical and transformation assays. In the first study, a small chemical library was screened using H-ras derived FRET-biosensors. Surprisingly from this screen, commonly used protein synthesis inhibitors (PSIs) were found to specifically increase H-ras nanoclustering and downstream signalling in a H-ras dependent manner. Using a representative PSI, increase in H-ras activity was shown to induce cancer stem cell (CSC)-enriched mammosphere formation and tumor growth of breast cancer cells. Moreover, PSIs do not increase K-ras nanoclustering, making this screening approach suitable for identifying Ras isoform-specific inhibitors. In the second study, a nanoncluster-directed screen using both H- and K-ras derived FRET biosensors identified CSC inhibitor salinomycin to specifically inhibit K-ras nanocluster organization and downstream signaling. A K-ras nanoclusteringassociated gene signature was established that predicts the drug sensitivity of cancer cells to CSC inhibitors. Interestingly, almost 8% of patient tumor samples in the The Cancer Genome Atlas (TCGA) database had the above gene signature and were associated with a significantly higher mortality. From this mechanistic insight, an additional microbial metabolite screen on H- and K-ras biosensors identified ophiobolin A and conglobatin A to specifically affect K-ras nanoclustering and to act as potential breast CSC inhibitors. In the third study, the Ras FRET-biosensor principle was used to investigate membrane anchorage and nanoclustering of myristoylated proteins such as heterotrimeric G-proteins, Yes- and Src-kinases. Furthermore, Yes-biosensor was validated to be a suitable platform for performing chemical and genetic screens to identify myristoylation inhibitors. The results of this thesis demonstrate the potential of the Ras-derived FRETbiosensor platform to differentiate and identify Ras-isoform specfic inhibitors. The results also highlight that most of the inhibitors identified predominantly perturb Ras subcellular distribution and membrane organization through some novel and yet unknown mechanisms. The results give new insights into the role of Ras nanoclusters as promising new molecular targets in cancer and in stem cells.