1
|
Cogno N, Axenie C, Bauer R, Vavourakis V. Agent-based modeling in cancer biomedicine: applications and tools for calibration and validation. Cancer Biol Ther 2024; 25:2344600. [PMID: 38678381 PMCID: PMC11057625 DOI: 10.1080/15384047.2024.2344600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 04/15/2024] [Indexed: 04/29/2024] Open
Abstract
Computational models are not just appealing because they can simulate and predict the development of biological phenomena across multiple spatial and temporal scales, but also because they can integrate information from well-established in vitro and in vivo models and test new hypotheses in cancer biomedicine. Agent-based models and simulations are especially interesting candidates among computational modeling procedures in cancer research due to the capability to, for instance, recapitulate the dynamics of neoplasia and tumor - host interactions. Yet, the absence of methods to validate the consistency of the results across scales can hinder adoption by turning fine-tuned models into black boxes. This review compiles relevant literature that explores strategies to leverage high-fidelity simulations of multi-scale, or multi-level, cancer models with a focus on verification approached as simulation calibration. We consolidate our review with an outline of modern approaches for agent-based models' validation and provide an ambitious outlook toward rigorous and reliable calibration.
Collapse
Affiliation(s)
- Nicolò Cogno
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
- Institute for Condensed Matter Physics, Technische Universit¨at Darmstadt, Darmstadt, Germany
| | - Cristian Axenie
- Computer Science Department and Center for Artificial Intelligence, Technische Hochschule Nürnberg Georg Simon Ohm, Nuremberg, Germany
| | - Roman Bauer
- Nature Inspired Computing and Engineering Research Group, Computer Science Research Centre, University of Surrey, Guildford, UK
| | - Vasileios Vavourakis
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
| |
Collapse
|
2
|
Pinton P. Machine Learning for Predicting Biologic Agent Efficacy in Ulcerative Colitis: An Analysis for Generalizability and Combination with Computational Models. Diagnostics (Basel) 2024; 14:1324. [PMID: 39001215 PMCID: PMC11240677 DOI: 10.3390/diagnostics14131324] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Revised: 05/24/2024] [Accepted: 06/20/2024] [Indexed: 07/16/2024] Open
Abstract
Machine learning (ML) has been applied to predict the efficacy of biologic agents in ulcerative colitis (UC). ML can offer precision, personalization, efficiency, and automation. Moreover, it can improve decision support in predicting clinical outcomes. However, it faces challenges related to data quality and quantity, overfitting, generalization, and interpretability. This paper comments on two recent ML models that predict the efficacy of vedolizumab and ustekinumab in UC. Models that consider multiple pathways, multiple ethnicities, and combinations of real-world and clinical trial data are required for optimal shared decision-making and precision medicine. This paper also highlights the potential of combining ML with computational models to enhance clinical outcomes and personalized healthcare. Key Insights: (1) ML offers precision, personalization, efficiency, and decision support for predicting the efficacy of biologic agents in UC. (2) Challenging aspects in ML prediction include data quality, overfitting, and interpretability. (3) Multiple pathways, multiple ethnicities, and combinations of real-world and clinical trial data should be considered in predictive models for optimal decision-making. (4) Combining ML with computational models may improve clinical outcomes and personalized healthcare.
Collapse
Affiliation(s)
- Philippe Pinton
- Clinical and Translational Sciences, International PharmaScience Center Ferring Pharmaceuticals, Amager Strandvej 405, 2770 Kastrup, Denmark
| |
Collapse
|
3
|
Ferdous J, Matthew Fricke G, Moses ME. More Is Faster: Why Population Size Matters in Biological Search. J Comput Biol 2024; 31:429-444. [PMID: 38754139 DOI: 10.1089/cmb.2023.0296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024] Open
Abstract
Many biological scenarios have multiple cooperating searchers, and the timing of the initial first contact between any one of those searchers and its target is critically important. However, we are unaware of biological models that predict how long it takes for the first of many searchers to discover a target. We present a novel mathematical model that predicts initial first contact times between searchers and targets distributed at random in a volume. We compare this model with the extreme first passage time approach in physics that assumes an infinite number of searchers all initially positioned at the same location. We explore how the number of searchers, the distribution of searchers and targets, and the initial distances between searchers and targets affect initial first contact times. Given a constant density of uniformly distributed searchers and targets, the initial first contact time decreases linearly with both search volume and the number of searchers. However, given only a single target and searchers placed at the same starting location, the relationship between the initial first contact time and the number of searchers shifts from a linear decrease to a logarithmic decrease as the number of searchers grows very large. More generally, we show that initial first contact times can be dramatically faster than the average first contact times and that the initial first contact times decrease with the number of searchers, while the average search times are independent of the number of searchers. We suggest that this is an underappreciated phenomenon in biology and other collective search problems.
Collapse
Affiliation(s)
- Jannatul Ferdous
- Department of Computer Science, The University of New Mexico, Albuquerque, New Mexico, USA
| | - George Matthew Fricke
- Department of Computer Science, The University of New Mexico, Albuquerque, New Mexico, USA
- Center for Advanced Research Computing, The University of New Mexico, Albuquerque, New Mexico, USA
| | - Melanie E Moses
- Department of Computer Science, The University of New Mexico, Albuquerque, New Mexico, USA
- Department of Biology, The University of New Mexico, Albuquerque, New Mexico, USA
- Santa Fe Institute, Santa Fe, New Mexico, USA
| |
Collapse
|
4
|
Cogno N, Bauer R, Durante M. Mechanistic model of radiotherapy-induced lung fibrosis using coupled 3D agent-based and Monte Carlo simulations. COMMUNICATIONS MEDICINE 2024; 4:16. [PMID: 38336802 PMCID: PMC10858213 DOI: 10.1038/s43856-024-00442-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Accepted: 01/22/2024] [Indexed: 02/12/2024] Open
Abstract
BACKGROUND Mechanistic modelling of normal tissue toxicities is unfolding as an alternative to the phenomenological normal tissue complication probability models. The latter, currently used in the clinics, rely exclusively on limited patient data and neglect spatial dose distribution information. Among the various approaches, agent-based models are appealing as they provide the means to include patient-specific parameters and simulate long-term effects in complex systems. However, Monte Carlo tools remain the state-of-the-art for modelling radiation transport and provide measurements of the delivered dose with unmatched precision. METHODS In this work, we develop and characterize a coupled 3D agent-based - Monte Carlo model that mechanistically simulates the onset of the radiation-induced lung fibrosis in an alveolar segment. To the best of our knowledge, this is the first such model. RESULTS Our model replicates extracellular matrix patterns, radiation-induced lung fibrosis severity indexes and functional subunits survivals that show qualitative agreement with experimental studies and are consistent with our past results. Moreover, in accordance with experimental results, higher functional subunits survival and lower radiation-induced lung fibrosis severity indexes are achieved when a 5-fractions treatment is simulated. Finally, the model shows increased sensitivity to more uniform protons dose distributions with respect to more heterogeneous ones from photon irradiation. CONCLUSIONS This study lays thus the groundwork for further investigating the effects of different radiotherapeutic treatments on the onset of radiation-induced lung fibrosis via mechanistic modelling.
Collapse
Affiliation(s)
- Nicolò Cogno
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291, Darmstadt, Germany
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289, Darmstadt, Germany
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Roman Bauer
- Department of Computer Science, University of Surrey, Guildford, GU2 7XH, UK
| | - Marco Durante
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291, Darmstadt, Germany.
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289, Darmstadt, Germany.
- Department of Physics "Ettore Pancini", University Federico II, Naples, Italy.
| |
Collapse
|
5
|
Cockx BJR, Foster T, Clegg RJ, Alden K, Arya S, Stekel DJ, Smets BF, Kreft JU. Is it selfish to be filamentous in biofilms? Individual-based modeling links microbial growth strategies with morphology using the new and modular iDynoMiCS 2.0. PLoS Comput Biol 2024; 20:e1011303. [PMID: 38422165 PMCID: PMC10947719 DOI: 10.1371/journal.pcbi.1011303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Revised: 03/18/2024] [Accepted: 02/01/2024] [Indexed: 03/02/2024] Open
Abstract
Microbial communities are found in all habitable environments and often occur in assemblages with self-organized spatial structures developing over time. This complexity can only be understood, predicted, and managed by combining experiments with mathematical modeling. Individual-based models are particularly suited if individual heterogeneity, local interactions, and adaptive behavior are of interest. Here we present the completely overhauled software platform, the individual-based Dynamics of Microbial Communities Simulator, iDynoMiCS 2.0, which enables researchers to specify a range of different models without having to program. Key new features and improvements are: (1) Substantially enhanced ease of use (graphical user interface, editor for model specification, unit conversions, data analysis and visualization and more). (2) Increased performance and scalability enabling simulations of up to 10 million agents in 3D biofilms. (3) Kinetics can be specified with any arithmetic function. (4) Agent properties can be assembled from orthogonal modules for pick and mix flexibility. (5) Force-based mechanical interaction framework enabling attractive forces and non-spherical agent morphologies as an alternative to the shoving algorithm. The new iDynoMiCS 2.0 has undergone intensive testing, from unit tests to a suite of increasingly complex numerical tests and the standard Benchmark 3 based on nitrifying biofilms. A second test case was based on the "biofilms promote altruism" study previously implemented in BacSim because competition outcomes are highly sensitive to the developing spatial structures due to positive feedback between cooperative individuals. We extended this case study by adding morphology to find that (i) filamentous bacteria outcompete spherical bacteria regardless of growth strategy and (ii) non-cooperating filaments outcompete cooperating filaments because filaments can escape the stronger competition between themselves. In conclusion, the new substantially improved iDynoMiCS 2.0 joins a growing number of platforms for individual-based modeling of microbial communities with specific advantages and disadvantages that we discuss, giving users a wider choice.
Collapse
Affiliation(s)
- Bastiaan J. R. Cockx
- Department of Environmental and Resource Engineering, Technical University of Demark, DTU Lyngby campus, Kgs. Lyngby, Denmark
| | - Tim Foster
- Centre for Computational Biology & Institute of Microbiology and Infection & School of Biosciences, University of Birmingham, Edgbaston, Birmingham, United Kingdom
| | - Robert J. Clegg
- Centre for Computational Biology & Institute of Microbiology and Infection & School of Biosciences, University of Birmingham, Edgbaston, Birmingham, United Kingdom
| | - Kieran Alden
- Centre for Computational Biology & Institute of Microbiology and Infection & School of Biosciences, University of Birmingham, Edgbaston, Birmingham, United Kingdom
| | - Sankalp Arya
- School of Biosciences, University of Nottingham, Sutton Bonington Campus, Loughborough, Leicestershire, United Kingdom
| | - Dov J. Stekel
- School of Biosciences, University of Nottingham, Sutton Bonington Campus, Loughborough, Leicestershire, United Kingdom
| | - Barth F. Smets
- Department of Environmental and Resource Engineering, Technical University of Demark, DTU Lyngby campus, Kgs. Lyngby, Denmark
| | - Jan-Ulrich Kreft
- Centre for Computational Biology & Institute of Microbiology and Infection & School of Biosciences, University of Birmingham, Edgbaston, Birmingham, United Kingdom
| |
Collapse
|
6
|
Johnson JA, Stein-O’Brien GL, Booth M, Heiland R, Kurtoglu F, Bergman DR, Bucher E, Deshpande A, Forjaz A, Getz M, Godet I, Lyman M, Metzcar J, Mitchell J, Raddatz A, Rocha H, Solorzano J, Sundus A, Wang Y, Gilkes D, Kagohara LT, Kiemen AL, Thompson ED, Wirtz D, Wu PH, Zaidi N, Zheng L, Zimmerman JW, Jaffee EM, Hwan Chang Y, Coussens LM, Gray JW, Heiser LM, Fertig EJ, Macklin P. Digitize your Biology! Modeling multicellular systems through interpretable cell behavior. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.17.557982. [PMID: 37745323 PMCID: PMC10516032 DOI: 10.1101/2023.09.17.557982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/26/2023]
Abstract
Cells are fundamental units of life, constantly interacting and evolving as dynamical systems. While recent spatial multi-omics can quantitate individual cells' characteristics and regulatory programs, forecasting their evolution ultimately requires mathematical modeling. We develop a conceptual framework-a cell behavior hypothesis grammar-that uses natural language statements (cell rules) to create mathematical models. This allows us to systematically integrate biological knowledge and multi-omics data to make them computable. We can then perform virtual "thought experiments" that challenge and extend our understanding of multicellular systems, and ultimately generate new testable hypotheses. In this paper, we motivate and describe the grammar, provide a reference implementation, and demonstrate its potential through a series of examples in tumor biology and immunotherapy. Altogether, this approach provides a bridge between biological, clinical, and systems biology researchers for mathematical modeling of biological systems at scale, allowing the community to extrapolate from single-cell characterization to emergent multicellular behavior.
Collapse
Affiliation(s)
- Jeanette A.I. Johnson
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Genevieve L. Stein-O’Brien
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
- Department of Neuroscience, Johns Hopkins University. Baltimore, MD USA
| | - Max Booth
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
| | - Randy Heiland
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Furkan Kurtoglu
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Daniel R. Bergman
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Elmar Bucher
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Atul Deshpande
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - André Forjaz
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University. Baltimore, MD USA
| | - Michael Getz
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Ines Godet
- Memorial Sloan Kettering Cancer Center. New York, NY USA
| | - Melissa Lyman
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - John Metzcar
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
- Department of Informatics, Indiana University. Bloomington, IN USA
| | - Jacob Mitchell
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
- Department of Human Genetics, Johns Hopkins University. Baltimore, MD USA
| | - Andrew Raddatz
- Department of Biomedical Engineering, Georgia Institute of Technology, Emory University. Atlanta, GA USA
| | - Heber Rocha
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Jacobo Solorzano
- Centre de Recherches en Cancerologie de Toulouse. Toulouse, France
| | - Aneequa Sundus
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Yafei Wang
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| | - Danielle Gilkes
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
| | - Luciane T. Kagohara
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Ashley L. Kiemen
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
- Department of Pathology, Johns Hopkins University. Baltimore, MD USA
| | | | - Denis Wirtz
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University. Baltimore, MD USA
- Department of Pathology, Johns Hopkins University. Baltimore, MD USA
- Department of Materials Science and Engineering, Johns Hopkins University. Baltimore, MD USA
| | - Pei-Hsun Wu
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University. Baltimore, MD USA
| | - Neeha Zaidi
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Lei Zheng
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Jacquelyn W. Zimmerman
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Elizabeth M. Jaffee
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
| | - Young Hwan Chang
- Department of Biomedical Engineering, Oregon Health & Science University. Portland, OR USA
| | - Lisa M. Coussens
- Department of Cell, Developmental and Cancer Biology, Oregon Health & Science University. Portland, OR USA
| | - Joe W. Gray
- Department of Biomedical Engineering, Oregon Health & Science University. Portland, OR USA
| | - Laura M. Heiser
- Department of Biomedical Engineering, Oregon Health & Science University. Portland, OR USA
| | - Elana J. Fertig
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University. Baltimore, MD USA
- Convergence Institute, Johns Hopkins University. Baltimore, MD USA
- Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, MD, USA
- Department of Biomedical Engineering, Johns Hopkins University. Baltimore, MD USA
| | - Paul Macklin
- Department of Intelligent Systems Engineering, Indiana University. Bloomington, IN USA
| |
Collapse
|
7
|
Palmer BJ, Almgren AS, Johnson CGM, Myers AT, Cannon WR. BMX: Biological modelling and interface exchange. Sci Rep 2023; 13:12235. [PMID: 37507417 PMCID: PMC10382537 DOI: 10.1038/s41598-023-39150-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 07/20/2023] [Indexed: 07/30/2023] Open
Abstract
High performance computing has a great potential to provide a range of significant benefits for investigating biological systems. These systems often present large modelling problems with many coupled subsystems, such as when studying colonies of bacteria cells. The aim to understand cell colonies has generated substantial interest as they can have strong economic and societal impacts through their roles in in industrial bioreactors and complex community structures, called biofilms, found in clinical settings. Investigating these communities through realistic models can rapidly exceed the capabilities of current serial software. Here, we introduce BMX, a software system developed for the high performance modelling of large cell communities by utilising GPU acceleration. BMX builds upon the AMRex adaptive mesh refinement package to efficiently model cell colony formation under realistic laboratory conditions. Using simple test scenarios with varying nutrient availability, we show that BMX is capable of correctly reproducing observed behavior of bacterial colonies on realistic time scales demonstrating a potential application of high performance computing to colony modelling. The open source software is available from the zenodo repository https://doi.org/10.5281/zenodo.8084270 under the BSD-2-Clause licence.
Collapse
Affiliation(s)
- Bruce J Palmer
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Washington, USA
| | - Ann S Almgren
- Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Connah G M Johnson
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Washington, USA.
| | - Andrew T Myers
- Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - William R Cannon
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Washington, USA
| |
Collapse
|
8
|
De Schutter E. Efficient simulation of neural development using shared memory parallelization. Front Neuroinform 2023; 17:1212384. [PMID: 37547492 PMCID: PMC10400717 DOI: 10.3389/fninf.2023.1212384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2023] [Accepted: 06/30/2023] [Indexed: 08/08/2023] Open
Abstract
The Neural Development Simulator, NeuroDevSim, is a Python module that simulates the most important aspects of brain development: morphological growth, migration, and pruning. It uses an agent-based modeling approach inherited from the NeuroMaC software. Each cycle has agents called fronts execute model-specific code. In the case of a growing dendritic or axonal front, this will be a choice between extension, branching, or growth termination. Somatic fronts can migrate to new positions and any front can be retracted to prune parts of neurons. Collision detection prevents new or migrating fronts from overlapping with existing ones. NeuroDevSim is a multi-core program that uses an innovative shared memory approach to achieve parallel processing without messaging. We demonstrate linear strong parallel scaling up to 96 cores for large models and have run these successfully on 128 cores. Most of the shared memory parallelism is achieved without memory locking. Instead, cores have only write privileges to private sections of arrays, while being able to read the entire shared array. Memory conflicts are avoided by a coding rule that allows only active fronts to use methods that need writing access. The exception is collision detection, which is needed to avoid the growth of physically overlapping structures. For collision detection, a memory-locking mechanism was necessary to control access to grid points that register the location of nearby fronts. A custom approach using a serialized lock broker was able to manage both read and write locking. NeuroDevSim allows easy modeling of most aspects of neural development for models simulating a few complex or thousands of simple neurons or a mixture of both. Code available at https://github.com/CNS-OIST/NeuroDevSim.
Collapse
Affiliation(s)
- Erik De Schutter
- Computational Neuroscience Unit, Okinawa Institute of Science and Technology, Okinawa, Japan
- Department of Biomedical Sciences, University of Antwerp, Antwerpen, Belgium
| |
Collapse
|
9
|
Du J, Zhou Y, Jin L, Sheng K. Gell: A GPU-powered 3D hybrid simulator for large-scale multicellular system. PLoS One 2023; 18:e0288721. [PMID: 37463167 DOI: 10.1371/journal.pone.0288721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2023] [Accepted: 07/03/2023] [Indexed: 07/20/2023] Open
Abstract
As a powerful but computationally intensive method, hybrid computational models study the dynamics of multicellular systems by evolving discrete cells in reacting and diffusing extracellular microenvironments. As the scale and complexity of studied biological systems continuously increase, the exploding computational cost starts to limit large-scale cell-based simulations. To facilitate the large-scale hybrid computational simulation and make it feasible on easily accessible computational devices, we develop Gell (GPU Cell), a fast and memory-efficient open-source GPU-based hybrid computational modeling platform for large-scale system modeling. We fully parallelize the simulations on GPU for high computational efficiency and propose a novel voxel sorting method to further accelerate the modeling of massive cell-cell mechanical interaction with negligible additional memory footprint. As a result, Gell efficiently handles simulations involving tens of millions of cells on a personal computer. We compare the performance of Gell with a state-of-the-art paralleled CPU-based simulator on a hanging droplet spheroid growth task and further demonstrate Gell with a ductal carcinoma in situ (DCIS) simulation. Gell affords ~150X acceleration over the paralleled CPU method with one-tenth of the memory requirement.
Collapse
Affiliation(s)
- Jiayi Du
- Department of Radiation Oncology, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Yu Zhou
- Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Lihua Jin
- Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Ke Sheng
- Department of Radiation Oncology, University of California, San Francisco, San Francisco, California, United States of America
| |
Collapse
|
10
|
de Montigny J, Sernagor E, Bauer R. Retinal self-organization: a model of retinal ganglion cells and starburst amacrine cells mosaic formation. Open Biol 2023; 13:220217. [PMID: 37015288 PMCID: PMC10072945 DOI: 10.1098/rsob.220217] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/06/2023] Open
Abstract
Individual retinal cell types exhibit semi-regular spatial patterns called retinal mosaics. Retinal ganglion cells (RGCs) and starburst amacrine cells (SACs) are known to exhibit such layouts. Mechanisms responsible for the formation of mosaics are not well understood but follow three main principles: (i) homotypic cells prevent nearby cells from adopting the same type, (ii) cell tangential migration and (iii) cell death. Alongside experiments in mouse, we use BioDynaMo, an agent-based simulation framework, to build a detailed and mechanistic model of mosaic formation. We investigate the implications of the three theories for RGC's mosaic formation. We report that the cell migration mechanism yields the most regular mosaics. In addition, we propose that low-density RGC type mosaics exhibit on average low regularities, and thus we question the relevance of regular spacing as a criterion for a group of RGCs to form a RGC type. We investigate SAC mosaics formation and interactions between the ganglion cell layer (GCL) and inner nuclear layer (INL) populations. We propose that homotypic interactions between the GCL and INL populations during mosaics creation are required to reproduce the observed SAC mosaics' characteristics. This suggests that the GCL and INL populations of SACs might not be independent during retinal development.
Collapse
Affiliation(s)
- Jean de Montigny
- Biosciences Institute, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | - Evelyne Sernagor
- Biosciences Institute, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | - Roman Bauer
- Department of Computer Science, University of Surrey, Guildford GU2 7XH, UK
| |
Collapse
|
11
|
Gonçalves IG, García-Aznar JM. Hybrid computational models of multicellular tumour growth considering glucose metabolism. Comput Struct Biotechnol J 2023; 21:1262-1271. [PMID: 36814723 PMCID: PMC9939553 DOI: 10.1016/j.csbj.2023.01.044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 01/30/2023] [Accepted: 01/30/2023] [Indexed: 02/04/2023] Open
Abstract
Cancer cells metabolize glucose through metabolic pathways that differ from those used by healthy and differentiated cells. In particular, tumours have been shown to consume more glucose than their healthy counterparts and to use anaerobic metabolic pathways, even under aerobic conditions. Nevertheless, scientists have still not been able to explain why cancer cells evolved to present an altered metabolism and what evolutionary advantage this might provide them. Experimental and computational models have been increasingly used in recent years to understand some of these biological questions. Multicellular tumour spheroids are effective experimental models as they replicate the initial stages of avascular solid tumour growth. Furthermore, these experiments generate data which can be used to calibrate and validate computational studies that aim to simulate tumour growth. Hybrid models are of particular relevance in this field of research because they model cells as individual agents while also incorporating continuum representations of the substances present in the surrounding microenvironment that may participate in intracellular metabolic networks as concentration or density distributions. Henceforth, in this review, we explore the potential of computational modelling to reveal the role of metabolic reprogramming in tumour growth.
Collapse
Key Words
- ABM, agent-based model
- ATP, adenosine triphosphate
- CA, cellular automata
- CPM, cellular Potts model
- ECM, extracellular matrix
- FBA, Flux Balance Analysis
- FDG-PET, [18F]-fluorodeoxyglucose-positron emission tomography
- MCTS, multicellular tumour spheroids
- ODEs, ordinary differential equations
- PDEs, partial differential equations
- SBML, Systems Biology Markup Language
- Warburg effect
- agent-based models
- glucose metabolism
- hybrid modelling
- multicellular simulations
Collapse
Affiliation(s)
- Inês G Gonçalves
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, Aragon Institute of Engineering Research (I3A), University of Zaragoza, Zaragoza 50018, Aragon, Spain
| | - José Manuel García-Aznar
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, Aragon Institute of Engineering Research (I3A), University of Zaragoza, Zaragoza 50018, Aragon, Spain
| |
Collapse
|
12
|
Cogno N, Bauer R, Durante M. An Agent-Based Model of Radiation-Induced Lung Fibrosis. Int J Mol Sci 2022; 23:13920. [PMID: 36430398 PMCID: PMC9693125 DOI: 10.3390/ijms232213920] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2022] [Revised: 11/03/2022] [Accepted: 11/05/2022] [Indexed: 11/16/2022] Open
Abstract
Early- and late-phase radiation-induced lung injuries, namely pneumonitis and lung fibrosis (RILF), severely constrain the maximum dose and irradiated volume in thoracic radiotherapy. As the most radiosensitive targets, epithelial cells respond to radiation either by undergoing apoptosis or switching to a senescent phenotype that triggers the immune system and damages surrounding healthy cells. Unresolved inflammation stimulates mesenchymal cells' proliferation and extracellular matrix (ECM) secretion, which irreversibly stiffens the alveolar walls and leads to respiratory failure. Although a thorough understanding is lacking, RILF and idiopathic pulmonary fibrosis share multiple pathways and would mutually benefit from further insights into disease progression. Furthermore, current normal tissue complication probability (NTCP) models rely on clinical experience to set tolerance doses for organs at risk and leave aside mechanistic interpretations of the undergoing processes. To these aims, we implemented a 3D agent-based model (ABM) of an alveolar duct that simulates cell dynamics and substance diffusion following radiation injury. Emphasis was placed on cell repopulation, senescent clearance, and intra/inter-alveolar bystander senescence while tracking ECM deposition. Our ABM successfully replicates early and late fibrotic response patterns reported in the literature along with the ECM sigmoidal dose-response curve. Moreover, surrogate measures of RILF severity via a custom indicator show qualitative agreement with published fibrosis indices. Finally, our ABM provides a fully mechanistic alveolar survival curve highlighting the need to include bystander damage in lung NTCP models.
Collapse
Affiliation(s)
- Nicolò Cogno
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291 Darmstadt, Germany
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289 Darmstadt, Germany
| | - Roman Bauer
- Department of Computer Science, University of Surrey, Guildford GU2 7XH, UK
| | - Marco Durante
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291 Darmstadt, Germany
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289 Darmstadt, Germany
| |
Collapse
|
13
|
Senthilkumar I, Howley E, McEvoy E. Thermodynamically-motivated chemo-mechanical models and multicellular simulation to provide new insight into active cell and tumour remodelling. Exp Cell Res 2022; 419:113317. [PMID: 36028058 DOI: 10.1016/j.yexcr.2022.113317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Revised: 07/19/2022] [Accepted: 08/08/2022] [Indexed: 11/25/2022]
Abstract
Computational models can shape our understanding of cell and tissue remodelling, from cell spreading, to active force generation, adhesion, and growth. In this mini-review, we discuss recent progress in modelling of chemo-mechanical cell behaviour and the evolution of multicellular systems. In particular, we highlight recent advances in (i) free-energy based single cell models that can provide new fundamental insight into cell spreading, cancer cell invasion, stem cell differentiation, and remodelling in disease, and (ii) mechanical agent-based models to simulate large numbers of discrete interacting cells in proliferative tumours. We describe how new biological understanding has emerged from such theoretical models, and the trade-offs and constraints associated with current approaches. Ultimately, we aim to make a case for why theory should be integrated with an experimental workflow to optimise new in-vitro studies, to predict feedback between cells and their microenvironment, and to deepen understanding of active cell behaviour.
Collapse
Affiliation(s)
- Irish Senthilkumar
- School of Computer Science, College of Science and Engineering, National University of Ireland Galway, Ireland; Biomedical Engineering, College of Science and Engineering, National University of Ireland Galway, Ireland
| | - Enda Howley
- School of Computer Science, College of Science and Engineering, National University of Ireland Galway, Ireland
| | - Eoin McEvoy
- Biomedical Engineering, College of Science and Engineering, National University of Ireland Galway, Ireland.
| |
Collapse
|
14
|
Arjoca S, Robu A, Neagu M, Neagu A. Mathematical and computational models in spheroid-based biofabrication. Acta Biomater 2022:S1742-7061(22)00418-4. [PMID: 35853599 DOI: 10.1016/j.actbio.2022.07.024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Revised: 06/25/2022] [Accepted: 07/12/2022] [Indexed: 11/01/2022]
Abstract
Ubiquitous in embryonic development, tissue fusion is of interest to tissue engineers who use tissue spheroids or organoids as building blocks of three-dimensional (3D) multicellular constructs. This review presents mathematical models and computer simulations of the fusion of tissue spheroids. The motivation of this study stems from the need to predict the post-printing evolution of 3D bioprinted constructs. First, we provide a brief overview of differential adhesion, the main morphogenetic mechanism involved in post-printing structure formation. It will be shown that clusters of cohesive cells behave as an incompressible viscous fluid on the time scale of hours. The discussion turns then to mathematical models based on the continuum hydrodynamics of highly viscous liquids and on statistical mechanics. Next, we analyze the validity and practical use of computational models of multicellular self-assembly in live constructs created by tissue spheroid bioprinting. Finally, we discuss the perspectives of the field as machine learning starts to reshape experimental design, and modular robotic workstations tend to alleviate the burden of repetitive tasks in biofabrication. STATEMENT OF SIGNIFICANCE: Bioprinted constructs are living systems, which evolve via morphogenetic mechanisms known from developmental biology. This review presents mathematical and computational tools devised for modeling post-printing structure formation. They help achieving a desirable outcome without expensive optimization experiments. While previous reviews mainly focused on assumptions, technical details, strengths, and limitations of computational models of multicellular self-assembly, this article discusses their validity and practical use in biofabrication. It also presents an overview of mathematical models that proved to be useful in the evaluation of experimental data on tissue spheroid fusion, and in the calibration of computational models. Finally, the perspectives of the field are discussed in the advent of robotic biofabrication platforms and bioprinting process optimization by machine learning.
Collapse
Affiliation(s)
- Stelian Arjoca
- Center for Modeling Biological Systems and Data Analysis, Department of Functional Sciences, Victor Babes University of Medicine and Pharmacy Timisoara, Piata Eftimie Murgu Nr. 2-4, Timisoara 300041, Romania
| | - Andreea Robu
- Department of Automation and Applied Informatics, Politehnica University of Timisoara, Timisoara 300006, Romania
| | - Monica Neagu
- Center for Modeling Biological Systems and Data Analysis, Department of Functional Sciences, Victor Babes University of Medicine and Pharmacy Timisoara, Piata Eftimie Murgu Nr. 2-4, Timisoara 300041, Romania
| | - Adrian Neagu
- Center for Modeling Biological Systems and Data Analysis, Department of Functional Sciences, Victor Babes University of Medicine and Pharmacy Timisoara, Piata Eftimie Murgu Nr. 2-4, Timisoara 300041, Romania; Department of Physics & Astronomy, University of Missouri-Columbia, Columbia, MO 65211, USA.
| |
Collapse
|
15
|
Ofițeru ID, Picioreanu C. No model is perfect, but some are useful. Science 2022; 376:914-916. [PMID: 35617381 DOI: 10.1126/science.abq0956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Agent-based model should inform the action plan to curb algal blooms in Lake Erie.
Collapse
Affiliation(s)
- Irina D Ofițeru
- School of Engineering, Newcastle University, Merz Court, Newcastle upon Tyne NE1 7RU, UK
| | - Cristian Picioreanu
- Water Desalination and Reuse Center, Biological and Environmental Sciences and Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
16
|
Demetriades M, Zivanovic M, Hadjicharalambous M, Ioannou E, Ljujic B, Vucicevic K, Ivosevic Z, Dagovic A, Milivojevic N, Kokkinos O, Bauer R, Vavourakis V. Interrogating and Quantifying In Vitro Cancer Drug Pharmacodynamics via Agent-Based and Bayesian Monte Carlo Modelling. Pharmaceutics 2022; 14:749. [PMID: 35456583 PMCID: PMC9029523 DOI: 10.3390/pharmaceutics14040749] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Revised: 03/25/2022] [Accepted: 03/28/2022] [Indexed: 02/01/2023] Open
Abstract
The effectiveness of chemotherapy in cancer cell regression is often limited by drug resistance, toxicity, and neoplasia heterogeneity. However, due to the significant complexities entailed by the many cancer growth processes, predicting the impact of interference and symmetry-breaking mechanisms is a difficult problem. To quantify and understand more about cancer drug pharmacodynamics, we combine in vitro with in silico cancer models. The anti-proliferative action of selected cytostatics is interrogated on human colorectal and breast adenocarcinoma cells, while an agent-based computational model is employed to reproduce experiments and shed light on the main therapeutic mechanisms of each chemotherapeutic agent. Multiple drug administration scenarios on each cancer cell line are simulated by varying the drug concentration, while a Bayesian-based method for model parameter optimisation is employed. Our proposed procedure of combining in vitro cancer drug screening with an in silico agent-based model successfully reproduces the impact of chemotherapeutic drugs in cancer growth behaviour, while the mechanisms of action of each drug are characterised through model-derived probabilities of cell apoptosis and division. We suggest that our approach could form the basis for the prospective generation of experimentally-derived and model-optimised pharmacological variables towards personalised cancer therapy.
Collapse
Affiliation(s)
- Marios Demetriades
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Marko Zivanovic
- Department of Science, Institute for Information Technologies Kragujevac, University of Kragujevac, 34000 Kragujevac, Serbia; (M.Z.); (N.M.)
| | - Myrianthi Hadjicharalambous
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Eleftherios Ioannou
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Biljana Ljujic
- Faculty of Medical Sciences, Human Genetics, University of Kragujevac, 34000 Kragujevac, Serbia; (B.L.); (Z.I.)
| | - Ksenija Vucicevic
- Department for Pharmaceutical Technologies, Faculty of Medical Sciences, University of Kragujevac, 34000 Kragujevac, Serbia;
| | - Zeljko Ivosevic
- Faculty of Medical Sciences, Human Genetics, University of Kragujevac, 34000 Kragujevac, Serbia; (B.L.); (Z.I.)
| | - Aleksandar Dagovic
- Oncology and Radiotherapy Centre, Faculty of Medical Sciences, University of Kragujevac, 34000 Kragujevac, Serbia;
| | - Nevena Milivojevic
- Department of Science, Institute for Information Technologies Kragujevac, University of Kragujevac, 34000 Kragujevac, Serbia; (M.Z.); (N.M.)
| | - Odysseas Kokkinos
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Roman Bauer
- Department of Computer Science, University of Surrey, Guilford GU2 7XH, UK
| | - Vasileios Vavourakis
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
- Department of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, UK
| |
Collapse
|
17
|
Abstract
Understanding the pathophysiology of lung fibrosis is of paramount importance to elaborate targeted and effective therapies. As it onsets, the randomly accumulating extracellular matrix (ECM) breaks the symmetry of the branching lung structure. Interestingly, similar pathways have been reported for both idiopathic pulmonary fibrosis and radiation-induced lung fibrosis (RILF). Individuals suffering from the disease, the worldwide incidence of which is growing, have poor prognosis and a short mean survival time. In this context, mathematical and computational models have the potential to shed light on key underlying pathological mechanisms, shorten the time needed for clinical trials, parallelize hypotheses testing, and improve personalized drug development. Agent-based modeling (ABM) has proven to be a reliable and versatile simulation tool, whose features make it a good candidate for recapitulating emergent behaviors in heterogeneous systems, such as those found at multiple scales in the human body. In this paper, we detail the implementation of a 3D agent-based model of lung fibrosis using a novel simulation platform, namely, BioDynaMo, and prove that it can qualitatively and quantitatively reproduce published results. Furthermore, we provide additional insights on late-fibrosis patterns through ECM density distribution histograms. The model recapitulates key intercellular mechanisms, while cell numbers and types are embodied by alveolar segments that act as agents and are spatially arranged by a custom algorithm. Finally, our model may hold potential for future applications in the context of lung disorders, ranging from RILF (by implementing radiation-induced cell damage mechanisms) to COVID-19 and inflammatory diseases (such as asthma or chronic obstructive pulmonary disease).
Collapse
|