Spatial and stochastic cellular modeling with the Smoldyn simulator

Particle-based MR modeling with diffusion, microstructure, and enzymatic reactions

PURPOSE

To develop a novel particle-based in silico MR model and demonstrate applications of this model to signal mechanisms which are affected by the spatial organization of particles, including metabolic reaction kinetics, microstructural effects on diffusion, and radiofrequency (RF) refocusing effects in gradient-echo sequences.

METHODS

The model was developed by integrating a forward solution of the Bloch equations with a Brownian dynamics simulator. Simulation configurations were then designed to model MR signal dynamics of interest, with a primary focus on hyperpolarized 13C MRI methods. Phantom scans and spectrophotometric assays were conducted to validate model results in vitro.

RESULTS

The model accurately reproduced the reaction kinetics of enzyme-mediated conversion of pyruvate to lactate. When varying proportions of restrictive structure were added to the reaction volume, nonlinear changes in the reaction rate measured in vitro were replicated in silico. Modeling of RF refocusing effects characterized the degree of diffusion-weighted contribution from preserved residual magnetization in nonspoiled gradient-echo sequences.

CONCLUSIONS

These results show accurate reproduction of a range of MR signal mechanisms, establishing the model's capability to investigate the multifactorial signal dynamics such as those underlying hyperpolarized 13C MRI data.

Theoretical and Practical Aspects of Multienzyme Organization and Encapsulation

The advent of biotechnology has enabled metabolic engineers to assemble heterologous pathways in cells to produce a variety of products of industrial relevance, often in a sustainable way. However, many pathways face challenges of low product yield. These pathways often suffer from issues that are difficult to optimize, such as low pathway flux and off-target pathway consumption of intermediates. These issues are exacerbated by the need to balance pathway flux with the health of the cell, particularly when a toxic intermediate builds up. Nature faces similar challenges and has evolved spatial organization strategies to increase metabolic pathway flux and efficiency. Inspired by these strategies, bioengineers have developed clever strategies to mimic spatial organization in nature. This review explores the use of spatial organization strategies, including protein scaffolding and protein encapsulation inside of proteinaceous shells, toward overcoming bottlenecks in metabolic engineering efforts.

Active Transport by Cytoplasmic Dynein Maintains the Localization of MAP2 in Developing Neurons

MAP2 has been widely used as a marker of neuronal dendrites because of its extensive restriction in the somatodendritic region of neurons. Despite that, how the precise localization of such a soluble protein is established and maintained against thermal forces and diffusion has been elusive and long remained a mystery in neuroscience. In this study, we aimed to uncover the mechanism behind how MAP2 is retained in the somatodendritic region. Using GFP-tagged MAP2 expressed in cultured hippocampal neurons, we discovered a crucial protein region responsible for the localization of MAP2, the serine/proline-rich (S/P) region. Our pulse-chase live-cell imaging revealed the slow but steady migration of MAP2 toward distal dendrites, which was not observed in a MAP2 mutant lacking the S/P region, indicating that S/P-dependent transport is vital for the proper localization of MAP2. Furthermore, our experiments using an inhibitor of cytoplasmic Dynein, ciliobrevin D, as well as Dynein knockdown, showed that cytoplasmic Dynein is involved in the transport of MAP2 in dendrites. We also found that Dynein complex binds to MAP2 through the S/P region in heterologous cells. Using mathematical modeling based on experimental data, we confirmed that an intermittent active transport mechanism is essential. Thus, we propose that the cytoplasmic Dynein recruits and transports free MAP2 toward distal dendrites, thereby maintaining the precise dendritic localization of MAP2 in neurons. Our findings shed light on the previously unknown mechanism behind MAP2 localization and provide a new direction for soluble protein trafficking research in the field of cell biology of neurons.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser).

SIGNIFICANCE

There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

The Role of in silico Research in Developing Nanoparticle-Based Therapeutics

Nanoparticles (NPs) hold great potential as therapeutics, particularly in the realm of drug delivery. They are effective at functional cargo delivery and offer a great degree of amenability that can be used to offset toxic side effects or to target drugs to specific regions in the body. However, there are many challenges associated with the development of NP-based drug formulations that hamper their successful clinical translation. Arguably, the most significant barrier in the way of efficacious NP-based drug delivery systems is the tedious and time-consuming nature of NP formulation—a process that needs to account for downstream effects, such as the onset of potential toxicity or immunogenicity, in vivo biodistribution and overall pharmacokinetic profiles, all while maintaining desirable therapeutic outcomes. Computational and AI-based approaches have shown promise in alleviating some of these restrictions. Via predictive modeling and deep learning, in silico approaches have shown the ability to accurately model NP-membrane interactions and cellular uptake based on minimal data, such as the physicochemical characteristics of a given NP. More importantly, machine learning allows computational models to predict how specific changes could be made to the physicochemical characteristics of a NP to improve functional aspects, such as drug retention or endocytosis. On a larger scale, they are also able to predict the in vivo pharmacokinetics of NP-encapsulated drugs, predicting aspects such as circulatory half-life, toxicity, and biodistribution. However, the convergence of nanomedicine and computational approaches is still in its infancy and limited in its applicability. The interactions between NPs, the encapsulated drug and the body form an intricate network of interactions that cannot be modeled with absolute certainty. Despite this, rapid advancements in the area promise to deliver increasingly powerful tools capable of accelerating the development of advanced nanoscale therapeutics. Here, we describe computational approaches that have been utilized in the field of nanomedicine, focusing on approaches for NP design and engineering.

