Parameter sensitivity analysis of stochastic models: Application to catalytic reaction networks

A general numerical methodology for parametric sensitivity analysis is proposed, which allows to determine the parameters exerting the greatest influence on the output of a stochastic computational model, especially when the knowledge about the actual value of a parameter is insufficient. An application of the procedure is performed on a model of protocell, in order to detect the kinetic rates mainly affecting the capability of a catalytic reaction network enclosed in a semi-permeable membrane to retain material from its environment and to generate a variety of molecular species within its boundaries. It is shown that the former capability is scarcely sensitive to variations in the model parameters, whereas a kinetic rate responsible for profound modifications of the latter can be identified and it depends on the specific reaction network. A faster uptaking of limited resources from the environment may have represented a significant advantage from an evolutionary point of view and this result is a first indication in order to decipher which kind of structures are more suitable to achieve a viable evolution.

[1]  Hua Wu,et al.  Parametric sensitivity in chemical systems , 1999 .

[2]  H Christopher Frey,et al.  OF SENSITIVITY ANALYSIS , 2001 .

[3]  S. Kauffman.  Autocatalytic sets of proteins. , 1986 .

[4]  M. Eigen,et al.  The Hypercycle: A principle of natural self-organization , 2009 .

[5]  M Villani,et al.  Sufficient conditions for emergent synchronization in protocell models. , 2008, Journal of theoretical biology.

[6]  John W. Pepper,et al.  Simple Models of Assortment through Environmental Feedback , 2007, Artificial Life.

[7]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[8]  J. Szostak,et al.  The origins of cellular life. , 2010, Cold Spring Harbor perspectives in biology.

[9]  Tom Lenaerts,et al.  Advances in Artificial Life, ECAL 2011: Proceedings of the Eleventh European Conference on the Synthesis and Simulation of Living Systems , 2011 .

[10]  Roberto Serra,et al.  A stochastic model of autocatalytic reaction networks , 2012, Theory in Biosciences.

[11]  Roberto Serra,et al.  The role of energy in a stochastic model of the emergence of autocatalytic sets , 2011, ECAL.

[12]  H. Nürnberg.  The Hypercycle. A Principle of Natural Self Organization. , 1981 .

[13]  Kunihiko Kaneko,et al.  Life: An Introduction to Complex Systems Biology (Understanding Complex Systems) , 2006 .

[14]  M. Kramer,et al.  Sensitivity Analysis in Chemical Kinetics , 1983 .

[15]  Roberto Serra,et al.  Non-linear protocell models: synchronization and chaos , 2010 .

[16]  D. Hamby.  A review of techniques for parameter sensitivity analysis of environmental models , 1994, Environmental monitoring and assessment.

[17]  Roberto Serra,et al.  A stochastic model of the emergence of autocatalytic cycles , 2011 .

[18]  Peter Schuster,et al.  A principle of natural self-organization , 1977, Naturwissenschaften.

[19]  Kim Dam-Johansen,et al.  Visualization Methods in Analysis of Detailed Chemical Kinetics Modelling , 2001, Comput. Chem..

[20]  Athel Cornish-Bowden,et al.  Self-organization at the origin of life. , 2008, Journal of theoretical biology.

[21]  S. Kauffman,et al.  Autocatalytic replication of polymers , 1986 .

[22]  Yang Cao,et al.  Sensitivity analysis of discrete stochastic systems. , 2005, Biophysical journal.

[23]  M. Eigen,et al.  The hypercycle. A principle of natural self-organization. Part A: Emergence of the hypercycle. , 1977, Die Naturwissenschaften.

[24]  Jérôme Morio,et al.  Global and local sensitivity analysis methods for a physical system , 2011 .

[25]  A. Graudenzi,et al.  The influence of different kinetic rates on the dynamics of a simple model of catalytic reaction network , 2012 .

[26]  E. Parzen.  On Estimation of a Probability Density Function and Mode , 1962 .

[27]  S. A. Kauffman,et al.  Autocatalytic sets of proteins , 1986, Origins of life and evolution of the biosphere.

[28]  T. Turányi.  Sensitivity analysis of complex kinetic systems. Tools and applications , 1990 .

[29]  Pasquale Stano,et al.  The Minimal Size of Liposome‐Based Model Cells Brings about a Remarkably Enhanced Entrapment and Protein Synthesis , 2009, Chembiochem : a European journal of chemical biology.

[30]  Sanjay Jain,et al.  Autocatalytic sets and the growth of complexity in an evolutionary model , 1998, adap-org/9809003.

[31]  Roberto Serra,et al.  Synchronization Phenomena in Surface-Reaction Models of Protocells , 2007, Artificial Life.

[32]  John H. Seinfeld,et al.  Stochastic sensitivity analysis in chemical kinetics , 1981 .

[33]  M. Eigen,et al.  The Hypercycle , 2004, Naturwissenschaften.

[34]  Michael P H Stumpf,et al.  Sensitivity, robustness, and identifiability in stochastic chemical kinetics models , 2011, Proceedings of the National Academy of Sciences.

[35]  Roberto Serra,et al.  A model of protocell based on the introduction of a semi-permeable membrane in a stochastic model of catalytic reaction networks , 2013, Wivace.

[36]  C. Keating,et al.  Polymer encapsulation within giant lipid vesicles. , 2007, Langmuir : the ACS journal of surfaces and colloids.

[37]  Stefano Tarantola,et al.  Sensitivity Analysis as an Ingredient of Modeling , 2000 .

[38]  A S Perelson,et al.  Modeling adaptive biological systems. , 1989, Bio Systems.

[39]  Pasquale Stano,et al.  Achievements and open questions in the self-reproduction of vesicles and synthetic minimal cells. , 2010, Chemical communications.

[40]  J. Szostak,et al.  Expanding roles for diverse physical phenomena during the origin of life. , 2010, Annual review of biophysics.

[41]  Martin Nilsson,et al.  Bridging Nonliving and Living Matter , 2003, Artificial Life.

[42]  A. Saltelli,et al.  Sensitivity Anaysis as an Ingredient of Modeling , 2000 .

[43]  Roberto Serra,et al.  SYNCHRONIZATION PHENOMENA IN PROTOCELL MODELS , 2008 .