A simple numerical model of calcium spark formation and detection in cardiac myocytes.

The elementary events of excitation-contraction coupling in heart muscle are Ca2+ sparks, which arise from one or more ryanodine receptors in the sarcoplasmic reticulum (SR). Here a simple numerical model is constructed to explore Ca2+ spark formation, detection, and interpretation in cardiac myocytes. This model includes Ca2+ release, cytosolic diffusion, resequestration by SR Ca2+-ATPases, and the association and dissociation of Ca2+ with endogenous Ca2+-binding sites and a diffusible indicator dye (fluo-3). Simulations in a homogeneous, isotropic cytosol reproduce the brightness and the time course of a typical cardiac Ca2+ spark, but underestimate its spatial size (approximately 1.1 micron vs. approximately 2.0 micron). Back-calculating [Ca2+]i by assuming equilibrium with indicator fails to provide a good estimate of the free Ca2+ concentration even when using blur-free fluorescence data. A parameter sensitivity study reveals that the mobility, kinetics, and concentration of the indicator are essential determinants of the shape of Ca2+ sparks, whereas the stationary buffers and pumps are less influential. Using a geometrically more complex version of the model, we show that the asymmetric shape of Ca2+ sparks is better explained by anisotropic diffusion of Ca2+ ions and indicator dye rather than by subsarcomeric inhomogeneities of the Ca2+ buffer and transport system. In addition, we examine the contribution of off-center confocal sampling to the variance of spark statistics.     

On the roles of Ca2+ diffusion, Ca2+ buffers, and the endoplasmic reticulum in IP3-induced Ca2+ waves.

We have investigated the effects of Ca2+ diffusion, mobile and stationary Ca2+ buffers in the cytosol, and Ca2+ handling by the endoplasmic reticulum on inositol 1,4,5-trisphosphate-induced Ca2+ wave propagation. Rapid equilibration of free and bound Ca2+ is used to describe Ca2+ sequestration by buffers in both the cytosol and endoplasmic reticulum (ER) lumen. Cytosolic Ca2+ regulation is based on a kinetic model of the inositol 1,4,5-trisphosphate (IP3) receptor of De Young and Keizer that includes activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane and SERCA Ca2+ pumps in the ER. Diffusion of Ca2+ in the cytosol and the ER and the breakdown and diffusion of IP3 are also included in our calculations. Although Ca2+ diffusion is severely limited because of buffering, when conditions are chosen just below the threshold for Ca2+ oscillations, a pulse of IP3 or Ca2+ results in a solitary trigger wave that requires diffusion of Ca2+ for its propagation. In the oscillatory regime repetitive wave trains are observed, but for this type of wave neither the wave shape nor the speed is strongly dependent on the diffusion of Ca2+. Local phase differences lead to waves that are predominately kinematic in nature, so that the wave speed (c) is related to the wavelength (lambda) and the period of the oscillations (tau) approximately by the formula c = lambda/tau. The period is determined by features that control the oscillations, including [IP3] and pump activity, which are related to recent experiments. Both solitary waves and wave trains are accompanied by a Ca2+ depletion wave in the ER lumen, similar to that observed in cortical preparations from sea urchin eggs. We explore the effect of endogenous and exogenous Ca2+ buffers on wave speed and wave shape, which can be explained in terms of three distinct effects of buffering, and show that exogenous buffers or Ca2+ dyes can have considerable influence on the amplitude and width of the waves.     

Minimal model for membrane oscillations in the pancreatic beta-cell.

Following the experimental findings of Atwater et al. (In Biochemistry Biophysics of the Pancreatic-beta-Cell, George Thieme Verlag, New York, 100-107), we have formulated a mathematical model for ionic and electrical events that take place in pancreatic-beta-cells. Our formulation incorporates a Hodgkin-Huxley type gating mechanism for Ca2+ and K+ channels, in addition to Ca2+ gated K+-channels. Consistent with the experimental observations, our model generates spikes and bursts in beta-cell membrane potentials and gives the correct responses to additions of glucose, quinine, and tetraethylammonium ions. The response of the oscillations to ouabain and changing concentrations of external K+ can be incorporated into the present model, although a more complete treatment would require inclusion of the Na+/K+ pump.     

ATP-sensitive potassium channel and bursting in the pancreatic beta cell. A theoretical study.

Based on the existence of ATP-sensitive potassium channels in the plasma membrane of pancreatic beta cells, we develop a quantitative explanation of the electrical activity observed in pancreatic islets. The proposed mechanism involves the voltage-dependent inward calcium and outward potassium currents described by Rorsman and Trube (1986), which are voltage-activated when an increase in the cytoplasmic ATP/ADP ratio decreases the conductance of the ATP-sensitive potassium channels. It is proposed that modulation of the ATP/ADP ratio occurs through calcium inhibition of oxidative phosphorylation. In this picture the mitochondria serve as a transducer of metabolic activity whose sensitivity is modulated by cytosolic calcium. Solution of the differential equations that describe this mechanism gives rise to both bursting and continuous spiking electrical activity similar to that observed experimentally. While the mechanism for bursting in this model involves the ATP/ADP ratio, the feedback is still provided by calcium, as originally proposed by Chay and Keizer (1983) using a Ca2+-activated potassium conductance. A mixed-model, which includes both ATP-sensitive and Ca2+-activated potassium conductances, also reproduces the experimentally observed electrical activity and may correspond more closely to the actual situation in vivo.     

First-generation transgenic plants and statistics

The statistical analyses of populations of first-generation transgenic plants are commonly based on mean and variance and generally require a test of normality. Since in many cases the assumptions of normality are not met, analyses can result in erroneous conclusions. Transformation of data to normality, the use of other distributions, or distribution-free statistical tests should then be used to obtain valid conclusions from these populations.