Mepivacaine reduces calcium transients in isolated murine ventricular cardiomyocytes.

BACKGROUND: The potential mechanism of mepivacaine's myocardial depressant effect observed in papillary muscle has not yet been investigated at cellular level. Therefore, we evaluated mepivacaine's effects on Ca2+ transient in isolated adult mouse cardiomyocytes. METHODS: Single ventricular myocytes were enzymatically isolated from wild-type C57Bl/6 mice and loaded with 10 µM fluorescent Ca2+ indicator Fluo-4-AM to record intracellular Ca2+ transients upon electrical stimulation. The mepivacaine effects at half-maximal inhibitory concentration (IC50) was determined on calibrated cardiomyocytes' Ca2+ transients by non-parametric statistical analyses on biophysical parameters. Combination of mepivacaine with NCX blockers ORM-10103 or NiCl2 were used to test a possible mechanism to explain mepivacaine-induced Ca2+ transients' reduction. RESULTS: A significant inhibition at mepivacaine's IC50 (50 µM) on Ca2+ transients was measured in biophysical parameters such as peak (control: 528.6 ± 73.61 nM vs mepivacaine: 130.9 ± 15.63 nM; p < 0.05), peak area (control: 401.7 ± 63.09 nM*s vs mepivacaine: 72.14 ± 10.46 nM*s; p < 0.05), slope (control: 7699 ± 1110 nM/s vs mepivacaine: 1686 ± 226.6 nM/s; p < 0.05), time to peak (control: 107.9 ± 8.967 ms vs mepivacaine: 83.61 ± 7.650 ms; p < 0.05) and D50 (control: 457.1 ± 47.16 ms vs mepivacaine: 284.5 ± 22.71 ms; p < 0.05). Combination of mepivacaine with NCX blockers ORM-10103 or NiCl2 showed a significant increase in the baseline of [Ca2+] and arrhythmic activity upon electrical stimulation. CONCLUSION: At cellular level, mepivacaine blocks Na+ channels, enhancing the reverse mode activity of NCX, leading to a significant reduction of Ca2+ transients. These results suggest a new mechanism for the mepivacaine-reduction contractility effect.

Simulation Strategies for Calcium Microdomains and Calcium Noise.

In this article, we present an overview of simulation strategies in the context of subcellular domains where calcium-dependent signaling plays an important role. The presentation follows the spatial and temporal scales involved and represented by each algorithm. As an exemplary cell type, we will mainly cite work done on striated muscle cells, i.e. skeletal and cardiac muscle. For these cells, a wealth of ultrastructural, biophysical and electrophysiological data is at hand. Moreover, these cells also express ubiquitous signaling pathways as they are found in many other cell types and thus, the generalization of the methods and results presented here is straightforward.The models considered comprise the basic calcium signaling machinery as found in most excitable cell types including Ca2+ ions, diffusible and stationary buffer systems, and calcium regulated calcium release channels. Simulation strategies can be differentiated in stochastic and deterministic algorithms. Historically, deterministic approaches based on the macroscopic reaction rate equations were the first models considered. As experimental methods elucidated highly localized Ca2+ signaling events occurring in femtoliter volumes, stochastic methods were increasingly considered. However, detailed simulations of single molecule trajectories are rarely performed as the computational cost implied is too large. On the mesoscopic level, Gillespie's algorithm is extensively used in the systems biology community and with increasing frequency also in models of microdomain calcium signaling. To increase computational speed, fast approximations were derived from Gillespie's exact algorithm, most notably the chemical Langevin equation and the τ-leap algorithm. Finally, in order to integrate deterministic and stochastic effects in multiscale simulations, hybrid algorithms are increasingly used. These include stochastic models of ion channels combined with deterministic descriptions of the calcium buffering and diffusion system on the one hand, and algorithms that switch between deterministic and stochastic simulation steps in a context-dependent manner on the other. The basic assumptions of the listed methods as well as implementation schemes are given in the text. We conclude with a perspective on possible future developments of the field.

S100A1 DNA-based Inotropic Therapy Protects Against Proarrhythmogenic Ryanodine Receptor 2 Dysfunction.

