An image-based model of calcium waves in differentiated neuroblastoma cells.

Calcium waves produced by bradykinin-induced inositol-1,4, 5-trisphosphate (InsP(3))-mediated release from endoplasmic reticulum (ER) have been imaged in N1E-115 neuroblastoma cells. A model of this process was built using the "virtual cell," a general computational system for integrating experimental image, biochemical, and electrophysiological data. The model geometry was based on a cell for which the calcium wave had been experimentally recorded. The distributions of the relevant cellular components [InsP(3) receptor (InsP(3)R)], sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) pumps, bradykinin receptors, and ER] were based on 3D confocal immunofluorescence images. Wherever possible, known biochemical and electrophysiological data were used to constrain the model. The simulation closely matched the spatial and temporal characteristics of the experimental calcium wave. Predictions on different patterns of calcium signals after InsP(3) uncaging or for different cell geometries were confirmed experimentally, thus helping to validate the model. Models in which the spatial distributions of key components are altered suggest that initiation of the wave in the center of the neurite derives from an interplay of soma-biased ER distribution and InsP(3) generation biased toward the neurite. Simulations demonstrate that mobile buffers (like the indicator fura-2) significantly delay initiation and lower the amplitude of the wave. Analysis of the role played by calcium diffusion indicated that the speed of the wave is only slightly dependent on the ability of calcium to diffuse to and activate neighboring InsP(3) receptor sites.

Morphological control of inositol-1,4,5-trisphosphate-dependent signals.

Inositol-1,4,5-trisphosphate (InsP(3))-mediated calcium signals represent an important mechanism for transmitting external stimuli to the cell. However, information about intracellular spatial patterns of InsP(3) itself is not generally available. In particular, it has not been determined how the interplay of InsP(3) generation, diffusion, and degradation within complex cellular geometries can control the patterns of InsP(3) signaling. Here, we explore the spatial and temporal characteristics of [InsP(3)](cyt) during a bradykinin-induced calcium wave in a neuroblastoma cell. This is achieved by using a unique image-based computer modeling system, Virtual Cell, to integrate experimental data on the rates and spatial distributions of the key molecular components of the process. We conclude that the characteristic calcium dynamics requires rapid, high-amplitude production of [InsP(3)](cyt) in the neurite. This requisite InsP(3) spatiotemporal profile is provided, in turn, as an intrinsic consequence of the cell's morphology, demonstrating how geometry can locally and dramatically intensify cytosolic signals that originate at the plasma membrane. In addition, the model predicts, and experiments confirm, that stimulation of just the neurite, but not the soma or growth cone, is sufficient to generate a calcium response throughout the cell.

Determination of time-dependent inositol-1,4,5-trisphosphate concentrations during calcium release in a smooth muscle cell.

The level of [InsP3]cyt required for calcium release in A7r5 cells, a smooth muscle cell line, was determined by a new set of procedures using quantitative confocal microscopy to measure release of InsP3 from cells microinjected with caged InsP3. From these experiments, the [InsP3]cyt required to evoke a half-maximal calcium response is 100 nM. Experiments with caged glycerophosphoryl-myo-inositol 4, 5-bisphosphate (GPIP2), a slowly metabolized analogue of InsP3, gave a much slower recovery and a half-maximal response of an order of magnitude greater than InsP3. Experimental data and highly constrained variables were used to construct a mathematical model of the InsP3-dependent [Ca2+]cyt changes; the resulting simulations show high fidelity to experiment. Among the elements considered in constructing this model were the mechanism of the InsP3-receptor, InsP3 degradation, calcium buffering in the cytosol, and refilling of the ER stores via sarcoplasmic endoplasmic reticulum ATPase (SERCA) pumps. The model predicts a time constant of 0.8 s for InsP3 degradation and 13 s for GPIP2. InsP3 degradation was found to be a prerequisite for [Ca2+]cyt recovery to baseline levels and is therefore critical to the pattern of the overall [Ca2+]cyt signal. Analysis of the features of this model provides insights into the individual factors controlling the amplitude and shape of the InsP3-mediated calcium signal.

A general computational framework for modeling cellular structure and function.

The "Virtual Cell" provides a general system for testing cell biological mechanisms and creates a framework for encapsulating the burgeoning knowledge base comprising the distribution and dynamics of intracellular biochemical processes. It approaches the problem by associating biochemical and electrophysiological data describing individual reactions with experimental microscopic image data describing their subcellular localizations. Individual processes are collected within a physical and computational infrastructure that accommodates any molecular mechanism expressible as rate equations or membrane fluxes. An illustration of the method is provided by a dynamic simulation of IP3-mediated Ca2+ release from endoplasmic reticulum in a neuronal cell. The results can be directly compared to experimental observations and provide insight into the role of experimentally inaccessible components of the overall mechanism.