Optogenetic manipulation of cardiac electrical dynamics using sub-threshold illumination: dissecting the role of cardiac alternans in terminating rapid rhythms.

Cardiac action potential (AP) shape and propagation are regulated by several key dynamic factors such as ion channel recovery and intracellular Ca2+ cycling. Experimental methods for manipulating AP electrical dynamics commonly use ion channel inhibitors that lack spatial and temporal specificity. In this work, we propose an approach based on optogenetics to manipulate cardiac electrical activity employing a light-modulated depolarizing current with intensities that are too low to elicit APs (sub-threshold illumination), but are sufficient to fine-tune AP electrical dynamics. We investigated the effects of sub-threshold illumination in isolated cardiomyocytes and whole hearts by using transgenic mice constitutively expressing a light-gated ion channel (channelrhodopsin-2, ChR2). We find that ChR2-mediated depolarizing current prolongs APs and reduces conduction velocity (CV) in a space-selective and reversible manner. Sub-threshold manipulation also affects the dynamics of cardiac electrical activity, increasing the magnitude of cardiac alternans. We used an optical system that uses real-time feedback control to generate re-entrant circuits with user-defined cycle lengths to explore the role of cardiac alternans in spontaneous termination of ventricular tachycardias (VTs). We demonstrate that VT stability significantly decreases during sub-threshold illumination primarily due to an increase in the amplitude of electrical oscillations, which implies that cardiac alternans may be beneficial in the context of self-termination of VT.

Arrhythmia susceptibility in a rat model of acute atrial dilation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia, associated with an increased risk of stroke and heart failure. Acute AF occurs in response to sudden increases of atrial hemodynamic load, leading to atrial stretch. The mechanisms of stretch-induced AF were investigated in large mammals with controversial results. We optimized an approach to monitor rat atrial electrical activity using a red-shifted voltage sensitive dye (VSD). The methodology includes cauterization of the main ventricular coronary arteries, allowing improved atrial staining by the VSD and appropriate atrial perfusion for long experiments. Next, we developed a rat model of acute biatrial dilation (ABD) through the insertion of latex balloons into both atria, which could be inflated with controlled volumes. A chronic model of atrial dilation (spontaneous hypertensive rats; SHR) was used for comparison. ABD was performed on atria from healthy Wistar-Kyoto (WKY) rats (WKY-ABD). The atria were characterized in terms of arrhythmias susceptibility, action potential duration and conduction velocity. The occurrence of arrhythmias in WKY-ABD was significantly higher compared to non-dilated WKY atria. In WKY-ABD we found a reduction of conduction velocity, similar to that observed in SHR atria, while action potential duration was unchanged. Low-dose caffeine was used to introduce a drop of CV in WKY atria (WKY-caff), quantitatively similar to the one observed after ABD, but no increased arrhythmia susceptibility was observed with caffeine only. In conclusion, CV decrease is not sufficient to promote arrhythmias; enlargement of atrial surface is essential to create a substrate for acute reentry-based arrhythmias.

Real-time optical manipulation of cardiac conduction in intact hearts.

