Novel Multiscale Modeling Tool Applied to Pseudomonas aeruginosa Biofilm Formation

Multiscale modeling is used to represent biological systems with increasing frequency and success. Multiscale models are often hybrids of different modeling frameworks and programming languages. We present the MATLAB-NetLogo extension (MatNet) as a novel tool for multiscale modeling. We demonstrate the utility of the tool with a multiscale model of Pseudomonas aeruginosa biofilm formation that incorporates both an agent-based model (ABM) and constraint-based metabolic modeling. The hybrid model correctly recapitulates oxygen-limited biofilm metabolic activity and predicts increased growth rate via anaerobic respiration with the addition of nitrate to the growth media. In addition, a genome-wide survey of metabolic mutants and biofilm formation exemplifies the powerful analyses that are enabled by this computational modeling tool.

[1]  F. Doyle,et al.  Dynamic flux balance analysis of diauxic growth in Escherichia coli. , 2002, Biophysical journal.

[2]  P. Stewart,et al.  Spatial Physiological Heterogeneity inPseudomonas aeruginosa Biofilm Is Determined by Oxygen Availability , 1998, Applied and Environmental Microbiology.

[3]  O. Wanner.  Modeling of Biofilms , 2003 .

[4]  Jason A. Papin,et al.  Metabolic Network Analysis of Pseudomonas aeruginosa during Chronic Cystic Fibrosis Lung Infection , 2010, Journal of bacteriology.

[5]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[6]  Claude Lobry,et al.  Bacteria can form interconnected microcolonies when a self-excreted product reduces their surface motility: evidence from individual-based model simulations , 2010, Theory in Biosciences.

[7]  Jay D. Humphrey,et al.  Ensuring Congruency in Multiscale Modeling: Towards Linking Agent Based and Continuum Biomechanical Models of Arterial Adaptation , 2011, Annals of Biomedical Engineering.

[8]  Dieter Jahn,et al.  Anaerobic adaptation in Pseudomonas aeruginosa: definition of the Anr and Dnr regulons. , 2009, Environmental microbiology.

[9]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[10]  Philip S. Stewart,et al.  Diffusion in Biofilms , 2003, Journal of bacteriology.

[11]  O. Lesouhaitier,et al.  TypA is involved in virulence, antimicrobial resistance and biofilm formation in Pseudomonas aeruginosa , 2013, BMC Microbiology.

[12]  Nicholas R. Magliocca,et al.  Using Pattern‐oriented Modeling (POM) to Cope with Uncertainty in Multi‐scale Agent‐based Models of Land Change , 2013, Trans. GIS.

[13]  N. Høiby,et al.  Pseudomonas aeruginosa biofilms in the respiratory tract of cystic fibrosis patients , 2009, Pediatric pulmonology.

[14]  Steven F. Railsback,et al.  Agent-Based and Individual-Based Modeling: A Practical Introduction , 2011 .

[15]  J. Chave.  The problem of pattern and scale in ecology: what have we learned in 20 years? , 2013, Ecology letters.

[16]  G. Deffuant,et al.  Viability and Resilience of a Bacterial Biofilm Individual-Based Model , 2011 .

[17]  Jason A. Papin,et al.  Multiscale computational analysis of Xenopus laevis morphogenesis reveals key insights of systems-level behavior , 2007, BMC Systems Biology.

[18]  C. Keel,et al.  Role of Gluconic Acid Production in the Regulation of Biocontrol Traits of Pseudomonas fluorescens CHA0 , 2009, Applied and Environmental Microbiology.

[19]  Tao Jia,et al.  An empirical study on human mobility and its agent-based modeling , 2012 .

[20]  Colin Harbour,et al.  Gene expression characteristics of a cystic fibrosis epidemic strain of Pseudomonas aeruginosa during biofilm and planktonic growth. , 2009, FEMS microbiology letters.

[21]  Garth D. Ehrlich,et al.  Oxygen Limitation Contributes to Antibiotic Tolerance of Pseudomonas aeruginosa in Biofilms , 2004, Antimicrobial Agents and Chemotherapy.

[22]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[23]  J. Seoane,et al.  Growth dependence of conjugation explains limited plasmid invasion in biofilms: an individual-based modelling study. , 2011, Environmental microbiology.

[24]  Jason A. Papin,et al.  * Corresponding authors , 2006 .

[25]  Erwin P. Gianchandani,et al.  The application of flux balance analysis in systems biology , 2010, Wiley interdisciplinary reviews. Systems biology and medicine.

[26]  Cristian Picioreanu,et al.  iDynoMiCS: next-generation individual-based modelling of biofilms. , 2011, Environmental microbiology.

[27]  Jay D. Humphrey,et al.  Toward a Multi-Scale Computational Model of Arterial Adaptation in Hypertension: Verification of a Multi-Cell Agent Based Model , 2011, Front. Physio..

[28]  Shayn M Peirce,et al.  Multiscale computational models of complex biological systems. , 2013, Annual review of biomedical engineering.

[29]  Volker Grimm,et al.  NetLogo meets R: Linking agent-based models with a toolbox for their analysis , 2010, Environ. Model. Softw..

[30]  D. Griffeath,et al.  Quantitative cellular automaton model for biofilms , 2001 .