7,415 research outputs found
Data-driven modelling of biological multi-scale processes
Biological processes involve a variety of spatial and temporal scales. A
holistic understanding of many biological processes therefore requires
multi-scale models which capture the relevant properties on all these scales.
In this manuscript we review mathematical modelling approaches used to describe
the individual spatial scales and how they are integrated into holistic models.
We discuss the relation between spatial and temporal scales and the implication
of that on multi-scale modelling. Based upon this overview over
state-of-the-art modelling approaches, we formulate key challenges in
mathematical and computational modelling of biological multi-scale and
multi-physics processes. In particular, we considered the availability of
analysis tools for multi-scale models and model-based multi-scale data
integration. We provide a compact review of methods for model-based data
integration and model-based hypothesis testing. Furthermore, novel approaches
and recent trends are discussed, including computation time reduction using
reduced order and surrogate models, which contribute to the solution of
inference problems. We conclude the manuscript by providing a few ideas for the
development of tailored multi-scale inference methods.Comment: This manuscript will appear in the Journal of Coupled Systems and
Multiscale Dynamics (American Scientific Publishers
MOLNs: A cloud platform for interactive, reproducible and scalable spatial stochastic computational experiments in systems biology using PyURDME
Computational experiments using spatial stochastic simulations have led to
important new biological insights, but they require specialized tools, a
complex software stack, as well as large and scalable compute and data analysis
resources due to the large computational cost associated with Monte Carlo
computational workflows. The complexity of setting up and managing a
large-scale distributed computation environment to support productive and
reproducible modeling can be prohibitive for practitioners in systems biology.
This results in a barrier to the adoption of spatial stochastic simulation
tools, effectively limiting the type of biological questions addressed by
quantitative modeling. In this paper, we present PyURDME, a new, user-friendly
spatial modeling and simulation package, and MOLNs, a cloud computing appliance
for distributed simulation of stochastic reaction-diffusion models. MOLNs is
based on IPython and provides an interactive programming platform for
development of sharable and reproducible distributed parallel computational
experiments
Emergent Properties of Tumor Microenvironment in a Real-life Model of Multicell Tumor Spheroids
Multicellular tumor spheroids are an important {\it in vitro} model of the
pre-vascular phase of solid tumors, for sizes well below the diagnostic limit:
therefore a biophysical model of spheroids has the ability to shed light on the
internal workings and organization of tumors at a critical phase of their
development. To this end, we have developed a computer program that integrates
the behavior of individual cells and their interactions with other cells and
the surrounding environment. It is based on a quantitative description of
metabolism, growth, proliferation and death of single tumor cells, and on
equations that model biochemical and mechanical cell-cell and cell-environment
interactions. The program reproduces existing experimental data on spheroids,
and yields unique views of their microenvironment. Simulations show complex
internal flows and motions of nutrients, metabolites and cells, that are
otherwise unobservable with current experimental techniques, and give novel
clues on tumor development and strong hints for future therapies.Comment: 20 pages, 10 figures. Accepted for publication in PLOS One. The
published version contains links to a supplementary text and three video
file
Efficient Reactive Brownian Dynamics
We develop a Split Reactive Brownian Dynamics (SRBD) algorithm for particle
simulations of reaction-diffusion systems based on the Doi or volume reactivity
model, in which pairs of particles react with a specified Poisson rate if they
are closer than a chosen reactive distance. In our Doi model, we ensure that
the microscopic reaction rules for various association and disassociation
reactions are consistent with detailed balance (time reversibility) at
thermodynamic equilibrium. The SRBD algorithm uses Strang splitting in time to
separate reaction and diffusion, and solves both the diffusion-only and
reaction-only subproblems exactly, even at high packing densities. To
efficiently process reactions without uncontrolled approximations, SRBD employs
an event-driven algorithm that processes reactions in a time-ordered sequence
over the duration of the time step. A grid of cells with size larger than all
of the reactive distances is used to schedule and process the reactions, but
unlike traditional grid-based methods such as Reaction-Diffusion Master
Equation (RDME) algorithms, the results of SRBD are statistically independent
of the size of the grid used to accelerate the processing of reactions. We use
the SRBD algorithm to compute the effective macroscopic reaction rate for both
reaction- and diffusion-limited irreversible association in three dimensions.
We also study long-time tails in the time correlation functions for reversible
association at thermodynamic equilibrium. Finally, we compare different
particle and continuum methods on a model exhibiting a Turing-like instability
and pattern formation. We find that for models in which particles diffuse off
lattice, such as the Doi model, reactions lead to a spurious enhancement of the
effective diffusion coefficients.Comment: To appear in J. Chem. Phy
A Software for Particle-Based Reaction-Diffusion Dynamics in Crowded Cellular Environments
We introduce the software package ReaDDy for simulation of detailed
spatiotemporal mechanisms of dynamical processes in the cell, based on
reaction-diffusion dynamics with particle resolution. In contrast to other
particle-based reaction kinetics programs, ReaDDy supports particle
interaction potentials. This permits effects such as space exclusion,
molecular crowding and aggregation to be modeled. The biomolecules simulated
can be represented as a sphere, or as a more complex geometry such as a domain
structure or polymer chain. ReaDDy bridges the gap between small-scale but
highly detailed molecular dynamics or Brownian dynamics simulations and large-
scale but little-detailed reaction kinetics simulations. ReaDDy has a modular
design that enables the exchange of the computing core by efficient platform-
specific implementations or dynamical models that are different from Brownian
dynamics
- …