Roy Kerckhoffs

Coupling of a 3D Finite Element Model of Cardiac Ventricular Mechanics to Lumped Systems Models of the Systemic and Pulmonic Circulation

In this study we present a novel, robust method to couple finite element (FE) models of cardiac mechanics to systems models of the circulation (CIRC), independent of cardiac phase. For each time step through a cardiac cycle, left and right ventricular pressures were calculated using ventricular compliances from the FE and CIRC models. These pressures served as boundary conditions in the FE and CIRC models. In succeeding steps, pressures were updated to minimize cavity volume error (FE minus CIRC volume) using Newton iterations. Coupling was achieved when a predefined criterion for the volume error was satisfied. Initial conditions for the multi-scale model were obtained by replacing the FE model with a varying elastance model, which takes into account direct ventricular interactions. Applying the coupling, a novel multi-scale model of the canine cardiovascular system was developed. Global hemodynamics and regional mechanics were calculated for multiple beats in two separate simulations with a left ventricular ischemic region and pulmonary artery constriction, respectively. After the interventions, global hemodynamics changed due to direct and indirect ventricular interactions, in agreement with previously published experimental results. The coupling method allows for simulations of multiple cardiac cycles for normal and pathophysiology, encompassing levels from cell to system.

Homogeneity of cardiac contraction despite physiological asynchrony of depolarization : a model study

AbstractThe use of mathematical models combining wave propagation and wall mechanics may provide new insights in the interpretation of cardiac deformation toward various forms of cardiac pathology. In the present study we investigated whether combining accepted mechanisms on propagation of the depolarization wave, time variant mechanical properties of cardiac tissue after depolarization, and hemodynamic load of the left ventricle (LV) by the aortic impedance in a three-dimensional finite element model results in a physiological pattern of cardiac contraction. We assumed that the delay between depolarization for all myocytes and the onset of crossbridge formation was constant. Two simulations were performed, one in which contraction was initiated according to the regular depolarization pattern (NORM simulation), and another in which contraction was initiated after synchronous depolarization (SYNC simulation). In the NORM simulation propagation of depolarization was physiological, but wall strain was unphysiologically inhomogeneous. When simulating LV mechanics with unphysiological synchronous depolarization (SYNC) myofiber strain was more homogeneous and more physiologic. Apparently, the assumption of a constant delay between depolarization and onset of crossbridge formation results in an unrealistic contraction pattern. The present finding may indicate that electromechanical delay times are heterogeneously distributed, such that a contraction in a normal heart is more synchronous than depolarization.

Patient-specific models of cardiac biomechanics

Patient-specific models of cardiac function have the potential to improve diagnosis and management of heart disease by integrating medical images with heterogeneous clinical measurements subject to constraints imposed by physical first principles and prior experimental knowledge. We describe new methods for creating three-dimensional patient-specific models of ventricular biomechanics in the failing heart. Three-dimensional bi-ventricular geometry is segmented from cardiac CT images at end-diastole from patients with heart failure. Human myofiber and sheet architecture is modeled using eigenvectors computed from diffusion tensor MR images from an isolated, fixed human organ-donor heart and transformed to the patient-specific geometric model using large deformation diffeomorphic mapping. Semi-automated methods were developed for optimizing the passive material properties while simultaneously computing the unloaded reference geometry of the ventricles for stress analysis. Material properties of active cardiac muscle contraction were optimized to match ventricular pressures measured by cardiac catheterization, and parameters of a lumped-parameter closed-loop model of the circulation were estimated with a circulatory adaptation algorithm making use of information derived from echocardiography. These components were then integrated to create a multi-scale model of the patient-specific heart. These methods were tested in five heart failure patients from the San Diego Veterans Affairs Medical Center who gave informed consent. The simulation results showed good agreement with measured echocardiographic and global functional parameters such as ejection fraction and peak cavity pressures.

Patient-Specific Modeling of Dyssynchronous Heart Failure: A Case Study

The development and clinical use of patient-specific models of the heart is now a feasible goal. Models have the potential to aid in diagnosis and support decision-making in clinical cardiology. Several groups are now working on developing multi-scale models of the heart for understanding therapeutic mechanisms and better predicting clinical outcomes of interventions such as cardiac resynchronization therapy. Here we describe the methodology for generating a patient-specific model of the failing heart with a myocardial infarct and left ventricular bundle branch block. We discuss some of the remaining challenges in developing reliable patient-specific models of cardiac electromechanical activity, and identify some of the main areas for focusing future research efforts. Key challenges include: efficiently generating accurate patient-specific geometric meshes and mapping regional myofiber architecture to them; modeling electrical activation patterns based on cellular alterations in human heart failure, and estimating regional tissue conductivities based on clinically available electrocardiographic recordings; estimating unloaded ventricular reference geometry and material properties for biomechanical simulations; and parameterizing systemic models of circulatory dynamics from available hemodynamic measurements.

Computational Methods for Cardiac Electromechanics