Restoring expression levels of the EF-hand calcium (Ca(2+)) sensor protein S100A1 has emerged as a key factor in reconstituting normal Ca(2+) handling in failing myocardium. Improved sarcoplasmic reticulum (SR) function with enhanced Ca(2+) resequestration appears critical for S100A1's cyclic adenosine monophosphate-independent inotropic effects but raises concerns about potential diastolic SR Ca(2+) leakage that might trigger fatal arrhythmias. This study shows for the first time a diminished interaction between S100A1 and ryanodine receptors (RyR2s) in experimental HF. Restoring this link in failing cardiomyocytes, engineered heart tissue and mouse hearts, respectively, by means of adenoviral and adeno-associated viral S100A1 cDNA delivery normalizes diastolic RyR2 function and protects against Ca(2+)- and ß-adrenergic receptor-triggered proarrhythmogenic SR Ca(2+) leakage in vitro and in vivo. S100A1 inhibits diastolic SR Ca(2+) leakage despite aberrant RyR2 phosphorylation via protein kinase A and calmodulin-dependent kinase II and stoichiometry with accessory modulators such as calmodulin, FKBP12.6 or sorcin. Our findings demonstrate that S100A1 is a regulator of diastolic RyR2 activity and beneficially modulates diastolic RyR2 dysfunction. S100A1 interaction with the RyR2 is sufficient to protect against basal and catecholamine-triggered arrhythmic SR Ca(2+) leak in HF, combining antiarrhythmic potency with chronic inotropic actions.

Colocalization properties of elementary Ca(2+) release signals with structures specific to the contractile filaments and the tubular system of intact mouse skeletal muscle fibers.

Ca(2+) regulates several important intracellular processes. We combined second harmonic generation (SHG) and two photon excited fluorescence microscopy (2PFM) to simultaneously record the SHG signal of the myosin filaments and localized elementary Ca(2+) release signals (LCSs). We found LCSs associated with Y-shaped structures of the myosin filament pattern (YMs), so called verniers, in intact mouse skeletal muscle fibers under hypertonic treatment. Ion channels crucial for the Ca(2+) regulation are located in the tubular system, a system that is important for Ca(2+) regulation and excitation-contraction coupling. We investigated the tubular system of intact, living mouse skeletal muscle fibers using 2PFM and the fluorescent Ca(2+) indicator Fluo-4 dissolved in the external solution or the membrane dye di-8-ANEPPS. We simultaneously measured the SHG signal from the myosin filaments of the skeletal muscle fibers. We found that at least a subset of the YMs observed in SHG images are closely juxtaposed with Y-shaped structures of the transverse tubules (YTs). The distances of corresponding YMs and YTs yield values between 1.3 µm and 4.1 µm including pixel uncertainty with a mean distance of 2.52±0.10 µm (S.E.M., n=41). Additionally, we observed that some of the linear-shaped areas in the tubular system are colocalized with linear-shaped areas in the SHG images.

MD simulations of the central pore of ryanodine receptors and sequence comparison with 2B protein from coxsackie virus.

The regulation of intracellular Ca(2+) triggers a multitude of vital processes in biological cells. Ca(2+) permeable ryanodine receptors (RyRs) are the biggest known ion channels and play a key role in the regulation of intracellular calcium concentrations, particularly in muscle cells. In this study, we construct a computational model of the pore region of the skeletal RyR and perform molecular dynamics (MD) simulations. The dynamics and distribution of Ca(2+) around the luminal pore entry of the RyR suggest that Ca(2+) ions are channeled to the pore entry due to the arrangement of (acidic) amino acids at the extramembrane surface of the protein. This efficient mechanism of Ca(2+) supply is thought to be part of the mechanism of Ca(2+) conductance of RyRs. Viral myocarditis is predominantly caused by coxsackie viruses that induce the expression of the protein 2B which is known to affect intracellular Ca(2+) homeostasis in infected cells. From our sequence comparison, it is hypothesized, that modulation of RyR could be due to replacement of its transmembrane domains (TMDs) by those domains of the viral channel forming protein 2B of coxsackie virus. This article is part of a Special Issue entitled: Viral Membrane Proteins - Channels for Cellular Networking.

Enabling the detection of UV signal in multimodal nonlinear microscopy with catalogue lens components.

Using an optical system made from fused silica catalogue optical components, third-order nonlinear microscopy has been enabled on conventional Ti:sapphire laser-based multiphoton microscopy setups. The optical system is designed using two lens groups with straightforward adaptation to other microscope stands when one of the lens groups is exchanged. Within the theoretical design, the optical system collects and transmits light with wavelengths between the near ultraviolet and the near infrared from an object field of at least 1 mm in diameter within a resulting numerical aperture of up to 0.56. The numerical aperture can be controlled with a variable aperture stop between the two lens groups of the condenser. We demonstrate this new detection capability in third harmonic generation imaging experiments at the harmonic wavelength of ∼300 nm and in multimodal nonlinear optical imaging experiments using third-order sum frequency generation and coherent anti-Stokes Raman scattering microscopy so that the wavelengths of the detected signals range from ∼300 nm to ∼660 nm.

