PhysiBoSS: a multi-scale agent-based modelling framework integrating physical dimension and cell signalling

Abstract Motivation Due to the complexity and heterogeneity of multicellular biological systems, mathematical models that take into account cell signalling, cell population behaviour and the extracellular environment are particularly helpful. We present PhysiBoSS, an open source software which combines intracellular signalling using Boolean modelling (MaBoSS) and multicellular behaviour using agent-based modelling (PhysiCell). Results PhysiBoSS provides a flexible and computationally efficient framework to explore the effect of environmental and genetic alterations of individual cells at the population level, bridging the critical gap from single-cell genotype to single-cell phenotype and emergent multicellular behaviour. PhysiBoSS thus becomes very useful when studying heterogeneous population response to treatment, mutation effects, different modes of invasion or isomorphic morphogenesis events. To concretely illustrate a potential use of PhysiBoSS, we studied heterogeneous cell fate decisions in response to TNF treatment. We explored the effect of different treatments and the behaviour of several resistant mutants. We highlighted the importance of spatial information on the population dynamics by considering the effect of competition for resources like oxygen. Availability and implementation PhysiBoSS is freely available on GitHub (https://github.com/sysbio-curie/PhysiBoSS), with a Docker image (https://hub.docker.com/r/gletort/physiboss/). It is distributed as open source under the BSD 3-clause license. Supplementary information Supplementary data are available at Bioinformatics online.

[1]  Joseph M. Mahaffy,et al.  Genetic control models with diffusion and delays , 1988 .

[2]  Jordi Garcia-Ojalvo,et al.  Modeling gene expression in time and space. , 2013, Annual review of biophysics.

[3]  A. Czirók,et al.  Collective cell motion in endothelial monolayers , 2010, Physical biology.

[4]  Alexander G. Fletcher,et al.  Comparing individual-based approaches to modelling the self-organization of multicellular tissues , 2016, bioRxiv.

[5]  Robert Noble,et al.  Spatial competition constrains resistance to targeted cancer therapy , 2017, Nature Communications.

[6]  Nadine Peyriéras,et al.  A cell-based computational model of early embryogenesis coupling mechanical behaviour and gene regulation , 2017, Nature Communications.

[7]  Manuel Théry,et al.  Cell shape and cell division. , 2006, Current opinion in cell biology.

[8]  Alexander G. Fletcher,et al.  Chaste: An Open Source C++ Library for Computational Physiology and Biology , 2013, PLoS Comput. Biol..

[9]  Ryan A. Kellogg,et al.  Digital signaling decouples activation probability and population heterogeneity , 2015, eLife.

[10]  Denis Thieffry,et al.  Mathematical Modelling of Cell-Fate Decision in Response to Death Receptor Engagement , 2010, PLoS Comput. Biol..

[11]  Andrew C. Oates,et al.  Delta-Notch signalling in segmentation , 2017, Arthropod structure & development.

[12]  Roeland M. H. Merks,et al.  Breakthrough Technologies VirtualLeaf : An Open-Source Framework for Cell-Based Modeling of Plant Tissue Growth and Development 1 , 2010 .

[13]  Justin M. Wozniak,et al.  High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow , 2017, BMC Bioinformatics.

[14]  Ignacio Ramis-Conde,et al.  Multi-scale modelling of cancer cell intravasation: the role of cadherins in metastasis , 2009, Physical biology.

[15]  Denis Wirtz,et al.  Modeling the two-way feedback between contractility and matrix realignment reveals a nonlinear mode of cancer cell invasion , 2017, Proceedings of the National Academy of Sciences.

[16]  Jean Clairambault,et al.  An age-and-cyclin-structured cell population model for healthy and tumoral tissues , 2008, Journal of mathematical biology.

[17]  Zev J. Gartner,et al.  A strategy for tissue self-organization that is robust to cellular heterogeneity and plasticity , 2015, Proceedings of the National Academy of Sciences.

[18]  Hiroshi Momiji,et al.  Oscillatory Notch-pathway activity in a delay model of neuronal differentiation. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  G. An,et al.  Agent‐based models in translational systems biology , 2009, Wiley interdisciplinary reviews. Systems biology and medicine.

[20]  H. Frieboes,et al.  Predictive oncology: A review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth , 2007, NeuroImage.

[21]  Randy Heiland,et al.  PhysiCell: An open source physics-based cell simulator for 3-D multicellular systems , 2017, bioRxiv.

[22]  Sotiris A. Prokopiou,et al.  IL-2 sensitivity and exogenous IL-2 concentration gradient tune the productive contact duration of CD8+ T cell-APC: a multiscale modeling study , 2016, BMC Systems Biology.

[23]  Ying Liu,et al.  Simultaneous detection of cell-secreted TNF-α and IFN-γ using micropatterned aptamer-modified electrodes. , 2012, Biomaterials.

[24]  Albert Goldbeter,et al.  An automaton model for the cell cycle , 2011, Interface Focus.

[25]  Malcolm S. Steinberg,et al.  Reconstruction of Tissues by Dissociated Cells , 1963 .

[26]  Claudine Chaouiya,et al.  EpiLog: A software for the logical modelling of epithelial dynamics , 2018, F1000Research.

[27]  Barbara Di Ventura,et al.  From in vivo to in silico biology and back , 2006, Nature.

[28]  V. Quaranta,et al.  Integrative mathematical oncology , 2008, Nature Reviews Cancer.

[29]  Vincent Hakim,et al.  Collective cell migration: a physics perspective , 2017, Reports on progress in physics. Physical Society.

[30]  Stefan Hoehme,et al.  Biomechanical and Nutrient Controls in the Growth of Mammalian Cell Populations , 2010 .

[31]  Vittorio Cristini,et al.  Patient-calibrated agent-based modelling of ductal carcinoma in situ (DCIS): from microscopic measurements to macroscopic predictions of clinical progression. , 2012, Journal of theoretical biology.

[32]  Thomas S Deisboeck,et al.  The effects of EGF-receptor density on multiscale tumor growth patterns. , 2005, Journal of theoretical biology.

[33]  A. Mogilner,et al.  Quantitative modeling in cell biology: what is it good for? , 2006, Developmental cell.

[34]  D. Drasdo,et al.  A single-cell-based model of tumor growth in vitro: monolayers and spheroids , 2005, Physical biology.

[35]  Weixiong Zhong,et al.  Reversible modulation of cell cycle kinetics in NIH/3T3 mouse fibroblasts by inducible overexpression of mitochondrial manganese superoxide dismutase. , 2004, Antioxidants & redox signaling.

[36]  Ulrich S. Schwarz,et al.  Dynamics of Cell Ensembles on Adhesive Micropatterns: Bridging the Gap between Single Cell Spreading and Collective Cell Migration , 2016, PLoS Comput. Biol..

[37]  Alexander G. Fletcher,et al.  Chaste: A test-driven approach to software development for biological modelling , 2009, Comput. Phys. Commun..

[38]  Simon Tanaka,et al.  Simulation Frameworks for Morphogenetic Problems , 2015, Comput..

[39]  Ignacio Ramis-Conde,et al.  Modeling the influence of the E-cadherin-beta-catenin pathway in cancer cell invasion: a multiscale approach. , 2008, Biophysical journal.

[40]  R. Sutherland.  Cell and environment interactions in tumor microregions: the multicell spheroid model. , 1988, Science.

[41]  Anna Haeger,et al.  Cell jamming: collective invasion of mesenchymal tumor cells imposed by tissue confinement. , 2014, Biochimica et biophysica acta.

[42]  J. Sharpe.  Computer modeling in developmental biology: growing today, essential tomorrow , 2017, Development.

[43]  Paul A. Bates,et al.  Matrix geometry determines optimal cancer cell migration strategy and modulates response to interventions , 2013, Nature Cell Biology.

[44]  J. Freyer,et al.  Regulation of growth saturation and development of necrosis in EMT6/Ro multicellular spheroids by the glucose and oxygen supply. , 1986, Cancer research.

[45]  Andre Levchenko,et al.  Transient IκB Kinase Activity Mediates Temporal NF-κB Dynamics in Response to a Wide Range of Tumor Necrosis Factor-α Doses* , 2006, Journal of Biological Chemistry.

[46]  Alexander R. A. Anderson,et al.  The Influence of Cellular Characteristics on the Evolution of Shape Homeostasis , 2017, Artificial Life.

[47]  Alissa M. Weaver,et al.  Tumor Morphology and Phenotypic Evolution Driven by Selective Pressure from the Microenvironment , 2006, Cell.

[48]  Eduardo Moreno,et al.  Tissue Crowding Induces Caspase-Dependent Competition for Space , 2016, Current Biology.

[49]  U. Schwarz,et al.  Cell organization in soft media due to active mechanosensing , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[50]  W. Rappel,et al.  Alignment of cellular motility forces with tissue flow as a mechanism for efficient wound healing , 2013, Proceedings of the National Academy of Sciences.

[51]  Jan Lammerding,et al.  Nuclear envelope rupture and repair during cancer cell migration , 2016, Science.

[52]  L. Kunz-Schughart,et al.  Multicellular tumor spheroids: an underestimated tool is catching up again. , 2010, Journal of biotechnology.

[53]  Vittorio Cristini,et al.  Multiscale cancer modeling. , 2010, Annual review of biomedical engineering.

[54]  Ryoichiro Kageyama,et al.  Oscillations in Notch Signaling Regulate Maintenance of Neural Progenitors , 2008, Neuron.

[55]  Samuel H. Friedman,et al.  BioFVM: an efficient, parallelized diffusive transport solver for 3-D biological simulations , 2015, Bioinform..

[56]  S. Hoehme,et al.  On the Role of Physics in the Growth and Pattern Formation of Multi-Cellular Systems: What can we Learn from Individual-Cell Based Models? , 2007 .

[57]  B. Frieden,et al.  Adaptive therapy. , 2009, Cancer research.

[58]  N. Novère.  Quantitative and logic modelling of molecular and gene networks , 2015, Nature Reviews Genetics.

[59]  Manuel Théry,et al.  The size-speed-force relationship governs migratory cell response to tumorigenic factors , 2017, Molecular biology of the cell.

[60]  H. Hirata,et al.  Oscillatory Expression of the bHLH Factor Hes1 Regulated by a Negative Feedback Loop , 2002, Science.

[61]  Ignacio Ramis-Conde,et al.  From genotypes to phenotypes: classification of the tumour profiles for different variants of the cadherin adhesion pathway , 2012, Physical biology.

[62]  Aurélien Naldi,et al.  MaBoSS 2.0: an environment for stochastic Boolean modeling , 2017, Bioinform..

[63]  Timothy K Lee,et al.  Single-cell NF-κB dynamics reveal digital activation and analogue information processing , 2010, Nature.

[64]  Dietmar W Hutmacher,et al.  A multiscale road map of cancer spheroids – incorporating experimental and mathematical modelling to understand cancer progression , 2013, Journal of Cell Science.

[65]  Julian Lewis.  Autoinhibition with Transcriptional Delay A Simple Mechanism for the Zebrafish Somitogenesis Oscillator , 2003, Current Biology.

[66]  Ignacio Ramis-Conde,et al.  Multi-scale modelling of the dynamics of cell colonies: insights into cell-adhesion forces and cancer invasion from in silico simulations , 2015, Journal of The Royal Society Interface.

[67]  J P Freyer,et al.  In situ oxygen consumption rates of cells in V‐79 multicellular spheroids during growth , 1984, Journal of cellular physiology.

[68]  Wouter-Jan Rappel,et al.  Velocity alignment leads to high persistence in confined cells. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[69]  Sotiris A. Prokopiou,et al.  Multiscale Modeling of the Early CD8 T-Cell Immune Response in Lymph Nodes: An Integrative Study , 2014, Comput..

[70]  M. S. Steinberg,et al.  Cadherin-mediated cell adhesion and tissue segregation: qualitative and quantitative determinants. , 2003, Developmental biology.

[71]  James A. Glazier,et al.  Virtual-tissue computer simulations define the roles of cell adhesion and proliferation in the onset of kidney cystic disease , 2016, Molecular biology of the cell.

[72]  A. Shaw,et al.  Tumour heterogeneity and resistance to cancer therapies , 2018, Nature Reviews Clinical Oncology.

[73]  William A. Thomas,et al.  Force measurements in E-cadherin–mediated cell doublets reveal rapid adhesion strengthened by actin cytoskeleton remodeling through Rac and Cdc42 , 2004, The Journal of cell biology.

[74]  Jacques Prost,et al.  Compressive stress inhibits proliferation in tumor spheroids through a volume limitation. , 2014, Biophysical journal.

[75]  P. Friedl,et al.  The Journal of Cell Biology , 2002 .

[76]  Emmanuel Barillot,et al.  Continuous time boolean modeling for biological signaling: application of Gillespie algorithm , 2012, BMC Systems Biology.

[77]  Bojana Gligorijevic,et al.  Multiparametric Classification Links Tumor Microenvironments with Tumor Cell Phenotype , 2014, PLoS biology.

[78]  M. Takeichi,et al.  Cadherin cell adhesion receptors as a morphogenetic regulator. , 1991, Science.

[79]  S Etienne-Manneville,et al.  Polarity proteins in migration and invasion , 2008, Oncogene.

[80]  Takuya Kato,et al.  A mechanically active heterotypic E-cadherin/N-cadherin adhesion enables fibroblasts to drive cancer cell invasion , 2017, Nature Cell Biology.

[81]  Andrei Zinovyev,et al.  Conceptual and computational framework for logical modelling of biological networks deregulated in diseases , 2019, Briefings Bioinform..

[82]  Mandar M. Inamdar,et al.  Proteolytic and non-proteolytic regulation of collective cell invasion: tuning by ECM density and organization , 2016, Scientific Reports.

[83]  W. Fiers,et al.  TNF-induced intracellular signaling leading to gene induction or to cytotoxicity by necrosis or by apoptosis. , 1995, Journal of inflammation.

[84]  Brian Seed,et al.  Fas triggers an alternative, caspase-8–independent cell death pathway using the kinase RIP as effector molecule , 2000, Nature Immunology.

[85]  C. Pao,et al.  Models of genetic control by repression with time delays and spatial effects , 1984, Journal of mathematical biology.

[86]  Alex Mogilner,et al.  Comparison of cell migration mechanical strategies in three-dimensional matrices: a computational study , 2016, Interface Focus.

[87]  Miquel Marin-Riera,et al.  Computational modeling of development by epithelia, mesenchyme and their interactions: a unified model , 2015, Bioinform..

[88]  Mohammad A. Qasaimeh,et al.  NF-κB signalling and cell fate decisions in response to a short pulse of tumour necrosis factor , 2016, Scientific Reports.

[89]  Scott A. Williams,et al.  Homo naledi, a new species of the genus Homo from the Dinaledi Chamber, South Africa , 2015, eLife.

[90]  Andreas Hellander,et al.  Orchestral: A Lightweight Framework for Parallel Simulations of Cell-Cell Communication , 2018, 2018 IEEE 14th International Conference on e-Science (e-Science).

[91]  Daniel A Fletcher,et al.  Tissue Geometry Determines Sites of Mammary Branching Morphogenesis in Organotypic Cultures , 2006, Science.

[92]  A. Anderson,et al.  Evolution of cell motility in an individual-based model of tumour growth. , 2009, Journal of theoretical biology.

[93]  François Nédélec.  preconfig: A Versatile Configuration File Generator for Varying Parameters , 2017 .

[94]  Tom Duke,et al.  A kinetic mechanism for cell sorting based on local variations in cell motility , 2014, Interface Focus.

[95]  Peter T. Cummings,et al.  An off-lattice hybrid discrete-continuum model of tumor growth and invasion. , 2010, Biophysical journal.

[96]  Justin M. Wozniak,et al.  High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow , 2017 .