Peter Bovendeerd

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.

Optimization of cardiac fiber orientation for homogeneous fiber strain during ejection

AbstractThe strain of muscle fibers in the heart is likely to be distributed uniformly over the cardiac walls during the ejection period of the cardiac cycle. Mathematical models of left ventricular (LV) wall mechanics have shown that the distribution of fiber strain during ejection is sensitive to the orientation of muscle fibers in the wall. In the present study, we tested the hypothesis that fiber orientation in the LV wall is such that fiber strain during ejection is as homogeneous as possible. A finite-element model of LV wall mechanics was set up to compute the distribution of fiber strain at the beginning (BE) and end (EE) of the ejection period of the cardiac cycle, with respect to a middiastolic reference state. The distribution of fiber orientation over the LV wall, quantified by three parameters, was systematically varied to minimize regional differences in fiber shortening during ejection and in the average of fiber strain at BE and EE. A well-defined optimum in the distribution of fiber orientation was found which was not significantly different from anatomical measurements. After optimization, the average of fiber strain at BE and EE was 0.025 ± 0.011 (mean ± standard deviation) and the difference in fiber strain during ejection was 0.214 ± 0.018. The results indicate that the LV structure is designed for maximum homogeneity of fiber strain during ejection.

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.

Dependence of Intramyocardial Pressure and Coronary Flow on Ventricular Loading and Contractility: A Model Study

The phasic coronary arterial inflow during the normal cardiac cycle has been explained with simple (waterfall, intramyocardial pump) models, emphasizing the role of ventricular pressure. To explain changes in isovolumic and low afterload beats, these models were extended with the effect of three-dimensional wall stress, nonlinear characteristics of the coronary bed, and extravascular fluid exchange. With the associated increase in the number of model parameters, a detailed parameter sensitivity analysis has become difficult. Therefore we investigated the primary relations between ventricular pressure and volume, wall stress, intramyocardial pressure and coronary blood flow, with a mathematical model with a limited number of parameters. The model replicates several experimental observations: the phasic character of coronary inflow is virtually independent of maximum ventricular pressure, the amplitude of the coronary flow signal varies about proportionally with cardiac contractility, and intramyocardial pressure in the ventricular wall may exceed ventricular pressure. A parameter sensitivity analysis shows that the normalized amplitude of coronary inflow is mainly determined by contractility, reflected in ventricular pressure and, at low ventricular volumes, radial wall stress. Normalized flow amplitude is less sensitive to myocardial coronary compliance and resistance, and to the relation between active fiber stress, time, and sarcomere shortening velocity.

Towards model-based analysis of cardiac MR tagging data: relation between left ventricular shear strain and myofiber orientation

Many cardiac pathologies are reflected in abnormal myocardial deformation, accessible through magnetic resonance tagging (MRT). Interpretation of the MRT data is difficult, since the relation between pathology and deformation is not straightforward. Mathematical models of cardiac mechanics could be used to translate measured abnormalities into the underlying pathology, but, so far, they even fail to correctly simulate myocardial deformation in the healthy heart. In this study we investigated to what extent (1) our previously published three-dimensional finite element model of cardiac mechanics [Kerckhoffs, R.C.P., Bovendeerd, P.H.M., Kotte, J.C.S., Prinzen, F.W., Smits, K., Arts, T., 2003. Homogeneity of cardiac contraction despite physiological asynchrony of depolarization: a model study. Ann. Biomed. Eng. 31, 536-547] can simulate measured cardiac deformation, and (2) discrepancies between strains in model and experiment are related to the choice of the myofiber orientation in the model. To this end, we measured midwall circumferential strain E(cc) and circumferential-radial shear strain E(cr) in three healthy subjects using MRT. E(cc) as computed in the model agreed well with measured E(cc). Computed E(cr) differed significantly from measured E(cr). The time course of E(cr) was found to be very sensitive to the choice of the myofiber orientation, in particular to the choice of the transverse angle. Discrepancies between circumferential-radial shear strain in model and experiment were reduced strongly by increasing the transverse angle in the original model by 25%.

Computational modeling of volumetric soft tissue growth: application to the cardiac left ventricle

As an initial step to investigate stimulus–response relations in growth and remodeling (G&R) of cardiac tissue, this study aims to develop a method to simulate 3D-inhomogeneous volumetric growth. Growth is regarded as a deformation that is decomposed into a plastic component which describes unconstrained growth and an elastic component to satisfy continuity of the tissue after growth. In current growth models, a single reference configuration is used that remains fixed throughout the entire growth process. However, considering continuous turnover to occur together with growth, such a fixed reference is unlikely to exist in reality. Therefore, we investigated the effect of tissue turnover on growth by incrementally updating the reference configuration. With both a fixed reference and an updated reference, strain-induced cardiac growth in magnitude of 30% could be simulated. However, with an updated reference, the amplitude of the stimulus for growth decreased over time, whereas with a fixed reference this amplitude increased. We conclude that, when modeling volumetric growth, the choice of the reference configuration is of great importance for the computed growth.

Computational analysis of the myocardial structure: Adaptation of cardiac myofiber orientations through deformation