KEY POINTS: Although optogenetics has clearly demonstrated the feasibility of cardiac manipulation, current optical stimulation strategies lack the capability to react acutely to ongoing cardiac wave dynamics. Here, we developed an all-optical platform to monitor and control electrical activity in real-time. The methodology was applied to restore normal electrical activity after atrioventricular block and to manipulate the intraventricular propagation of the electrical wavefront. The closed-loop approach was also applied to simulate a re-entrant circuit across the ventricle. The development of this innovative optical methodology provides the first proof-of-concept that a real-time all-optical stimulation can control cardiac rhythm in normal and abnormal conditions. ABSTRACT: Optogenetics has provided new insights in cardiovascular research, leading to new methods for cardiac pacing, resynchronization therapy and cardioversion. Although these interventions have clearly demonstrated the feasibility of cardiac manipulation, current optical stimulation strategies do not take into account cardiac wave dynamics in real time. Here, we developed an all-optical platform complemented by integrated, newly developed software to monitor and control electrical activity in intact mouse hearts. The system combined a wide-field mesoscope with a digital projector for optogenetic activation. Cardiac functionality could be manipulated either in free-run mode with submillisecond temporal resolution or in a closed-loop fashion: a tailored hardware and software platform allowed real-time intervention capable of reacting within 2 ms. The methodology was applied to restore normal electrical activity after atrioventricular block, by triggering the ventricle in response to optically mapped atrial activity with appropriate timing. Real-time intraventricular manipulation of the propagating electrical wavefront was also demonstrated, opening the prospect for real-time resynchronization therapy and cardiac defibrillation. Furthermore, the closed-loop approach was applied to simulate a re-entrant circuit across the ventricle demonstrating the capability of our system to manipulate heart conduction with high versatility even in arrhythmogenic conditions. The development of this innovative optical methodology provides the first proof-of-concept that a real-time optically based stimulation can control cardiac rhythm in normal and abnormal conditions, promising a new approach for the investigation of the (patho)physiology of the heart.

Novel insights on the relationship between T-tubular defects and contractile dysfunction in a mouse model of hypertrophic cardiomyopathy.

Abnormalities of cardiomyocyte Ca(2+) homeostasis and excitation-contraction (E-C) coupling are early events in the pathogenesis of hypertrophic cardiomyopathy (HCM) and concomitant determinants of the diastolic dysfunction and arrhythmias typical of the disease. T-tubule remodelling has been reported to occur in HCM but little is known about its role in the E-C coupling alterations of HCM. Here, the role of T-tubule remodelling in the electro-mechanical dysfunction associated to HCM is investigated in the Δ160E cTnT mouse model that expresses a clinically-relevant HCM mutation. Contractile function of intact ventricular trabeculae is assessed in Δ160E mice and wild-type siblings. As compared with wild-type, Δ160E trabeculae show prolonged kinetics of force development and relaxation, blunted force-frequency response with reduced active tension at high stimulation frequency, and increased occurrence of spontaneous contractions. Consistently, prolonged Ca(2+) transient in terms of rise and duration are also observed in Δ160E trabeculae and isolated cardiomyocytes. Confocal imaging in cells isolated from Δ160E mice reveals significant, though modest, remodelling of T-tubular architecture. A two-photon random access microscope is employed to dissect the spatio-temporal relationship between T-tubular electrical activity and local Ca(2+) release in isolated cardiomyocytes. In Δ160E cardiomyocytes, a significant number of T-tubules (>20%) fails to propagate action potentials, with consequent delay of local Ca(2+) release. At variance with wild-type, we also observe significantly increased variability of local Ca(2+) transient rise as well as higher Ca(2+)-spark frequency. Although T-tubule structural remodelling in Δ160E myocytes is modest, T-tubule functional defects determine non-homogeneous Ca(2+) release and delayed myofilament activation that significantly contribute to mechanical dysfunction.

CaMKII activation and dynamics are independent of the holoenzyme structure: an infinite subunit holoenzyme approximation.

The combinatorial explosion produced by the multi-state, multi-subunit character of CaMKII has made analysis and modeling of this key signaling protein a significant challenge. Using rule-based and particle-based approaches, we construct exact models of CaMKII holoenzyme dynamics and study these models as a function of the number of subunits per holoenzyme, N. Without phosphatases the dynamics of activation are independent of the holoenzyme structure unless phosphorylation significantly alters the kinase activity of a subunit. With phosphatases the model is independent of holoenzyme size for N > 6. We introduce an infinite subunit holoenzyme approximation (ISHA), which simplifies the modeling by eliminating the combinatorial complexities encountered in any finite holoenzyme model. The ISHA is an excellent approximation to the full system over a broad range of physiologically relevant parameters. Finally, we demonstrate that the ISHA reproduces the behavior of exact models during synaptic plasticity protocols, which justifies its use as a module in large models of synaptic plasticity.

