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.

Propagation through cutting technique of species ocurring in the Lower São Francisco River in Sergipe State with different concentrations of indolbutiric acid

The objective of this work was to evaluate the feasibility of vegetative propagation through cutting technique of seven tree species with strong occurrence in the riparian forest of the Lower São Francisco River in Sergipe State, under different concentrations of indolbutiric acid at 0, 2500, and 5000 mg.L-1, for potentialization of its use in soil bioengineering technique. It was used a complete random block design with three replicates, and a total of twenty-one treatments. The evaluation period was 120 days for each species, and the data collection was made in intervals of fifteen days, in a total of eight evaluations for each species. The evaluated parameters were: Survival Rate, callus formation, and Root Dry matter Weight. Among the studied species, Schinus terebinthifolius Raddi presented the best results related to cutting technique mainly under the indolbutiric acid concentration of 2500 mg.L-1.

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.

Vegetative propagation of bauhinia x blakeana, an ornamental sterile tree

Bauhinia x blakeana (B. purpurea x B. variegata) is a natural hybrid that has been cultivated in gardens, streets and parks. Due to its sterility, it must be vegetatively propagated. The objective of this work was to evaluate the viability of cuttings and grafting on its propagation. Semi-woody cuttings were collected during four seasons and treated with 0; 1,000; 2,000; and 3,000 mg L-1 of IBA. The experimental design was entirely randomized and the treatments were arranged in a 4x4 factorial scheme (four collecting times x four IBA concentrations) and five replications with 10 cuttings each, per collecting time and per IBA concentration. Characteristics of roots and shoots were evaluated after 90 days. The data means were compared by the Tukey test and submitted to the polynomial regression analysis. For the grafting experiment, B. variegata and B. variegata var. candida plants of six and 12 months were used as rootstocks and the splice graft and T-budding methods were tested. The experimental design was entirely randomized and the treatments were arranged in a 2x2x2 factorial scheme (two rootstock species x two grafting methods x two rootstock ages) and four replications with five plants each, per rootstock species, per grafting method and per rootstock age. Characteristics of shoots were evaluated after 90 days and the data means were compared by the Tukey test. B. x blakeana can be propagated by semi-woody cuttings collected in spring, without IBA application, or in summer, with the application of 3,000 mg L-1 of IBA. The tested grafting methods were not effective.

VEGETATIVE RESCUE AND CUTTINGS PROPAGATION OF Araucaria angustifolia (Bertol.) Kuntze

Brazilian pine or araucaria (Araucaria angustifolia) is a coniferous tree with great economic, social and environmental importance in southern Brazil, being exploited for both wood production and for its edible pine nuts. However, no efficient cloning techniques are available and, therefore, the purpose of this study was to evaluate the effectiveness of vegetative rescue methods for cuttings propagation of the species. Shoots/cuttings were generated in two ways: 26 years old trees underwent coppicing and 20 years old trees had the primary branches on the upper third of crown pruned at 2, 20 and 50 cm from the main trunk. Orthotropic shoots were rooted after application of indole-3-butyric acid (IBA) at 0, 2, 4 and 6 g.L-1. Coppicing produced 47 cuttings per plant with 90% orthotropic shoots, while pruning resulted in 182 cuttings per plant with 44% orthotropic shoots. Rooting success indexes were low with no influence of IBA, although they are slightly superior to the ones available in the literature for the species, ranging from 12 to 30% for the coppice shoots and from 0 to 28% for the branches shoots. We conclude that both vegetative rescue techniques are viable and have potentially important applications. Coppicing is recommended for the propagation aiming the production of wood, while shoots derived from the side branches of the crown are more appropriate for seeds orchards formation.

CLONAL PROPAGATION OF NEEM (Azadirachta indica A. Juss.) VIA DIRECT AND INDIRECT IN VITRO REGENERATION

