Modeling an Excitable Biosynthetic Tissue with Inherent Variability for Paired Computational-Experimental Studies

A probabilistic modeling framework for the prediction of spontaneous premature beats and reentry initiation

BACKGROUND

Spontaneously occurring life-threatening reentrant arrhythmias result when a propagating premature beat encounters a region with significant dispersion of refractoriness. Although localized structural tissue heterogeneities and prescribed cell functional gradients have been incorporated into computational electrophysiologic models, a quantitative framework for the evolution from normal to abnormal behavior that occurs by disease is lacking.

OBJECTIVE

The purpose of this study was to develop a probabilistic modeling framework representing the complex interplay of cell function and tissue structure in health and disease that predicts the emergence of premature beats and the initiation of reentry.

METHODS

An action potential model of the rabbit was developed with data-driven uncertainty characterization as done previously. A novel tissue model using the discrete-cell monodomain equations was developed by implementing cellular uncertainty as a random spatial field.

RESULTS

Cellular action potentials exhibited a wide range of duration and even a variety of behaviors, with 67% exhibiting normal repolarization, 27% displaying early afterdepolarizations, and 6% showing repolarization failure. Nevertheless, simulations in tissue resulted in localized synchronized repolarization. Thus, cellular variability provided "tissue-level robustness," and premature beats and reentry induction were never observed even with abnormalities in cell function (IKr block) or tissue structure (increased tissue resistance). Alterations of both cell function and tissue structure were necessary for the generation of premature beats and arrhythmia initiation.

CONCLUSION

Once extended to whole hearts and validated for a specific context, this modeling framework provides a means to predict the probability of the initiation of life-threatening arrhythmias.

Prokaryotic voltage-gated sodium channels are more effective than endogenous Nav1.5 channels in rescuing cardiac action potential conduction: an in silico study

Methods to augment Na+ current in cardiomyocytes hold potential for the treatment of various cardiac arrhythmias involving conduction slowing. Because the gene coding cardiac Na+ channel (Nav1.5) is too large to fit in a single adeno-associated virus (AAV) vector, new gene therapies are being developed to enhance endogenous Nav1.5 current (by overexpression of chaperon molecules or use of multiple AAV vectors) or to exogenously introduce prokaryotic voltage-gated Na+ channels (BacNav) whose gene size is significantly smaller than that of the Nav1.5. In this study, based on experimental measurements in heterologous expression systems, we developed an improved computational model of the BacNav channel, NavSheP D60A. We then compared in silico how NavSheP D60A expression vs. Nav1.5 augmentation affects the electrophysiology of cardiac tissue. We found that the incorporation of BacNav channels in both adult guinea pig and human cardiomyocyte models increased their excitability and reduced action potential duration. When compared with equivalent augmentation of Nav1.5 current in simulated settings of reduced tissue excitability, the addition of the BacNav current was superior in improving the safety of conduction under conditions of current source-load mismatch, reducing the vulnerability to unidirectional conduction block during premature pacing, preventing the instability and breakup of spiral waves, and normalizing the conduction and ECG in Brugada syndrome tissues with mutated Nav1.5. Overall, our studies show that compared with a potential enhancement of the endogenous Nav1.5 current, expression of the BacNav channels with their slower inactivation kinetics can provide greater anti-arrhythmic benefits in hearts with compromised action potential conduction.NEW & NOTEWORTHY Slow action potential conduction is a common cause of various cardiac arrhythmias; yet, current pharmacotherapies cannot augment cardiac conduction. This in silico study compared the efficacy of recently proposed antiarrhythmic gene therapy approaches that increase peak sodium current in cardiomyocytes. When compared with the augmentation of endogenous sodium current, expression of slower-inactivating bacterial sodium channels was superior in preventing conduction block and arrhythmia induction. These results further the promise of antiarrhythmic gene therapies targeting sodium channels.

On-demand serum-free media formulations for human hematopoietic cell expansion using a high dimensional search algorithm

Substitution of serum and other clinically incompatible reagents is requisite for controlling product quality in a therapeutic cell manufacturing process. However, substitution with chemically defined compounds creates a complex, large-scale optimization problem due to the large number of possible factors and dose levels, making conventional process optimization methods ineffective. We present a framework for high-dimensional optimization of serum-free formulations for the expansion of human hematopoietic cells. Our model-free approach utilizes evolutionary computing principles to drive an experiment-based feedback control platform. We validate this method by optimizing serum-free formulations for first, TF-1 cells and second, primary T-cells. For each cell type, we successfully identify a set of serum-free formulations that support cell expansions similar to the serum-containing conditions commonly used to culture these cells, by experimentally testing less than 1 × 10-5 % of the total search space. We also demonstrate how this iterative search process can provide insights into factor interactions that contribute to supporting cell expansion.

