Diffusion on a Curved Surface Coupled to Diffusion in the Volume: Application to Cell Biology

Midpoint projection algorithm for stochastic differential equations on manifolds

Stochastic differential equations projected onto manifolds occur in physics, chemistry, biology, engineering, nanotechnology, and optimization, with interdisciplinary applications. Intrinsic coordinate stochastic equations on the manifold are sometimes computationally impractical, and numerical projections are therefore useful in many cases. In this paper a combined midpoint projection algorithm is proposed that uses a midpoint projection onto a tangent space, combined with a subsequent normal projection to satisfy the constraints. We also show that the Stratonovich form of stochastic calculus is generally obtained with finite bandwidth noise in the presence of a strong enough external potential that constrains the resulting physical motion to a manifold. Numerical examples are given for a wide range of manifolds, including circular, spheroidal, hyperboloidal, and catenoidal cases, higher-order polynomial constraints that give a quasicubical surface, and a ten-dimensional hypersphere. In all cases the combined midpoint method has greatly reduced errors compared to other methods used for comparison, namely, a combined Euler projection approach and a tangential projection algorithm. We derive intrinsic stochastic equations for spheroidal and hyperboloidal surfaces for comparison purposes to verify the results. Our technique can handle multiple constraints, which allows manifolds that embody several conserved quantities. The algorithm is accurate, simple, and efficient. A reduction of an order of magnitude in the diffusion distance error is found compared to the other methods and an up to several orders of magnitude reduction in constraint function errors.

Thermodynamically Consistent Models for Coupled Bulk and Surface Dynamics

We present a constructive paradigm to derive thermodynamically consistent models coupling the bulk and surface dynamics hierarchically following the generalized Onsager principle. In the model, the bulk and surface thermodynamical variables are allowed to be different and the free energy of the model comprises the bulk, surface, and coupling energy, which can be weakly or strongly non-local. We illustrate the paradigm using a phase field model for binary materials and show that the model includes the existing thermodynamically consistent ones for the binary material system in the literature as special cases. In addition, we present a set of such phase field models for a few selected mobility operators and free energies to show how boundary dynamics impart changes to bulk dynamics and vice verse. As an example, we show numerically how reactive transport on the boundary impacts the dynamics in the bulk using a reactive transport model for binary reactive fluids by adopting a structure-preserving algorithm to solve the model equations in a rectangular domain.

Reaction-diffusion models in weighted and directed connectomes

Connectomes represent comprehensive descriptions of neural connections in a nervous system to better understand and model central brain function and peripheral processing of afferent and efferent neural signals. Connectomes can be considered as a distinctive and necessary structural component alongside glial, vascular, neurochemical, and metabolic networks of the nervous systems of higher organisms that are required for the control of body functions and interaction with the environment. They are carriers of functional epiphenomena such as planning behavior and cognition, which are based on the processing of highly dynamic neural signaling patterns. In this study, we examine more detailed connectomes with edge weighting and orientation properties, in which reciprocal neuronal connections are also considered. Diffusion processes are a further necessary condition for generating dynamic bioelectric patterns in connectomes. Based on our high-precision connectome data, we investigate different diffusion-reaction models to study the propagation of dynamic concentration patterns in control and lesioned connectomes. Therefore, differential equations for modeling diffusion were combined with well-known reaction terms to allow the use of connection weights, connectivity orientation and spatial distances.

Three reaction-diffusion systems Gray-Scott, Gierer-Meinhardt and Mimura-Murray were investigated. For this purpose, implicit solvers were implemented in a numerically stable reaction-diffusion system within the framework of neuroVIISAS. The implemented reaction-diffusion systems were applied to a subconnectome which shapes the mechanosensitive pathway that is strongly affected in the multiple sclerosis demyelination disease. It was found that demyelination modeling by connectivity weight modulation changes the oscillations of the target region, i.e. the primary somatosensory cortex, of the mechanosensitive pathway.

In conclusion, a new application of reaction-diffusion systems to weighted and directed connectomes has been realized. Because the implementation were performed in the neuroVIISAS framework many possibilities for the study of dynamic reaction-diffusion processes in empirical connectomes as well as specific randomized network models are available now.

Fission Yeast Polarization: Modeling Cdc42 Oscillations, Symmetry Breaking, and Zones of Activation and Inhibition