Python interfaces for the Smoldyn simulator

MOTIVATION

Smoldyn is a particle-based biochemical simulator that is frequently used for systems biology and biophysics research. Previously, users could only define models using text-based input or a C/C++ application programming interface (API), which were convenient, but limited extensibility.

RESULTS

We added a Python API to Smoldyn to improve integration with other software tools, such as Jupyter notebooks, other Python code libraries and other simulators. It includes low-level functions that closely mimic the existing C/C++ API and higher-level functions that are more convenient to use. These latter functions follow modern object-oriented Python conventions.

AVAILABILITY AND IMPLEMENTATION

Smoldyn is open source and free, available at http://www.smoldyn.org and can be installed with the Python package manager pip. It runs on Mac, Windows and Linux.

Documentation is available at http://www.smoldyn.org/SmoldynManual.pdf and https://smoldyn.readthedocs.io/en/latest/python/api.html.

Tools for computational analysis of moving boundary problems in cellular mechanobiology

A cell's ability to change shape is one of the most fundamental biological processes and is essential for maintaining healthy organisms. When the ability to control shape goes awry, it often results in a diseased system. As such, it is important to understand the mechanisms that allow a cell to sense and respond to its environment so as to maintain cellular shape homeostasis. Because of the inherent complexity of the system, computational models that are based on sound theoretical understanding of the biochemistry and biomechanics and that use experimentally measured parameters are an essential tool. These models involve an inherent feedback, whereby shape is determined by the action of regulatory signals whose spatial distribution depends on the shape. To carry out computational simulations of these moving boundary problems requires special computational techniques. A variety of alternative approaches, depending on the type and scale of question being asked, have been used to simulate various biological processes, including cell motility, division, mechanosensation, and cell engulfment. In general, these models consider the forces that act on the system (both internally generated, or externally imposed) and the mechanical properties of the cell that resist these forces. Moving forward, making these techniques more accessible to the non-expert will help improve interdisciplinary research thereby providing new insight into important biological processes that affect human health. This article is categorized under: Cancer > Cancer>Computational Models Cancer > Cancer>Molecular and Cellular Physiology.

Rule-Based Modeling Using Wildcards in the Smoldyn Simulator

Many biological molecules exist in multiple variants, such as proteins with different posttranslational modifications, DNAs with different sequences, and phospholipids with different chain lengths. Representing these variants as distinct species, as most biochemical simulators do, leads to the problem that the number of species, and chemical reactions that interconvert them, typically increase combinatorially with the number of ways that the molecules can vary. This can be alleviated by "rule-based modeling methods," in which software generates the chemical reaction network from relatively simple "rules." This chapter presents a new approach to rule-based modeling. It is based on wildcards that match to species names, much as wildcards can match to file names in computer operating systems. It is much simpler to use than the formal rule-based modeling approaches developed previously but can lead to unintended consequences if not used carefully. This chapter demonstrates rule-based modeling with wildcards through examples for signaling systems, protein complexation, polymerization, nucleic acid sequence copying and mutation, the "SMILES" chemical notation, and others. The method is implemented in Smoldyn, a spatial and stochastic biochemical simulator, for both generate-first and on-the-fly expansion, meaning whether the reaction network is generated before or during the simulation.

Using NEURON for Reaction-Diffusion Modeling of Extracellular Dynamics

