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.

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.

Modelling of combustion and sorbent reactions in three-dimensional flow environment of a circulating fluidized bed furnace

This thesis presents a three-dimensional, semi-empirical, steady state model for simulating the combustion, gasification, and formation of emissions in circulating fluidized bed (CFB) processes. In a large-scale CFB furnace, the local feeding of fuel, air, and other input materials, as well as the limited mixing rate of different reactants produce inhomogeneous process conditions. To simulate the real conditions, the furnace should be modelled three-dimensionally or the three-dimensional effects should be taken into account. The only available methods for simulating the large CFB furnaces three-dimensionally are semi-empirical models, which apply a relatively coarse calculation mesh and a combination of fundamental conservation equations, theoretical models and empirical correlations. The number of such models is extremely small. The main objective of this work was to achieve a model which can be applied to calculating industrial scale CFB boilers and which can simulate all the essential sub-phenomena: fluid dynamics, reactions, the attrition of particles, and heat transfer. The core of the work was to develop the model frame and the required sub-models for determining the combustion and sorbent reactions. The objective was reached, and the developed model was successfully used for studying various industrial scale CFB boilers combusting different types of fuel. The model for sorbent reactions, which includes the main reactions for calcitic limestones, was applied for studying the new possible phenomena occurring in the oxygen-fired combustion. The presented combustion and sorbent models and principles can be utilized in other model approaches as well, including other empirical and semi-empirical model approaches, and CFD based simulations. The main achievement is the overall model frame which can be utilized for the further development and testing of new sub-models and theories, and for concentrating the knowledge gathered from the experimental work carried out at bench scale, pilot scale and industrial scale apparatus, and from the computational work performed by other modelling methods.

Physiological evaluation of heat stress in gestating sows under different housing systems in bedding and concrete floor

The purpose of this paper was to observe the use of bedding (wood shavings) in physiological variables that indicate thermal stress in gestating sows. The experiment was conducted in order to evaluate the effect of two types of floor (concrete and wood shavings). Worse microclimatic conditions were observed in bedding systems (P<0.05), with an increase in temperature and enthalpy of 1.14 ºC and 2.37 kJ.kg dry air-1, respectively. The floor temperature at the dirty area was higher in the bedding presence in comparison to its absence. In spite of the worse microclimatic conditions in the bedding, the rectal temperature did not differ significantly (P>0.05) but the skin surface temperature was higher in the bedding systems. The same occurred with the respiratory rates. The physical characteristics of the floor material influenced the rate of heat loss by conductance. Estimated values were 35.04 and 7.99 W m-2 for the conductive heat loss between the animal and floor for treatments with or without bedding, respectively. The use of bedding in sow rearing has a negative impact on microclimatic conditions, what implies in thermoregulatory damages.

Heat loss in one day old pullets inside a hatchery

After hatching, pullets are transported to brooding area and vaccinated. One day old chicks have not already developed thermoregulation ability; thus, brooding temperature variations may affect pullet quality leading to broiler meat production losses. This research aimed to calculate sensible heat loss in one day old pullets in hatching area and vaccination room. Ten one day old pullets were randomly selected from hatching area of a commercial hatchery. Infrared images were used to calculate bird surface temperature. Exposure areas for the two conditions were quantified, and both air temperature and wind speed was recorded. Total sensible heat loss was calculated as heat loss by radiation plus heat loss by convection. It was found that heat transfer occurs in different ways at different bird body parts. Total heat loss found for hatching baskets was equivalent to 0.81 J s-1 while for vaccination room was 1.16 J s-1. Pullet nutrition is based on energy loss from brooding to farm, and the overall pullet heat loss from hatchery to farm accepted is 13.95 J s-1. Thus, "starter feed" has relevant excess of energy input. These findings indicate that less energy can be used in initial feed, once heat loss is lower than assumed nowadays. Improved knowledge on these conditions may enhance broiler farm feeding strategies and economics during first rearing week.

The impact of evaporative cooling on the thermoregulation and sensible heat loss of sows during farrowing

