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.

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.

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.

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.

Electric field-directed fibroblast locomotion involves cell surface molecular reorganization and is calcium independent.

Directional cellular locomotion is thought to involve localized intracellular calcium changes and the lateral transport of cell surface molecules. We have examined the roles of both calcium and cell surface glycoprotein redistribution in the directional migration of two murine fibroblastic cell lines, NIH 3T3 and SV101. These cell types exhibit persistent, cathode directed motility when exposed to direct current electric fields. Using time lapse phase contrast microscopy and image analysis, we have determined that electric field-directed locomotion in each cell type is a calcium independent process. Both exhibit cathode directed motility in the absence of extracellular calcium, and electric fields cause no detectable elevations or gradients of cytosolic free calcium. We find evidence suggesting that galvanotaxis in these cells involves the lateral redistribution of plasma membrane glycoproteins. Electric fields cause the lateral migration of plasma membrane concanavalin A receptors toward the cathode in both NIH 3T3 and SV101 fibroblasts. Exposure of directionally migrating cells to Con A inhibits the normal change of cell direction following a reversal of electric field polarity. Additionally, when cells are plated on Con A-coated substrata so that Con A receptors mediate cell-substratum adhesion, cathode-directed locomotion and a cathodal accumulation of Con A receptors are observed. Immunofluorescent labeling of the fibronectin receptor in NIH 3T3 fibroblasts suggests the recruitment of integrins from large clusters to form a more diffuse distribution toward the cathode in field-treated cells. Our results indicate that the mechanism of electric field directed locomotion in NIH 3T3 and SV101 fibroblasts involves the lateral redistribution of plasma membrane glycoproteins involved in cell-substratum adhesion.

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.

Localized membrane depolarizations and localized calcium influx during electric field-guided neurite growth.

Our study explores the mechanisms behind neurite galvanotropism. Using phase, differential interference contrast and ratiometric fluorescence microscopy, we reveal four responses of N1E-115 mouse neuroblastoma cells to 0.1-1.0 mV/microns uniform DC electric fields: cathode-directed neurite initiation and elongation, cathode-biased growth cone filopodial protrusions, transient cathode-localized calcium increases, and persistent cathode-localized membrane depolarizations. These newly demonstrated events are temporally and spatially correlated, suggesting that they are causally related. The calcium increases are prevented by calcium channel blockers and by the removal of extracellular calcium. We therefore propose that the observed field-induced membrane depolarizations activate voltage-dependent calcium channels, resulting in cathode-localized calcium influx. This, in turn, may initiate the observed cathode-biased growth cone filopodial protrusions, followed by the cathode-directed neurite elongation.

Distinct electric potentials in soma and neurite membranes.

Structurally similar voltage-dependent ion channels may behave differently in different locations along the surface of a neuron. A possible reason could be that channels experience nonuniform electrical potentials along the plasmalemma. Here, we map the electrical potentials along the membrane of differentiated N1E-115 neuroblastoma cells with a potential-sensitive dye. We find that the intramembrane potential gradient is indeed more positive in the membranes of neurites than in the membranes of somata. This is not attributable to differences in ion conductances or surface charge densities between the membranes of neurites and somata; instead, it can be explained by differences in lipid composition. The spatial variation in intramembrane electrical potential may help account for electrophysiological and functional differences between neurites and somata.

Insertion of amphiphilic molecules into membranes is catalyzed by a high molecular weight non-ionic surfactant.

We have employed an amphiphilic fluorescent probe to elucidate the mechanism by which a class of oxyethylene-oxypropylene copolymers catalyzes the insertion of hydrophobic or amphiphilic molecules into membranes. The rate of binding can be accelerated by over two orders of magnitude in the presence of the catalyst which does not itself disrupt the lipid bilayer. The rate of probe binding to lipid vesicles does not depend on the lipid concentration in the presence or absence of catalyst but is linearly related to the concentration of the catalyst. Probe binding to the polyol surfactant appears to be a component of the catalytic mechanism and equilibrium binding parameters can be determined; these are used to indirectly establish quantitative binding parameters for the probe to the vesicle membrane. The polyol surfactant is also shown to catalyze insertion of the probe into the outer leaflet of a hemispherical lipid bilayer and the plasma membrane of HeLa cells. The latter were also stained by catalyzed transfer of a fluorescent lipid from lipid vesicles. The permeability of the cell membrane is not significantly altered under any of the catalytic conditions. These data, taken together, suggest that the polyol surfactant extracts a monomeric substrate molecule from its aggregate or microcrystal and passes it to the membrane via a loose and transient contact.

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.

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.

Isolation, characterization and partial purification of a transferable membrane channel (amoebapore) produced by Entamoeba histolytica.

Entamoeba histolytica kills cells by contact dependent cytolysis. The mechanism underlying this process must be of rapid onset because target cells round up and show marked zeiosis within 15 min of contact. In earlier work, we identified a remarkable ion-channel forming protein which we named amoebapore, that may contribute to the amoeba-induced target cell killing. Within the amoeba it exists as part of a supramolecular aggregate together with other proteins of unknown function. In this work we report the purification of a solubilized form of the amoebapore. Amoebapore was found to exist as an apparent dimer of the previously reported protein whose molecular weight had been determined under denaturing conditions. Two isoforms of this dimer, with pI values of 6.8 and 5.3 present at a ratio of 7 to 1, were identified and purified. Both isoforms demonstrate ion-channel forming activity in planar lipid membranes. These channels show a unit conductance of 5-20 pS and remain open for less than 1 s. Upon lateral aggregation, opening becomes concerted to a greater degree with channel conductance are observed. The isolated particulate form of amoebapore depolarizes cells.

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.

Dye screening and signal-to-noise ratio for retrogradely transported voltage-sensitive dyes.

Using a novel method for retrogradely labeling specific neuronal populations, we tested different styryl dyes in attempt to find dyes whose staining would be specific, rapid, and lead to large activity dependent signals. The dyes were injected into the ventral roots of the isolated chick spinal cord from embryos at days E9-E12. The voltage-sensitive dye signals were recorded from synaptically activated motoneurons using a 464 element photodiode array. The best labeling and optical signals were obtained using the relatively hydrophobic dyes di-8-ANEPPQ and di-12-ANEPEQ. Over the 24 h period we examined, these dyes bound specifically to the cells with axons in the ventral roots. The dyes responded with an increase in fluorescence of 1-3% (delta F/F) in response to synaptic depolarization of the motoneurons. The signal-to-noise ratio obtained in a single trial from a detector that received light from a 14 x 14 microns2 area of the motoneuron population was about 10:1. Nonetheless, signals on neighboring diodes were similar, suggesting that we were not detecting the activity of individual neurons. Retrograde labeling and optical recording with voltage-sensitive dyes provides a means for monitoring the activity of identified neurons in situations where microelectrode recordings are not feasible.

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.

Design and characterization of electrochromic membrane probes.

Electrochronic membrane probes display a spectroscopic response to membrane potential by a direct electronic mechanism. This allows such probes to be designed a priori via quantum-chemical techniques. The detailed behavior of potentiometric optical probes can be elucidated with an apparatus based on phase-sensitive detection from a hemispherical lipid bilayer; several different types of response spectra can be obtained with this apparatus allowing distinction between the electrochromic mechanism and the more common molecular-motion based mechanisms. The development of 'fast' potentiometric dyes has now reached a stage where practical and exciting applications are rapidly appearing. It is anticipated that the emergence of a complementary set of electrochromic probes will lead to new applications; in particular, it may be possible to elucidate the molecular events which underlie biological or physiological phenomena.

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.