Development of credible clinically-relevant brain simulations has been slowed due to a focus on electrophysiology in computational neuroscience, neglecting the multiscale whole-tissue modeling approach used for simulation in most other organ systems. We have now begun to extend the NEURON simulation platform in this direction by adding extracellular modeling. The extracellular medium of neural tissue is an active medium of neuromodulators, ions, inflammatory cells, oxygen, NO and other gases, with additional physiological, pharmacological and pathological agents. These extracellular agents influence, and are influenced by, cellular electrophysiology, and cellular chemophysiology-the complex internal cellular milieu of second-messenger signaling and cascades. NEURON's extracellular reaction-diffusion is supported by an intuitive Python-based where/who/what command sequence, derived from that used for intracellular reaction diffusion, to support coarse-grained macroscopic extracellular models. This simulation specification separates the expression of the conceptual model and parameters from the underlying numerical methods. In the volume-averaging approach used, the macroscopic model of tissue is characterized by free volume fraction-the proportion of space in which species are able to diffuse, and tortuosity-the average increase in path length due to obstacles. These tissue characteristics can be defined within particular spatial regions, enabling the modeler to account for regional differences, due either to intrinsic organization, particularly gray vs. white matter, or to pathology such as edema. We illustrate simulation development using spreading depression, a pathological phenomenon thought to play roles in migraine, epilepsy and stroke. Simulation results were verified against analytic results and against the extracellular portion of the simulation run under FiPy. The creation of this NEURON interface provides a pathway for interoperability that can be used to automatically export this class of models into complex intracellular/extracellular simulations and future cross-simulator standardization.

Graphics processing units in bioinformatics, computational biology and systems biology

Several studies in Bioinformatics, Computational Biology and Systems Biology rely on the definition of physico-chemical or mathematical models of biological systems at different scales and levels of complexity, ranging from the interaction of atoms in single molecules up to genome-wide interaction networks. Traditional computational methods and software tools developed in these research fields share a common trait: they can be computationally demanding on Central Processing Units (CPUs), therefore limiting their applicability in many circumstances. To overcome this issue, general-purpose Graphics Processing Units (GPUs) are gaining an increasing attention by the scientific community, as they can considerably reduce the running time required by standard CPU-based software, and allow more intensive investigations of biological systems. In this review, we present a collection of GPU tools recently developed to perform computational analyses in life science disciplines, emphasizing the advantages and the drawbacks in the use of these parallel architectures. The complete list of GPU-powered tools here reviewed is available at http://bit.ly/gputools.

Smoldyn: particle-based simulation with rule-based modeling, improved molecular interaction and a library interface

Motivation

Smoldyn is a spatial and stochastic biochemical simulator. It treats each molecule of interest as an individual particle in continuous space, simulating molecular diffusion, molecule-membrane interactions and chemical reactions, all with good accuracy. This article presents several new features.

Results

Smoldyn now supports two types of rule-based modeling. These are a wildcard method, which is very convenient, and the BioNetGen package with extensions for spatial simulation, which is better for complicated models. Smoldyn also includes new algorithms for simulating the diffusion of surface-bound molecules and molecules with excluded volume. Both are exact in the limit of short time steps and reasonably good with longer steps. In addition, Smoldyn supports single-molecule tracking simulations. Finally, the Smoldyn source code can be accessed through a C/C ++ language library interface.

Availability and Implementation

Smoldyn software, documentation, code, and examples are at http://www.smoldyn.org .

Contact

steven.s.andrews@gmail.com.

Dimensionality of Motion and Binding Valency Govern Receptor-Ligand Kinetics As Revealed by Agent-Based Modeling

Mathematical modeling and computer simulations have become an integral part of modern biological research. The strength of theoretical approaches is in the simplification of complex biological systems. We here consider the general problem of receptor–ligand binding in the context of antibody–antigen binding. On the one hand, we establish a quantitative mapping between macroscopic binding rates of a deterministic differential equation model and their microscopic equivalents as obtained from simulating the spatiotemporal binding kinetics by stochastic agent-based models. On the other hand, we investigate the impact of various properties of B cell-derived receptors—such as their dimensionality of motion, morphology, and binding valency—on the receptor–ligand binding kinetics. To this end, we implemented an algorithm that simulates antigen binding by B cell-derived receptors with a Y-shaped morphology that can move in different dimensionalities, i.e., either as membrane-anchored receptors or as soluble receptors. The mapping of the macroscopic and microscopic binding rates allowed us to quantitatively compare different agent-based model variants for the different types of B cell-derived receptors. Our results indicate that the dimensionality of motion governs the binding kinetics and that this predominant impact is quantitatively compensated by the bivalency of these receptors.

ML-Space: Hybrid Spatial Gillespie and Particle Simulation of Multi-Level Rule-Based Models in Cell Biology

Spatio-temporal dynamics of cellular processes can be simulated at different levels of detail, from (deterministic) partial differential equations via the spatial Stochastic Simulation algorithm to tracking Brownian trajectories of individual particles. We present a spatial simulation approach for multi-level rule-based models, which includes dynamically hierarchically nested cellular compartments and entities. Our approach ML-Space combines discrete compartmental dynamics, stochastic spatial approaches in discrete space, and particles moving in continuous space. The rule-based specification language of ML-Space supports concise and compact descriptions of models and to adapt the spatial resolution of models easily.

Monte Carlo Analysis of Molecule Absorption Probabilities in Diffusion-Based Nanoscale Communication Systems with Multiple Receivers