Cardiac and respiratory dysfunction in Duchenne muscular dystrophy and the role of second messengers.

Duchenne muscular dystrophy (DMD) affects young boys and is characterized by the absence of dystrophin, a large cytoskeletal protein present in skeletal and cardiac muscle cells and neurons. The heart and diaphragm become necrotic in DMD patients and animal models of DMD, resulting in cardiorespiratory failure as the leading cause of death. The major consequences of the absence of dystrophin are high levels of intracellular Ca(2+) and the unbalanced production of NO that can finally trigger protein degradation and cell death. Cytoplasmic increase in Ca(2+) concentration directly and indirectly triggers different processes such as necrosis, fibrosis, and activation of macrophages. The absence of the neuronal isoform of nitric oxide synthase (nNOS) and the overproduction of NO by the inducible isoform (iNOS) further increase the intracellular Ca(2+) via a hypernitrosylation of the ryanodine receptor. NO overproduction, which further induces the expression of iNOS but decreases the expression of the endothelial isoform (eNOS), deregulates the muscle tissue blood flow creating an ischemic situation. The high levels of Ca(2+) in dystrophic muscles and the ischemic state of the muscle tissue would culminate in a positive feedback loop. While efforts continue toward optimizing cardiac and respiratory care of DMD patients, both Ca(2+) and NO in cardiac and respiratory muscle pathways have been shown to be important to the etiology of the disease. Understanding the mechanisms behind the fine regulation of Ca(2+) -NO may be important for a noninterventional and noninvasive supportive approach to treat DMD patients, improving the quality of life and natural history of DMD patients.

Molecular dynamics simulations and conductance studies of the interaction of VP1 N-terminus from Polio virus and gp41 fusion peptide from HIV-1 with lipid membranes.

The icosahedral Polio virus capsid consists of 60 copies of each of the coat proteins VP1, VP2, VP3 and myristolyated VP4 (myrVP4). Catalyzed by the host cell receptor the Polio virus enters the host cell via externalization of myrVP4 and the N terminal part of VP1. There are several assumptions about the individual role of both of the proteins in the mechanism of membrane attachment and genome injection. We use the first 32 N terminal amino acids of VP1 and applied molecular dynamics simulations to assess its mechanism of function when attached and inserted into hydrated lipid membranes (POPC). Helical models are placed in various positions in regard to the lipid membrane to start with. As a comparison, the first 33 amino acids of the fusion peptide of gp41 of HIV-1 are simulated under identical conditions. Computational data support the idea that VP1 is not penetrating into the membrane to form a pore; it rather lays on the membrane surface and only perturbs the membrane. Furthermore, this idea is strengthened by channel recordings of both peptides showing irregular openings.

Simulation strategies for calcium microdomains and calcium-regulated calcium channels.

In this article, we present an overview of simulation strategies in the context of subcellular domains where calcium-dependent signaling plays an important role. The presentation follows the spatial and temporal scales involved and represented by each algorithm. As an exemplary cell type, we will mainly cite work done on striated muscle cells, i.e. skeletal and cardiac muscle. For these cells, a wealth of ultrastructural, biophysical and electrophysiological data is at hand. Moreover, these cells also express ubiquitous signaling pathways as they are found in many other cell types and thus, the generalization of the methods and results presented here is straightforward.The models considered comprise the basic calcium signaling machinery as found in most excitable cell types including Ca(2+) ions, diffusible and stationary buffer systems, and calcium regulated calcium release channels. Simulation strategies can be differentiated in stochastic and deterministic algorithms. Historically, deterministic approaches based on the macroscopic reaction rate equations were the first models considered. As experimental methods elucidated highly localized Ca(2+) signaling events occurring in femtoliter volumes, stochastic methods were increasingly considered. However, detailed simulations of single molecule trajectories are rarely performed as the computational cost implied is too large. On the mesoscopic level, Gillespie's algorithm is extensively used in the systems biology community and with increasing frequency also in models of microdomain calcium signaling. To increase computational speed, fast approximations were derived from Gillespie's exact algorithm, most notably the chemical Langevin equation and the τ-leap algorithm. Finally, in order to integrate deterministic and stochastic effects in multiscale simulations, hybrid algorithms are increasingly used. These include stochastic models of ion channels combined with deterministic descriptions of the calcium buffering and diffusion system on the one hand, and algorithms that switch between deterministic and stochastic simulation steps in a context-dependent manner on the other. The basic assumptions of the listed methods as well as implementation schemes are given in the text. We conclude with a perspective on possible future developments of the field.