Microheterogeneity-induced conduction slowing and wavefront collisions govern macroscopic conduction behavior: A computational and experimental study

The incidence of cardiac arrhythmias is known to be associated with tissue heterogeneities including fibrosis. However, the impact of microscopic structural heterogeneities on conduction in excitable tissues remains poorly understood. In this study, we investigated how acellular microheterogeneities affect macroscopic conduction under conditions of normal and reduced excitability by utilizing a novel platform of paired in vitro and in silico studies to examine the mechanisms of conduction. Regular patterns of nonconductive micro-obstacles were created in confluent monolayers of the previously described engineered-excitable Ex293 cell line. Increasing the relative ratio of obstacle size to intra-obstacle strand width resulted in significant conduction slowing up to 23.6% and a significant increase in wavefront curvature anisotropy, a measure of spatial variation in wavefront shape. Changes in bulk electrical conductivity and in path tortuosity were insufficient to explain these observed macroscopic changes. Rather, microscale behaviors including local conduction slowing due to microscale branching, and conduction acceleration due to wavefront merging were shown to contribute to macroscopic phenomena. Conditions of reduced excitability led to further conduction slowing and a reversal of wavefront curvature anisotropy due to spatially non-uniform effects on microscopic slowing and acceleration. This unique experimental and computation platform provided critical mechanistic insights in the impact of microscopic heterogeneities on macroscopic conduction, pertinent to settings of fibrotic heart disease.

It is well known that perturbations in the heart structure are associated with the initiation and maintenance of clinically significant cardiac arrhythmia. While previous studies have examined how single structural perturbations affect local electrical conduction, our understanding of how numerous microscopic heterogeneities act in aggregate to alter macroscopic electrical behavior is limited. In this study, we utilized simplified engineered excitable cells that contain the minimal machinery of excitability and can be directly computationally modeled. By pairing experimental and computational studies, we showed that the microscopic branching and collisions of electrical waves slow and speed conduction, respectively, resulting in macroscopic changes in the speed and pattern of electrical activation. These microscale behaviors are significantly altered under reduced excitability, resulting in exaggerated collision effects. Overall, this study helps improve our understanding of how microscopic structural heterogeneities in excitable tissue lead to abnormal action potential propagation, conducive to arrhythmias.

Generation and customization of biosynthetic excitable tissues for electrophysiological studies and cell-based therapies

We describe a two-stage protocol to generate electrically excitable and actively conducting cell networks with stable and customizable electrophysiological phenotypes. Using this method, we have engineered monoclonally derived excitable tissues as a robust and reproducible platform to investigate how specific ion channels and mutations affect action potential (AP) shape and conduction. In the first stage of the protocol, we combine computational modeling, site-directed mutagenesis, and electrophysiological techniques to derive optimal sets of mammalian and/or prokaryotic ion channels that produce specific AP shape and conduction characteristics. In the second stage of the protocol, selected ion channels are stably expressed in unexcitable human cells by means of viral or nonviral delivery, followed by flow cytometry or antibiotic selection to purify the desired phenotype. This protocol can be used with traditional heterologous expression systems or primary excitable cells, and application of this method to primary fibroblasts may enable an alternative approach to cardiac cell therapy. Compared with existing methods, this protocol generates a well-defined, relatively homogeneous electrophysiological phenotype of excitable cells that facilitates experimental and computational studies of AP conduction and can decrease arrhythmogenic risk upon cell transplantation. Although basic cell culture and molecular biology techniques are sufficient to generate excitable tissues using the described protocol, experience with patch-clamp techniques is required to characterize and optimize derived cell populations.

Virtual cardiac monolayers for electrical wave propagation

The complex structure of cardiac tissue is considered to be one of the main determinants of an arrhythmogenic substrate. This study is aimed at developing the first mathematical model to describe the formation of cardiac tissue, using a joint in silico-in vitro approach. First, we performed experiments under various conditions to carefully characterise the morphology of cardiac tissue in a culture of neonatal rat ventricular cells. We considered two cell types, namely, cardiomyocytes and fibroblasts. Next, we proposed a mathematical model, based on the Glazier-Graner-Hogeweg model, which is widely used in tissue growth studies. The resultant tissue morphology was coupled to the detailed electrophysiological Korhonen-Majumder model for neonatal rat ventricular cardiomyocytes, in order to study wave propagation. The simulated waves had the same anisotropy ratio and wavefront complexity as those in the experiment. Thus, we conclude that our approach allows us to reproduce the morphological and physiological properties of cardiac tissue.