For biomedical applications of nanonetworks, employing molecular communication for information transport is advantageous over nano-electromagnetic communication: molecular communication is potentially biocompatible and inherently energy-efficient. Recently, several studies have modeled receivers in diffusion-based molecular communication systems as "perfectly monitoring" or "perfectly absorbing" spheres based on idealized descriptions of chemoreception. In this paper, we focus on perfectly absorbing receivers and present methods to improve the accuracy of simulation procedures that are used to analyze these receivers. We employ schemes available from the chemical physics and biophysics literature and outline a Monte Carlo simulation algorithm that accounts for the possibility of molecule absorption during discrete time steps, leading to a more accurate analysis of absorption probabilities. Unlike most existing studies that consider a single receiver, this paper analyzes absorption probabilities for multiple receivers deterministically or randomly deployed in a region. For random deployments, the ultimate absorption probabilities as a function of transmitter-receiver distance are shown to fit well to power laws; the exponents derived become more negative as the number of receivers increases up to a limit beyond which no additional receivers can be "packed" in the deployment region. This paper is expected to impact the design of molecular nanonetworks with multiple absorbing receivers.

Analysis of Ligand-Receptor Association and Intermediate Transfer Rates in Multienzyme Nanostructures with All-Atom Brownian Dynamics Simulations

We present the second-generation GeomBD Brownian dynamics software for determining interenzyme intermediate transfer rates and substrate association rates in biomolecular complexes. Substrate and intermediate association rates for a series of enzymes or biomolecules can be compared between the freely diffusing disorganized configuration and various colocalized or complexed arrangements for kinetic investigation of enhanced intermediate transfer. In addition, enzyme engineering techniques, such as synthetic protein conjugation, can be computationally modeled and analyzed to better understand changes in substrate association relative to native enzymes. Tools are provided to determine nonspecific ligand-receptor association residence times, and to visualize common sites of nonspecific association of substrates on receptor surfaces. To demonstrate features of the software, interenzyme intermediate substrate transfer rate constants are calculated and compared for all-atom models of DNA origami scaffold-bound bienzyme systems of glucose oxidase and horseradish peroxidase. Also, a DNA conjugated horseradish peroxidase enzyme was analyzed for its propensity to increase substrate association rates and substrate local residence times relative to the unmodified enzyme. We also demonstrate the rapid determination and visualization of common sites of nonspecific ligand-receptor association by using HIV-1 protease and an inhibitor, XK263. GeomBD2 accelerates simulations by precomputing van der Waals potential energy grids and electrostatic potential grid maps, and has a flexible and extensible support for all-atom and coarse-grained force fields. Simulation software is written in C++ and utilizes modern parallelization techniques for potential grid preparation and Brownian dynamics simulation processes. Analysis scripts, written in the Python scripting language, are provided for quantitative simulation analysis. GeomBD2 is applicable to the fields of biophysics, bioengineering, and enzymology in both predictive and explanatory roles.

Simulated single molecule microscopy with SMeagol

Summary: SMeagol is a software tool to simulate highly realistic microscopy data based on spatial systems biology models, in order to facilitate development, validation and optimization of advanced analysis methods for live cell single molecule microscopy data.

Availability and implementation: SMeagol runs on Matlab R2014 and later, and uses compiled binaries in C for reaction–diffusion simulations. Documentation, source code and binaries for Mac OS, Windows and Ubuntu Linux can be downloaded from http://smeagol.sourceforge.net.

Contact:johan.elf@icm.uu.se

Supplementary information: Supplementary data are available at Bioinformatics online.

Modeling Excitable Dynamics of Chemotactic Networks

The study of chemotaxis has benefited greatly from computational models that describe the response of cells to chemoattractant stimuli. These models must keep track of spatially and temporally varying distributions of numerous intracellular species. Moreover, recent evidence suggests that these are not deterministic interactions, but also include the effect of stochastic variations that trigger an excitable network. In this chapter we illustrate how to create simulations of excitable networks using the Virtual Cell modeling environment.

Systems Biology Approaches for Understanding Genome Architecture

The linear and three-dimensional arrangement and composition of chromatin in eukaryotic genomes underlies the mechanisms directing gene regulation. Understanding this organization requires the integration of many data types and experimental results. Here we describe the approach of integrating genome-wide protein-DNA binding data to determine chromatin states. To investigate spatial aspects of genome organization, we present a detailed description of how to run stochastic simulations of protein movements within a simulated nucleus in 3D. This systems level approach enables the development of novel questions aimed at understanding the basic mechanisms that regulate genome dynamics.

Discrete-continuous reaction-diffusion model with mobile point-like sources and sinks

