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.

Orthogonal function method in modelling paper web quality variation

This work is devoted to the analysis of signal variation of the Cross-Direction and Machine-Direction measurements from paper web. The data that we possess comes from the real paper machine. Goal of the work is to reconstruct the basis weight structure of the paper and to predict its behaviour to the future. The resulting synthetic data is needed for simulation of paper web. The main idea that we used for describing the basis weight variation in the Cross-Direction is Empirical Orthogonal Functions (EOF) algorithm, which is closely related to Principal Component Analysis (PCA) method. Signal forecasting in time is based on Time-Series analysis. Two principal mathematical procedures that we used in the work are Autoregressive-Moving Average (ARMA) modelling and Ornstein–Uhlenbeck (OU) process.

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.

Effects of thermal environment and supplementation levels on the physiological parameters of moxotó goats in confined and semi-confined rising systems

The objective of this work was to determine the effect of environmental variables and supplementation levels on physiological parameters of Moxotó goats in confined and semi-confined rising systems, in the Brazilian semi-arid region. The semi-confined individuals were kept on a grass based diet during the day and arrested in the end of the afternoon. The confined animals were kept in a management center, receiving two diets composed by forage cactus and maniçoba hay into two different levels (0.5 and 1.5% of the body weight). Inside the management center and in the external environment the environmental comfort parameters were set high during the afternoon period characterizing a situation of thermal discomfort for the animals. During the morning the semi-confined animals presented an average respiratory frequency (69.5 mov min-1) and rectal temperature (39.5 ºC) higher than the confined ones (62.6 mov min-1 and 39.0 ºC, respectively). The confined and semi-confined animals were able to maintain their rectal temperature within normal limits, with increase in the cardiac beatings rate and respiratory frequency. The greater percentage of the used supplementations (1.5%) seemed to increase rectal temperature in the two studied rising systems.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Thermal comfort indices in individual shelters for dairy calves with different types of roofs

Among the building materials used in rural facilities, roofs are noteworthy for being largely responsible for thermal comfort, influencing the thermal balance within the shelter. This study aimed to evaluate the influence of roof on the Enthalpy (H), Thermal Load of Radiation (TLR), and Black Globe Temperature and Humidity Index (BGHI) in individual shelters for dairy calves. The design was completely randomized with three treatments: Z - zinc tile, AC - asbestos-cement tile and ACW - asbestos-cement tile painted white on the upper side. The averages were compared by the Scott Knott test at 1% probability. The results showed no statistical difference between treatments (P<0.01) and the external environment for H. For TLR, there was statistical difference among all treatments, where ACW showed the lowest TLR, 489.28 W m-2, followed by AC with 506.72 W m-2 and Z with the highest TLR, 523.55 W m-2. For BGHI, the lowest values were observed for ACW (76.8) and AC (77.4), differing significantly from Z, which obtained the highest value (81.6). The tiles with white paint on the upper side promoted the lowest TLR and the lowest BGHI, favoring the thermal environment in the shelter.

Thermal performance and concentration of gases in facilities for pigs in semiarid region from State of Paraiba - Brazil

The study was conducted in a facility for pigs during the nursery and finishing in the town of 'Montadas', in the semiarid of the state of Paraiba, Brazil, in the rainy and dry season, aiming to evaluate the concentration of oxygen, methane, carbon monoxide and ammonia, and the bioclimatic indexes: ambient temperature (AT), relative humidity (RH) and the index of black globe temperature and humidity (IBGTH). These indexes differed significantly (P>0.05) between the periods and times. The AT in the rainy season was in the thermal comfort zone(TCZ) in most of the times in the nursery; for the finishing phase, thermal discomfort occurred; during the dry season, there was thermal comfort in the nursery phase; in the finishing phase the thermal discomfort occurred at all times. In the rainy season, the IBGTH was in TCZ; in the dry season, it was above the TCZ. The RH in the rainy period was in the TCZ; in the dry season, in most of the times, below the range of the TCZ. The concentration of gases showed no differences (P > 0.05) between periods and between the times, and the carbon monoxide, hydrogen sulfide and methane were below 1.0 ppm, and the ammonia showed a mean of 5.2 ppm. None of the analyzed gases exceeded the limits established by Brazilian and international standards for animals and workers.

Estimation of transpiration of the 'Valencia' orange young plant using thermal dissipation probe method

Most studies on measures of transpiration of plants, especially woody fruit, relies on methods of heat supply in the trunk. This study aimed to calibrate the Thermal Dissipation Probe Method (TDP) to estimate the transpiration, study the effects of natural thermal gradients and determine the relation between outside diameter and area of xylem in 'Valencia' orange young plants. TDP were installed in 40 orange plants of 15 months old, planted in boxes of 500 L, in a greenhouse. It was tested the correction of the natural thermal differences (DTN) for the estimation based on two unheated probes. The area of the conductive section was related to the outside diameter of the stem by means of polynomial regression. The equation for estimation of sap flow was calibrated having as standard lysimeter measures of a representative plant. The angular coefficient of the equation for estimating sap flow was adjusted by minimizing the absolute deviation between the sap flow and daily transpiration measured by lysimeter. Based on these results, it was concluded that the method of TDP, adjusting the original calibration and correction of the DTN, was effective in transpiration assessment.

