John C. Williams

A comprehensive multiscale framework for simulating optogenetics in the heart

Optogenetics has emerged as an alternative method for electrical control of the heart, where illumination is used to elicit a bioelectric response in tissue modified to express photosensitive proteins (opsins). This technology promises to enable evocation of spatiotemporally precise responses in targeted cells or tissues, thus creating new possibilities for safe and effective therapeutic approaches to ameliorate cardiac function. Here, we present a comprehensive framework for multi-scale modelling of cardiac optogenetics, allowing both mechanistic examination of optical control and exploration of potential therapeutic applications. The framework incorporates accurate representations of opsin channel kinetics and delivery modes, spatial distribution of photosensitive cells, and tissue illumination constraints, making possible the prediction of emergent behaviour resulting from interactions at sub-organ scales. We apply this framework to explore how optogenetic delivery characteristics determine energy requirements for optical stimulation and to identify cardiac structures that are potential pacemaking targets with low optical excitation threshold.

Computational Optogenetics: Empirically-Derived Voltage- and Light-Sensitive Channelrhodopsin-2 Model

Channelrhodospin-2 (ChR2), a light-sensitive ion channel, and its variants have emerged as new excitatory optogenetic tools not only in neuroscience, but also in other areas, including cardiac electrophysiology. An accurate quantitative model of ChR2 is necessary for in silico prediction of the response to optical stimulation in realistic tissue/organ settings. Such a model can guide the rational design of new ion channel functionality tailored to different cell types/tissues. Focusing on one of the most widely used ChR2 mutants (H134R) with enhanced current, we collected a comprehensive experimental data set of the response of this ion channel to different irradiances and voltages, and used these data to develop a model of ChR2 with empirically-derived voltage- and irradiance- dependence, where parameters were fine-tuned via simulated annealing optimization. This ChR2 model offers: 1) accurate inward rectification in the current-voltage response across irradiances; 2) empirically-derived voltage- and light-dependent kinetics (activation, deactivation and recovery from inactivation); and 3) accurate amplitude and morphology of the response across voltage and irradiance settings. Temperature-scaling factors (Q10) were derived and model kinetics was adjusted to physiological temperatures. Using optical action potential clamp, we experimentally validated model-predicted ChR2 behavior in guinea pig ventricular myocytes. The model was then incorporated in a variety of cardiac myocytes, including human ventricular, atrial and Purkinje cell models. We demonstrate the ability of ChR2 to trigger action potentials in human cardiomyocytes at relatively low light levels, as well as the differential response of these cells to light, with the Purkinje cells being most easily excitable and ventricular cells requiring the highest irradiance at all pulse durations. This new experimentally-validated ChR2 model will facilitate virtual experimentation in neural and cardiac optogenetics at the cell and organ level and provide guidance for the development of in vivo tools.

Effect of Ventricular Extrasystoles on Closure of Mitral Valve

Mitral valve function was assessed by roentgen videodensitometry. Mitral reflux was rare when single ventricular ectopic systoles were produced by electronic stimulation of the right or left ventricle at various times in the cardiac cycle. It was also rare during the compensatory pause after the ectopic systole or with the following postectopic systole. Recurrent ventricular ectopic systoles interposed once per cycle were associated with minor reflux when introduced in mid-cycle. Such extrasystoles occurred late enough in the cardiac cycle for the ventricle to relax after the primary systole and for the mitral valve to open before the extrasystole. The extrasystolic contractions were, however, weak and incapable of opening the aortic valve. When interposed early in the cardiac cycle, extrasystolic potentiation of the primary ventricular contraction occurred, and no or minimal mitral reflux was observed.

Channelrhodopsin2 Current During the Action Potential: “Optical AP Clamp” and Approximation

The most widely used optogenetic tool, Channelrhodopsin2 (ChR2), is both light- and voltage-sensitive. A light-triggered action potential or light-driven perturbations of ongoing electrical activity provide instant voltage feedback, shaping ChR2 current. Therefore, depending on the cell type and the light pulse duration, the typically reported voltage-clamp-measured ChR2 current traces are often not a good surrogate for the ChR2 current during optically-triggered action potentials. We discuss two experimental methods to reveal ChR2 current during an action potential: an “optical AP clamp” and its approximation employing measured current-voltage curve for ChR2. The methods are applicable to voltage- and light-sensitive ion currents operating in excitable cells, e.g. cardiomyocytes or neurons.

Relations between cortical thickness, serotonin 1A receptor binding, and structural connectivity: A multimodal imaging study.

