This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.104.052126v1 88/4/2541    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Zhou, B. Articles by Liao, J. C. Search for Related Content PubMed PubMed Citation Articles by Zhou, B. Articles by Liao, J. C.

Markov Chain Modeling of Pyelonephritis-Associated Pili Expression in Uropathogenic Escherichia coli

Department of Chemical Engineering, University of California at Los Angeles, California

Correspondence: Address reprint requests to James C. Liao, 5531 Boelter Hall, Dept. of Chemical Engineering, University of California at Los Angeles, Los Angeles, CA 90095. Tel.: 310-825-1656; Fax: 310-206-4107; E-mail: liaoj{at}ucla.edu.

Pyelonephritis-associated pili (Pap) expression in uropathogenic Escherichia coli is regulated by a complex phase variation mechanism involving the competition between leucine-responsive regulatory protein (Lrp) and DNA adenine methylase (Dam). Population dynamics of pap gene expression has been studied extensively and the detailed molecular mechanism has been largely elucidated, providing sufficient information for mathematical modeling. Although the Gillespie algorithm is suited for modeling of stochastic systems such as the pap operon, it becomes computationally expensive when detailed molecular steps are explicitly modeled in a population. Here we developed a Markov Chain model to simplify the computation. Our model is analytically derived from the molecular mechanism. The model presented here is able to reproduce results presented using the Gillespie method, but since the regulatory information is incorporated before simulation, our model runs more efficiently and allows investigation of additional regulatory features. The model predictions are consistent with experimental data obtained in this work and in the literature. The results show that pap expression in uropathogenic E. coli is initial-state-dependent, as previously reported. However, without environment stimuli, the pap-expressing fraction in a population will reach an equilibrium level after 50–100 generations. The transient time before reaching equilibrium is determined by PapI stability and Lrp and Dam copy numbers per cell. This work demonstrates that the Markov Chain model captures the essence of the complex molecular mechanism and greatly simplifies the computation.

This article has been cited by other articles:

				X. Fang, W. E. Bentley, and E. Zafiriou Stochastic Modeling of Gene Positive Autoregulation Networks Involving Signal Molecules Biophys. J., October 1, 2008; 95(7): 3137 - 3145. [Abstract] [Full Text] [PDF]