In many applications in soft and biological physics, there are multiple time and length scales involved but often with a distinct separation between them. For instance, in enzyme kinetics, enzymes are relatively large, move slowly and their copy numbers are typically small, while the metabolites (being transformed by these enzymes) are often present in abundance, are small in size and diffuse fast. It seems thus natural to apply different techniques to different time and length levels and couple them. Here we explore this possibility by constructing a stochastic-deterministic discrete-continuous reaction-diffusion model with mobile sources and sinks. Such an approach allows in particular to separate different sources of stochasticity. We demonstrate its application by modelling enzyme-catalysed reactions with freely diffusing enzymes and a heterogeneous source of metabolites. Our calculations suggest that using a higher amount of less active enzymes, as compared to fewer more active enzymes, reduces the metabolite pool size and correspondingly the lag time, giving rise to a faster response to external stimuli. The methodology presented can be extended to more complex systems and offers exciting possibilities for studying problems where spatial heterogeneities, stochasticity or discreteness play a role.

Modeling for (physical) biologists: an introduction to the rule-based approach

Models that capture the chemical kinetics of cellular regulatory networks can be specified in terms of rules for biomolecular interactions. A rule defines a generalized reaction, meaning a reaction that permits multiple reactants, each capable of participating in a characteristic transformation and each possessing certain, specified properties, which may be local, such as the state of a particular site or domain of a protein. In other words, a rule defines a transformation and the properties that reactants must possess to participate in the transformation. A rule also provides a rate law. A rule-based approach to modeling enables consideration of mechanistic details at the level of functional sites of biomolecules and provides a facile and visual means for constructing computational models, which can be analyzed to study how system-level behaviors emerge from component interactions.

Modeling biomolecular site dynamics in immunoreceptor signaling systems

The immune system plays a central role in human health. The activities of immune cells, whether defending an organism from disease or triggering a pathological condition such as autoimmunity, are driven by the molecular machinery of cellular signaling systems. Decades of experimentation have elucidated many of the biomolecules and interactions involved in immune signaling and regulation, and recently developed technologies have led to new types of quantitative, systems-level data. To integrate such information and develop nontrivial insights into the immune system, computational modeling is needed, and it is essential for modeling methods to keep pace with experimental advances. In this chapter, we focus on the dynamic, site-specific, and context-dependent nature of interactions in immunoreceptor signaling (i.e., the biomolecular site dynamics of immunoreceptor signaling), the challenges associated with capturing these details in computational models, and how these challenges have been met through use of rule-based modeling approaches.

Multiscale reaction-diffusion simulations with Smoldyn

Summary: Smoldyn is a software package for stochastic modelling of spatial biochemical networks and intracellular systems. It was originally developed with an accurate off-lattice particle-based model at its core. This has recently been enhanced with the addition of a computationally efficient on-lattice model, which can be run stand-alone or coupled together for multiscale simulations using both models in regions where they are most required, increasing the applicability of Smoldyn to larger molecule numbers and spatial domains. Simulations can switch between models with only small additions to their configuration file, enabling users with existing Smoldyn configuration files to run the new on-lattice model with any reaction, species or surface descriptions they might already have.

Availability and Implementation: Source code and binaries freely available for download at www.smoldyn.org, implemented in C/C++ and supported on Linux, Mac OSX and MS Windows.

Contact:martin.robinson@maths.ox.ac.uk

Supplementary Information: Supplementary data are available at Bioinformatics online and include additional details on model specification and modelling of surfaces, as well as the Smoldyn configuration file used to generate Figure 1.

Adaptive two-regime method: application to front propagation

The Adaptive Two-Regime Method (ATRM) is developed for hybrid (multiscale) stochastic simulation of reaction-diffusion problems. It efficiently couples detailed Brownian dynamics simulations with coarser lattice-based models. The ATRM is a generalization of the previously developed Two-Regime Method [Flegg et al., J. R. Soc., Interface 9, 859 (2012)] to multiscale problems which require a dynamic selection of regions where detailed Brownian dynamics simulation is used. Typical applications include a front propagation or spatio-temporal oscillations. In this paper, the ATRM is used for an in-depth study of front propagation in a stochastic reaction-diffusion system which has its mean-field model given in terms of the Fisher equation [R. Fisher, Ann. Eugen. 7, 355 (1937)]. It exhibits a travelling reaction front which is sensitive to stochastic fluctuations at the leading edge of the wavefront. Previous studies into stochastic effects on the Fisher wave propagation speed have focused on lattice-based models, but there has been limited progress using off-lattice (Brownian dynamics) models, which suffer due to their high computational cost, particularly at the high molecular numbers that are necessary to approach the Fisher mean-field model. By modelling only the wavefront itself with the off-lattice model, it is shown that the ATRM leads to the same Fisher wave results as purely off-lattice models, but at a fraction of the computational cost. The error analysis of the ATRM is also presented for a morphogen gradient model.