Computational modeling provides a potentially powerful way to integrate structural properties measured in vitro to physiological functions measured in vivo. Focusing on the various scales (cell-tissue-organ-system), we give an overview of the importance and applications of numerical models of ventricular anatomy, electrophysiology, mechanics, and circulatory models. The integration of these models in one multiscale model of cardiac electromechanics is discussed in the light of applications to hypothesis generation, diagnosis, surgery(planning, training, and outcome of interventions), and therapies. Special attention is paid to practical use in terms of computational demand. Because of growing computer power and the development of efficient algorithms, we expect that real-time simulations with multiscale models of cardiac electromechanics become feasible in 2008 (despite the increasing complexity of models due to data accumulation on molecular and cellular mechanisms).

Current progress in patient-specific modeling

We present a survey of recent advancements in the emerging field of patient-specific modeling (PSM). Researchers in this field are currently simulating a wide variety of tissue and organ dynamics to address challenges in various clinical domains. The majority of this research employs three-dimensional, image-based modeling techniques. Recent PSM publications mostly represent feasibility or preliminary validation studies on modeling technologies, and these systems will require further clinical validation and usability testing before they can become a standard of care. We anticipate that with further testing and research, PSM-derived technologies will eventually become valuable, versatile clinical tools.

Cardiac resynchronization: Insight from experimental and computational models

Cardiac resynchronization therapy (CRT) is a promising therapy for heart failure patients with a conduction disturbance, such as left bundle branch block. The aim of CRT is to resynchronize contraction between and within ventricles. However, about 30% of patients do not respond to this therapy. Therefore, a better understanding is needed for the relation between electrical and mechanical activation. In this paper, we focus on to what extent animal experiments and mathematical models can help in order to understand the pathophysiology of asynchrony to further improve CRT.

Effects of biventricular pacing and scar size in a computational model of the failing heart with left bundle branch block

OBJECTIVES To study the impact of biventricular pacing (BiV) and scar size on left ventricular (LV) regional and global function using a detailed finite element model of three-dimensional electromechanics in the failing canine heart. BACKGROUND Cardiac resynchronization therapy (CRT) clinical trials have demonstrated that up to 30% of patients may be classified as non-responders. The presence of a scar appears to contribute to those that do not respond to CRT. A recent study in patients with myocardial scar showed that LV dyssynchrony was the sole independent predictor of reverse remodeling, and not scar location or size. METHODS Two activation sequences were simulated: left bundle branch block (LBBB) and acute simultaneous BiV (with leads in the left and right ventricle) in hearts with chronic scars of various sizes. The dependence of regional function (mean fiber ejection strain, variance of fiber isovolumic strain and fraction of tissue stretched during ejection) and global function (left ventricular dP/dt(max), ejection fraction, stroke work) on scar size and pacing protocol was tested. RESULTS Global function and regional function averaged over the whole LV during LBBB and BiV decreased as a function of scar size. In the non-scarred regions, however, regional function was largely independent of scar size for a fixed pacing site. CONCLUSIONS The model results suggest that uniformity of mechanical contraction in non-scarred regions in the failing heart during biventricular pacing is independent of scar size for a fixed pacing site.

Timing of depolarization and contraction in the paced canine left ventricle : model and experiment

Introduction: For efficient pump function, contraction of the heart should be as synchronous as possible. Ventricular pacing induces asynchrony of depolarization and contraction. The degree of asynchrony depends on the position of the pacing electrode. The aim of this study was to extend an existing numerical model of electromechanics in the left ventricle (LV) to the application of ventricular pacing. With the model, the relation between pacing site and patterns of depolarization and contraction was investigated.

Ventricular Dilation and Electrical Dyssynchrony Synergistically Increase Regional Mechanical Nonuniformity But Not Mechanical Dyssynchrony A Computational Model

Background—Heart failure (HF) in combination with mechanical dyssynchrony is associated with a high mortality rate. To quantify contractile dysfunction in patients with HF, investigators have proposed several indices of mechanical dyssynchrony, including percentile range of time to peak shortening (WTpeak), circumferential uniformity ratio estimate (CURE), and internal stretch fraction (ISF). The goal of this study was to compare the sensitivity of these indices to 4 major abnormalities responsible for cardiac dysfunction in dyssynchronous HF: dilation, negative inotropy, negative lusitropy, and dyssynchronous activation. Methods and Results—All combinations of these 4 major abnormalities were included in 3D computational models of ventricular electromechanics. Compared with a nonfailing heart model, ventricles were dilated, inotropy was reduced, twitch duration was prolonged, and activation sequence was changed from normal to left bundle branch block. In the nonfailing heart, CURE, ISF, and WTpeak were 0.97±0.004, 0.010±0.002, and 78±1 milliseconds, respectively. With dilation alone, CURE decreased 2.0±0.07%, ISF increased 58±47%, and WTpeak increased 31±3%. With dyssynchronous activation alone, CURE decreased 15±0.6%, ISF increased 14-fold (±3), and WTpeak increased 121±4%. With the combination of dilation and dyssynchronous activation, CURE decreased 23±0.8%, ISF increased 20-fold (±5), and WTpeak increased 147±5%. Conclusions—Dilation and left bundle branch block combined synergistically decreased regional cardiac function. CURE and ISF were sensitive to this combination, but WTpeak was not. CURE and ISF also reflected the relative nonuniform distribution of regional work better than WTpeak. These findings might explain why CURE and ISF are better predictors of reverse remodeling in cardiac resynchronization therapy.