English   español  
Please use this identifier to cite or link to this item: http://hdl.handle.net/10261/158922
logo share SHARE   Add this article to your Mendeley library MendeleyBASE
Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL
Exportar a otros formatos:


Mathematical modelling, analysis and numerical simulation for a general class of gene regulatory networks

AuthorsPájaro, Manuel
AdvisorAlonso, Antonio A. ; Vázquez, Carlos
Issue Date2017
AbstractThe research work developed in this thesis is mainly oriented to the mathematical modelling of biological systems, the behaviour of which is inherently stochastic, as it is the case of gene regulatory networks. Their relevance emerges from the fact that all necessary information for life cycle is encoded in the DNA. Consequently, the study of DNA expression, transcription into messenger RNA and translation into proteins, together with their regulation becomes essential to predict cells response to environmental signals. The inherent stochastic nature of gene expression makes these systems to be far away from the classical kinetic limit where the (macroscopic) deterministic methods are valid. In modelling these systems, we need to employ microscopic methods which incorporate the underlying stochastic behaviour. The Chemical Master Equation (CME) remains at the basis for the modelling of these phenomena. However, a closed form solution of the CME is unavailable in general, due to the large number (eventually in nity) of coupled equations. A widespread technique to approximate the CME solution is the Stochastic Simulation Algorithm (SSA), a computationally involved Monte Carlo type method. Although many numerical approximations emerge to reduce the complexity of the CME, we will focus on the Partial Integro-Di erential Equation (PIDE) or Friedman model, which represents a continuous approximation of the CME. For the one dimensional version (self-regulation), the PIDE model has an analytic solution for its steady state. This fact will allow us to characterize the regions in the space of parameters in which the system changes its behaviour (unimodal, bimodal). Also we have carried out an stability analysis by means of entropy methods
Moreover, we obtain a multidimensional version of the Friedman model to handle more complex gene regulatory networks with more than one gene. The mathematical properties of the corresponding equation will be exhaustively analyzed, also with special emphasis on stability of the solution using entropy methods. In addition, we propose two semi-Lagrangian methods for the numerical solution of the multidimensional model. The rst method results very e cient and scalable to higher dimensions, as the numerical results illustrate, although in practice exhibits rst order convergence in time and space. Solutions provided by the proposed method are compared with those obtained by SSA to assess e ciency, accuracy and computational costs. For the second semi-Lagrangian method we develop the theoretical numerical analysis, thus proving second order convergence in time and space. This is clearly illustrated by a numerical example. However, the computational cost of this second approach results much higher, so that the scalability to higher dimensions seems a di cult task. All the numerical techniques have been implemented on a user friendly toolbox (SELANSI) which is detailed in the Appendix
Description327 pages
Publisher version (URL)http://www.investigo.biblioteca.uvigo.es/xmlui/handle/11093/917
Appears in Collections:(IIM) Tesis
Files in This Item:
File Description SizeFormat 
Tesis_Manuel_Pajaro.pdf16,69 MBAdobe PDFThumbnail
Show full item record
Review this work

WARNING: Items in Digital.CSIC are protected by copyright, with all rights reserved, unless otherwise indicated.