Role and regulatory mechanisms of heat shock factor 2

Protein homeostasis is essential for cells to prosper and survive. Various forms of stress, such as elevated temperatures, oxidative stress, heavy metals or bacterial infections cause protein damage, which might lead to improper folding and formation of toxic protein aggregates. Protein aggregation is associated with serious pathological conditions such as Alzheimer’s and Huntington’s disease. The heat shock response is a defense mechanism that protects the cell against protein-damaging stress. Its ancient origin and high conservation among eukaryotes suggest that the response is crucial for survival. The main regulator of the heat shock response is the transcription factor heat shock factor 1 (HSF1), which induces transcription of genes encoding protective molecular chaperones. In vertebrates, a family of four HSFs exists (HSF1-4), with versatile functions not only in coping with acute stress, but also in development, longevity and cancer. Thus, knowledge of the HSFs will aid in our understanding on how cells survive suboptimal circumstances, but will also provide insights into normal physiological processes as well as diseaseassociated conditions. In this study, the function and regulation of HSF2 have been investigated. Earlier gene inactivation experiments in mice have revealed roles for HSF2 in development, particularly in corticogenesis and spermatogenesis. Here, we demonstrate that HSF2 holds a role also in the heat shock response and influences stress-induced expression of heat shock proteins. Intriguingly, DNA-binding activity of HSF2 upon stress was dependent on the presence of intact HSF1, suggesting functional interplay between HSF1 and HSF2. The underlying mechanism for this phenomenon could be configuration of heterotrimers between the two factors, a possibility that was experimentally verified. By changing the levels of HSF2, the expression of HSF1-HSF2 heterotrimer target genes was altered, implementing HSF2 as a modulator of HSF-mediated transcription. The results further indicate that HSF2 activity is dependent on its concentration, which led us to ask the question of how accurate HSF2 levels are achieved. Using mouse spermatogenesis as a model system, HSF2 was found to be under direct control of miR-18, a miRNA belonging to the miR-17~92 cluster/Oncomir-1 and whose physiological function had remained unclear. Investigations on spermatogenesis are severely hampered by the lack of cell systems that would mimic the complex differentiation processes that constitute male germ cell development. Therefore, to verify that HSF2 is regulated by miR-18 in spermatogenesis, a novel method named T-GIST (Transfection of Germ cells in Intact Seminiferous Tubules) was developed. Employing this method, the functional consequences of miR-18-mediated regulation in vivo were demonstrated; inhibition of miR- 18 led to increased expression of HSF2 and altered the expression of HSF2 target genes Ssty2 and Speer4a. Consequently, the results link miR-18 to HSF2-mediated processes such as germ cell maturation and quality control and provide miR-18 with a physiological role in gene expression during spermatogenesis.Taken together, this study presents compelling evidence that HSF2 is a transcriptional regulator in the heat shock response and establishes the concept of physical interplay between HSF2 and HSF1 and functional consequences thereof. This is also the first study describing miRNA-mediated regulation of an HSF.

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.

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.

Regulation of HSF2 and its function in mitosis