Pigs are more sensitive to high environmental temperatures explained by the inability of sweating and panting properly when compared to other species of farmed livestock. The evaporative cooling system might favor the thermal comfort of animals during exposure to extreme environmental heat and reduce the harmful effects of heat stress. The purpose of this study was to assess the sensible heat loss and thermoregulation parameters from lactating sows during summer submitted to two different acclimatization systems: natural and evaporative cooling. The experiment was carried out in a commercial farm with 72 lactating sows. The ambient variables (temperature, relative humidity and air velocity) and sows physiological parameters (rectal temperature, surface temperature and respiratory rate) were monitored and then the sensible heat loss at 21days lactation was calculated. The results of rectal temperature did not differ between treatments. However, the evaporative cooling led to a significant reduction in surface temperature and respiratory rate and a significant increase in the sow's sensible heat loss. It was concluded that the use of evaporative cooling system was essential to increase sensible heat loss; thus, it should reduce the negative effects of heat on the sows' thermoregulation during summer.

Mathematical simulation parameters for drying of cassava starch pellets

Based on experimental tests, it was obtained the equations for drying, equilibrium moisture content, latent heat of vaporization of water contained in the product and the equation of specific heat of cassava starch pellets, essential parameters for realizing modeling and mathematical simulation of mechanical drying of cassava starch for a new technique proposed, consisting of preformed by pelleting and subsequent artificial drying of starch pellets. Drying tests were conducted in an experimental chamber by varying the air temperature, relative humidity, air velocity and product load. The specific heat of starch was determined by differential scanning calorimetry. The generated equations were validated through regression analysis, finding an appropriate correlation of the data, which indicates that by using these equations, can accurately model and simulate the drying process of cassava starch pellets.

Low-flow estimates in regions of extrapolation of the regionalization equations: a new concept

ABSTRACT Knowledge of natural water availability, which is characterized by low flows, is essential for planning and management of water resources. One of the most widely used hydrological techniques to determine streamflow is regionalization, but the extrapolation of regionalization equations beyond the limits of sample data is not recommended. This paper proposes a new method for reducing overestimation errors associated with the extrapolation of regionalization equations for low flows. The method is based on the use of a threshold value for the maximum specific low flow discharge estimated at the gauging sites that are used in the regionalization. When a specific low flow, which has been estimated using the regionalization equation, exceeds the threshold value, the low flow can be obtained by multiplying the drainage area by the threshold value. This restriction imposes a physical limit to the low flow, which reduces the error of overestimating flows in regions of extrapolation. A case study was done in the Urucuia river basin, in Brazil, and the results showed the regionalization equation to perform positively in reducing the risk of extrapolation.

Fuzzy modeling of broiler performance, raised from 1 to 21 days, subject to heat stress

ABSTRACT Given the need to obtain systems to better control broiler production environment, we performed an experiment with broilers from 1 to 21 days, which were submitted to different intensities and air temperature durations in conditioned wind tunnels and the results were used for validation of afuzzy model. The model was developed using as input variables: duration of heat stress (days), dry bulb air temperature (°C) and as output variable: feed intake (g) weight gain (g) and feed conversion (g.g-1). The inference method used was Mamdani, 20 rules have been prepared and the defuzzification technique used was the Center of Gravity. A satisfactory efficiency in determining productive responses is evidenced in the results obtained in the model simulation, when compared with the experimental data, where R2 values ​​calculated for feed intake, weight gain and feed conversion were 0.998, 0.981 and 0.980, respectively.

Effect of a simulated heat wave in thermal and aerial environment broiler-rearing environment