Exact and approximate stochastic simulation of intracellular calcium dynamics.

In simulations of chemical systems, the main task is to find an exact or approximate solution of the chemical master equation (CME) that satisfies certain constraints with respect to computation time and accuracy. While Brownian motion simulations of single molecules are often too time consuming to represent the mesoscopic level, the classical Gillespie algorithm is a stochastically exact algorithm that provides satisfying results in the representation of calcium microdomains. Gillespie's algorithm can be approximated via the tau-leap method and the chemical Langevin equation (CLE). Both methods lead to a substantial acceleration in computation time and a relatively small decrease in accuracy. Elimination of the noise terms leads to the classical, deterministic reaction rate equations (RRE). For complex multiscale systems, hybrid simulations are increasingly proposed to combine the advantages of stochastic and deterministic algorithms. An often used exemplary cell type in this context are striated muscle cells (e.g., cardiac and skeletal muscle cells). The properties of these cells are well described and they express many common calcium-dependent signaling pathways. The purpose of the present paper is to provide an overview of the aforementioned simulation approaches and their mutual relationships in the spectrum ranging from stochastic to deterministic algorithms.

Second harmonic generation microscopy probes different states of motor protein interaction in myofibrils.

The second harmonic generation (SHG) signal intensity sourced from skeletal muscle myosin II strongly depends on the polarization of the incident laser beam relative to the muscle fiber axis. This dependence is related to the second-order susceptibility χ((2)), which can be described by a single component ratio γ under generally assumed symmetries. We precisely extracted γ from SHG polarization dependence curves with an extended focal field model. In murine myofibrillar preparations, we have found two distinct polarization dependencies: With the actomyosin system in the rigor state, γ(rig) has a mean value of γ(rig) = 0.52 (SD = 0.04, n = 55); in a relaxed state where myosin is not bound to actin, γ(rel) has a mean value of γ(rel) = 0.24 (SD = 0.07, n = 70). We observed a similar value in an activated state where the myosin power stroke was pharmacologically inhibited using N-benzyl-p-toluene sulfonamide. In summary, different actomyosin states can be visualized noninvasively with SHG microscopy. Specifically, SHG even allows us to distinguish different actin-bound states of myosin II using γ as a parameter.

Modulating the activity of the channel-forming segment of Vpr protein from HIV-1.

Viral protein of regulation (Vpr) encoded by human immunodeficiency virus type 1 (HIV-1) is a short auxiliary protein that is 96 amino acids in length. During the viral life cycle, Vpr is released into the blood serum and is able to enter cellular membranes of noninfected cells. In this study a short peptide, Vpr(55-83), was shown to exhibit ion-channel-like activity when reconstituted into (1) planar lipid bilayers and (2) lipid bilayers held at the tip of a glass pipette. The two set-ups led to differences in the oligomerization state of the peptide, which was reflected in differences in the conductance levels. Experiments under applied hydrostatic pressure affect the dynamics of the protein within the membrane.

Determination of the Maximum Velocity of Filaments in the in vitro Motility Assay.

The in vitro motility assay (IVMA) is a powerful tool commonly used in basic muscle research and for drug screenings with the aim to find treatment options for neuromuscular disorders. In brief, the sliding movement of fluorescence-labeled actin filaments on myosin motor proteins is recorded, and the sliding velocity is analyzed via image analysis methods. Due to low signal-to-noise ratios and large variability in the velocity signal, accurate determination of the maximum sliding velocity is challenging. We introduce a new method and software program named Actin Phase Velocity (ActiPHV). The method extracts the maximum velocity from filament tracking data. Based on simulated and real reference data we show that our method yields a higher accuracy compared to previous methods. As a result, our method enables enhancing the sensitivity of the IVMA to better exploit its full potential.

Augmentation of myocardial If dysregulates calcium homeostasis and causes adverse cardiac remodeling.

