Computational Modeling of Cardiac Excitation-Contraction Coupling in Physiological and Pathological Conditions

The cardiomyocyte is a complex biological system where many mechanisms interact non-linearly to regulate the coupling between electrical excitation and mechanical contraction. For this reason, the development of mathematical models is fundamental in the field of cardiac electrophysiology, where the use of computational tools has become complementary to the classical experimentation. My doctoral research has been focusing on the development of such models for investigating the regulation of ventricular excitation-contraction coupling at the single cell level. In particular, the following researches are presented in this thesis: 1) Study of the unexpected deleterious effect of a Na channel blocker on a long QT syndrome type 3 patient. Experimental results were used to tune a Na current model that recapitulates the effect of the mutation and the treatment, in order to investigate how these influence the human action potential. Our research suggested that the analysis of the clinical phenotype is not sufficient for recommending drugs to patients carrying mutations with undefined electrophysiological properties. 2) Development of a model of L-type Ca channel inactivation in rabbit myocytes to faithfully reproduce the relative roles of voltage- and Ca-dependent inactivation. The model was applied to the analysis of Ca current inactivation kinetics during normal and abnormal repolarization, and predicts arrhythmogenic activity when inhibiting Ca-dependent inactivation, which is the predominant mechanism in physiological conditions. 3) Analysis of the arrhythmogenic consequences of the crosstalk between β-adrenergic and Ca-calmodulin dependent protein kinase signaling pathways. The descriptions of the two regulatory mechanisms, both enhanced in heart failure, were integrated into a novel murine action potential model to investigate how they concur to the development of cardiac arrhythmias. These studies show how mathematical modeling is suitable to provide new insights into the mechanisms underlying cardiac excitation-contraction coupling and arrhythmogenesis.

Closed-loop investigation of the dynamical properties of cardiomyocyte electrophysiology by means of Dynamic clamp and Mathematical modelling

The research field of my PhD concerns mathematical modeling and numerical simulation, applied to the cardiac electrophysiology analysis at a single cell level. This is possible thanks to the development of mathematical descriptions of single cellular components, ionic channels, pumps, exchangers and subcellular compartments. Due to the difficulties of vivo experiments on human cells, most of the measurements are acquired in vitro using animal models (e.g. guinea pig, dog, rabbit). Moreover, to study the cardiac action potential and all its features, it is necessary to acquire more specific knowledge about single ionic currents that contribute to the cardiac activity. Electrophysiological models of the heart have become very accurate in recent years giving rise to extremely complicated systems of differential equations. Although describing the behavior of cardiac cells quite well, the models are computationally demanding for numerical simulations and are very difficult to analyze from a mathematical (dynamical-systems) viewpoint. Simplified mathematical models that capture the underlying dynamics to a certain extent are therefore frequently used. The results presented in this thesis have confirmed that a close integration of computational modeling and experimental recordings in real myocytes, as performed by dynamic clamp, is a useful tool in enhancing our understanding of various components of normal cardiac electrophysiology, but also arrhythmogenic mechanisms in a pathological condition, especially when fully integrated with experimental data.

Correlations between cardiac magnetic resonance and histologic findings in Anderson Fabry disease

ABSTRACT Background Cardiac magnetic resonance (CMR) has been shown as promising diagnostic tool in Anderson-Fabry disease (AFD) cardiomyopathy due to its ability to detect fat deposits through lower native T1 values. However no histological validation has been provided to date. Objectives To correlate CMR and histologic findings in different cardiac stages of AFD focusing on T1 mapping. Methods Fifteen AFD patients (49 years [IQR 39-63], 60% females) undergoing CMR (cines, native T1 and T2 mapping, LGE and post-contrast T1 imaging) and endomyocardial biopsy (EMB, n=11) or septal myectomy (n=4), were retrospectively evaluated. Tissue specimens were analyzed with light/electron microscopy and vacuolization amount calculated as percentages of vacuolated myocytes and vacuolated myocyte area (%VMA) through a quantitative histomorphometric color-based analysis. Results In patients without increased indexed left ventricular mass (LVMi) at CMR (67%), T1 fell as %VMA increased (r= -0.883; p<0.001), whereas no clear relationship was evident once increased LVMi occurred (r= -0.501; p=0.389). At least 45% of vacuolized myocytes and 10% of VMA were needed for low T1 to occur. %VMA positively correlate with maximal wall thickness (MWT, r=0.860, p<0.0001) and LVMi (r= 0.762; p<0.001). Increased MWT and LVMi were present with at least 45% and 80% of vacuolated myocytes, respectively, and 18% and 22% of VMA. Conclusions This study demonstrated an inverse correlation between native T1 and the vacuolization amount in patients without increased LVMi at CMR, providing a histological validation of low native T1 in AFD. Importantly, a significant vacuolization burden was needed before low T1 and left ventricle hypertrophy occurred.