ABSTRACT Global warming increases the occurrence of events such as extreme heat waves. Research on thermal and air conditions affecting broiler-rearing environment are important to evaluate the animal welfare under extreme heat aiming mitigation measures. This study aimed at evaluating the effect of a simulated heat wave, in a climatic chamber, on the thermal and air environment of 42-day-old broilers. One hundred and sixty broilers were housed and reared for 42 days in a climatic chamber; the animals were divided into eight pens. Heat wave simulation was performed on the 42nd day, the period of great impact and data sampling. The analyzed variables were room and litter temperatures, relative humidity, concentrations of oxygen, carbon monoxide and ammonia at each pen. These variables were assessed each two hours, starting at 8 am, simulating a day heating up to 4 pm, when it is reached the maximum temperature. By the results, we concluded that increasing room temperatures promoted a proportional raise in litter temperatures, contributing to ammonia volatilization. In addition, oxygen concentrations decreased with increasing temperatures; and the carbon monoxide was only observed at temperatures above 27.0 °C, relative humidity higher than 88.4% and litter temperatures superior to 30.3 °C.

Word Equations and Related Topics. Independence, Decidability and Characterizations

The three main topics of this work are independent systems and chains of word equations, parametric solutions of word equations on three unknowns, and unique decipherability in the monoid of regular languages. The most important result about independent systems is a new method giving an upper bound for their sizes in the case of three unknowns. The bound depends on the length of the shortest equation. This result has generalizations for decreasing chains and for more than three unknowns. The method also leads to shorter proofs and generalizations of some old results. Hmelevksii’s theorem states that every word equation on three unknowns has a parametric solution. We give a significantly simplified proof for this theorem. As a new result we estimate the lengths of parametric solutions and get a bound for the length of the minimal nontrivial solution and for the complexity of deciding whether such a solution exists. The unique decipherability problem asks whether given elements of some monoid form a code, that is, whether they satisfy a nontrivial equation. We give characterizations for when a collection of unary regular languages is a code. We also prove that it is undecidable whether a collection of binary regular languages is a code.

Laskentasovelluksen kehittäminen modulaariselle vesiputkikattilalle

Modulaarisella vesiputkikattilalla tarkoitetaan täysin konepajalla valmistettavissa olevaa kattilaa, joka voidaan kuljettaa yhtenä tai muutamana suurena moduulina työmaalle. Tässä diplomityössä käsiteltiin modulaarisen vesiputkikattilan laskentasovelluksen ke-hittämistä KPA Unicon Oy:lle. Työn tavoitteena oli tarkastella modulaarisen vesiputkikattilan lämpöteknistä mitoitusta ja suunnittelua sekä kehittää laskentasovellus, jonka avulla voidaan arvioida kattilan mittoja ja painoa. Laskentasovellus laadittiin Microsoftin Excel-alustalle, josta se on myöhemmin mahdollista siirtää muille ohjelma-alustoille. Laskentasovelluksessa hyödynnetään lämmönsiirron ja virtaustekniikan laskentayhtälöitä sekä -menetelmiä. Sovellukseen valitut laskentayhtälöt sekä -menetelmät ovat yleisesti hyväksyttyjä ja käytännössä testattuja. Diplomityön tuloksena valmistui laskentasovellus, joka kykenee modulaarisen vesiputkikattilan lämpötekniseen mitoittamiseen. Sovelluksen avulla voidaan mitoittaa kattilan tulipesä, tulistimet, höyrystinpinnat sekä ekonomaiseri. Laskentasovellusta on tarkoitus hyödyntää yrityksen tarjousvaiheen projekteissa sekä mahdollisesti kattiloiden esisuunnittelussa. Laskentasovelluksen laatimista varten ei ollut mahdollista hyödyntää toiminnassa olevien kattiloiden prosessitietoja, koska ensimmäiset tämän tyyppiset kattilat ovat asenteilla. Sen sijaan sovelluksen antamia laskentatuloksia verrattiin toisen mitoitusohjelman antamiin tuloksiin, joiden perusteella laskentasovelluksen voidaan olettaa antavan oikeita tuloksia.

Fundamentals of heat transfer

