Mode switching of Inositol 1,4,5-trisphosphate receptor channel shapes the Spatiotemporal scales of Ca2+ signals

TraceSpecks: A Software for Automated Idealization of Noisy Patch-Clamp and Imaging Data

Experimental records of single molecules or ion channels from fluorescence microscopy and patch-clamp electrophysiology often include high-frequency noise and baseline fluctuations that are not generated by the system under investigation and have to be removed. Moreover, multiple channels or conductance levels can be present at a time in the data that need to be quantified to accurately understand the behavior of the system. Manual procedures for removing these fluctuations and extracting conducting states or multiple channels are laborious, prone to subjective bias, and likely to hinder the processing of often very large data sets. We introduce a maximal likelihood formalism for separating signal from a noisy and drifting background such as fluorescence traces from imaging of elementary Ca²⁺ release events called puffs arising from clusters of channels, and patch-clamp recordings of ion channels. Parameters such as the number of open channels or conducting states, noise level, and background signal can all be optimized using the expectation-maximization algorithm. We implement our algorithm following the Baum-Welch approach to expectation-maximization in the portable Java language with a user-friendly graphical interface and test the algorithm on both synthetic and experimental data from the patch-clamp electrophysiology of Ca²⁺ channels and fluorescence microscopy of a cluster of Ca²⁺ channels and Ca²⁺ channels with multiple conductance levels. The resulting software is accurate, fast, and provides detailed information usually not available through manual analysis. Options for visual inspection of the raw and processed data with key parameters are provided, in addition to a range of statistics such as the mean open probabilities, mean open times, mean close times, dwell-time distributions for different number of channels open or conductance levels, amplitude distribution of all opening events, and number of transitions between different number of open channels or conducting levels in asci format with a single click.

Data-driven modeling of mitochondrial dysfunction in Alzheimer's disease

Intracellular accumulation of oligomeric forms of β amyloid (Aβ) are now believed to play a key role in the earliest phase of Alzheimer's disease (AD) as their rise correlates well with the early symptoms of the disease. Extensive evidence points to impaired neuronal Ca²⁺ homeostasis as a direct consequence of the intracellular Aβ oligomers. However, little is known about the downstream effects of the resulting Ca²⁺ rise on the many intracellular Ca²⁺-dependent pathways. Here we use multiscale modeling in conjunction with patch-clamp electrophysiology of single inositol 1,4,5-trisphosphate (IP₃) receptor (IP₃R) and fluorescence imaging of whole-cell Ca²⁺ response, induced by exogenously applied intracellular Aβ₄₂ oligomers to show that Aβ₄₂ inflicts cytotoxicity by impairing mitochondrial function. Driven by patch-clamp experiments, we first model the kinetics of IP₃R, which is then extended to build a model for the whole-cell Ca²⁺ signals. The whole-cell model is then fitted to fluorescence signals to quantify the overall Ca²⁺ release from the endoplasmic reticulum by intracellular Aβ₄₂ oligomers through G-protein-mediated stimulation of IP₃ production. The estimated IP₃ concentration as a function of intracellular Aβ₄₂ content together with the whole-cell model allows us to show that Aβ₄₂ oligomers impair mitochondrial function through pathological Ca²⁺ uptake and the resulting reduced mitochondrial inner membrane potential, leading to an overall lower ATP and increased production of reactive oxygen species and H₂O₂. We further show that mitochondrial function can be restored by the addition of Ca²⁺ buffer EGTA, in accordance with the observed abrogation of Aβ₄₂ cytotoxicity by EGTA in our live cells experiments.

High expression of type I inositol 1,4,5-trisphosphate receptor in the kidney of rats with hepatorenal syndrome

AIM

To detect the expression of type I inositol 1,4,5-trisphosphate receptor (IP3RI) in the kidney of rats with hepatorenal syndrome (HRS).

METHODS

One hundred and twenty-five Sprague-Dawley rats were randomly divided into four groups to receive an intravenous injection of D-galactosamine (D-GalN) plus lipopolysaccharide (LPS; group G/L, n = 50), D-GalN alone (group G, n = 25), LPS alone (group L, n = 25), and normal saline (group NS, n = 25), respectively. At 3, 6, 9, 12, and 24 h after injection, blood, liver, and kidney samples were collected. Hematoxylin-eosin staining of liver tissue was performed to assess hepatocyte necrosis. Electron microscopy was used to observe ultrastructural changes in the kidney. Western blot analysis and real-time PCR were performed to detect the expression of IP3RI protein and mRNA in the kidney, respectively.

RESULTS

Hepatocyte necrosis was aggravated gradually, which was most significant at 12 h after treatment with D-galactosamine/lipopolysaccharide, and was characterized by massive hepatocyte necrosis. At the same time, serum levels of biochemical indicators including liver and kidney function indexes were all significantly changed. The structure of the renal glomerulus and tubules was normal at all time points. Western blot analysis indicated that IP3RI protein expression began to rise at 3 h (P < 0.05) and peaked at 12 h (P < 0.01). Real-time PCR demonstrated that IP3RI mRNA expression began to rise at 3 h (P < 0.05) and peaked at 9 h (P < 0.01).

CONCLUSION

IP3RI protein expression is increased in the kidney of HRS rats, and may be regulated at the transcriptional level.