HCN channels underlie the depolarizing funny current (If) that contributes importantly to cardiac pacemaking. If is upregulated in failing and infarcted hearts, but its implication in disease mechanisms remained unresolved. We generated transgenic mice (HCN4tg/wt) to assess functional consequences of HCN4 overexpression-mediated If increase in cardiomyocytes to levels observed in human heart failure. HCN4tg/wt animals exhibit a dilated cardiomyopathy phenotype with increased cellular arrhythmogenicity but unchanged heart rate and conduction parameters. If augmentation induces a diastolic Na+ influx shifting the Na+/Ca2+ exchanger equilibrium towards 'reverse mode' leading to increased [Ca2+]i. Changed Ca2+ homeostasis results in significantly higher systolic [Ca2+]i transients and stimulates apoptosis. Pharmacological inhibition of If prevents the rise of [Ca2+]i and protects from ventricular remodeling. Here we report that augmented myocardial If alters intracellular Ca2+ homeostasis leading to structural cardiac changes and increased arrhythmogenicity. Inhibition of myocardial If per se may constitute a therapeutic mechanism to prevent cardiomyopathy.

Alterations in intracellular Ca2+-homeostasis of skeletal muscle fibers during sepsis.

OBJECTIVE: To investigate changes in intracellular Ca2+-regulation and Ca2+-sensitivity of the contractile apparatus in murine skeletal muscle fibers during sepsis. DESIGN AND SETTING: Animal study in a university-based research laboratory. SUBJECTS: Isolated muscle fibers (M. extensor digitorum longus) of septic mice. INTERVENTIONS: In one group, sepsis was induced in "black six" mice using cecal ligation and puncture (CLP). In a second group, laparotomy (SHAM), and in a third group, general anesthesia (GA) was performed. Saponin-skinned skeletal muscle fibers were examined 2, 3, 5, and 7 days after treatment, and caffeine-induced Ca2+-release from the sarcoplasmic reticulum (SR) as well as Ca2+-sensitivity of the contractile apparatus were assessed. MEASUREMENTS AND RESULTS: In the CLP group, Ca2+-release significantly decreased over 5 days and increased again after 7 days. In the SHAM group, Ca2+-release decreased at days 2 and 3, whereas no changes were observed in the GA group. Ca2+-sensitivity significantly increased over 5 days in the CLP group and decreased again at day 7. In the SHAM group, Ca2+-sensitivity increased at days 2 and 3, and no changes were seen in the GA group. CONCLUSIONS: In murine skeletal muscle fibers, Ca2+-release from the SR decreases during sepsis, with effects being most pronounced 2-3 days after CLP. In parallel, Ca2+-sensitivity of the contractile apparatus is increased, and all changes are reversible. Thus, these effects might be involved in skeletal muscle dysfunction during sepsis as corresponding changes are less pronounced or absent in control groups.

Velocity distributions of single F-actin trajectories from a fluorescence image series using trajectory reconstruction and optical flow mapping.

We present an approach for the computation of single-object velocity statistics in a noisy fluorescence image series. The algorithm is applied to molecular imaging data from an in vitro actin-myosin motility assay. We compare the relative efficiency of wavelet and curvelet transform denoising in terms of noise reduction and object restoration. It is shown that while both algorithms reduce background noise efficiently, curvelet denoising restores the curved edges of actin filaments more reliably. Noncrossing spatiotemporal actin trajectories are unambiguously identified using a novel segmentation scheme that locally combines the information of 2-D and 3-D segmentation. Finally, the optical flow vector field for the image sequence is computed via the 3-D structure tensor and mapped to the segmented trajectories. Using single-trajectory statistics, the global velocity distribution extracted from an image sequence is decomposed into the contributions of individual trajectories. The technique is further used to analyze the distribution of the x and y components of the velocity vectors separately, and it is shown that directed actin motion is found in myosin extracts from single skeletal muscle fibers. The presented approach may prove helpful to identify actin filament subpopulations and to analyze actin-myosin interaction kinetics under biochemical regulation.

In silico assessment of the conduction mechanism of the Ryanodine Receptor 1 reveals previously unknown exit pathways.