Following over 170+ pages and additional appendixes are formed based on content of Course: Fundamentals of Heat Transfer. Mainly this summarizes relevant parts on Book of Fundamentals of Heat and Mass Transfer (Incropera), but also other references introducing the same concepts are included. Student’s point of view has been consideredwith following highlights: (1) Relevant topics are presented in a nutshell to provide fast digestion of principles of heat transfer. (2) Appendixes include terminology dictionary. (3) Totally 22 illustrating examples are connecting theory to practical applications and quantifying heat transfer to understandable forms as: temperatures, heat transfer rates, heat fluxes, resistances and etc. (4) Most important Learning outcomes are presented for each topic separately. The Book, Fundamentals of Heat and Mass Transfer (Incropera), is certainly recommended for those going beyond basic knowledge of heat transfer. Lecture Notes consists of four primary content-wise objectives: (1) Give understanding to physical mechanisms of heat transfer, (2)Present basic concepts and terminology relevant for conduction, convection and radiation (3) Introduce thermal performance analysis methods for steady state and transient conduction systems. (4) Provide fast-to-digest phenomenological understanding required for basic design of thermal models

Validation of a 3D model for oxygen combustion in a circulating fluidized bed

In this thesis, a model called CFB3D is validated for oxygen combustion in circulating fluidized bed boiler. The first part of the work consists of literature review in which circulating fluidized bed and oxygen combustion technologies are studied. In addition, the modeling of circulating fluidized bed furnaces is discussed and currently available industrial scale three-dimensional furnace models are presented. The main features of CFB3D model are presented along with the theories and equations related to the model parameters used in this work. The second part of this work consists of the actual research and modeling work including measurements, model setup, and modeling results. The objectives of this thesis is to study how well CFB3D model works with oxygen combustion compared to air combustion in circulating fluidized bed boiler and what model parameters need to be adjusted when changing from air to oxygen combustion. The study is performed by modeling two air combustion cases and two oxygen combustion cases with comparable boiler loads. The cases are measured at Ciuden 30 MWth Flexi-Burn demonstration plant in April 2012. The modeled furnace temperatures match with the measurements as well in oxygen combustion cases as in air combustion cases but the modeled gas concentrations differ from the measurements clearly more in oxygen combustion cases. However, the same model parameters are optimal for both air and oxygen combustion cases. When the boiler load is changed, some combustion and heat transfer related model parameters need to be adjusted. To improve the accuracy of modeling results, better flow dynamics model should be developed in the CFB3D model. Additionally, more measurements are needed from the lower furnace to find the best model parameters for each case. The validation work needs to be continued in order to improve the modeling results and model predictability.

Limestone Reactions in a Fluidized Bed Heat Exchanger of an Oxy-fuel Circulating Fluidized Bed Boiler

Oxy-fuel combustion in a circulating fluidized bed (CFB) boiler appears to be a promising option for capturing CO2 in power plants. Oxy-fuel combustion is based on burning of fuel in the mixture of oxygen and re-circulated flue gas instead of air. Limestone (CaCO3) is typically used for capturing of SO2 in CFB boilers where limestone calcines to calcium oxide (CaO). Because of high CO2 concentration in oxy-fuel combustion, calcination reaction may be hindered or carbonation, the reverse reaction of calcination, may occur. Carbonation of CaO particles can cause problems especially in the circulation loop of a CFB boiler where temperature level is lower than in the furnace. The aim of the thesis was to examine carbonation of CaO in a fluidized bed heat exchanger of a CFB boiler featuring oxy-fuel combustion. The calculations and analyzing were based on measurement data from an oxy-fuel pilot plant and on 0-dimensional (0D) gas balance of a fluidized bed heat exchanger. Additionally, the objective was to develop a 1-dimensional (1D) model of a fluidized bed heat exchanger by searching a suitable pre-exponential factor for a carbonation rate constant. On the basis of gas measurement data and the 0D gas balance, it was found that the amount of fluidization gas decreased as it flew through the fluidized bed heat exchanger. Most likely the reason for this was carbonation of CaO. It was discovered that temperature has a promoting effect on the reaction rate of carbonation. With the 1D model, a suitable pre-exponential factor for the equation of carbonation rate constant was found. However, during measurements there were several uncertainties, and in the calculations plenty of assumptions were made. Besides, the temperature level in the fluidized bed heat exchanger was relatively low during the measurements. Carbonation should be considered when fluidized bed heat exchangers and the capacity of related fans are designed for a CFB boiler with oxy-fuel combustion.