Integrated luminal and cytosolic aspects of the calcium release control

Luminal Ca(2+) dynamics during IP3R mediated signals

The role of cytosolic Ca(2+) on the kinetics of Inositol 1,4,5-triphosphate receptors (IP3Rs) and on the dynamics of IP3R-mediated Ca(2+) signals has been studied at large both experimentally and by modeling. The role of luminal Ca(2+) has not been investigated with that much detail although it has been found that it is relevant for signal termination in the case of Ca(2+) release through ryanodine receptors. In this work we present the results of observing the dynamics of luminal and cytosolic Ca(2+) simultaneously in Xenopus laevis oocytes. Combining observations and modeling we conclude that there is a rapid mechanism that guarantees the availability of free Ca(2+) in the lumen even when a relatively large Ca(2+) release is evoked. Comparing the dynamics of cytosolic and luminal Ca(2+) during a release, we estimate that they are consistent with a 80% of luminal Ca(2+) being buffered. The rapid availability of free luminal Ca(2+) correlates with the observation that the lumen occupies a considerable volume in several regions across the images.

An allosteric model of the inositol trisphosphate receptor with nonequilibrium binding

The inositol trisphosphate receptor (IPR) is a crucial ion channel that regulates the Ca(2+) influx from the endoplasmic reticulum (ER) to the cytoplasm. A thorough study of the IPR channel contributes to a better understanding of calcium oscillations and waves. It has long been observed that the IPR channel is a typical biological system which performs adaptation. However, recent advances on the physical essence of adaptation show that adaptation systems with a negative feedback mechanism, such as the IPR channel, must break detailed balance and always operate out of equilibrium with energy dissipation. Almost all previous IPR models are equilibrium models assuming detailed balance and thus violate the dissipative nature of adaptation. In this article, we constructed a nonequilibrium allosteric model of single IPR channels based on the patch-clamp experimental data obtained from the IPR in the outer membranes of isolated nuclei of the Xenopus oocyte. It turns out that our model reproduces the patch-clamp experimental data reasonably well and produces both the correct steady-state and dynamic properties of the channel. Particularly, our model successfully describes the complicated bimodal [Ca(2+)] dependence of the mean open duration at high [IP3], a steady-state behavior which fails to be correctly described in previous IPR models. Finally, we used the patch-clamp experimental data to validate that the IPR channel indeed breaks detailed balance and thus is a nonequilibrium system which consumes energy.

Inositol trisphosphate receptor and ion channel models based on single-channel data

The inositol trisphosphate receptor (IPR) plays an important role in controlling the dynamics of intracellular Ca(2+). Single-channel patch-clamp recordings are a typical way to study these receptors as well as other ion channels. Methods for analyzing and using this type of data have been developed to fit Markov models of the receptor. The usual method of parameter fitting is based on maximum-likelihood techniques. However, Bayesian inference and Markov chain Monte Carlo techniques are becoming more popular. We describe the application of the Bayesian methods to real experimental single-channel data in three ion channels: the ryanodine receptor, the K(+) channel, and the IPR. One of the main aims of all three studies was that of model selection with different approaches taken. We also discuss the modeling implications for single-channel data that display different levels of channel activity within one recording.

A model-based method for estimating Ca2+ release fluxes from linescan images in Xenopus oocytes

We propose a model-based method of interpreting linescan images observed in Xenopus oocytes with the use of Oregon Green-1 as a fluorescent dye. We use a detailed modeling formalism based on numerical simulations that incorporate physical barriers for local diffusion, and, by assuming a Gaussian distribution of release durations, we derive the distributions of release Ca(2+) amounts and currents, fluorescence amplitudes, and puff widths. We analyze a wide set of available data collected from 857 and 281 events observed in the animal and the vegetal hemispheres of the oocyte, respectively. A relatively small fraction of events appear to involve coupling of two or three adjacent clusters of Ca(2+) releasing channels. In the animal hemisphere, the distribution of release currents with a mean of 1.4 pA presents a maximum at 1.0 pA and a rather long tail extending up to 5 pA. The overall distribution of liberated Ca(2+) amounts exhibits a dominant peak at 120 fC, a smaller peak at 375 fC, and an average of 166 fC. Ca(2+) amounts and release fluxes in the vegetal hemisphere appear to be 3.6 and 1.6 times smaller than in the animal hemisphere, respectively. Predicted diameters of elemental release sites are approximately 1.0 microm in the animal and approximately 0.5 microm in the vegetal hemisphere, but the side-to-side separation between adjacent sites appears to be identical (approximately 0.4 microm). By fitting the model to individual puffs we can estimate the quantity of liberated calcium, the release current, the orientation of the scan line, and the dimension of the corresponding release site.