Cells polarize for growth, motion, or mating through regulation of membrane-bound small GTPases between active GTP-bound and inactive GDP-bound forms. Activators (GEFs, GTP exchange factors) and inhibitors (GAPs, GTPase activating proteins) provide positive and negative feedbacks. We show that a reaction-diffusion model on a curved surface accounts for key features of polarization of model organism fission yeast. The model implements Cdc42 membrane diffusion using measured values for diffusion coefficients and dissociation rates and assumes a limiting GEF pool (proteins Gef1 and Scd1), as in prior models for budding yeast. The model includes two types of GAPs, one representing tip-localized GAPs, such as Rga3; and one representing side-localized GAPs, such as Rga4 and Rga6, that we assume switch between fast and slow diffusing states. After adjustment of unknown rate constants, the model reproduces active Cdc42 zones at cell tips and the pattern of GEF and GAP localization at cell tips and sides. The model reproduces observed tip-to-tip oscillations with periods of the order of several minutes, as well as asymmetric to symmetric oscillations transitions (corresponding to NETO "new end take off"), assuming the limiting GEF amount increases with cell size.

VoroCrust: Voronoi Meshing Without Clipping

Polyhedral meshes are increasingly becoming an attractive option with particular advantages over traditional meshes for certain applications. What has been missing is a robust polyhedral meshing algorithm that can handle broad classes of domains exhibiting arbitrary curved boundaries and sharp features. In addition, the power of primal-dual mesh pairs, exemplified by Voronoi-Delaunay meshes, has been recognized as an important ingredient in numerous formulations. The VoroCrust algorithm is the first provably correct algorithm for conforming Voronoi meshing for non-convex and possibly non-manifold domains with guarantees on the quality of both surface and volume elements. A robust refinement process estimates a suitable sizing field that enables the careful placement of Voronoi seeds across the surface circumventing the need for clipping and avoiding its many drawbacks. The algorithm has the flexibility of filling the interior by either structured or random samples, while all sharp features are preserved in the output mesh. We demonstrate the capabilities of the algorithm on a variety of models and compare against state-of-the-art polyhedral meshing methods based on clipped Voronoi cells establishing the clear advantage of VoroCrust output.

Exploration and stabilization of Ras1 mating zone: A mechanism with positive and negative feedbacks

In mating fission yeast cells, sensing and response to extracellular pheromone concentrations occurs through an exploratory Cdc42 patch that stochastically samples the cell cortex before stabilizing towards a mating partner. Active Ras1 (Ras1-GTP), an upstream regulator of Cdc42, and Gap1, the GTPase-activating protein for Ras1, localize at the patch. We developed a reaction-diffusion model of Ras1 patch appearance and disappearance with a positive feedback by a Guanine nucleotide Exchange Factor (GEF) and Gap1 inhibition. The model is based on new estimates of Ras1-GDP, Ras1-GTP and Gap1 diffusion coefficients and rates of cytoplasmic exchange studied by FRAP. The model reproduces exploratory patch behavior and lack of Ras1 patch in cells lacking Gap1. Transition to a stable patch can occur by change of Gap1 rates constants or local increase of the positive feedback rate constants. The model predicts that the patch size and number of patches depend on the strength of positive and negative feedbacks. Measurements of Ras1 patch size and number in cells overexpressing the Ras1 GEF or Gap1 are consistent with the model.

Unicellular fission yeasts mate by fusing with partners of the opposite mating type. Each pair member grows towards its selected partner that signals its presence through secreted pheromone. The process of partner selection occurs through an exploratory patch (containing activated signaling protein Cdc42 and upstream regulator Ras1) that assembles and disassembles on the cell cortex, stabilizing in regions of higher opposite pheromone concentration. We present a computational model of the molecular mechanisms driving the dynamical pattern of patch exploration and stabilization. The model is based on reaction and diffusion along the curved cell membrane, with diffusion coefficients measured experimentally. In the model, a positive Ras1 activation feedback loop generates a patch containing most of the activating protein (Ras1 GEF). The fast diffusing inhibitor Gap1 that is recruited locally from the cytoplasm spreads on the cell membrane, limiting patch size and causing its decay. Spontaneous reinitiation of Ras1 activation elsewhere on the cortex provides a mechanism for exploration. Transition of the system’s behavior to that of a single stable patch is possible upon simulated pheromone sensing. The computational model provides predictions for the number of patches and patch size dependence on parameters that we tested experimentally.

Numerical Approach to Spatial Deterministic-Stochastic Models Arising in Cell Biology