The ryanodine receptor 1 is a large calcium ion channel found in mammalian skeletal muscle. The ion channel gained a lot of attention recently, after multiple independent authors published near-atomic cryo electron microscopy data. Taking advantage of the unprecedented quality of structural data, we performed molecular dynamics simulations on the entire ion channel as well as on a reduced model. We calculated potentials of mean force for Ba2+, Ca2+, Mg2+, K+, Na+ and Cl- ions using umbrella sampling to identify the key residues involved in ion permeation. We found two main binding sites for the cations, whereas the channel is strongly repulsive for chloride ions. Furthermore, the data is consistent with the model that the receptor achieves its ion selectivity by over-affinity for divalent cations in a calcium-block-like fashion. We reproduced the experimental conductance for potassium ions in permeation simulations with applied voltage. The analysis of the permeation paths shows that ions exit the pore via multiple pathways, which we suggest to be related to the experimental observation of different subconducting states.

S100A1 decreases calcium spark frequency and alters their spatial characteristics in permeabilized adult ventricular cardiomyocytes.

S100A1, a Ca2+-sensor protein of the EF-hand type, exerts positive inotropic effects in the heart via enhanced cardiac ryanodine receptor (RyR2) activity. Here we report that S100A1 protein (0.1microM) interacts with the RyR2 in resting permeabilized cardiomyocytes at free Ca2+-levels comparable to diastolic Ca2+-concentrations ( approximately 150nM). Alterations of RyR2 function due to S100A1 binding was assessed via analysis of Ca2+-spark characteristics. Ca2+-spark frequency, amplitude and duration were all reduced upon perfusion with 0.1microM S100A1 protein by 38%, 14% and 18%, respectively. Most likely, these effects were conveyed through the S100A1 C-terminus (S100A1-ct; amino acids 75-94) as the corresponding S100A1-ct peptide (0.1microM) inhibited S100A1 protein binding to the RyR2 and similarly attenuated frequency, amplitude and duration of Ca2+-sparks by 52%, 8% and 26%, respectively. Accordingly, the sarcoplasmic reticulum (SR) Ca2+-content was slightly increased but the stoichiometry of other accessory RyR2 modulators (sorcin/FKBP12.6) remained unaltered by S100A1. Hence, we propose S100A1 as a novel inhibitory modulator of RyR2 function at diastolic Ca2+-concentrations in rabbit ventricular cardiomyocytes.

The neuroactive steroids alphaxalone and pregnanolone increase the conductance of single GABAA channels in newborn rat hippocampal neurons.

The effects of the neuroactive steroids alphaxalone and pregnanolone on single GABA(A) receptor channels were tested in cell-attached and inside-out patches from cultured newborn rat hippocampal neurons. The conductance of these single channels ranged between 10 and 80 pS when exposed to low (0.5-3 microM) GABA concentrations. These GABA concentrations activated low-conducting channels (<40 pS) in 78% of the patches, 22% of patches had channels with a maximum conductance above 40 pS. Alphaxalone at concentrations above 1 microM, and pregnanolone at concentrations above 0.1 microM, significantly increased the conductance of initially low-conducting single channels activated by GABA up to seven-fold and at all concentrations tested, both drugs increased open probability and mean open time and decreased closed probability and mean closed time of channels. Both steroids at higher concentrations could directly activate high conductance (>40 pS) chloride channels. Both the directly activated channels and those channels that had been previously affected by alphaxalone were modulated by diazepam, a benzodiazepine drug that is known to specifically modulate GABA(A) channels. The present study is the first one to show that neurosteroids can significantly increase single GABA(A) channel conductance, thus enlarging our current knowledge on the molecular mechanism of action of these compounds.

Fast XYT imaging of elementary calcium release events in muscle with multifocal multiphoton microscopy and wavelet denoising and detection.

We used multifocal multiphoton microscopy to image fast, localized elevations of the cytosolic Ca2+ concentration in two spatial dimensions plus time (XYT). This technique extends the common spatially 1-D XT imaging and allows the acquisition of more than ten times longer time series (>500 images) and ten times larger areas of interest than for previously used confocal XYT imaging techniques due to lower phototoxicity and fast multifocal scanning. We recorded spontaneously occurring elementary Ca2+ release events in chemically permeabilized adult mammalian skeletal muscle fibers using two-photon excitation of the fluorescent dye Fluo-4. The resulting time series were analyzed with an automated denoising and detection algorithm based on the à trous implementation of the discrete wavelet transform. Wavelet coefficient hard-thresholding is used for denoising and event detection is performed across several wavelet scales. The spatiotemporal characteristics of the detected Ca2+ release events are followed throughout the XYT stack and are parametrized using a biophysically valid anisotropic Gaussian event model. The proposed method allows a detailed spatiotemporal analysis of elementary Ca2+ release events underlying the excitation-contraction coupling process in muscle.