From qualitative data to quantitative models: analysis of the phage shock protein stress response in Escherichia coli

BackgroundBacteria have evolved a rich set of mechanisms for sensing and adapting to adverse conditions in their environment. These are crucial for their survival, which requires them to react to extracellular stresses such as heat shock, ethanol treatment or phage infection. Here we focus on studying the phage shock protein (Psp) stress response in Escherichia coli induced by a phage infection or other damage to the bacterial membrane. This system has not yet been theoretically modelled or analysed in silico.ResultsWe develop a model of the Psp response system, and illustrate how such models can be constructed and analyzed in light of available sparse and qualitative information in order to generate novel biological hypotheses about their dynamical behaviour. We analyze this model using tools from Petri-net theory and study its dynamical range that is consistent with currently available knowledge by conditioning model parameters on the available data in an approximate Bayesian computation (ABC) framework. Within this ABC approach we analyze stochastic and deterministic dynamics. This analysis allows us to identify different types of behaviour and these mechanistic insights can in turn be used to design new, more detailed and time-resolved experiments.ConclusionsWe have developed the first mechanistic model of the Psp response in E. coli. This model allows us to predict the possible qualitative stochastic and deterministic dynamic behaviours of key molecular players in the stress response. Our inferential approach can be applied to stress response and signalling systems more generally: in the ABC framework we can condition mathematical models on qualitative data in order to delimit e.g. parameter ranges or the qualitative system dynamics in light of available end-point or qualitative information.

[1]  Michael P. H. Stumpf,et al.  Simulation-based model selection for dynamical systems in systems and population biology , 2009, Bioinform..

[2]  Chris J. Myers,et al.  Engineering Genetic Circuits , 2009 .

[3]  A. J. Darwin,et al.  The phage‐shock‐protein response , 2005, Molecular microbiology.

[4]  Attila Csikász-Nagy,et al.  Stochastic Petri Net extension of a yeast cell cycle model. , 2008, Journal of theoretical biology.

[5]  Michael P H Stumpf,et al.  Induction and Function of the Phage Shock Protein Extracytoplasmic Stress Response in Escherichia coli* , 2006, Journal of Biological Chemistry.

[6]  Claudine Chaouiya,et al.  Petri net modelling of biological networks , 2007, Briefings Bioinform..

[7]  J. Kormanec,et al.  Pushing the envelope: extracytoplasmic stress responses in bacterial pathogens , 2006, Nature Reviews Microbiology.

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

[9]  K. Liberek,et al.  [Regulation of Escherichia coli heat shock response]. , 1995, Postepy biochemii.

[10]  Tadao Murata,et al.  Petri nets: Properties, analysis and applications , 1989, Proc. IEEE.

[11]  Gong-Xin Yu,et al.  Ruleminer: a Knowledge System for Supporting High-throughput Protein Function Annotations , 2004, J. Bioinform. Comput. Biol..

[12]  M. Stumpf,et al.  The Evolution of the Phage Shock Protein Response System: Interplay between Protein Function, Genomic Organization, and System Function , 2010, Molecular biology and evolution.

[13]  Liming Ying,et al.  In vivo localizations of membrane stress controllers PspA and PspG in Escherichia coli , 2009, Molecular microbiology.

[14]  G. Jovanovic,et al.  Identification, nucleotide sequence, and characterization of PspF, the transcriptional activator of the Escherichia coli stress-induced psp operon , 1996, Journal of bacteriology.

[15]  Michael P. H. Stumpf,et al.  Tutorial on ABC rejection and ABC SMC for parameter estimation and model selection , 2009, 0910.4472.

[16]  J. Dworkin,et al.  The Escherichia coli phage‐shock‐protein (psp) operon , 1997, Molecular microbiology.

[17]  Erika Cule,et al.  ABC-SysBio—approximate Bayesian computation in Python with GPU support , 2010, Bioinform..

[18]  Monika Heiner,et al.  Model Validation of Biological Pathways Using Petri Nets - Demonstrated for Apoptosis , 2003, CMSB.

[19]  References , 1971 .

[20]  Alexander van Oudenaarden,et al.  Variability in gene expression underlies incomplete penetrance , 2009, Nature.

[21]  Manuel Silva,et al.  A Simple and Fast Algorithm to Obtain All Invariants of a Generalized Petri Net , 1980, Selected Papers from the First and the Second European Workshop on Application and Theory of Petri Nets.

[22]  M. Heiner,et al.  Petri nets in biology, chemistry, and medicine : bibliography , 2002 .

[23]  P J Goss,et al.  Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Maciej Koutny,et al.  Application and Theory of Petri Nets 2001, 22nd International Conference, ICATPN 2001, Newcastle upon Tyne, UK, June 25-29, 2001, Proceedings , 2001, ICATPN.

[25]  Monika Heiner,et al.  From Petri Nets to Differential Equations - An Integrative Approach for Biochemical Network Analysis , 2006, ICATPN.

[26]  Tina Toni,et al.  Model-based evolutionary analysis: the natural history of phage-shock stress response. , 2009, Biochemical Society transactions.

[27]  Goran Jovanovic,et al.  Physical, functional and conditional interactions between ArcAB and phage shock proteins upon secretin-induced stress in Escherichia coli , 2009, Molecular microbiology.

[28]  D R Westhead,et al.  Petri Net representations in systems biology. , 2003, Biochemical Society transactions.

[29]  M. Stumpf,et al.  Managing membrane stress: the phage shock protein (Psp) response, from molecular mechanisms to physiology. , 2010, FEMS microbiology reviews.

[30]  P. Model,et al.  Phage shock protein, a stress protein of Escherichia coli. , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[31]  Mark M. Tanaka,et al.  Sequential Monte Carlo without likelihoods , 2007, Proceedings of the National Academy of Sciences.

[32]  M. Feldman,et al.  Population growth of human Y chromosomes: a study of Y chromosome microsatellites. , 1999, Molecular biology and evolution.

[33]  Pierre N. Robillard,et al.  Modeling and Simulation of Molecular Biology Systems Using Petri Nets: Modeling Goals of Various Approaches , 2004, J. Bioinform. Comput. Biol..

[34]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[35]  Falko Bause,et al.  Stochastic Petri Nets , 1996 .

[36]  Darren J. Wilkinson Stochastic Modelling for Systems Biology , 2006 .

[37]  T. Silhavy,et al.  Quality control in the bacterial periplasm. , 2004, Biochimica et biophysica acta.

[38]  C. Petri Kommunikation mit Automaten , 1962 .

[39]  Michael Hucka,et al.  A Correction to the Review Titled "Rules for Modeling Signal-Transduction Systems" by W. S. Hlavacek et al. , 2006, Science's STKE.

[40]  Masasuke Yoshida,et al.  Escherichia coli phage‐shock protein A (PspA) binds to membrane phospholipids and repairs proton leakage of the damaged membranes , 2007, Molecular microbiology.

[41]  Tina Toni,et al.  The ABC of reverse engineering biological signalling systems. , 2009, Molecular bioSystems.