Hybrid deterministic-stochastic methods provide an efficient alternative to a fully stochastic treatment of models which include components with disparate levels of stochasticity. However, general-purpose hybrid solvers for spatially resolved simulations of reaction-diffusion systems are not widely available. Here we describe fundamentals of a general-purpose spatial hybrid method. The method generates realizations of a spatially inhomogeneous hybrid system by appropriately integrating capabilities of a deterministic partial differential equation solver with a popular particle-based stochastic simulator, Smoldyn. Rigorous validation of the algorithm is detailed, using a simple model of calcium 'sparks' as a testbed. The solver is then applied to a deterministic-stochastic model of spontaneous emergence of cell polarity. The approach is general enough to be implemented within biologist-friendly software frameworks such as Virtual Cell.

A computational method for the coupled solution of reaction-diffusion equations on evolving domains and manifolds: Application to a model of cell migration and chemotaxis

In this paper, we devise a moving mesh finite element method for the approximate solution of coupled bulk-surface reaction-diffusion equations on an evolving two dimensional domain. Fundamental to the success of the method is the robust generation of bulk and surface meshes. For this purpose, we use a novel moving mesh partial differential equation (MMPDE) approach. The developed method is applied to model problems with known analytical solutions; these experiments indicate second-order spatial and temporal accuracy. Coupled bulk-surface problems occur frequently in many areas; in particular, in the modelling of eukaryotic cell migration and chemotaxis. We apply the method to a model of the two-way interaction of a migrating cell in a chemotactic field, where the bulk region corresponds to the extracellular region and the surface to the cell membrane.

Stability analysis and simulations of coupled bulk-surface reaction-diffusion systems

In this article, we formulate new models for coupled systems of bulk-surface reaction-diffusion equations on stationary volumes. The bulk reaction-diffusion equations are coupled to the surface reaction-diffusion equations through linear Robin-type boundary conditions. We then state and prove the necessary conditions for diffusion-driven instability for the coupled system. Owing to the nature of the coupling between bulk and surface dynamics, we are able to decouple the stability analysis of the bulk and surface dynamics. Under a suitable choice of model parameter values, the bulk reaction-diffusion system can induce patterning on the surface independent of whether the surface reaction-diffusion system produces or not, patterning. On the other hand, the surface reaction-diffusion system cannot generate patterns everywhere in the bulk in the absence of patterning from the bulk reaction-diffusion system. For this case, patterns can be induced only in regions close to the surface membrane. Various numerical experiments are presented to support our theoretical findings. Our most revealing numerical result is that, Robin-type boundary conditions seem to introduce a boundary layer coupling the bulk and surface dynamics.

A conservative algorithm for parabolic problems in domains with moving boundaries

We describe a novel conservative algorithm for parabolic problems in domains with moving boundaries developed for modeling in cell biology. The spatial discretization is accomplished by applying Voronoi decomposition to a fixed rectangular grid. In the vicinity of the boundary, the procedure generates irregular Voronoi cells that conform to the domain shape and merge seamlessly with regular control volumes in the domain interior. Consequently, our algorithm is free of the CFL stability issue due to moving interfaces and does not involve cell-merging or mass redistribution. Local mass conservation is ensured by finite-volume discretization and natural-neighbor interpolation. Numerical experiments with two-dimensional geometries demonstrate exact mass conservation and indicate an order of convergence in space between one and two. The use of standard meshing techniques makes extension of the method to three dimensions conceptually straightforward.

Simple computation of reaction-diffusion processes on point clouds

The study of reaction-diffusion processes is much more complicated on general curved surfaces than on standard Cartesian coordinate spaces. Here we show how to formulate and solve systems of reaction-diffusion equations on surfaces in an extremely simple way, using only the standard Cartesian form of differential operators, and a discrete unorganized point set to represent the surface. Our method decouples surface geometry from the underlying differential operators. As a consequence, it becomes possible to formulate and solve rather general reaction-diffusion equations on general surfaces without having to consider the complexities of differential geometry or sophisticated numerical analysis. To illustrate the generality of the method, computations for surface diffusion, pattern formation, excitable media, and bulk-surface coupling are provided for a variety of complex point cloud surfaces.

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.

Computational analysis of Rho GTPase cycling