An investigation of models of the IP3R channel in Xenopus oocyte

We consider different models of inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R) channels in order to fit nuclear membrane patch clamp data of the stationary open probability, mean open time, and mean close time of channels in the Xenopus oocyte. Our results indicate that rather than to treat the tetrameric IP(3)R as four independent and identical subunits, one should assume sequential binding-unbinding processes of Ca(2+) ions and IP(3) messengers. Our simulations also favor the assumption that a channel opens through a conformational transition from a close state to an active state.

Modulation of calcium signals by fluorescent dyes in the presence of tubular endoplasmic reticulum: a modelling approach

The complex network of Ca(2+) signals uses local events as building blocks for generating global calcium signals with different shapes. However, the nature of the large time- and space-scales of local calcium signals observed in Xenopus oocytes has remained unclear. By numeric simulations that include optical blurring of the image and the geometrical restrictions imposed by tubules of the endoplasmic reticulum or other cell structures, we investigate how the fluorescent dye affect the observed features of calcium events, such as rate of signal decay, spatial size, fluorescence amplitude, or the apparent diffusion like from a point source in a spherically symmetric space. We add more evidence that, irrespective of the dye properties, local calcium signals produced in the presence of tubular cellular structures are consistently wider than expected in a homogeneous environment. Moreover, the spatial dimension and the decay time of the event increase with the quantity of liberated Ca(2+). Our results also indicate that a fast binding Ca(2+) indicator that does not bind to cytosolic proteins yields fast signals when the event is observed in the front of the release site, and slow signals when the event is viewed from the opposite side of the tubule. We propose several ways to test our model by various experimental procedures.

A kinetic model of single and clustered IP3 receptors in the absence of Ca2+ feedback

Ca2+ liberation through inositol 1,4,5-trisphosphate receptor (IP3R) channels generates complex patterns of spatiotemporal cellular Ca2+ signals owing to the biphasic modulation of channel gating by Ca2+ itself. These processes have been extensively studied in Xenopus oocytes, where imaging studies have revealed local Ca2+ signals ("puffs") arising from clusters of IP3R, and patch-clamp studies on isolated oocyte nuclei have yielded extensive data on IP3R gating kinetics. To bridge these two levels of experimental data, we developed an IP3R model and applied stochastic simulation and transition matrix theory to predict the behavior of individual and clustered IP3R channels. The channel model consists of four identical, independent subunits, each of which has an IP3-binding site together with one activating and one inactivating Ca2+-binding site. The channel opens when at least three subunits undergo a conformational change to an "active" state after binding IP3 and Ca2+. The model successfully reproduces patch-clamp data; including the dependence of open probability, mean open duration, and mean closed duration on [IP3] and [Ca2+]. Notably, the biexponential distribution of open-time duration and the dependence of mean open time on [Ca2+] are explained by populations of openings involving either three or four active subunits. As a first step toward applying the single IP3R model to describe cellular responses, we then simulated measurements of puff latency after step increases of [IP3]. Assuming that stochastic opening of a single IP3R at basal cytosolic [Ca2+] and any given [IP3] has a high probability of rapidly triggering neighboring channels by calcium-induced calcium release to evoke a puff, optimal correspondence with experimental data of puff latencies after photorelease of IP3 was obtained when the cluster contained a total of 40-70 IP3Rs.

Hybrid stochastic and deterministic simulations of calcium blips

