1
|
Montefusco A, Helfmann L, Okunola T, Winkelmann S, Schütte C. Partial mean-field model for neurotransmission dynamics. Math Biosci 2024; 369:109143. [PMID: 38220067 DOI: 10.1016/j.mbs.2024.109143] [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: 07/04/2023] [Revised: 12/07/2023] [Accepted: 01/09/2024] [Indexed: 01/16/2024]
Abstract
This article addresses reaction networks in which spatial and stochastic effects are of crucial importance. For such systems, particle-based models allow us to describe all microscopic details with high accuracy. However, they suffer from computational inefficiency if particle numbers and density get too large. Alternative coarse-grained-resolution models reduce computational effort tremendously, e.g., by replacing the particle distribution by a continuous concentration field governed by reaction-diffusion PDEs. We demonstrate how models on the different resolution levels can be combined into hybrid models that seamlessly combine the best of both worlds, describing molecular species with large copy numbers by macroscopic equations with spatial resolution while keeping the spatial-stochastic particle-based resolution level for the species with low copy numbers. To this end, we introduce a simple particle-based model for the binding dynamics of ions and vesicles at the heart of the neurotransmission process. Within this framework, we derive a novel hybrid model and present results from numerical experiments which demonstrate that the hybrid model allows for an accurate approximation of the full particle-based model in realistic scenarios.
Collapse
Affiliation(s)
- Alberto Montefusco
- Mathematics of Complex Systems, Zuse-Institut Berlin, Takustraße 7, Berlin, 14195, Germany
| | - Luzie Helfmann
- Mathematics of Complex Systems, Zuse-Institut Berlin, Takustraße 7, Berlin, 14195, Germany
| | - Toluwani Okunola
- Mathematics of Complex Systems, Zuse-Institut Berlin, Takustraße 7, Berlin, 14195, Germany; Institute Of Mathematics, Technische Universität Berlin, Straße des 17. Juni 136, Berlin, 10623, Germany
| | - Stefanie Winkelmann
- Mathematics of Complex Systems, Zuse-Institut Berlin, Takustraße 7, Berlin, 14195, Germany.
| | - Christof Schütte
- Mathematics of Complex Systems, Zuse-Institut Berlin, Takustraße 7, Berlin, 14195, Germany; Institute of Mathematics, Freie Universität Berlin, Arnimallee 6, Berlin, 14195, Germany
| |
Collapse
|
2
|
Smith ER, Theodorakis PE. Multiscale simulation of fluids: coupling molecular and continuum. Phys Chem Chem Phys 2024; 26:724-744. [PMID: 38113114 DOI: 10.1039/d3cp03579d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.
Collapse
Affiliation(s)
- Edward R Smith
- Department of Mechanical and Aerospace Engineering, Brunel University London, Uxbridge, Middlesex UB8 3PH, UK.
| | | |
Collapse
|
3
|
Liang Q, Yang J, Fan WTL, Lo WC. Patch formation driven by stochastic effects of interaction between viruses and defective interfering particles. PLoS Comput Biol 2023; 19:e1011513. [PMID: 37782667 PMCID: PMC10569632 DOI: 10.1371/journal.pcbi.1011513] [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: 02/07/2023] [Revised: 10/12/2023] [Accepted: 09/12/2023] [Indexed: 10/04/2023] Open
Abstract
Defective interfering particles (DIPs) are virus-like particles that occur naturally during virus infections. These particles are defective, lacking essential genetic materials for replication, but they can interact with the wild-type virus and potentially be used as therapeutic agents. However, the effect of DIPs on infection spread is still unclear due to complicated stochastic effects and nonlinear spatial dynamics. In this work, we develop a model with a new hybrid method to study the spatial-temporal dynamics of viruses and DIPs co-infections within hosts. We present two different scenarios of virus production and compare the results from deterministic and stochastic models to demonstrate how the stochastic effect is involved in the spatial dynamics of virus transmission. We compare the spread features of the virus in simulations and experiments, including the formation and the speed of virus spread and the emergence of stochastic patchy patterns of virus distribution. Our simulations simultaneously capture observed spatial spread features in the experimental data, including the spread rate of the virus and its patchiness. The results demonstrate that DIPs can slow down the growth of virus particles and make the spread of the virus more patchy.
Collapse
Affiliation(s)
- Qiantong Liang
- Department of Mathematics, City University of Hong Kong, Hong Kong, China
| | - Johnny Yang
- Department of Mathematics, Indiana University, Bloomington, Indiana, United States of America
| | - Wai-Tong Louis Fan
- Department of Mathematics, Indiana University, Bloomington, Indiana, United States of America
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
| | - Wing-Cheong Lo
- Department of Mathematics, City University of Hong Kong, Hong Kong, China
| |
Collapse
|
4
|
Jørgensen ACS, Hill CS, Sturrock M, Tang W, Karamched SR, Gorup D, Lythgoe MF, Parrinello S, Marguerat S, Shahrezaei V. Data-driven spatio-temporal modelling of glioblastoma. ROYAL SOCIETY OPEN SCIENCE 2023; 10:221444. [PMID: 36968241 PMCID: PMC10031411 DOI: 10.1098/rsos.221444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 02/23/2023] [Indexed: 06/18/2023]
Abstract
Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.
Collapse
Affiliation(s)
| | - Ciaran Scott Hill
- Department of Neurosurgery, The National Hospital for Neurology and Neurosurgery, London WC1N 3BG, UK
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Marc Sturrock
- Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, Dublin D02 YN77, Ireland
| | - Wenhao Tang
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| | - Saketh R. Karamched
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Dunja Gorup
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Mark F. Lythgoe
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Simona Parrinello
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Samuel Marguerat
- Genomics Translational Technology Platform, UCL Cancer Institute, University College London, London WC1E 6DD, UK
| | - Vahid Shahrezaei
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| |
Collapse
|
5
|
Ruiz-Martinez A, Gong C, Wang H, Sové RJ, Mi H, Kimko H, Popel AS. Simulations of tumor growth and response to immunotherapy by coupling a spatial agent-based model with a whole-patient quantitative systems pharmacology model. PLoS Comput Biol 2022; 18:e1010254. [PMID: 35867773 PMCID: PMC9348712 DOI: 10.1371/journal.pcbi.1010254] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [Revised: 08/03/2022] [Accepted: 05/26/2022] [Indexed: 12/23/2022] Open
Abstract
Quantitative systems pharmacology (QSP) models and spatial agent-based models (ABM) are powerful and efficient approaches for the analysis of biological systems and for clinical applications. Although QSP models are becoming essential in discovering predictive biomarkers and developing combination therapies through in silico virtual trials, they are inadequate to capture the spatial heterogeneity and randomness that characterize complex biological systems, and specifically the tumor microenvironment. Here, we extend our recently developed spatial QSP (spQSP) model to analyze tumor growth dynamics and its response to immunotherapy at different spatio-temporal scales. In the model, the tumor spatial dynamics is governed by the ABM, coupled to the QSP model, which includes the following compartments: central (blood system), tumor, tumor-draining lymph node, and peripheral (the rest of the organs and tissues). A dynamic recruitment of T cells and myeloid-derived suppressor cells (MDSC) from the QSP central compartment has been implemented as a function of the spatial distribution of cancer cells. The proposed QSP-ABM coupling methodology enables the spQSP model to perform as a coarse-grained model at the whole-tumor scale and as an agent-based model at the regions of interest (ROIs) scale. Thus, we exploit the spQSP model potential to characterize tumor growth, identify T cell hotspots, and perform qualitative and quantitative descriptions of cell density profiles at the invasive front of the tumor. Additionally, we analyze the effects of immunotherapy at both whole-tumor and ROI scales under different tumor growth and immune response conditions. A digital pathology computational analysis of triple-negative breast cancer specimens is used as a guide for modeling the immuno-architecture of the invasive front.
Collapse
Affiliation(s)
- Alvaro Ruiz-Martinez
- Department of Biomedical Engineering, Johns Hopkins, University School of Medicine, Baltimore, Maryland, United States of America
| | - Chang Gong
- Department of Biomedical Engineering, Johns Hopkins, University School of Medicine, Baltimore, Maryland, United States of America
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins, University School of Medicine, Baltimore, Maryland, United States of America
| | - Richard J. Sové
- Department of Biomedical Engineering, Johns Hopkins, University School of Medicine, Baltimore, Maryland, United States of America
| | - Haoyang Mi
- Department of Biomedical Engineering, Johns Hopkins, University School of Medicine, Baltimore, Maryland, United States of America
| | - Holly Kimko
- Clinical Pharmacology & Quantitative Pharmacology, AstraZeneca, Gaithersburg, Maryland, United States of America
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins, University School of Medicine, Baltimore, Maryland, United States of America
- Department of Oncology and Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University, Baltimore, Maryland, United States of America
| |
Collapse
|
6
|
Lee D, Green A, Wu H, Kwon JS. Hybrid
PDE‐kMC
modeling approach to simulate multivalent lectin‐glycan binding process. AIChE J 2021. [DOI: 10.1002/aic.17453] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Dongheon Lee
- Department of Biomedical Engineering Duke University Durham North Carolina USA
| | - Aaron Green
- Artie McFerrin Department of Chemical Engineering Texas A&M University Texas USA
| | - Hung‐Jen Wu
- Artie McFerrin Department of Chemical Engineering Texas A&M University Texas USA
| | - Joseph Sang‐Il Kwon
- Artie McFerrin Department of Chemical Engineering Texas A&M University Texas USA
| |
Collapse
|
7
|
Sego TJ, Aponte-Serrano JO, Gianlupi JF, Glazier JA. Generation of multicellular spatiotemporal models of population dynamics from ordinary differential equations, with applications in viral infection. BMC Biol 2021; 19:196. [PMID: 34496857 PMCID: PMC8424622 DOI: 10.1186/s12915-021-01115-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2021] [Accepted: 08/02/2021] [Indexed: 01/10/2023] Open
Abstract
BACKGROUND The biophysics of an organism span multiple scales from subcellular to organismal and include processes characterized by spatial properties, such as the diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales, through the generation of representative models. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent the spatial and stochastic features of a biological system, limiting their insights and applications. However, spatial models describing biological systems with spatial information are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them and highlights the need for simpler methods able to model the spatial features of biological systems. RESULTS In this work, we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models, the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates of extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters. We additionally investigate objects and processes implicitly represented by ODE model terms and parameters and improve the reproducibility of spatial, stochastic models. CONCLUSION We developed and demonstrate a method for generating spatiotemporal, multicellular models from non-spatial population dynamics models of multicellular systems. We envision employing our method to generate new ODE model terms from spatiotemporal and multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes.
Collapse
Affiliation(s)
- T J Sego
- Department of Intelligent Systems Engineering and Biocomplexity Institute, Indiana University, Bloomington, IN, USA.
| | - Josua O Aponte-Serrano
- Department of Intelligent Systems Engineering and Biocomplexity Institute, Indiana University, Bloomington, IN, USA
| | - Juliano F Gianlupi
- Department of Intelligent Systems Engineering and Biocomplexity Institute, Indiana University, Bloomington, IN, USA
| | - James A Glazier
- Department of Intelligent Systems Engineering and Biocomplexity Institute, Indiana University, Bloomington, IN, USA
| |
Collapse
|
8
|
Murphy RJ, Buenzli PR, Tambyah TA, Thompson EW, Hugo HJ, Baker RE, Simpson MJ. The role of mechanical interactions in EMT. Phys Biol 2021; 18. [PMID: 33789261 DOI: 10.1088/1478-3975/abf425] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 03/31/2021] [Indexed: 02/07/2023]
Abstract
The detachment of cells from the boundary of an epithelial tissue and the subsequent invasion of these cells into surrounding tissues is important for cancer development and wound healing, and is strongly associated with the epithelial-mesenchymal transition (EMT). Chemical signals, such as TGF-β, produced by surrounding tissue can be uptaken by cells and induce EMT. In this work, we present a novel cell-based discrete mathematical model of mechanical cellular relaxation, cell proliferation, and cell detachment driven by chemically-dependent EMT in an epithelial tissue. A continuum description of the model is then derived in the form of a novel nonlinear free boundary problem. Using the discrete and continuum models we explore how the coupling of chemical transport and mechanical interactions influences EMT, and postulate how this could be used to help control EMT in pathological situations.
Collapse
Affiliation(s)
- Ryan J Murphy
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| | - Pascal R Buenzli
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| | - Tamara A Tambyah
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| | - Erik W Thompson
- Queensland University of Technology, Institute of Health and Biomedical Innovation, Brisbane, Australia.,Queensland University of Technology, School of Biomedical Sciences, Faculty of Health, Brisbane, Australia.,Translational Research Institute, Brisbane, Australia
| | - Honor J Hugo
- Queensland University of Technology, Institute of Health and Biomedical Innovation, Brisbane, Australia.,Queensland University of Technology, School of Biomedical Sciences, Faculty of Health, Brisbane, Australia.,Translational Research Institute, Brisbane, Australia
| | - Ruth E Baker
- University of Oxford, Mathematical Institute, Oxford, United Kingdom
| | - Matthew J Simpson
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| |
Collapse
|
9
|
Smith CA, Yates CA. Incorporating domain growth into hybrid methods for reaction-diffusion systems. J R Soc Interface 2021; 18:20201047. [PMID: 33849339 DOI: 10.1098/rsif.2020.1047] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Reaction-diffusion mechanisms are a robust paradigm that can be used to represent many biological and physical phenomena over multiple spatial scales. Applications include intracellular dynamics, the migration of cells and the patterns formed by vegetation in semi-arid landscapes. Moreover, domain growth is an important process for embryonic growth and wound healing. There are many numerical modelling frameworks capable of simulating such systems on growing domains; however, each of these may be well suited to different spatial scales and particle numbers. Recently, spatially extended hybrid methods on static domains have been produced to bridge the gap between these different modelling paradigms in order to represent multi-scale phenomena. However, such methods have not been developed with domain growth in mind. In this paper, we develop three hybrid methods on growing domains, extending three of the prominent static-domain hybrid methods. We also provide detailed algorithms to allow others to employ them. We demonstrate that the methods are able to accurately model three representative reaction-diffusion systems accurately and without bias.
Collapse
Affiliation(s)
- Cameron A Smith
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| | - Christian A Yates
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| |
Collapse
|
10
|
Karagöz Z, Rijns L, Dankers PY, van Griensven M, Carlier A. Towards understanding the messengers of extracellular space: Computational models of outside-in integrin reaction networks. Comput Struct Biotechnol J 2020; 19:303-314. [PMID: 33425258 PMCID: PMC7779863 DOI: 10.1016/j.csbj.2020.12.025] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 12/15/2020] [Accepted: 12/16/2020] [Indexed: 02/06/2023] Open
Abstract
The interactions between cells and their extracellular matrix (ECM) are critically important for homeostatic control of cell growth, proliferation, differentiation and apoptosis. Transmembrane integrin molecules facilitate the communication between ECM and the cell. Since the characterization of integrins in the late 1980s, there has been great advancement in understanding the function of integrins at different subcellular levels. However, the versatility in molecular pathways integrins are involved in, the high diversity in their interaction partners both outside and inside the cell as well as on the cell membrane and the short lifetime of events happening at the cell-ECM interface make it difficult to elucidate all the details regarding integrin function experimentally. To overcome the experimental challenges and advance the understanding of integrin biology, computational modeling tools have been used extensively. In this review, we summarize the computational models of integrin signaling while we explain the function of integrins at three main subcellular levels (outside the cell, cell membrane, cytosol). We also discuss how these computational modeling efforts can be helpful in other disciplines such as biomaterial design. As such, this review is a didactic modeling summary for biomaterial researchers interested in complementing their experimental work with computational tools or for seasoned computational scientists that would like to advance current in silico integrin models.
Collapse
Affiliation(s)
- Zeynep Karagöz
- Department of Cell Biology-Inspired Tissue Engineering, MERLN Institute for Technology-Inspired Regenerative Medicine, Maastricht University, Universiteitssingel 40, 6229 ER Maastricht, the Netherlands
| | - Laura Rijns
- Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven, the Netherlands
| | - Patricia Y.W. Dankers
- Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven, the Netherlands
| | - Martijn van Griensven
- Department of Cell Biology-Inspired Tissue Engineering, MERLN Institute for Technology-Inspired Regenerative Medicine, Maastricht University, Universiteitssingel 40, 6229 ER Maastricht, the Netherlands
| | - Aurélie Carlier
- Department of Cell Biology-Inspired Tissue Engineering, MERLN Institute for Technology-Inspired Regenerative Medicine, Maastricht University, Universiteitssingel 40, 6229 ER Maastricht, the Netherlands
| |
Collapse
|
11
|
Yates CA, George A, Jordana A, Smith CA, Duncan AB, Zygalakis KC. The blending region hybrid framework for the simulation of stochastic reaction-diffusion processes. J R Soc Interface 2020; 17:20200563. [PMID: 33081647 PMCID: PMC7653393 DOI: 10.1098/rsif.2020.0563] [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] [Indexed: 11/19/2022] Open
Abstract
The simulation of stochastic reaction–diffusion systems using fine-grained representations can become computationally prohibitive when particle numbers become large. If particle numbers are sufficiently high then it may be possible to ignore stochastic fluctuations and use a more efficient coarse-grained simulation approach. Nevertheless, for multiscale systems which exhibit significant spatial variation in concentration, a coarse-grained approach may not be appropriate throughout the simulation domain. Such scenarios suggest a hybrid paradigm in which a computationally cheap, coarse-grained model is coupled to a more expensive, but more detailed fine-grained model, enabling the accurate simulation of the fine-scale dynamics at a reasonable computational cost. In this paper, in order to couple two representations of reaction–diffusion at distinct spatial scales, we allow them to overlap in a ‘blending region’. Both modelling paradigms provide a valid representation of the particle density in this region. From one end of the blending region to the other, control of the implementation of diffusion is passed from one modelling paradigm to another through the use of complementary ‘blending functions’ which scale up or down the contribution of each model to the overall diffusion. We establish the reliability of our novel hybrid paradigm by demonstrating its simulation on four exemplar reaction–diffusion scenarios.
Collapse
Affiliation(s)
- Christian A Yates
- Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| | - Adam George
- Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| | - Armand Jordana
- Centre de Mathématiques et de Leurs Applications, CNRS, ENS Paris-Saclay, Université Paris-Saclay, 94235 Cachan cedex, France
| | - Cameron A Smith
- Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| | - Andrew B Duncan
- Department of Mathematics, Imperial College London, London SW7 2AZ, UK
| | - Konstantinos C Zygalakis
- School of Mathematics, University of Edinburgh, James Clerk Maxwell Building, The King's Buildings, Peter Guthrie Tait Road, Edinburgh EH9 3FD, UK
| |
Collapse
|
12
|
Abstract
Domain growth is a key process in many areas of biology, including embryonic development, the growth of tissue, and limb regeneration. As a result, mechanisms for incorporating it into traditional models for cell movement, interaction, and proliferation are of great importance. A previously well-used method to incorporate domain growth into on-lattice reaction-diffusion models causes a buildup of particles on the boundaries of the domain, which is particularly evident when diffusion is low in comparison to the rate of domain growth. Here we present an alternative method which addresses this unphysical buildup of particles at the boundaries and demonstrate that it is accurate for scenarios in which the previous method fails. Further, we discuss for which parameter regimes it is feasible to continue using the original method due to diffusion dominating the domain growth mechanism.
Collapse
Affiliation(s)
- Cameron A Smith
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, United Kingdom
| | - Cécile Mailler
- Probability Laboratory, Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, United Kingdom
| | - Christian A Yates
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, United Kingdom
| |
Collapse
|
13
|
Ruiz-Martínez Á, Bartol TM, Sejnowski TJ, Tartakovsky DM. Stochastic self-tuning hybrid algorithm for reaction-diffusion systems. J Chem Phys 2019; 151:244117. [PMID: 31893874 PMCID: PMC7341680 DOI: 10.1063/1.5125022] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2019] [Accepted: 12/01/2019] [Indexed: 02/06/2023] Open
Abstract
Many biochemical phenomena involve reactants with vastly different concentrations, some of which are amenable to continuum-level descriptions, while the others are not. We present a hybrid self-tuning algorithm to model such systems. The method combines microscopic (Brownian) dynamics for diffusion with mesoscopic (Gillespie-type) methods for reactions and remains efficient in a wide range of regimes and scenarios with large variations of concentrations. Its accuracy, robustness, and versatility are balanced by redefining propensities and optimizing the mesh size and time step. We use a bimolecular reaction to demonstrate the potential of our method in a broad spectrum of scenarios: from almost completely reaction-dominated systems to cases where reactions rarely occur or take place very slowly. The simulation results show that the number of particles present in the system does not degrade the performance of our method. This makes it an accurate and computationally efficient tool to model complex multireaction systems.
Collapse
Affiliation(s)
- Á Ruiz-Martínez
- Department of Mechanical and Aerospace Engineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
| | - T M Bartol
- Computational Neurobiology Laboratory, Salk Institute for Biological Studies, 10010 North Torrey Pines Road, La Jolla, California 92037, USA
| | - T J Sejnowski
- Computational Neurobiology Laboratory, Salk Institute for Biological Studies, 10010 North Torrey Pines Road, La Jolla, California 92037, USA
| | - D M Tartakovsky
- Department of Energy Resources Engineering, Stanford University, 367 Panama Street, Stanford, California 94305, USA
| |
Collapse
|
14
|
Surendran A, Plank MJ, Simpson MJ. Spatial structure arising from chase-escape interactions with crowding. Sci Rep 2019; 9:14988. [PMID: 31628421 PMCID: PMC6800429 DOI: 10.1038/s41598-019-51565-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 10/03/2019] [Indexed: 12/17/2022] Open
Abstract
Movement of individuals, mediated by localised interactions, plays a key role in numerous processes including cell biology and ecology. In this work, we investigate an individual-based model accounting for various intraspecies and interspecies interactions in a community consisting of two distinct species. In this framework we consider one species to be chasers and the other species to be escapees, and we focus on chase-escape dynamics where the chasers are biased to move towards the escapees, and the escapees are biased to move away from the chasers. This framework allows us to explore how individual-level directional interactions scale up to influence spatial structure at the macroscale. To focus exclusively on the role of motility and directional bias in determining spatial structure, we consider conservative communities where the number of individuals in each species remains constant. To provide additional information about the individual-based model, we also present a mathematically tractable deterministic approximation based on describing the evolution of the spatial moments. We explore how different features of interactions including interaction strength, spatial extent of interaction, and relative density of species influence the formation of the macroscale spatial patterns.
Collapse
Affiliation(s)
- Anudeep Surendran
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - Michael J Plank
- School of Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand.,Te Pūnaha Matatini, A New Zealand Centre of Research Excellence, Auckland, New Zealand
| | - Matthew J Simpson
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia.
| |
Collapse
|
15
|
Fu Y, Yogurtcu ON, Kothari R, Thorkelsdottir G, Sodt AJ, Johnson ME. An implicit lipid model for efficient reaction-diffusion simulations of protein binding to surfaces of arbitrary topology. J Chem Phys 2019; 151:124115. [PMID: 31575182 DOI: 10.1063/1.5120516] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Localization of proteins to a membrane is an essential step in a broad range of biological processes such as signaling, virion formation, and clathrin-mediated endocytosis. The strength and specificity of proteins binding to a membrane depend on the lipid composition. Single-particle reaction-diffusion methods offer a powerful tool for capturing lipid-specific binding to membrane surfaces by treating lipids explicitly as individual diffusible binding sites. However, modeling lipid particle populations is expensive. Here, we present an algorithm for reversible binding of proteins to continuum surfaces with implicit lipids, providing dramatic speed-ups to many body simulations. Our algorithm can be readily integrated into most reaction-diffusion software packages. We characterize changes to kinetics that emerge from explicit vs implicit lipids as well as surface adsorption models, showing excellent agreement between our method and the full explicit lipid model. Compared to models of surface adsorption, which couple together binding affinity and lipid concentration, our implicit lipid model decouples them to provide more flexibility for controlling surface binding properties and lipid inhomogeneity, thus reproducing binding kinetics and equilibria. Crucially, we demonstrate our method's application to membranes of arbitrary curvature and topology, modeled via a subdivision limit surface, again showing excellent agreement with explicit lipid simulations. Unlike adsorption models, our method retains the ability to bind lipids after proteins are localized to the surface (through, e.g., a protein-protein interaction), which can greatly increase the stability of multiprotein complexes on the surface. Our method will enable efficient cell-scale simulations involving proteins localizing to realistic membrane models, which is a critical step for predictive modeling and quantification of in vitro and in vivo dynamics.
Collapse
Affiliation(s)
- Yiben Fu
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, 3400 N. Charles St., Baltimore, Maryland 21218, USA
| | - Osman N Yogurtcu
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, 3400 N. Charles St., Baltimore, Maryland 21218, USA
| | - Ruchita Kothari
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, USA
| | - Gudrun Thorkelsdottir
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, USA
| | - Alexander J Sodt
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, USA
| | - Margaret E Johnson
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, 3400 N. Charles St., Baltimore, Maryland 21218, USA
| |
Collapse
|
16
|
Kang HW, Erban R. Multiscale Stochastic Reaction-Diffusion Algorithms Combining Markov Chain Models with Stochastic Partial Differential Equations. Bull Math Biol 2019; 81:3185-3213. [PMID: 31165406 PMCID: PMC6677718 DOI: 10.1007/s11538-019-00613-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2018] [Accepted: 05/09/2019] [Indexed: 12/21/2022]
Abstract
Two multiscale algorithms for stochastic simulations of reaction-diffusion processes are analysed. They are applicable to systems which include regions with significantly different concentrations of molecules. In both methods, a domain of interest is divided into two subsets where continuous-time Markov chain models and stochastic partial differential equations (SPDEs) are used, respectively. In the first algorithm, Markov chain (compartment-based) models are coupled with reaction-diffusion SPDEs by considering a pseudo-compartment (also called an overlap or handshaking region) in the SPDE part of the computational domain right next to the interface. In the second algorithm, no overlap region is used. Further extensions of both schemes are presented, including the case of an adaptively chosen boundary between different modelling approaches.
Collapse
Affiliation(s)
- Hye-Won Kang
- Department of Mathematics and Statistics, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250 USA
| | - Radek Erban
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG UK
| |
Collapse
|
17
|
Gunaratne RS, Wilson DB, Flegg MB, Erban R. Multi-resolution dimer models in heat baths with short-range and long-range interactions. Interface Focus 2019; 9:20180070. [PMID: 31065341 PMCID: PMC6501348 DOI: 10.1098/rsfs.2018.0070] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/11/2019] [Indexed: 11/16/2022] Open
Abstract
This work investigates multi-resolution methodologies for simulating dimer models. The solvent particles which make up the heat bath interact with the monomers of the dimer either through direct collisions (short-range) or through harmonic springs (long-range). Two types of multi-resolution methodologies are considered in detail: (a) describing parts of the solvent far away from the dimer by a coarser approach; (b) describing each monomer of the dimer by using a model with different level of resolution. These methodologies are then used to investigate the effect of a shared heat bath versus two uncoupled heat baths, one for each monomer. Furthermore, the validity of the multi-resolution methods is discussed by comparison to dynamics of macroscopic Langevin equations.
Collapse
Affiliation(s)
- Ravinda S. Gunaratne
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Daniel B. Wilson
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Mark B. Flegg
- School of Mathematical Sciences, Monash University, 9 Rainforest walk, Clayton campus, Victoria 3168, Australia
| | - Radek Erban
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| |
Collapse
|
18
|
A hybrid method for micro-mesoscopic stochastic simulation of reaction-diffusion systems. Math Biosci 2019; 312:23-32. [PMID: 30998936 DOI: 10.1016/j.mbs.2019.04.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2018] [Revised: 04/13/2019] [Accepted: 04/14/2019] [Indexed: 12/19/2022]
Abstract
The present paper introduces a new micro-meso hybrid algorithm based on the Ghost Cell Method concept in which the microscopic subdomain is governed by the Reactive Multi-Particle Collision (RMPC) dynamics. The mesoscopic subdomain is modeled using the Reaction-Diffusion Master Equation (RDME). The RDME is solved by means of the Inhomogeneous Stochastic Simulation Algorithm. No hybrid algorithm has hitherto used the RMPC dynamics for modeling reactions and the trajectories of each individual particle. The RMPC is faster than other molecular based methods and has the advantage of conserving mass, energy and momentum in the collision and free streaming steps. The new algorithm is tested on three reaction-diffusion systems. In all the systems studied, very good agreement with the deterministic solutions of the corresponding differential equations is obtained. In addition, it has been shown that proper discretization of the computational domain results in significant speed-ups in comparison with the full RMPC algorithm.
Collapse
|
19
|
Smith CA, Yates CA. The auxiliary region method: a hybrid method for coupling PDE- and Brownian-based dynamics for reaction-diffusion systems. ROYAL SOCIETY OPEN SCIENCE 2018; 5:180920. [PMID: 30225082 PMCID: PMC6124063 DOI: 10.1098/rsos.180920] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/18/2018] [Accepted: 06/28/2018] [Indexed: 06/08/2023]
Abstract
Reaction-diffusion systems are used to represent many biological and physical phenomena. They model the random motion of particles (diffusion) and interactions between them (reactions). Such systems can be modelled at multiple scales with varying degrees of accuracy and computational efficiency. When representing genuinely multiscale phenomena, fine-scale models can be prohibitively expensive, whereas coarser models, although cheaper, often lack sufficient detail to accurately represent the phenomenon at hand. Spatial hybrid methods couple two or more of these representations in order to improve efficiency without compromising accuracy. In this paper, we present a novel spatial hybrid method, which we call the auxiliary region method (ARM), which couples PDE- and Brownian-based representations of reaction-diffusion systems. Numerical PDE solutions on one side of an interface are coupled to Brownian-based dynamics on the other side using compartment-based 'auxiliary regions'. We demonstrate that the hybrid method is able to simulate reaction-diffusion dynamics for a number of different test problems with high accuracy. Furthermore, we undertake error analysis on the ARM which demonstrates that it is robust to changes in the free parameters in the model, where previous coupling algorithms are not. In particular, we envisage that the method will be applicable for a wide range of spatial multi-scales problems including filopodial dynamics, intracellular signalling, embryogenesis and travelling wave phenomena.
Collapse
Affiliation(s)
- Cameron A. Smith
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| | | |
Collapse
|