The Rho family of GTPases control actin organization during diverse cellular responses (migration, cytokinesis and endocytosis). Although the primary members of this family (RhoA, Rac and Cdc42) have different downstream effects on actin remodeling, the basic mechanism involves targeting to the plasma membrane and activation by GTP binding. Our hypothesis is that the details of GTPase cycling between membrane and cytosol are key to the differential upstream regulation of these biochemical switches. Accordingly, we developed a modeling framework to analyze experimental data for these systems. This analysis can reveal details of GDI-mediated cycling and help distinguish between GDI-dependent and -independent mechanisms, including vesicle trafficking and direct association-dissociation of GTPase with membrane molecules. Analysis of experimental data for Rac membrane cycling reveals that the lower apparent affinity of GDI for RacGTP compared to RacGDP can be fully explained by the faster dissociation of the latter from the membrane. Non-dimensional steady-state solutions for membrane fraction of GTPase are presented in multidimensional charts. This methodology is then used to analyze glucose stimulated Rac cycling in pancreatic β-cells. The charts are used to illustrate the effects of GEFs/GAPs and regulated affinities between GTPases and membrane and/or GDI on the amount of membrane bound GTPase. In a similar fashion, the charts can be used as a guide in assessing how targeted modifications may compensate for altered GTPase-GDI balance in disease scenarios.

Among the functions of the small GTPases Rac, RhoA and Cdc42 are the regulation of protein traffic, insulin secretion, cell shape, survival and motility. The last two are important steps for tumor growth and metastasis. The function of these proteins relies on their expression levels, proper membrane localization and activation. In addition, all three proteins compete for the same protein ‘GDI’, which modulates their cycling. These proteins are ubiquitous in mammalian cells, but also studied in simpler in vitro systems and cultured yeast. Here we show, using a series of computational analyses, that for each of these experimental systems the dominant pathway for membrane cycling of GTPases seems to differ. This means that the researcher interested in the physiological function of any of those proteins must make sure that the experimental system is appropriate. We present a methodology to identify the dominant pathways by measuring the apparent membrane dissociation rate of the protein as a function of GDI concentration. We provide charts generated from parametric scans. This analysis is then applied to the Rac-dependent insulin secretion pathway in pancreatic ß-cells, revealing that direct signaling between Rac and the membrane is an essential mechanism that emerges from the data.

Computational modeling of cellular signaling processes embedded into dynamic spatial contexts

Cellular signaling processes depend on specific spatiotemporal distributions of their molecular components. Multi-color high-resolution microscopy now permits detailed assessment of such distributions, providing the input for fine-grained computational models that explore the mechanisms governing dynamic assembly of multi-molecular complexes and their role in shaping cellular behavior. However, incorporating into such models both complex molecular reaction cascades and the spatial localization of signaling components within dynamic cellular morphologies presents substantial challenges. Here we introduce an approach that addresses these challenges by automatically generating computational representations of complex reaction networks based on simple bi-molecular interaction rules embedded into detailed, adaptive models of cellular morphology. Using examples of receptor-mediated cellular adhesion and signal-induced localized MAPK activation in yeast, we illustrate the capacity of this simulation technique to provide insights into cell biological processes. The modeling algorithms, implemented in a version of the Simmune tool set, are accessible through intuitive graphical interfaces as well as programming libraries.

Spatial modeling of cell signaling networks

The shape of a cell, the sizes of subcellular compartments, and the spatial distribution of molecules within the cytoplasm can all control how molecules interact to produce a cellular behavior. This chapter describes how these spatial features can be included in mechanistic mathematical models of cell signaling. The Virtual Cell computational modeling and simulation software is used to illustrate the considerations required to build a spatial model. An explanation of how to appropriately choose between physical formulations that implicitly or explicitly account for cell geometry and between deterministic versus stochastic formulations for molecular dynamics is provided, along with a discussion of their respective strengths and weaknesses. As a first step toward constructing a spatial model, the geometry needs to be specified and associated with the molecules, reactions, and membrane flux processes of the network. Initial conditions, diffusion coefficients, velocities, and boundary conditions complete the specifications required to define the mathematics of the model. The numerical methods used to solve reaction-diffusion problems both deterministically and stochastically are then described and some guidance is provided in how to set up and run simulations. A study of cAMP signaling in neurons ends the chapter, providing an example of the insights that can be gained in interpreting experimental results through the application of spatial modeling.

Virtual Cell: computational tools for modeling in cell biology

The Virtual Cell (VCell) is a general computational framework for modeling physicochemical and electrophysiological processes in living cells. Developed by the National Resource for Cell Analysis and Modeling at the University of Connecticut Health Center, it provides automated tools for simulating a wide range of cellular phenomena in space and time, both deterministically and stochastically. These computational tools allow one to couple electrophysiology and reaction kinetics with transport mechanisms, such as diffusion and directed transport, and map them onto spatial domains of various shapes, including irregular three-dimensional geometries derived from experimental images. In this article, we review new robust computational tools recently deployed in VCell for treating spatially resolved models.