Image analysis for assessing broiler breeder behavior response to thermal environment

The research proposes a methodology for assessing broiler breeder response to changes in rearing thermal environment. The continuous video recording of a flock analyzed may offer compelling evidences of thermal comfort, as well as other indications of welfare. An algorithm for classifying specific broiler breeder behavior was developed. Videos were recorded over three boxes where 30 breeders were reared. The boxes were mounted inside an environmental chamber were ambient temperature varied from cold to hot. Digital images were processed based on the number of pixels, according to their light intensity variation and binary contrast allowing a sequence of behaviors related to welfare. The system used the default of x, y coordinates, where x represents the horizontal distance from the top left of the work area to the point P, and y is the vertical distance. The video images were observed, and a grid was developed for identifying the area the birds stayed and the time they spent at that place. The sequence was analyzed frame by frame confronting the data with specific adopted thermal neutral rearing standards. The grid mask overlapped the real bird image. The resulting image allows the visualization of clusters, as birds in flock behave in certain patterns. An algorithm indicating the breeder response to thermal environment was developed.

Nocturnal thermal comfort in facilities for growing swines

In most of Brazilian pig farms, the environmental acclimatization systems run manually. For night and early morning periods, this practice isn't appropriate, because, in general, there are not employees available to run these manual systems. This research aimed to compare the bioclimatic profile of two differently constructed facilities to the external environment, considering the period from 6 p.m. to 6 a.m. during the spring, in order to show that night and early morning temperatures do not coincides with growing pig's thermoneutral zone. For this reason, acclimatization must be also carried out at these periods. It was analyzed the dry bulb temperature, relative air humidity, temperature-humidity index (THI) and enthalpy data of the sheds and external areas. Under the studied conditions, it was possible to conclude that the constructively appropriate shed appeared to be less influenced by the external environment, allowing better thermal control for growing pigs. Further research must be conducted to verify if automatic cooling systems is needed during night and early morning.

Thermal comfort in reduced models of broilers' houses, under different types of roofing materials

The research was developed to evaluate the use of different types of roofing materials regularly used in poultry houses. Measurements of thermal comfort were made through the use of techniques such as the Black Globe and Humidity Index (BGHI), the Thermal Heat Load (THL) and Enthalpy (H). Conducted in the State University of Goiás, during the months of April and May, 2011, the experiment was composed of five different treatments: AC - Asbestos cement tiles, BA -Bamboo tiles, BAP - Bamboo tiles painted in white, FB -Vegetable fiber tiles and bitumen, FBP -Vegetable fiber tiles and bitumen painted in white. The experiment consisted in 15 repetitions, which were considered the different days of measurements taken. Throughout the studied period, the time of the day considered the least comfortable was the one observed at 2:00pm, and the coverage of vegetable fiber and bitumen showed the highest value of BGHI (84.1) when compared to other types of coverage, characterizing a situation of lower thermal comfort, and no difference was found for THL and H on treatments in the studied region.

Thermal comfort temperature update for broiler chickens up to 21 days of age

Due to changes in genetics and nutrition, as well as in acclimatization of broiler chickens to the Brazilian climate, temperature values currently accepted as optimal may be outdated. The objective of this research was to update the environment temperatures that characterize the thermal comfort for broilers chickens from one to 21 days of age, under Brazilian production conditions. This research was conducted with 600 COBB birds, which were distributed in five growth chambers maintained at different temperatures during the first three weeks of age. During the experimental period, temperature values were progressively reduced, consisting in five treatments: T2724/21, T30/27/24, T33/30/27, T36/33/30 and T39/36/33. It was observed that the birds maintained in the T30(27-24) treatment presented better performance compared to other environment conditions. Based on the obtained regression models, the environment temperature values that provide greater weighing gain for the broiler chicken growth in the initial period were 31.3, 25.5 and 21.8 ºC, respectively for the first, second and third week of age.

Uncertainties of the 137Cs technique for validation of soil redistribution modelling in a semiarid meso-scale watershed

Scarcity of long-term series of sediment-related variables has led watershed managers to apply mathematical models to simulate sediment fluxes. Due to the high efforts for installation and maintenance of sedimentological gauges, tracers have been pointed out as an alternative to validate soil redistribution modelling. In this study, the 137Cs technique was used to assess the WASA-SED model performance at the Benguê watershed (933 km²), in the Brazilian semiarid. Qualitatively, good agreement was found among the 137Cs technique and the WASA-SED model results. Nonetheless, quantitatively great differences, up to two orders of magnitude, were found between the two methods. Among the uncertainties inherent to the 137Cs technique, definition of the reference inventory seems to be a major source of imprecision. In addition, estimations of water and sediment fluxes with mathematical models usually also present high uncertainty, contributing to the quantitative differences of the soil redistribution estimates with the two methods.