Dynamical modelling of feedback gene regulatory networks

Living cells are made up of networks of interacting genes, proteins and other bio-molecules. Simple interactions between network components in forms of feedback regulations can lead to complex collective dynamics. A key task in cell biology is to gain a thorough understanding of the dynamics of intracellular systems and processes. In this thesis, a combined approach of mathematical modelling, computational simulation and analytical techniques, has been used to obtain a deeper insight into the dynamical aspects of a variety of feedback systems commonly encountered in cells. These systems range from model system with detailed available molecular knowledge to general regulatory motifs with varying network structures. Deterministic as well as stochastic modelling techniques have been employed, depending primarily on the specific questions asked. The first part of the thesis focuses on dissecting the principles behind the regulatory design of the Tryptophan Operon system in Escherichia coli. It has evolved three negative feedback loops, namely repression, attenuation and enzyme inhibition, as core regulator mechanisms to control the intracellular level of tryptophan amino acid, which is taken up for protein synthesis. Despite extensive experimental knowledge, the roles of these seemingly redundant loops remain unclear from a dynamical point of view. We aim to understand why three loops, rather than one, have evolved. Using a large-scale perturbation/response analysis through modelling and simulations and novel metrics for transient dynamics quantification, it has been revealed that the multiple negative feedback loops employed by the tryptophan operon are not redundant. In fact, they have evolved to concertedly give rise to a much more efficient, adaptive and stable system, than any single mechanism would provide. Since even the full topology of feedback interactions within a network is insufficient to determine its behavioural dynamics, other factors underlying feedback loops must be characterised to better predict system dynamics. In the second part of the thesis, we aim to derive these factors and explore how they shape system dynamics. We develop an analytical approach for stability and bifurcation analysis and apply it to class of feedback systems commonly encountered in cells. Our analysis showed that the strength and the Hill coefficient of a feedback loop play key role in determining the dynamics of the system carrying the loop. Not only that, the position of the loop was also found to be crucial in this decision. The analytical method we developed also facilitates parameter sensitivity analysis in which we investigate how the production and degradation rates affect system dynamics. We find that these rates are quite different in the way they shape up system behaviour, with the degradation rates exhibiting a more intricate manner. We demonstrated that coupled-loop systems display greater complexity and a richer repertoire of behaviours in comparison with single-loop ones. Different combinations of the feedback strengths of individual loops give rise to different dynamical regimes. The final part of the thesis aims to understand the effects of molecular noise on dynamics of specific systems, in this case the Tryptophan Operon. We developed two stochastic models for the system and compared their predictions to those given by the deterministic model. By means of simulations, we have shown that noise can induce oscillatory behaviour. On the other hand, incorporating noise in an oscillatory system can alter the characteristics of oscillation by shifting the bifurcation point of certain parameters by a substantial amount. Measurement of fluctuations reveals that that noise at the transcript level is most significant while noise at the enzyme level is smallest. This study highlights that noise should not be neglected if we want to obtain a complete understanding of the dynamic behaviour of cells