Modeling capping protein FRAP and CALI experiments reveals in vivo regulation of actin dynamics

To gain insights on cellular mechanisms regulating actin polymerization, we used the Virtual Cell to model fluorescence recovery after photobleaching (FRAP) and chromophore-assisted laser inactivation (CALI) experiments on EGFP-capping protein (EGFP-CP). Modeling the FRAP kinetics demonstrated that the in vivo rate for the dissociation of CP from actin filaments is much faster (approximately 0.1 s(-1)) than that measured in vitro (0.01-0.0004 s(-1)). The CALI simulation revealed that in order to induce sustainable changes in cell morphology after CP inactivation, the cells should exhibit anticapping ability. We included the VASP protein as the anticapping agent in the modeling scheme. The model predicts that VASP affinity for barbed ends has a cooperative dependence on the concentration of VASP-barbed end complexes. This dependence produces a positive feedback that stabilizes the complexes and allows sustained growth at clustered filament tips. We analyzed the range of laser intensities that are sufficient to induce changes in cell morphology. This analysis demonstrates that FRAP experiments with EGFP-CP can be performed safely without changes in cell morphology, because, the intensity of the photobleaching beam is not high enough to produce the critical concentration of free barbed ends that will induce filament growth before diffusional replacement of EGFP-CP occurs.

The Moving Boundary Node Method: A level set-based, finite volume algorithm with applications to cell motility

Eukaryotic cell crawling is a highly complex biophysical and biochemical process, where deformation and motion of a cell are driven by internal, biochemical regulation of a poroelastic cytoskeleton. One challenge to building quantitative models that describe crawling cells is solving the reaction-diffusion-advection dynamics for the biochemical and cytoskeletal components of the cell inside its moving and deforming geometry. Here we develop an algorithm that uses the level set method to move the cell boundary and uses information stored in the distance map to construct a finite volume representation of the cell. Our method preserves Cartesian connectivity of nodes in the finite volume representation while resolving the distorted cell geometry. Derivatives approximated using a Taylor series expansion at finite volume interfaces lead to second order accuracy even on highly distorted quadrilateral elements. A modified, Laplacian-based interpolation scheme is developed that conserves mass while interpolating values onto nodes that join the cell interior as the boundary moves. An implicit time-stepping algorithm is used to maintain stability. We use the algoirthm to simulate two simple models for cellular crawling. The first model uses depolymerization of the cytoskeleton to drive cell motility and suggests that the shape of a steady crawling cell is strongly dependent on the adhesion between the cell and the substrate. In the second model, we use a model for chemical signalling during chemotaxis to determine the shape of a crawling cell in a constant gradient and to show cellular response upon gradient reversal.

Use of virtual cell in studies of cellular dynamics

The Virtual Cell (VCell) is a unique computational environment for modeling and simulation of cell biology. It has been specifically designed to be a tool for a wide range of scientists, from experimental cell biologists to theoretical biophysicists. The models created with VCell can range from the simple, to evaluate hypotheses or to interpret experimental data, to complex multilayered models used to probe the predicted behavior of spatially resolved, highly nonlinear systems. In this chapter, we discuss modeling capabilities of VCell and demonstrate representative examples of the models published by the VCell users.

Diffusion in cytoplasm: effects of excluded volume due to internal membranes and cytoskeletal structures

The intricate geometry of cytoskeletal networks and internal membranes causes the space available for diffusion in cytoplasm to be convoluted, thereby affecting macromolecule diffusivity. We present a first systematic computational study of this effect by approximating intracellular structures as mixtures of random overlapping obstacles of various shapes. Effective diffusion coefficients are computed using a fast homogenization technique. It is found that a simple two-parameter power law provides a remarkably accurate description of effective diffusion over the entire range of volume fractions and for any given composition of structures. This universality allows for fast computation of diffusion coefficients, once the obstacle shapes and volume fractions are specified. We demonstrate that the excluded volume effect alone can account for a four-to-sixfold reduction in diffusive transport in cells, relative to diffusion in vitro. The study lays the foundation for an accurate coarse-grain formulation that would account for cytoplasm heterogeneity on a micron scale and binding of tracers to intracellular structures.

An open model of actin dendritic nucleation