Virtual Cell modelling and simulation software environment.

The Virtual Cell (VCell; http://vcell.org/) is a problem solving environment, built on a central database, for analysis, modelling and simulation of cell biological processes. VCell integrates a growing range of molecular mechanisms, including reaction kinetics, diffusion, flow, membrane transport, lateral membrane diffusion and electrophysiology, and can associate these with geometries derived from experimental microscope images. It has been developed and deployed as a web-based, distributed, client-server system, with more than a thousand world-wide users. VCell provides a separation of layers (core technologies and abstractions) representing biological models, physical mechanisms, geometry, mathematical models and numerical methods. This separation clarifies the impact of modelling decisions, assumptions and approximations. The result is a physically consistent, mathematically rigorous, spatial modelling and simulation framework. Users create biological models and VCell will automatically (i) generate the appropriate mathematical encoding for running a simulation and (ii) generate and compile the appropriate computer code. Both deterministic and stochastic algorithms are supported for describing and running non-spatial simulations; a full partial differential equation solver using the finite volume numerical algorithm is available for reaction-diffusion-advection simulations in complex cell geometries including 3D geometries derived from microscope images. Using the VCell database, models and model components can be reused and updated, as well as privately shared among collaborating groups, or published. Exchange of models with other tools is possible via import/export of SBML, CellML and MatLab formats. Furthermore, curation of models is facilitated by external database binding mechanisms for unique identification of components and by standardised annotations compliant with the MIRIAM standard. VCell is now open source, with its native model encoding language (VCML) being a public specification, which stands as the basis for a new generation of more customised, experiment-centric modelling tools using a new plug-in based platform.

STAT module can function as a biphasic amplitude filter.

Signal transducer and actuator of transcription (STATs) are a family of transcription factors activated by various cytokines, growth factors and hormones. They are important mediators of immune responses and growth and differentiation of various cell types. The STAT signalling system represents a defined functional module with a pattern of signalling that is conserved from flies to mammals. In order to probe and gain insights into the signalling properties of the STAT module by computational means, we developed a simple non-linear ordinary differential equations model within the 'Virtual Cell' framework. Our results demonstrate that the STAT module can operate as a 'biphasic amplitude filter' with an ability to amplify input signals within a specific intermediate range. We show that dimerisation of phosphorylated STAT is crucial for signal amplification and the amplitude filtering function. We also demonstrate that maximal amplification at intermediate levels of STAT activation is a moderately robust property of STAT module. We propose that these observations can be extrapolated to the analogous SMAD signalling module.

Wavelength- and time-dependence of potentiometric non-linear optical signals from styryl dyes.

Second harmonic generation (SHG) imaging microscopy is an important emerging technique for biological research, complementing existing one- and two-photon fluorescence (2PF) methods. A non-linear phenomenon employing light from mode-locked Ti:sapphire or fiber-based lasers, SHG results in intrinsic optical sectioning without the need for a confocal aperture. Furthermore, as a second-order process SHG is confined to loci lacking a center of symmetry, a constraint that is readily satisfied by lipid membranes with only one leaflet stained by a dye. Of particular interest is "resonance-enhanced" SHG from styryl dyes in cellular membranes and the possibility that SHG is sensitive to transmembrane potential. We have previously confirmed this, using simultaneous voltage-clamping and non-linear imaging of cells to find that SHG is up to four times more sensitive to potential than fluorescence. In this work, we have extended these results in two directions. First, with a range of wavelengths available from a mode-locked Ti:sapphire laser and a fiber-based laser, we have more fully investigated SHG and 2PF voltage-sensitivity from ANEP and ASTAP chromophores, obtaining SHG sensitivity spectra that are consistent with resonance enhancements. Second, we have modified our system to coordinate the application of voltage-clamp steps with non-linear image acquisition to more precisely characterize the time dependence of SHG and 2PF voltage sensitivity, finding that, at least for some dyes, SHG responds more slowly than fluorescence to changes in transmembrane potential.

Initiation of sodium spikelets in basal dendrites of neocortical pyramidal neurons.

Cortical information processing relies critically on the processing of electrical signals in pyramidal neurons. Electrical transients mainly arise when excitatory synaptic inputs impinge upon distal dendritic regions. To study the dendritic aspect of synaptic integration one must record electrical signals in distal dendrites. Since thin dendritic branches, such as oblique and basal dendrites, do not support routine glass electrode measurements, we turned our effort towards voltage-sensitive dye recordings. Using the optical imaging approach we found and reported previously that basal dendrites of neocortical pyramidal neurons show an elaborate repertoire of electrical signals, including backpropagating action potentials and glutamate-evoked plateau potentials. Here we report a novel form of electrical signal, qualitatively and quantitatively different from backpropagating action potentials and dendritic plateau potentials. Strong glutamatergic stimulation of an individual basal dendrite is capable of triggering a fast spike, which precedes the dendritic plateau potential. The amplitude of the fast initial spikelet was actually smaller that the amplitude of the backpropagating action potential in the same dendritic segment. Therefore, the fast initial spike was dubbed "spikelet". Both the basal spikelet and plateau potential propagate decrementally towards the cell body, where they are reflected in the somatic whole-cell recordings. The low incidence of basal spikelets in the somatic intracellular recordings and the impact of basal spikelets on soma-axon action potential initiation are discussed.

Novel naphthylstyryl-pyridium potentiometric dyes offer advantages for neural network analysis.

The submucous plexus of the guinea pig intestine is a quasi-two-dimensional mammalian neural network that is particularly amenable to study using multiple site optical recording of transmembrane voltage (MSORTV) [Biol. Bull. 183 (1992) 344; J. Neurosci. 19 (1999) 3073]. For several years the potentiometric dye of choice for monitoring the electrical activity of its individual neurons has been di-8-ANEPPS [Neuron 9 (1992) 393], a naphthylstyryl-pyridinium dye with a propylsulfonate headgroup that provides relatively large fluorescence changes during action potentials and synaptic potentials. Limitations to the use of this dye, however, have been its phototoxicity and its low water solubility which requires the presence of DMSO and Pluronic F-127 in the staining solution. In searching for less toxic and more soluble dyes exhibiting larger fluorescence signals, we first tried the dienylstyryl-pyridinium dye RH795 [J. Neurosci. 14 (1994) 2545] which is highly soluble in water. This dye yielded relatively large signals, but it was internalized quickly by the submucosal neurons resulting in rapid degradation of the signal-to-noise ratio. We decided to synthesize a series of naphthylstyryl-pyridinium dyes (di-n-ANEPPDHQ) having the same chromophore as di-8-ANEPPS and the quaternary ammonium headgroup (DHQ) of RH795 (resulting in two positive charges versus the neutral propylsulfonate-ring nitrogen combination), and we tested the di-methyl (JPW3039), di-ethyl (JPW2081), di-propyl (JPW3031), di-butyl (JPW5029), and di-octyl (JPW5037) analogues, all of them soluble in ethanol. We found that the di-propyl (di-3-ANEPPDHQ) and the di-butyl (di-4-ANEPPDHQ) forms yielded the best combination of signal-to-noise ratio, moderate phototoxicity and absence of dye internalization.

The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models.

MOTIVATION: Molecular biotechnology now makes it possible to build elaborate systems models, but the systems biology community needs information standards if models are to be shared, evaluated and developed cooperatively. RESULTS: We summarize the Systems Biology Markup Language (SBML) Level 1, a free, open, XML-based format for representing biochemical reaction networks. SBML is a software-independent language for describing models common to research in many areas of computational biology, including cell signaling pathways, metabolic pathways, gene regulation, and others. AVAILABILITY: The specification of SBML Level 1 is freely available from http://www.sbml.org/

The Virtual Cell: a software environment for computational cell biology.

The newly emerging field of computational cell biology requires software tools that address the needs of a broad community of scientists. Cell biological processes are controlled by an interacting set of biochemical and electrophysiological events that are distributed within complex cellular structures. Computational modeling is familiar to researchers in fields such as molecular structure, neurobiology and metabolic pathway engineering, and is rapidly emerging in the area of gene expression. Although some of these established modeling approaches can be adapted to address problems of interest to cell biologists, relatively few software development efforts have been directed at the field as a whole. The Virtual Cell is a computational environment designed for cell biologists as well as for mathematical biologists and bioengineers. It serves to aid the construction of cell biological models and the generation of simulations from them. The system enables the formulation of both compartmental and spatial models, the latter with either idealized or experimentally derived geometries of one, two or three dimensions.

Second-harmonic imaging microscopy of living cells.

Second harmonic generation (SHG) has been developed in our laboratories as a high-resolution nonlinear optical imaging microscopy for cellular membranes and intact tissues. SHG shares many of the advantageous features for microscopy of another more established nonlinear optical technique: two-photon excited fluorescence (TPEF). Both are capable of optical sectioning to produce three-dimensional images of thick specimens and both result in less photodamage to living tissue than confocal microscopy. SHG is complementary to TPEF in that it uses a different contrast mechanism and is most easily detected in the transmitted light optical path. It can be used to image membrane probes with high membrane specificity and displays extraordinary sensitivity in reporting membrane potential; it also has the ability to image highly ordered structural proteins without any exogenous labels.

Computational cell biologists snowed in at Cranwell.

Cell biology is being inundated by an avalanche of data from the genomics and proteomics enterprises. The complexity and sheer volume of information threaten to overwhelm the ability of traditional cell biologists to grasp its implications and develop experimentally testable hypotheses. For this reason, some have begun to explore computational approaches towards organizing complex data into quantitative models. This requires communication and collaboration between the biological science community and and the physical and mathematical sciences communities. A recent meeting [The First International Symposium on Computational Cell Biology, Cranwell Resort, Lenox, MA, USA; 4-6 March 2001. Organizers: J.H. Carson, A. Cowan, and L.M. Loew (www.nrcam.uchc.edu/conference).] made a first attempt to bring these two communities together. Three feet of new snow fell during the meeting, but the 125 attendees, an unusual mixture of cell biologists, computer scientists, mathematicians, physicists, and engineers, were having too much fun defining the new field of computational cell biology to notice that they were literally snowed in.

Topology of the mitochondrial inner membrane: dynamics and bioenergetic implications.

Electron tomography indicates that the mitochondrial inner membrane is not normally comprised of baffle-like folds as depicted in textbooks. In actuality, this membrane is pleomorphic, with narrow tubular regions connecting the internal compartments (cristae) to each other and to the membrane periphery. The membrane topologies observed in condensed (matrix contracted) and orthodox (matrix expanded) mitochondria cannot be interconverted by passive folding and unfolding. Instead, transitions between these morphological states likely involve membrane fusion and fission. Formation of tubular junctions in the inner membrane appears to be energetically favored, because they form spontaneously in yeast mitochondria following large-amplitude swelling and recontraction. However, aberrant, unattached, vesicular cristae are also observed in these mitochondria, suggesting that formation of cristae junctions depends on factors (such as the distribution of key proteins and/or lipids) that are disrupted during extreme swelling. Computer modeling studies using the "Virtual Cell" program suggest that the shape of the inner membrane can influence mitochondrial function. Simulations indicate that narrow cristae junctions restrict diffusion between intracristal and external compartments, causing depletion of ADP and decreased ATP output inside the cristae.

Physiological modeling with virtual cell framework.

This article describes a computational framework for cell biological modeling and simulation that is based on the mapping of experimental biochemical and electrophysiological data onto experimental images. The framework is designed to enable the construction of complex general models that encompass the general class of problems coupling reaction and diffusion.

An image-based model of calcium waves in differentiated neuroblastoma cells.

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.

High-resolution nonlinear optical imaging of live cells by second harmonic generation.

By adapting a laser scanning microscope with a titanium sapphire femtosecond pulsed laser and transmission optics, we are able to produce live cell images based on the nonlinear optical phenomenon of second harmonic generation (SHG). Second harmonic imaging (SHIM) is an ideal method for probing membranes of living cells because it offers the high resolution of nonlinear optical microscopy with the potential for near-total avoidance of photobleaching and phototoxicity. The technique has been implemented on three cell lines labeled with membrane-staining dyes that have large nonlinear optical coefficients. The images can be obtained within physiologically relevant time scales. Both achiral and chiral dyes were used to compare image formation for the case of single- and double-leaflet staining, and it was found that chirality plays a significant role in the mechanism of contrast generation. It is also shown that SHIM is highly sensitive to membrane potential, with a depolarization of 25 mV resulting in an approximately twofold loss of signal intensity.

Morphological control of inositol-1,4,5-trisphosphate-dependent signals.

Inositol-1,4,5-trisphosphate (InsP(3))-mediated calcium signals represent an important mechanism for transmitting external stimuli to the cell. However, information about intracellular spatial patterns of InsP(3) itself is not generally available. In particular, it has not been determined how the interplay of InsP(3) generation, diffusion, and degradation within complex cellular geometries can control the patterns of InsP(3) signaling. Here, we explore the spatial and temporal characteristics of [InsP(3)](cyt) during a bradykinin-induced calcium wave in a neuroblastoma cell. This is achieved by using a unique image-based computer modeling system, Virtual Cell, to integrate experimental data on the rates and spatial distributions of the key molecular components of the process. We conclude that the characteristic calcium dynamics requires rapid, high-amplitude production of [InsP(3)](cyt) in the neurite. This requisite InsP(3) spatiotemporal profile is provided, in turn, as an intrinsic consequence of the cell's morphology, demonstrating how geometry can locally and dramatically intensify cytosolic signals that originate at the plasma membrane. In addition, the model predicts, and experiments confirm, that stimulation of just the neurite, but not the soma or growth cone, is sufficient to generate a calcium response throughout the cell.

Determination of time-dependent inositol-1,4,5-trisphosphate concentrations during calcium release in a smooth muscle cell.

The level of [InsP3]cyt required for calcium release in A7r5 cells, a smooth muscle cell line, was determined by a new set of procedures using quantitative confocal microscopy to measure release of InsP3 from cells microinjected with caged InsP3. From these experiments, the [InsP3]cyt required to evoke a half-maximal calcium response is 100 nM. Experiments with caged glycerophosphoryl-myo-inositol 4, 5-bisphosphate (GPIP2), a slowly metabolized analogue of InsP3, gave a much slower recovery and a half-maximal response of an order of magnitude greater than InsP3. Experimental data and highly constrained variables were used to construct a mathematical model of the InsP3-dependent [Ca2+]cyt changes; the resulting simulations show high fidelity to experiment. Among the elements considered in constructing this model were the mechanism of the InsP3-receptor, InsP3 degradation, calcium buffering in the cytosol, and refilling of the ER stores via sarcoplasmic endoplasmic reticulum ATPase (SERCA) pumps. The model predicts a time constant of 0.8 s for InsP3 degradation and 13 s for GPIP2. InsP3 degradation was found to be a prerequisite for [Ca2+]cyt recovery to baseline levels and is therefore critical to the pattern of the overall [Ca2+]cyt signal. Analysis of the features of this model provides insights into the individual factors controlling the amplitude and shape of the InsP3-mediated calcium signal.