Biocellion: accelerating computer simulation of multicellular biological system models

MOTIVATION

Biological system behaviors are often the outcome of complex interactions among a large number of cells and their biotic and abiotic environment. Computational biologists attempt to understand, predict and manipulate biological system behavior through mathematical modeling and computer simulation. Discrete agent-based modeling (in combination with high-resolution grids to model the extracellular environment) is a popular approach for building biological system models. However, the computational complexity of this approach forces computational biologists to resort to coarser resolution approaches to simulate large biological systems. High-performance parallel computers have the potential to address the computing challenge, but writing efficient software for parallel computers is difficult and time-consuming.

RESULTS

We have developed Biocellion, a high-performance software framework, to solve this computing challenge using parallel computers. To support a wide range of multicellular biological system models, Biocellion asks users to provide their model specifics by filling the function body of pre-defined model routines. Using Biocellion, modelers without parallel computing expertise can efficiently exploit parallel computers with less effort than writing sequential programs from scratch. We simulate cell sorting, microbial patterning and a bacterial system in soil aggregate as case studies.

AVAILABILITY AND IMPLEMENTATION

Biocellion runs on x86 compatible systems with the 64 bit Linux operating system and is freely available for academic use. Visit http://biocellion.com for additional information.

The long and viscous road: uncovering nuclear diffusion barriers in closed mitosis

Diffusion barriers are effective means for constraining protein lateral exchange in cellular membranes. In Saccharomyces cerevisiae, they have been shown to sustain parental identity through asymmetric segregation of ageing factors during closed mitosis. Even though barriers have been extensively studied in the plasma membrane, their identity and organization within the nucleus remains poorly understood. Based on different lines of experimental evidence, we present a model of the composition and structural organization of a nuclear diffusion barrier during anaphase. By means of spatial stochastic simulations, we propose how specialised lipid domains, protein rings, and morphological changes of the nucleus may coordinate to restrict protein exchange between mother and daughter nuclear lobes. We explore distinct, plausible configurations of these diffusion barriers and offer testable predictions regarding their protein exclusion properties and the diffusion regimes they generate. Our model predicts that, while a specialised lipid domain and an immobile protein ring at the bud neck can compartmentalize the nucleus during early anaphase; a specialised lipid domain spanning the elongated bridge between lobes would be entirely sufficient during late anaphase. Our work shows how complex nuclear diffusion barriers in closed mitosis may arise from simple nanoscale biophysical interactions.

Spatial segregation of molecular contents is often necessary for an accurate, timely accomplishment of cellular functions, such as signal transduction and cell-fate decisions. For instance, budding yeast division requires the asymmetric segregation of proteins to distinguish a newborn cell from its parent. However, the strategies to achieve this parental identity are poorly understood. This holds especially true for key proteins and molecular complexes involved in mitosis that diffuse within the nuclear envelope. In fact, segregation within the nuclear envelope has been experimentally verified, but both the nature and configuration of any plausible diffusion barrier remain unknown. In this work, we built virtual models of the nucleus and carried out simulations testing the plausibility of specialised lipid domains and protein rings constituting the diffusion barrier. Moreover, we explored distinct barrier configurations in early and late stages of cell division, and verified our simulation results match experimental observations. Our work shows that the biophysical properties of these molecules, coordinated with morphological changes in the nucleus, make them suitable components of the nuclear diffusion barrier. Importantly, our research approach offers a novel avenue to study diffusion barriers in other biological membranes.

Parallel solutions for voxel-based simulations of reaction-diffusion systems

There is an increasing awareness of the pivotal role of noise in biochemical processes and of the effect of molecular crowding on the dynamics of biochemical systems. This necessity has given rise to a strong need for suitable and sophisticated algorithms for the simulation of biological phenomena taking into account both spatial effects and noise. However, the high computational effort characterizing simulation approaches, coupled with the necessity to simulate the models several times to achieve statistically relevant information on the model behaviours, makes such kind of algorithms very time-consuming for studying real systems. So far, different parallelization approaches have been deployed to reduce the computational time required to simulate the temporal dynamics of biochemical systems using stochastic algorithms. In this work we discuss these aspects for the spatial TAU-leaping in crowded compartments (STAUCC) simulator, a voxel-based method for the stochastic simulation of reaction-diffusion processes which relies on the Sτ-DPP algorithm. In particular we present how the characteristics of the algorithm can be exploited for an effective parallelization on the present heterogeneous HPC architectures.

Python-based geometry preparation and simulation visualization toolkits for STEPS