The availability of quantitative experimental data on the kinetics of actin assembly has enabled the construction of many mathematical models focused on explaining specific behaviors of this complex system. However these ad hoc models are generally not reusable or accessible by the large community of actin biologists. In this work, we present a comprehensive model that integrates and unifies much of the in vitro data on the components of the dendritic nucleation mechanism for actin dynamics. More than 300 simulations have been run based on compartmental and three-dimensional spatial versions of this model. Several key findings are highlighted, including an explanation for the sharp boundary between actin assembly and disassembly in the lamellipodia of migrating cells. Because this model, with the simulation results, is "open source", in the sense that it is publicly available and editable through the Virtual Cell database (http://vcell.org), it can be accessed, analyzed, modified, and extended.

SOLVING PDES IN COMPLEX GEOMETRIES: A DIFFUSE DOMAIN APPROACH

We extend previous work and present a general approach for solving partial differential equations in complex, stationary, or moving geometries with Dirichlet, Neumann, and Robin boundary conditions. Using an implicit representation of the geometry through an auxilliary phase field function, which replaces the sharp boundary of the domain with a diffuse layer (e.g. diffuse domain), the equation is reformulated on a larger regular domain. The resulting partial differential equation is of the same order as the original equation, with additional lower order terms to approximate the boundary conditions. The reformulated equation can be solved by standard numerical techniques. We use the method of matched asymptotic expansions to show that solutions of the re-formulated equations converge to those of the original equations. We provide numerical simulations which confirm this analysis. We also present applications of the method to growing domains and complex three-dimensional structures and we discuss applications to cell biology and heteroepitaxy.

Virtual Cell modelling and simulation software environment

The Virtual Cell (VCell; http://vcell.org/) is a problem solving environment, built on a central database, for analysis, modelling and simulation of cell biological processes. VCell integrates a growing range of molecular mechanisms, including reaction kinetics, diffusion, flow, membrane transport, lateral membrane diffusion and electrophysiology, and can associate these with geometries derived from experimental microscope images. It has been developed and deployed as a web-based, distributed, client-server system, with more than a thousand world-wide users. VCell provides a separation of layers (core technologies and abstractions) representing biological models, physical mechanisms, geometry, mathematical models and numerical methods. This separation clarifies the impact of modelling decisions, assumptions and approximations. The result is a physically consistent, mathematically rigorous, spatial modelling and simulation framework. Users create biological models and VCell will automatically (i) generate the appropriate mathematical encoding for running a simulation and (ii) generate and compile the appropriate computer code. Both deterministic and stochastic algorithms are supported for describing and running non-spatial simulations; a full partial differential equation solver using the finite volume numerical algorithm is available for reaction-diffusion-advection simulations in complex cell geometries including 3D geometries derived from microscope images. Using the VCell database, models and model components can be reused and updated, as well as privately shared among collaborating groups, or published. Exchange of models with other tools is possible via import/export of SBML, CellML and MatLab formats. Furthermore, curation of models is facilitated by external database binding mechanisms for unique identification of components and by standardised annotations compliant with the MIRIAM standard. VCell is now open source, with its native model encoding language (VCML) being a public specification, which stands as the basis for a new generation of more customised, experiment-centric modelling tools using a new plug-in based platform.

Analysis of phosphatidylinositol-4,5-bisphosphate signaling in cerebellar Purkinje spines

A 3D model was developed and used to explore dynamics of phosphatidylinositol-4,5-bisphosphate (PIP2) signaling in cerebellar Purkinje neurons. Long-term depression in Purkinje neurons depends on coincidence detection of climbing fiber stimulus evoking extracellular calcium flux into the cell and parallel fiber stimulus evoking inositol-1,4,5-trisphosphate (IP3)-meditated calcium release from the endoplasmic reticulum. Experimental evidence shows that large concentrations of IP3 are required for calcium release. This study uses computational analysis to explore how the Purkinje cell provides sufficient PIP2 to produce large amounts of IP3. Results indicate that baseline PIP2 concentration levels in the plasma membrane are inadequate, even if the model allows for PIP2 replenishment by lateral diffusion from neighboring dendrite membrane. Lateral diffusion analysis indicates apparent anomalous diffusion of PIP2 in the spiny dendrite membrane, due to restricted diffusion through spine necks. Stimulated PIP2 synthesis and elevated spine PIP2 mediated by a local sequestering protein were explored as candidate mechanisms to supply sufficient PIP2. Stimulated synthesis can indeed lead to high IP3 amplitude of long duration; local sequestration produces high IP3 amplitude, but of short duration. Simulation results indicate that local sequestration could explain the experimentally observed finely tuned timing between parallel fiber and climbing fiber activation.