The cell is continuously subjected to various forms of external and intrinsic proteindamaging stresses, including hyperthermia, pathophysiological states, as well as cell differentiation and proliferation. Proteindamaging stresses result in denaturation and improper folding of proteins, leading to the formation of toxic aggregates that are detrimental for various pathological conditions, including Alzheimer’s and Huntington’s diseases. In order to maintain protein homeostasis, cells have developed different cytoprotective mechanisms, one of which is the evolutionary well-conserved heat shock response. The heat shock response results in the expression of heat shock proteins (Hsps), which act as molecular chaperones that bind to misfolded proteins, facilitate their refolding and prevent the formation of protein aggregates. Stress-induced expression of Hsps is mediated by a family of transcription factors, the heat shock factors, HSFs. Of the four HSFs found in vertebrates, HSF1-4, HSF1 is the major stress-responsive factor that is required for the induction of the heat shock response. HSF2 cannot alone induce Hsps, but modulates the heat shock response by forming heterotrimers with HSF1. HSFs are not only involved in the heat shock response, but they have also been found to have a function in development, neurodegenerative disorders, cancer, and longevity. Therefore, insight into how HSFs are regulated is important for the understanding of both normal physiological and disease processes. The activity of HSF1 is mainly regulated by intricate post-translational modifications, whereas the activity of HSF2 is concentrationdependent. However, there is only limited understanding of how the abundance of HSF2 is regulated. This study describes two different means of how HSF2 levels are regulated. In the first study it was shown that microRNA miR-18, a member of the miR-17~92 cluster, directly regulates Hsf2 mRNA stability and thus protein levels. HSF2 has earlier been shown to play a profound role in the regulation of male germ cell maturation during the spermatogenesis. The effect on miR-18 on HSF2 was examined in vivo by transfecting intact seminiferous tubules, and it was found that inhibition of miR-18 resulted in increased HSF2 levels and modified expression of the HSF2 targets Ssty2 and Speer4a. HSF2 has earlier been reported to modulate the heat shock response by forming heterotrimers with HSF1. In the second study, it was shown that HSF2 is cleared off the Hsp70 promoter and degraded by the ubiquitinproteasome pathway upon acute stress. By silencing components of the anaphase promoting complex/cyclosome (APC/C), including the co-activators Cdc20 and Cdh1, it was shown that APC/C mediates the heatinduced ubiquitylation of HSF2. Furthermore, down-regulation of Cdc20 was shown to alter the expression of heat shock-responsive genes. Next, we studied if APC/C-Cdc20, which controls cell cycle progression, also regulates HSF2 during the cell cycle. We found that both HSF2 mRNA and protein levels decreased during mitosis in several but not all human cell lines, indicating that HSF2 has a function in mitotic cells. Interestingly, although transcription is globally repressed during mitosis, mainly due to the displacement of RNA polymerase II and transcription factors, including HSF1, from the mitotic chromatin, HSF2 is capable of binding DNA during mitosis. Thus, during mitosis the heat shock response is impaired, leaving mitotic cells vulnerable to proteotoxic stress. However, in HSF2-deficient mitotic cells the Hsp70 promoter is accessible to both HSF1 and RNA polymerase II, allowing for stress-inducible Hsp expression to occur. As a consequence HSF2-deficient mitotic cells have a survival advantage upon acute heat stress. The results, presented in this thesis contribute to the understanding of the regulatory mechanisms of HSF2 and its function in the heat shock response in both interphase and mitotic cells.

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.

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.

Isentropic Exergy and Pressure of the Shock Wave Caused by the Explosion of a Pressure Vessel

An accidental burst of a pressure vessel is an uncontrollable and explosion-like batch process. In this study it is called an explosion. The destructive effectof a pressure vessel explosion is relative to the amount of energy released in it. However, in the field of pressure vessel safety, a mutual understanding concerning the definition of explosion energy has not yet been achieved. In this study the definition of isentropic exergy is presented. Isentropic exergy is the greatest possible destructive energy which can be obtained from a pressure vessel explosion when its state changes in an isentropic way from the initial to the final state. Finally, after the change process, the gas has similar pressure and flow velocity as the environment. Isentropic exergy differs from common exergy inthat the process is assumed to be isentropic and the final gas temperature usually differs from the ambient temperature. The explosion process is so fast that there is no time for the significant heat exchange needed for the common exergy.Therefore an explosion is better characterized by isentropic exergy. Isentropicexergy is a characteristic of a pressure vessel and it is simple to calculate. Isentropic exergy can be defined also for any thermodynamic system, such as the shock wave system developing around an exploding pressure vessel. At the beginning of the explosion process the shock wave system has the same isentropic exergyas the pressure vessel. When the system expands to the environment, its isentropic exergy decreases because of the increase of entropy in the shock wave. The shock wave system contains the pressure vessel gas and a growing amount of ambient gas. The destructive effect of the shock wave on the ambient structures decreases when its distance from the starting point increases. This arises firstly from the fact that the shock wave system is distributed to a larger space. Secondly, the increase of entropy in the shock waves reduces the amount of isentropic exergy. Equations concerning the change of isentropic exergy in shock waves are derived. By means of isentropic exergy and the known flow theories, equations illustrating the pressure of the shock wave as a function of distance are derived. Amethod is proposed as an application of the equations. The method is applicablefor all shapes of pressure vessels in general use, such as spheres, cylinders and tubes. The results of this method are compared to measurements made by various researchers and to accident reports on pressure vessel explosions. The test measurements are found to be analogous with the proposed method and the findings in the accident reports are not controversial to it.