Intracellular calcium release is a prime example for the role of stochastic effects in cellular systems. Recent models consist of deterministic reaction-diffusion equations coupled to stochastic transitions of calcium channels. The resulting dynamics is of multiple time and spatial scales, which complicates far-reaching computer simulations. In this article, we introduce a novel hybrid scheme that is especially tailored to accurately trace events with essential stochastic variations, while deterministic concentration variables are efficiently and accurately traced at the same time. We use finite elements to efficiently resolve the extreme spatial gradients of concentration variables close to a channel. We describe the algorithmic approach and we demonstrate its efficiency compared to conventional methods. Our single-channel model matches experimental data and results in intriguing dynamics if calcium is used as charge carrier. Random openings of the channel accumulate in bursts of calcium blips that may be central for the understanding of cellular calcium dynamics.

Characterization of local calcium signals in tubular networks of endoplasmic reticulum

To explain the large time and space scales of elementary calcium events in Xenopus oocytes it is assumed that the Ca2+ source is located on tubules of the endoplasmic reticulum, which provide local barriers for diffusion. The event duration, width and signal mass dependence on the total quantity of released Ca2+ is determined at different orientations of the scan line and different ionic currents. Excellent agreement with published data is obtained with on- and off-rate constants of the fluorescent indicator of 15 microM(-1) s(-1) and 2.55 s(-1), respectively. It is found that one signal mass unit, calculated with the classical method that assumes spherical symmetry of the cytosolic space surrounding the release site, corresponds to 0.189 fC of released Ca2+ in the presence of a tubular network. It is estimated that release Ca2+ currents and amounts are randomly distributed, with averages of 0.165 pA and 3.66 fC per event and average release duration of 22.2 ms. The total quantity of liberated Ca2+ and the release current amplitude in the presence of endoplasmic reticulum tubules is predicted to be about one order of magnitude lower than estimated within the isotropic diffusion formalism. This could have implications in muscle cell Ca2+ imaging as well.

Gating mechanisms of the type-1 inositol trisphosphate receptor

A large amount of data and observations on inositol 1,4,5-trisphosphate (IP(3)) binding to the IP(3) receptor/Ca(2+) channel, the steady-state activity of the channel, and its inactivation by IP(3) can be explained by assuming one activation and one inhibition module, both allosterically operated by Ca(2+), IP(3), and ATP, and one adaptation element, driven by IP(3), Ca(2+), and the interconversion between two possible conformations of the receptor. The adaptation module becomes completely insensitive to a second IP(3) pulse within 80 s. Observed kinetic responses are well reproduced if, in addition, two module open states are rendered inactive by the current charge carrier Mn(2+). The inactivation time constants are 59 s in the activation, and 0.75 s in the adaptation module. The in vivo open probability of the channel is predicted to be almost in coincidence with the behavior in lipid bilayers for IP(3) levels of 0.2 and 2 microM and one-order-higher at 0.02 microM IP(3), whereas at 180 microM IP(3) the maximal in vivo activity may be 2.5-orders higher than in bilayers and restricted to a narrower Ca(2+) domain (approximately 10 microM-wide versus approximately 100 microM-wide). IP(3) is likely to inhibit channel activity at < or =120 nM Ca(2+) in vivo.

Wolframin expression induces novel ion channel activity in endoplasmic reticulum membranes and increases intracellular calcium

Wolfram syndrome is an autosomal recessive neuro-degenerative disorder associated with juvenile onset non-autoimmune diabetes mellitus and progressive optic atrophy. The disease has been attributed to mutations in the WFS1 gene, which codes for a protein predicted to possess 9-10 transmembrane segments. Little is known concerning the function of the WFS1 protein (wolframin). Endoglycosidase H digestion, immunocytochemistry, and subcellular fractionation studies all indicated that wolframin is localized to the endoplasmic reticulum in rat brain hippocampus and rat pancreatic islet beta-cells, and after ectopic expression in Xenopus oocytes. Reconstitution of wolframin from oocyte membranes into planar lipid bilayers demonstrated that the protein induced a large cation-selective ion channel that was blocked by Mg2+ or Ca2+. Inositol triphosphate was capable of activating channels in the fused bilayers that were similar to channel components induced by wolframin expression. Expression of wolframin also increased cytosolic calcium levels in oocytes. Wolframin thus appears to be important in the regulation of intracellular Ca2+ homeostasis. Disruption of this function may place cells at risk to suffer inappropriate death decisions, thus accounting for the progressive beta-cell loss and neuronal degeneration associated with the disease.