ABSTRACTThe study was conducted with shoot tip explants of neem (Azadirachta indica A. Juss) to identify a viable regenerative process. Shoot tips were obtained from neem embryos cultured alternatingly in DKW medium supplemented with BAP and medium without hormones. Initial shoot development was influenced by cotyledon presence. Basal callus, excised from in vitro stem base, also presented organogenic potential. In some cases, plant lines, obtained from each seed, presented different characteristics. The most common characteristic observed in vitro was callus formation at the stem base. However, the rarest characteristics were stem callus formation and leaf senescence. The regenerated shoot tips were further subculture and rooted on a medium supplemented with IBA so that complete plants could be obtained. The rooted plants were transplanted to a greenhouse and successfully acclimatized. No significant differences in in vivo development were observed between neem plants from callus and from shoot tip propagation.

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.

Spatio-temporal modelling of the duration of the cotton cycle in the State Of Goiás, Brazil

The objective of this work was to develop and validate a mathematical model to estimate the duration of cotton (Gossypium hirsutum L. r. latifolium hutch) cycle in the State of Goiás, Brazil, by applying the method of growing degree-days (GD), and considering, simultaneously, its time-space variation. The model was developed as a linear combination of elevation, latitude, longitude, and Fourier series of time variation. The model parameters were adjusted by using multiple-linear regression to the observed GD accumulated with air temperature in the range of 15°C to 40°C. The minimum and maximum temperature records used to calculate the GD were obtained from 21 meteorological stations, considering data varying from 8 to 20 years of observation. The coefficient of determination, resulting from the comparison between the estimated and calculated GD along the year was 0.84. Model validation was done by comparing estimated and measured crop cycle in the period from cotton germination to the stage when 90 percent of bolls were opened in commercial crop fields. Comparative results showed that the model performed very well, as indicated by the Pearson correlation coefficient of 0.90 and Willmott agreement index of 0.94, resulting in a performance index of 0.85.

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.

Modelling head loss along emitting pipes using dimensional analysis

Local head losses must be considered in estimating properly the maximum length of drip irrigation laterals. The aim of this work was to develop a model based on dimensional analysis for calculating head loss along laterals accounting for in-line drippers. Several measurements were performed with 12 models of emitters to obtain the experimental data required for developing and assessing the model. Based on the Camargo & Sentelhas coefficient, the model presented an excellent result in terms of precision and accuracy on estimating head loss. The deviation between estimated and observed values of head loss increased according to the head loss and the maximum deviation reached 0.17 m. The maximum relative error was 33.75% and only 15% of the data set presented relative errors higher than 20%. Neglecting local head losses incurred a higher than estimated maximum lateral length of 19.48% for pressure-compensating drippers and 16.48% for non pressure-compensating drippers.

Software for design of experiments and response modelling of cake ltration applications

Filtration is a widely used unit operation in chemical engineering. The huge variation in the properties of materials to be ltered makes the study of ltration a challenging task. One of the objectives of this thesis was to show that conventional ltration theories are di cult to use when the system to be modelled contains all of the stages and features that are present in a complete solid/liquid separation process. Furthermore, most of the ltration theories require experimental work to be performed in order to obtain critical parameters required by the theoretical models. Creating a good overall understanding of how the variables a ect the nal product in ltration is somewhat impossible on a purely theoretical basis. The complexity of solid/liquid separation processes require experimental work and when tests are needed, it is advisable to use experimental design techniques so that the goals can be achieved. The statistical design of experiments provides the necessary tools for recognising the e ects of variables. It also helps to perform experimental work more economically. Design of experiments is a prerequisite for creating empirical models that can describe how the measured response is related to the changes in the values of the variable. A software package was developed that provides a ltration practitioner with experimental designs and calculates the parameters for linear regression models, along with the graphical representation of the responses. The developed software consists of two software modules. These modules are LTDoE and LTRead. The LTDoE module is used to create experimental designs for di erent lter types. The lter types considered in the software are automatic vertical pressure lter, double-sided vertical pressure lter, horizontal membrane lter press, vacuum belt lter and ceramic capillary action disc lter. It is also possible to create experimental designs for those cases where the variables are totally user de ned, say for a customized ltration cycle or di erent piece of equipment. The LTRead-module is used to read the experimental data gathered from the experiments, to analyse the data and to create models for each of the measured responses. Introducing the structure of the software more in detail and showing some of the practical applications is the main part of this thesis. This approach to the study of cake ltration processes, as presented in this thesis, has been shown to have good practical value when making ltration tests.