Deformation and structure of the cardiac wall can be assessed non-invasively by imaging techniques such as magnetic resonance imaging. Understanding the (patho-)physiology that underlies the observed deformation and structure is critical for clinical diagnosis. However, much about the genesis of deformation and structure is unknown. In the present computational model study, we hypothesize that myofibers locally adapt their orientation to achieve minimal fiber-cross fiber shear strain during the cardiac cycle. This hypothesis was tested in a 3D finite element model of left ventricular (LV) mechanics by computation of tissue deformations and subsequent adaptation of initial myofiber orientations towards those in the deformed tissue. As a consequence of adaptation, local tissue peak stress, strain during ejection and stroke work density were all found to increase by at least 10%, as well as to become 50% more homogeneous throughout the wall. Global LV work (peak systolic pressure, stroke volume and stroke work) increased significantly as well (>9%). The model-predicted myofiber orientations were found to be similar to those in experiments. To the best of our knowledge the presented model is the first that is able to simultaneously predict a realistic myocardial structure as well as to account for the experimentally observed homogeneity in local mechanics.

Optimization of cardiac fiber orientation for homogeneous fiber strain at beginning of ejection

Mathematical models of left ventricular (LV) wall mechanics show that fiber stress depends heavily on the choice of muscle fiber orientation in the wall. This finding brought us to the hypothesis that fiber orientation may be such that mechanical load in the wall is homogeneous. Aim of this study was to use the hypothesis to compute a distribution of fiber orientation within the wall. In a finite element model of LV wall mechanics, fiber stresses and strains were calculated at beginning of ejection (BE). Local fiber orientation was quantified by helix (HA) and transverse (TA) fiber angles using a coordinate system with local r-, c-, and l-directions perpendicular to the wall, along the circumference and along the meridian, respectively. The angle between the c-direction and the projection of the fiber direction on the cl-plane (HA) varied linearly with transmural position in the wall. The angle between the c-direction and the projection of the fiber direction on the cr-plane (TA) was zero at the epicardial and endocardial surfaces. Midwall TA increased with distance from the equator. Fiber orientation was optimized so that fiber strains at BE were as homogeneous as possible. By optimization with TA = 0 degree, HA was found to vary from 81.0 degrees at the endocardium to -35.8 degrees at the epicardium. Inclusion of TA in the optimization changed these angles to respectively 90.1 degrees and -48.2 degrees while maximum TA was 15.3 degrees. Then the standard deviation of fiber strain (epsilon f) at BE decreased from +/- 12.5% of mean epsilon f to +/- 9.5%. The root mean square (RMS) difference between computed HA and experimental data reported in literature was 15.0 degrees compared to an RMS difference of 11.6 degrees for a linear regression line through the latter data.

Determinants of left ventricular shear strain

Mathematical models of cardiac mechanics can potentially be used to relate abnormal cardiac deformation, as measured noninvasively by ultrasound strain rate imaging or magnetic resonance tagging (MRT), to the underlying pathology. However, with current models, the correct prediction of wall shear strain has proven to be difficult, even for the normal healthy heart. Discrepancies between simulated and measured strains have been attributed to 1) inadequate modeling of passive tissue behavior, 2) neglecting active stress development perpendicular to the myofiber direction, or 3) neglecting crossover of myofibers in between subendocardial and subepicardial layers. In this study, we used a finite-element model of left ventricular (LV) mechanics to investigate the sensitivity of midwall circumferential-radial shear strain (E(cr)) to settings of parameters determining passive shear stiffness, cross-fiber active stress development, and transmural crossover of myofibers. Simulated time courses of midwall LV E(cr) were compared with time courses obtained in three healthy volunteers using MRT. E(cr) as measured in the volunteers during the cardiac cycle was characterized by an amplitude of approximately 0.1. In the simulations, a realistic amplitude of the E(cr) signal could be obtained by tuning either of the three model components mentioned above. However, a realistic time course of E(cr), with virtually no change of E(cr) during isovolumic contraction and a correct base-to-apex gradient of E(cr) during ejection, could only be obtained by including transmural crossover of myofibers. Thus, accounting for this crossover seems to be essential for a realistic model of LV wall mechanics.

A method for the quantification of the pressure dependent 3D collagen configuration in the arterial adventitia

Collagen plays an important role in the response of the arterial wall to mechanical loading and presumably has a load-bearing function preventing overdistension. Collagen configuration is important for understanding this role, in particular in mathematical models of arterial wall mechanics. In this study a new method is presented to image and quantify this configuration. Collagen in the arterial adventitia is stained with CNA35, and imaged in situ at high resolution with confocal microscopy at luminal pressures from 0 to 140mm Hg. The images are processed with a new automatic approach, utilizing techniques intended for MRI-DTI data. Collagen configuration is quantified through three parameters: the waviness, the transmural angle and the helical angle. The method is demonstrated for the case of carotid arteries of the white New Zealand rabbit. The waviness indicated a gradual straightening between 40 and 80mm Hg. The transmural angle was about zero indicating that the fibers stayed within an axial-circumferential plane at all pressures. The helical angle was characterized by a symmetrical distribution around the axial direction, indicating a double symmetrical helix. The method is the first to combine high resolution imaging with a new automatic image processing approach to quantify the 3D configuration of collagen in the adventitia as a function of pressure.