Ionotropic and Metabotropic Mechanisms of Allosteric Modulation of α7 Nicotinic Receptor Intracellular Calcium

The pharmacological targeting of the α7 nicotinic acetylcholine receptor (α7) is a promising strategy in the development of new drugs for neurologic diseases. Because α7 receptors regulate cellular calcium, we investigated how the prototypical type II-positive allosteric modulator PNU120596 affects α7-mediated calcium signaling. Live imaging experiments show that PNU120596 augments ryanodine receptor-driven calcium-induced calcium release (CICR), inositol-induced calcium release (IICR), and phospholipase C activation by the α7 receptor. Both influx of calcium through the α7 nicotinic acetylcholine receptor (nAChR) channel as well as the binding of intracellular G proteins were involved in the effect of PNU120596 on intracellular calcium. This is evidenced by the findings that chelation of extracellular calcium, expression of α7D44A or α7345-348A mutant subunits, or blockade of calcium store release compromised the ability of PNU120596 to increase intracellular calcium transients generated by α7 ligand activation. Spatiotemporal stochastic modeling of calcium transient responses corroborates these results and indicates that α7 receptor activation enables calcium microdomains locally and to lesser extent in the distant cytosol. From the model, allosteric modulation of the receptor activates CICR locally via ryanodine receptors and augments IICR through enhanced calcium influx due to prolonged α7 nAChR opening. These findings provide a new mechanistic framework for understanding the effect of α7 receptor allosteric modulation on both local and global calcium dynamics.

Analyzing optical imaging of Ca2+ signals via TIRF microscopy: The limits on resolution due to chemical rates and depth of the channels

High resolution total internal reflection (TIRF) microscopy (TIRFM) together with detailed computational modeling provides a powerful approach towards the understanding of a wide range of Ca²⁺ signals mediated by the ubiquitous inositol 1,4,5-trisphosphate (IP₃) receptor (IP₃R) channel. Exploiting this fruitful collaboration further requires close agreement between the models and observations. However, elementary Ca²⁺ release events, puffs, imaged through TIRFM do not show the rapid single-channel openings and closings during and between puffs as are present in simulated puffs using data-driven single channel models. TIRFM also shows a rapid equilibration of 10ms after a channel opens or closes which is not achievable in simulation using standard Ca²⁺ diffusion coefficients and reaction rates between indicator dye and Ca²⁺. Furthermore, TIRFM imaging cannot decipher the depth of the channel with respect to the microscope, which will affect the change in fluorescence that the microscope detects, thereby affecting its sensitivity to fast single-channel activity. Using the widely used Ca²⁺ diffusion coefficients and reaction rates, our simulations show equilibration rates that are eight times slower than TIRFM imaging. We show that to get equilibrium rates consistent with observed values, the diffusion coefficients and reaction rates have to be significantly higher than the values reported in the literature, and predict the channel depth to be 200-250nm. Finally, we show that with the addition of noise, short events due to 1-2ms opening and closing of channels that are observed in computational models can be missed in TIRFM.

Mapping Interpuff Interval Distribution to the Properties of Inositol Trisphosphate Receptors

Tightly clustered inositol trisphosphate receptors (IP₃Rs) control localized Ca²⁺ liberation from the endoplasmic reticulum to generate repetitive Ca²⁺ puffs. Distributions of the interpuff interval (IPI), i.e., the waiting time between successive puffs, are found to be well characterized by a probability density function involving only two parameters, λ and ξ, which represent the basal rate of puff generation and the recovery rate from refractoriness, respectively. However, how the two parameters depend on the kinetic parameters of single IP₃Rs in a cluster is still unclear. In this article, using a stochastic puff model and a single-channel data-based IP₃R model, we establish the dependencies of λ and ξ on two important IP₃R model parameters, IP₃ concentration ([IP₃]) and the recovery rate from Ca²⁺ inhibition (r_low). By varying [IP₃] and r_low in physiologically plausible ranges, we find that the ξ-λ plane is comprised of only two disjoint regions, a biologically impermissible region and a region where each parameter set (ξ, λ) can be caused by using two different combinations of [IP₃] and r_low. The two combinations utilize very different mechanisms to maintain the same IPI distribution, and the mechanistic difference provides a way of identifying IP₃R kinetic parameters by observing properties of the IPI.

Modelling modal gating of ion channels with hierarchical Markov models

Many ion channels spontaneously switch between different levels of activity. Although this behaviour known as modal gating has been observed for a long time it is currently not well understood. Despite the fact that appropriately representing activity changes is essential for accurately capturing time course data from ion channels, systematic approaches for modelling modal gating are currently not available. In this paper, we develop a modular approach for building such a model in an iterative process. First, stochastic switching between modes and stochastic opening and closing within modes are represented in separate aggregated Markov models. Second, the continuous-time hierarchical Markov model, a new modelling framework proposed here, then enables us to combine these components so that in the integrated model both mode switching as well as the kinetics within modes are appropriately represented. A mathematical analysis reveals that the behaviour of the hierarchical Markov model naturally depends on the properties of its components. We also demonstrate how a hierarchical Markov model can be parametrized using experimental data and show that it provides a better representation than a previous model of the same dataset. Because evidence is increasing that modal gating reflects underlying molecular properties of the channel protein, it is likely that biophysical processes are better captured by our new approach than in earlier models.