Modeling and simulation of a synthetic genetic circuit that implements a multicelled behavior in a growing microcolony of E. coli

Synthetic Biology is a relatively new discipline, born at the beginning of the New Millennium, that brings the typical engineering approach (abstraction, modularity and standardization) to biotechnology. These principles aim to tame the extreme complexity of the various components and aid the construction of artificial biological systems with specific functions, usually by means of synthetic genetic circuits implemented in bacteria or simple eukaryotes like yeast. The cell becomes a programmable machine and its low-level programming language is made of strings of DNA. This work was performed in collaboration with researchers of the Department of Electrical Engineering of the University of Washington in Seattle and also with a student of the Corso di Laurea Magistrale in Ingegneria Biomedica at the University of Bologna: Marilisa Cortesi. During the collaboration I contributed to a Synthetic Biology project already started in the Klavins Laboratory. In particular, I modeled and subsequently simulated a synthetic genetic circuit that was ideated for the implementation of a multicelled behavior in a growing bacterial microcolony. In the first chapter the foundations of molecular biology are introduced: structure of the nucleic acids, transcription, translation and methods to regulate gene expression. An introduction to Synthetic Biology completes the section. In the second chapter is described the synthetic genetic circuit that was conceived to make spontaneously emerge, from an isogenic microcolony of bacteria, two different groups of cells, termed leaders and followers. The circuit exploits the intrinsic stochasticity of gene expression and intercellular communication via small molecules to break the symmetry in the phenotype of the microcolony. The four modules of the circuit (coin flipper, sender, receiver and follower) and their interactions are then illustrated. In the third chapter is derived the mathematical representation of the various components of the circuit and the several simplifying assumptions are made explicit. Transcription and translation are modeled as a single step and gene expression is function of the intracellular concentration of the various transcription factors that act on the different promoters of the circuit. A list of the various parameters and a justification for their value closes the chapter. In the fourth chapter are described the main characteristics of the gro simulation environment, developed by the Self Organizing Systems Laboratory of the University of Washington. Then, a sensitivity analysis performed to pinpoint the desirable characteristics of the various genetic components is detailed. The sensitivity analysis makes use of a cost function that is based on the fraction of cells in each one of the different possible states at the end of the simulation and the wanted outcome. Thanks to a particular kind of scatter plot, the parameters are ranked. Starting from an initial condition in which all the parameters assume their nominal value, the ranking suggest which parameter to tune in order to reach the goal. Obtaining a microcolony in which almost all the cells are in the follower state and only a few in the leader state seems to be the most difficult task. A small number of leader cells struggle to produce enough signal to turn the rest of the microcolony in the follower state. It is possible to obtain a microcolony in which the majority of cells are followers by increasing as much as possible the production of signal. Reaching the goal of a microcolony that is split in half between leaders and followers is comparatively easy. The best strategy seems to be increasing slightly the production of the enzyme. To end up with a majority of leaders, instead, it is advisable to increase the basal expression of the coin flipper module. At the end of the chapter, a possible future application of the leader election circuit, the spontaneous formation of spatial patterns in a microcolony, is modeled with the finite state machine formalism. The gro simulations provide insights into the genetic components that are needed to implement the behavior. In particular, since both the examples of pattern formation rely on a local version of Leader Election, a short-range communication system is essential. Moreover, new synthetic components that allow to reliably downregulate the growth rate in specific cells without side effects need to be developed. In the appendix are listed the gro code utilized to simulate the model of the circuit, a script in the Python programming language that was used to split the simulations on a Linux cluster and the Matlab code developed to analyze the data.

Sea level trends in the Mediterranean from tide gauges and satellite altimetry

Sea level variation is one of the parameters directly related to climate change. Monitoring sea level rise is an important scientific issue since many populated areas of the world and megacities are located in low-lying regions. At present, sea level is measured by means of two techniques: the tide gauges and the satellite radar altimetry. Tide gauges measure sea-level relatively to a ground benchmark, hence, their measurements are directly affected by vertical ground motions. Satellite radar altimetry measures sea-level relative to a geocentric reference and are not affected by vertical land motions. In this study, the linear relative sea level trends of 35 tide gauge stations distributed across the Mediterranean Sea have been computed over the period 1993-2014. In order to extract the real sea-level variation, the vertical land motion has been estimated using the observations of available GPS stations and removed from the tide gauges records. These GPS-corrected trends have then been compared with satellite altimetry measurements over the same time interval (AVISO data set). A further comparison has been performed, over the period 1993-2013, using the CCI satellite altimetry data set which has been generated using an updated modeling. The absolute sea level trends obtained from satellite altimetry and GPS-corrected tide gauge data are mostly consistent, meaning that GPS data have provided reliable corrections for most of the sites. The trend values range between +2.5 and +4 mm/yr almost everywhere in the Mediterranean area, the largest trends were found in the Northern Adriatic Sea and in the Aegean. These results are in agreement with estimates of the global mean sea level rise over the last two decades. Where GPS data were not available, information on the vertical land motion deduced from the differences between absolute and relative trends are in agreement with the results of other studies.

Modeling and Model Predictive Control of an Industrial Hydraulic System

This thesis aims to illustrate the construction of a mathematical model of a hydraulic system, oriented to the design of a model predictive control (MPC) algorithm. The modeling procedure starts with the basic formulation of a piston-servovalve system. The latter is a complex non linear system with some unknown and not measurable effects that constitute a challenging problem for the modeling procedure. The first level of approximation for system parameters is obtained basing on datasheet informations, provided workbench tests and other data from the company. Then, to validate and refine the model, open-loop simulations have been made for data matching with the characteristics obtained from real acquisitions. The final developed set of ODEs captures all the main peculiarities of the system despite some characteristics due to highly varying and unknown hydraulic effects, like the unmodeled resistive elements of the pipes. After an accurate analysis, since the model presents many internal complexities, a simplified version is presented. The latter is used to linearize and discretize correctly the non linear model. Basing on that, a MPC algorithm for reference tracking with linear constraints is implemented. The results obtained show the potential of MPC in this kind of industrial applications, thus a high quality tracking performances while satisfying state and input constraints. The increased robustness and flexibility are evident with respect to the standard control techniques, such as PID controllers, adopted for these systems. The simulations for model validation and the controlled system have been carried out in a Python code environment.

Numerical modeling on the short-term response of fiber-reinforced concrete round panels

Fiber-reinforced concrete is a composite material consisting of discrete, discontinuous, and uniformly distributed fibers in plain concrete primarily used to enhance the tensile properties of the concrete. FRC performance depends upon the fiber, interface, and matrix properties. The use of fiber-reinforced concrete has been increasing substantially in the past few years in different fields of the construction industry such as ground-level application in sidewalks and building floors, tunnel lining, aircraft parking, runways, slope stabilization, etc. Many experiments have been performed to observe the short-term and long-term mechanical behavior of fiber-reinforced concrete in the last decade and numerous numerical models have been formulated to accurately capture the response of fiber-reinforced concrete. The main purpose of this dissertation is to numerically calibrate the short-term response of the concrete and fiber parameters in mesoscale for the three-point bending test and cube compression test in the MARS framework which is based on the lattice discrete particle model (LDPM) and later validate the same parameters for the round panels. LDPM is the most validated theory in mesoscale theories for concrete. Different seeds representing the different orientations of concrete and fiber particles are simulated to produce the mean numerical response. The result of numerical simulation shows that the lattice discrete particle model for fiber-reinforced concrete can capture results of experimental tests on the behavior of fiber-reinforced concrete to a great extent.