STEPS is a stochastic reaction-diffusion simulation engine that implements a spatial extension of Gillespie's Stochastic Simulation Algorithm (SSA) in complex tetrahedral geometries. An extensive Python-based interface is provided to STEPS so that it can interact with the large number of scientific packages in Python. However, a gap existed between the interfaces of these packages and the STEPS user interface, where supporting toolkits could reduce the amount of scripting required for research projects. This paper introduces two new supporting toolkits that support geometry preparation and visualization for STEPS simulations.

Cellular and molecular structure as a unifying framework for whole-cell modeling

Whole-cell modeling has the potential to play a major role in revolutionizing our understanding of cellular biology over the next few decades. A computational model of the entire cell would allow cellular biologists to integrate data from many disparate sources in a single consistent framework. Such a comprehensive model would be useful both for hypothesis testing and in the discovery of new behaviors that emerge from complex biological networks. Cellular and molecular structure can and should be a key organizing principle in a whole-cell model, connecting models across time and length scales in a multiscale approach. Here I present a summary of recent research centered around using molecular and cellular structure to model the behavior of cells.

Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems

Rule-based modeling was developed to address the limitations of traditional approaches for modeling chemical kinetics in cell signaling systems. These systems consist of multiple interacting biomolecules (e.g., proteins), which themselves consist of multiple parts (e.g., domains, linear motifs, and sites of phosphorylation). Consequently, biomolecules that mediate information processing generally have the potential to interact in multiple ways, with the number of possible complexes and posttranslational modification states tending to grow exponentially with the number of binary interactions considered. As a result, only large reaction networks capture all possible consequences of the molecular interactions that occur in a cell signaling system, which is problematic because traditional modeling approaches for chemical kinetics (e.g., ordinary differential equations) require explicit network specification. This problem is circumvented through representation of interactions in terms of local rules. With this approach, network specification is implicit and model specification is concise. Concise representation results in a coarse graining of chemical kinetics, which is introduced because all reactions implied by a rule inherit the rate law associated with that rule. Coarse graining can be appropriate if interactions are modular, and the coarseness of a model can be adjusted as needed. Rules can be specified using specialized model-specification languages, and recently developed tools designed for specification of rule-based models allow one to leverage powerful software engineering capabilities. A rule-based model comprises a set of rules, which can be processed by general-purpose simulation and analysis tools to achieve different objectives (e.g., to perform either a deterministic or stochastic simulation).

Multiscale modeling and synaptic plasticity

Synaptic plasticity is a major convergence point for theory and computation, and the process of plasticity engages physiology, cell, and molecular biology. In its many manifestations, plasticity is at the hub of basic neuroscience questions about memory and development, as well as more medically themed questions of neural damage and recovery. As an important cellular locus of memory, synaptic plasticity has received a huge amount of experimental and theoretical attention. If computational models have tended to pick specific aspects of plasticity, such as STDP, and reduce them to an equation, some experimental studies are equally guilty of oversimplification each time they identify a new molecule and declare it to be the last word in plasticity and learning. Multiscale modeling begins with the acknowledgment that synaptic function spans many levels of signaling, and these are so tightly coupled that we risk losing essential features of plasticity if we focus exclusively on any one level. Despite the technical challenges and gaps in data for model specification, an increasing number of multiscale modeling studies have taken on key questions in plasticity. These have provided new insights, but importantly, they have opened new avenues for questioning. This review discusses a wide range of multiscale models in plasticity, including their technical landscape and their implications.

Efficient stochastic simulation of chemical kinetics networks using a weighted ensemble of trajectories

We apply the "weighted ensemble" (WE) simulation strategy, previously employed in the context of molecular dynamics simulations, to a series of systems-biology models that range in complexity from a one-dimensional system to a system with 354 species and 3680 reactions. WE is relatively easy to implement, does not require extensive hand-tuning of parameters, does not depend on the details of the simulation algorithm, and can facilitate the simulation of extremely rare events. For the coupled stochastic reaction systems we study, WE is able to produce accurate and efficient approximations of the joint probability distribution for all chemical species for all time t. WE is also able to efficiently extract mean first passage times for the systems, via the construction of a steady-state condition with feedback. In all cases studied here, WE results agree with independent "brute-force" calculations, but significantly enhance the precision with which rare or slow processes can be characterized. Speedups over "brute-force" in sampling rare events via the Gillespie direct Stochastic Simulation Algorithm range from ~10(12) to ~10(18) for characterizing rare states in a distribution, and ~10(2) to ~10(4) for finding mean first passage times.

The autoregulation of a eukaryotic DNA transposon

How do DNA transposons live in harmony with their hosts? Bacteria provide the only documented mechanisms for autoregulation, but these are incompatible with eukaryotic cell biology. Here we show that autoregulation of Hsmar1 operates during assembly of the transpososome and arises from the multimeric state of the transposase, mediated by a competition for binding sites. We explore the dynamics of a genomic invasion using a computer model, supported by in vitro and in vivo experiments, and show that amplification accelerates at first but then achieves a constant rate. The rate is proportional to the genome size and inversely proportional to transposase expression and its affinity for the transposon ends. Mariner transposons may therefore resist post-transcriptional silencing. Because regulation is an emergent property of the reaction it is resistant to selfish exploitation. The behavior of distantly related eukaryotic transposons is consistent with the same mechanism, which may therefore be widely applicable.