Serotonin 1A (5‐HT1A) receptors play a direct role in neuronal development, cell proliferation, and dendritic branching. We hypothesized that variability in 5‐HT1A binding can affect cortical thickness, and may account for a subtype of major depressive disorder (MDD) in which both are altered. To evaluate this, we measured cortical thickness from structural magnetic resonance imaging (MRI) and 5‐HT1A binding by positron emission tomography (PET) in an exploratory study. To examine a range of 5‐HT1A binding and cortical thickness values, we recruited 25 healthy controls and 19 patients with MDD. We hypothesized increased 5‐HT1A binding in the raphe nucleus (RN) would be negatively associated with cortical thickness due to reduced serotonergic transmission. Contrary to our hypothesis, raphe 5‐HT1A binding was positively correlated with cortical thickness in right posterior cingulate cortex (PCC), a region implicated in the default mode network. Cortical thickness was also positively correlated with 5‐HT1A in each cortical region. We further hypothesized that the strength of 5‐HT1A‐cortical thickness correlation depends on the number of axons between the raphe nucleus and each region. To explore this we related 5‐HT1A–cortical thickness correlation coefficients to the number of tracts connecting that region and the raphe, as measured by diffusion tensor imaging (DTI) in an independent sample. The 5‐HT1A–cortical thickness association correlated significantly with the number of tracts to each region, supporting our hypothesis. We posit a defect in the raphe may affect the PCC within the default mode network in MDD through serotonergic fibers, resulting in increased ruminative processing.

OptoGap: an optogenetics-enabled assay for quantification of cell-cell coupling in multicellular cardiac tissue

Intercellular electrical coupling is an essential means of communication between cells. It is important to obtain quantitative knowledge of such coupling between cardiomyocytes and nonexcitable cells when, for example, pathological electrical coupling between myofibroblasts and cardiomyocytes yields increased arrhythmia risk or during the integration of donor (e.g. cardiac progenitor) cells with native cardiomyocytes in cell-therapy approaches. Currently, there is no direct method for assessing heterocellular coupling within multicellular tissue. Here we demonstrate experimentally and computationally a new contactless assay for electrical coupling, OptoGap, based on selective illumination of inexcitable cells that express optogenetic actuators and optical sensing of the response of coupled excitable cells, e.g. cardiomyocytes, that are light-insensitive. Cell-cell coupling is quantified by the energy required to elicit an action potential via junctional current from the light-stimulated cell(s). The proposed technique is experimentally validated against the standard indirect approach, GapFRAP, using light-sensitive cardiac fibroblasts and non-transformed cardiomyocytes in a two-dimensional setting. It’s potential applicability to the complex three-dimensional setting of the native heart is corroborated by computational modeling and proper calibration. Intercellular coupling is a fundamental form of communication between cells, essential for the synchronization of physiological processes in different organs. Pathologically altered coupling or the emergence of de novo coupling between native and donor cells are problems of interest in many cardiac applications, e.g. during cell delivery and cell integration for cardiac repair therapy1,2. In particular, interactions between cardiomyocytes and fibroblasts are of interest, especially the pro-arrhythmic increase in coupling as the latter transition to myofibroblasts3-6. Electrical coupling in cardiac tissue is mediated primarily by low-resistance paths formed by gap-junctional proteins (connexins), that can link cardiomyocytes (CMs) to each other and to non-cardiomyocytes (nCMs), such as fibroblasts. Qualitative and quantitative methods, e.g. immunofluorescence, messenger RNA and Western blots, are often used to assay connexin expression levels as a surrogate measure of coupling, but they do not provide functional information. A method for direct quantification of cell-cell coupling within the multicellular tissue context is highly desirable.

OptoDyCE: Automated system for high-throughput all-optical dynamic cardiac electrophysiology

In the last two decades, <30% of drugs withdrawals from the market were due to cardiac toxicity, where unintended interactions with ion channels disrupt the heart’s normal electrical function. Consequently, all new drugs must undergo preclinical testing for cardiac liability, adding to an already expensive and lengthy process. Recognition that proarrhythmic effects often result from drug action on multiple ion channels demonstrates a need for integrative and comprehensive measurements. Additionally, patient-specific therapies relying on emerging technologies employing stem-cell derived cardiomyocytes (e.g. induced pluripotent stem-cell-derived cardiomyocytes, iPSC-CMs) require better screening methods to become practical. However, a high-throughput, cost-effective approach for cellular cardiac electrophysiology has not been feasible. Optical techniques for manipulation and recording provide a contactless means of dynamic, high-throughput testing of cells and tissues. Here, we consider the requirements for all-optical electrophysiology for drug testing, and we implement and validate OptoDyCE, a fully automated system for all-optical cardiac electrophysiology. We demonstrate the high-throughput capabilities using multicellular samples in 96-well format by combining optogenetic actuation with simultaneous fast high-resolution optical sensing of voltage or intracellular calcium. The system can also be implemented using iPSC-CMs and other cell-types by delivery of optogenetic drivers, or through the modular use of dedicated light-sensitive somatic cells in conjunction with non-modified cells. OptoDyCE provides a truly modular and dynamic screening system, capable of fully-automated acquisition of high-content information integral for improved discovery and development of new drugs and biologics, as well as providing a means of better understanding of electrical disturbances in the heart.

Disease Modeling in Human Induced Pluripotent Stem Cell Derived Cardiomyocytes Using High-Throughput All-Optical Dynamic Cardiac Electrophysiology

We present an all-optical high-throughput system for phenotyping and monitoring iPSC-CMs, with capabilities for performing personalized cardiotoxicity screening. We demonstrate the system’s utility for characterizing a new disease model in iPSC-CMs.