Abstract
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.