DOI:http://dx.doi.org/10.7554/eLife.00668.001

Transposons are regions of mobile DNA that can jump from one location in the genome to another. This represents a genetic burden to the host because there is always the risk that the transposon will inactivate a cellular gene. However, a greater problem is that transposition is accompanied by an increase in the number of copies of the transposon. Since each new copy will be a source of further new copies, amplification of transposons is necessarily exponential. The fact that eukaryotic cells are able to tolerate DNA transposons suggests the existence of regulatory mechanisms to defuse the inevitable genomic melt-down. Host-mediated epigenetic modifications and RNA interference will provide some level of protection. However, they are by no means completely effective and a well-adapted genomic parasite, such as a transposon, might be expected to have its own mechanism of regulation.

Now, Claeys Bouuaert, Lipkow and colleagues have used a computer model in combination with in vivo and in vitro experiments to search for this mechanism. Their experiments reveal how a DNA transposon is down-regulated by its own transposase. The transposase is the enzyme that catalyzes the ‘jump’ or transposition. It binds to specific sites at either end of the transposon and brings these together to make up a nucleoprotein complex called the transpososome. It is within this complex that the chemical steps of the reaction take place. When the number of transposons increases, so does the concentration of transposase. Claeys Bouuaert et al. show that the binding sites become saturated at a relatively low transposase concentration and that negative regulation arises from the resulting competition. Thus, the rate of transposition decreases as the number of transposons increases. They further use the computer model to explore how the amplification of the transposon is affected by transposon-specific and cellular-specific factors.

Claeys Bouuaert, Lipkow and colleagues based their study predominantly on a resurrected copy of the Hsmar1 transposon, which was active in the human genome 50 million years ago. However, they also tested two distantly related eukaryotic transposons and observed that their behavior was similar, which suggests that this could be a general mechanism that controls the activity of jumping genes. They also note that their competition mechanism is conceptually similar to the immunological ‘prozone effect’. This is a recurrent theme in protein chemistry and demonstrates once again that less is in fact sometimes more.

DOI:http://dx.doi.org/10.7554/eLife.00668.002

Modeling and simulation of biological systems from image data

This essay provides an introduction to the terminology, concepts, methods, and challenges of image-based modeling in biology. Image-based modeling and simulation aims at using systematic, quantitative image data to build predictive models of biological systems that can be simulated with a computer. This allows one to disentangle molecular mechanisms from effects of shape and geometry. Questions like "what is the functional role of shape" or "how are biological shapes generated and regulated" can be addressed in the framework of image-based systems biology. The combination of image quantification, model building, and computer simulation is illustrated here using the example of diffusion in the endoplasmic reticulum.

Modeling cellular compartmentation in one-carbon metabolism

Folate-mediated one-carbon metabolism (FOCM) is associated with risk for numerous pathological states including birth defects, cancers, and chronic diseases. Although the enzymes that constitute the biological pathways have been well described and their interdependency through the shared use of folate cofactors appreciated, the biological mechanisms underlying disease etiologies remain elusive. The FOCM network is highly sensitive to nutritional status of several B-vitamins and numerous penetrant gene variants that alter network outputs, but current computational approaches do not fully capture the dynamics and stochastic noise of the system. Combining the stochastic approach with a rule-based representation will help model the intrinsic noise displayed by FOCM, address the limited flexibility of standard simulation methods for coarse-graining the FOCM-associated biochemical processes, and manage the combinatorial complexity emerging from reactions within FOCM that would otherwise be intractable.

Detailed simulations of cell biology with Smoldyn 2.1

Most cellular processes depend on intracellular locations and random collisions of individual protein molecules. To model these processes, we developed algorithms to simulate the diffusion, membrane interactions, and reactions of individual molecules, and implemented these in the Smoldyn program. Compared to the popular MCell and ChemCell simulators, we found that Smoldyn was in many cases more accurate, more computationally efficient, and easier to use. Using Smoldyn, we modeled pheromone response system signaling among yeast cells of opposite mating type. This model showed that secreted Bar1 protease might help a cell identify the fittest mating partner by sharpening the pheromone concentration gradient. This model involved about 200,000 protein molecules, about 7000 cubic microns of volume, and about 75 minutes of simulated time; it took about 10 hours to run. Over the next several years, as faster computers become available, Smoldyn will allow researchers to model and explore systems the size of entire bacterial and smaller eukaryotic cells.