Mathematical modeling of left ventricular dimensional changes in mice during aging
© Yang et al.; licensee BioMed Central Ltd. 2012
Published: 17 December 2012
Skip to main content
© Yang et al.; licensee BioMed Central Ltd. 2012
Published: 17 December 2012
Cardiac aging is characterized by diastolic dysfunction of the left ventricle (LV), which is due in part to increased LV wall stiffness. In the diastolic phase, myocytes are relaxed and extracellular matrix (ECM) is a critical determinant to the changes of LV wall stiffness. To evaluate the effects of ECM composition on cardiac aging, we developed a mathematical model to predict LV dimension and wall stiffness changes in aging mice by integrating mechanical laws and our experimental results. We measured LV dimension, wall thickness, LV mass, and collagen content for wild type (WT) C57/BL6J mice of ages ranging from 7.3 months to those of 34.0 months. The model was established using the thick wall theory and stretch-induced tissue growth to an isotropic and homogeneous elastic composite with mixed constituents. The initial conditions of the simulation were set based on the data from the young mice. Matlab simulations of this mathematical model demonstrated that the model captured the major features of LV remodeling with age and closely approximated experimental results. Specifically, the temporal progression of the LV interior and exterior dimensions demonstrated the same trend and order-of-magnitude change as our experimental results. In conclusion, we present here a validated mathematical model of cardiac aging that applies the thick-wall theory and stretch-induced tissue growth to LV remodeling with age.
Over 70% of 50 million Americans over 60 years of age have cardiovascular disease (CVD) [1, 2]. CVD prevalence increases with age, and outcomes in older patients with acute coronary syndromes are poor [1, 3]. It is important to note, however, that aging itself, even in the absence of CVD, alters LV structure and function and impairs the ability of the LV to respond to stress and injury. Thus, understanding the mechanisms of cardiac aging has significant clinical relevance.
While most studies focus on the myocyte contribution to cardiac systolic function, indices for systolic function, such as ejection fraction, systolic velocities, and systolic isovolumic acceleration rate, have been shown to have little relation with age in both clinical and animal studies [4, 5]. In contrast, echocardiographic indices of diastolic function including early (E) and late (Atrial-A) diastolic peak filling velocities and the E and A velocity ratio (E/A ratio), mitral deceleration time (the time from the peak to the end of the Doppler E-wave), isovolumic relaxation time (the time between the closing artifact of the aortic valve), and the earliest detection of trans-mitral blood flow have been demonstrated to decline with age in clinical studies [6–8]. With physiological aging, the LV undergoes monotonic structural changes that include increased wall thickness, chamber diameter, and mass [6–8].
We have shown diastolic dysfunction at the organ level in mice during cardiac aging , and diastolic dysfunction is caused by increased myocardial stiffness at the tissue level . The myocardium is composed of myocytes (muscle) surrounded by the ECM environment. Accordingly, myocardial stiffness is determined by the volume ratio and the combined mechanical property of the myocytes and ECM. In the diastolic phase, myocytes are relaxed and ECM is a critical determinant to the change of LV wall stiffness. About 90% of cardiac ECM composition in the young LV is collagen I and III. We have shown that collagen content in senescent mouse hearts doubles compared to young hearts [9, 11]. Since collagen has a magnitude increase in stiffness over myocytes, the age-related increase of collagen content in the LV shifts myocardial mechanical properties from a myocyte-based stiffness to one influenced by collagen-based stiffness . Therefore, the goal of this study was to evaluate the effects of ECM composition on LV remodeling with aging using a mathematical model developed by integrating cardiac mechanics and our experimental results in mice.
Different LV wall stress models, such as the Laplace law based thin-wall models, thick-wall shell models, and finite element models, have been established to describe LV mechanics and compute stress [13, 14]. Most current thin-wall and thick-wall models of LV remodeling were established assuming idealistic spherical, spheroidal, or ellipsoidal geometries. While finite element models allow some flexibility in the LV geometry, they require high computational power. Currently, these models focus on the stress calculation taken from a particular snapshot of the cardiac cycle or for the entire cycle. However, the interplay between LV stress and strain and between remodeling and geometric evolution on the life-time scale has not been established.
Phenomenological models on geometric remodeling with aging have been developed to apply on arteries under hypertensive conditions [15, 16]. The arterial wall is considered as a thick-wall vessel, and the remodeling equation on either inner or outer artery surface has been postulated as a function of strain and stress at the corresponding location. While these models do not consider the intrinsic relationships between phenomenological assumptions, they provide a possible methodology to model the temporal progression of tissue remodeling with aging.
There are very few computational models available to study LV geometric adaptation with aging. Recently, we established a computational model of LV aging that incorporates Laplace law based stress model . This model captures the overall trend of LV radius change with age. However, this thin-wall model assumes a constant thickness and no stress variation across the wall, LV thickening with age was not addressed. Here, we improve this model by adopting a more sophisticated thick-wall theory and stretch-induced tissue growth theory for the model used in this study. The novelty of this model lies in the integration of both computational and experimental approaches. The wall radii remodeling model established in this study was subject to the temporal function of the total mass which was measured in our experiments. Additionally, the predictions of LV geometry and end diastolic pressure-volume relation from the mathematical model were compared with our experimental measurements to validate predictions of the mathematical model.
The LV was modeled as a concentric spherical shell with interior and exterior radii and , respectively, at the no-load state. Due to the transmural pressure difference on the inner and outer surfaces, the deformed inner and outer radii were and , respectively. The ventricular wall was modeled as an isotropic and linear elastic material of two constituents, cardiac myocytes and ECM which is represented by collagen.
Internal and external radii, ie., R i and R o , as well as wall thickness R o - R i , all have non-decreasing trend with age.
It's concluded that Poisson's ratios of all stable isotropic materials falls in the region (-1, 0.5) and only auxetic materials with honeycomb structures and networks have been found to have negative values . Apparently, heart tissue is not auxetic. Therefore, for LV, its Poisson's ratios falls in the region 0 < v < 0.5. As a consequence, it is apparent that the loaded wall thickness, i.e. r o - r i is smaller than the zero-load thickness, i.e. R o - R i . The LV transmural pressure P has been shown to alter with aging. We adopted the temporal profile of pressure difference reported in the literature [11, 19].
In equation (7), V c (t)/V m (t) is the volume occupied by collagen/muscle at time t, and M c (t)/M m (t) is the mass of collagen/muscle at time t. Parameters ρ c and ρ m are mass densities of collagen and myocytes. The density of collagen is an intrinsic property of the material and therefore does not change with time unless the collagen molecular structure undergoes significant changes. The density of collagen is adopted as ρ c = 1.70 g/ml . For muscle, due to the fact that approximately 80% of the mass is water as has been measured in experiments for multiple mammalian muscle (either cardiac or non-cardiac), the density is close to that of water, with variation contributed by tissue and tissue solids. We adopt myocyte density as ρ m = 1.06 g/ml [23, 24].
Since we assumed that the material is isotropic and homogeneous in our model (assumption 3), the spatial remodeling can only occur in the radial direction. As a result, the spherical geometry is always maintained. Therefore, only radii, e.g., the free radii R i and R o alter with aging.
in which, f (t) is the function that characterized the rate of change for the outer radius, parameter τ R is the remodeling characteristic time for R o (t), variable is the ratio of the change between deformed and free radii raised to the power D on the exterior surface, parameter τ f is the decay constant for function f which guarantees a bounded outer dimension of the LV. The parameters τ R , D and τ f are constant parameters to be determined.
Equations (5-7) and (12-14), together with the definitions of α(t) and λ(t), formed a mathematical model to characterize the temporal progression of the LV geometry by applying thick-wall theory and stretch-induced tissue growth. Given the temporal functions M c (t), M m (t) and P(t) with aging, we solved the mathematical model numerically.
Constant parameters and initial values adopted in numerical calculation and the source of the data.
Initial value R i (t = 7.5mons) (mm)
Our experiments on over 148 mice
Muscle density ρ m (gcm-3)
Collagen density ρ c (gcm-3)
Muscle elastic modulus E m (Pa)
Collagen elastic modulus E c (Pa)
Poisson' ratio v
Time step h (hrs)
Further, from equations (2-3), the pressure difference is the only external load that deforms the LV wall. With retreat of this pressure, or establishment of balance between inner and outer surfaces, any stress or strain disappears. Now if we compare the elastic modulus and the pressure difference, the highest pressure, i.e. 7.1 mmHg, is equivalent to 944 Pa, less than 2 orders of magnitude lower the magnitude of Young's modulus. Therefore, without detailed computation, we can see that there should exist none noteworthy deformation (< 1%) and the strained dimensions should be extremely close to the free ones. This was proved by detailed computation shown later.
The radii and were real positive numbers subject to the condition d(t) > c 3(t).
Read initial values R i (t 0), R o (t 0), calculate r i (t 0) and r o (t 0) and , note that α(t 0) = 1. Iterate steps 2-5 until the final step N.
The numerical calculations were performed using in-house Matlab codes. In all simulations, parameter τ R = 0.22 month, D = 2.5 and τ f = 5.6 month, initial values used in our model were listed in Table-1. The initial radii were determined by minimizing the error from experimental measurements considering both end diastolic dimension/diameter (EDD) and wall thickness.
As is seen in Eq. (2), and its contribution is a force toward inner direction. The magnitude of all these averaged stresses went up with time, which was a natural result of enhanced pressure P(t) (Figure 3A). We also showed that the magnitude of the circumferential or meridional stress was higher than that the radial (comparing Figures 4E and 4F).
And the above inequalities do not rely on a concrete system. Both sides are known functions without parameter dependence. Therefore, it is a universal condition for the selection of the type of rate function. f(t) Though this is a stringent condition, functions fall in the proposed region will guarantee the stability of the system.
In our simulations, parameters τ R = 0.22 months, D = 2.5 and τ f = 5.6 months were chosen by minimizing the deviations of computed EDD and wall thickness from measurements at senescent age, i.e. (R i - R i , exp )2 + (WT - WT exp )2, in which WT means wall thickness. The growth function f(t) determined by these parameters, the dotted cyan curve, falls in the bounded region. The analysis of the boundary of the growth function provided a standard procedure for us to determine the possible sets for parameters with or without simulation results, which should be able to be applied in the modeling of age-dependent remodeling process of different species given experimental data, e.g. rat, dog or even human. In addition, if a growth function f(t) of this form is given, we can also use this boundary condition to predict whether LV geometry and wall thickness will monotonically increase or not.
In this study, we established the first analytical mathematical model to quantify long-term LV geometric remodeling dynamics with aging by applying thick-wall theory and stretch-induced tissue growth postulate. This model addressed the temporal progression of cardiac aging and is an advancement from most current models that analyze the static equilibrium of LV remodeling and our previous thin-wall model of LV remodeling. In addition, we analytically determined the boundary for the tissue growth rate to guarantee the stability of the remodeling. This boundary might provide a reference for experimental measurement to examine the remodeling outcomes. Further, the parameters in the mathematical model were determined by using our experimental results from over 140 C57/BL6 mice. Some parameters and initial conditions were selected by minimizing the error between our computational predictions and experimental measurements. In addition, the model was validated by comparing the predictions of internal and external radii and wall thickness in diastole to experimental results. The predicted LV geometry trends were consistent with the experimental results. Thus, this study is a real integration of computational and experimental approaches for model establishment and validation using mice data. This approach and the proposed model can be applied to establish a cardiac aging model for human in the future.
The P-V curve in Figure 5 shows that both pressure and end diastolic volume escalate with age. This has been observed in a LV P-V relationship experiment on rat aging study by Pacher and colleagues . In their study, LV end diastolic volume, LV end diastolic pressure, and the slope of EDPVR were greater in senescent than in young rats. In another experimental aging study on C57BL/6 female mice , the pressure was augmented in 16-month-old mice compared to 6-month-old group. In our experimental study, we observed increased LV volume with age. These experimental results confirmed our computational predictions.
Though our computational results follow the trend of LV geometry changes, the proposed model has some limitations. First, we assumed the LV to be a spherical shell. With this assumption, all directions grow at the same pace to maintain the spherical symmetry. Sphere is a special geometry that requires minimum linear dimension for a constant volume. However, LV does not have a spherical geometry in reality. In addition, we assumed that LV was composed of elastic, isotropic, and homogeneous materials and applied the linear mixture theory to calculate the Young's modulus of the myocardium. Though the assumptions simplify the mathematical model and the computational analysis, such simplification might lead to prediction errors. From our computational results, the radii growth calculated are 3.43% for external, 2.37% for internal radius, and 5.72% for wall thickness collected from C57/BL6 mice. Some of these values were small compared to our echo experiment on C57/BL6 mice, in which EDD grows 2.60% and wall thickness is up 13.20% from young to senescent. In addition, our previous experimental work on CB6F1 mice reported EDD and wall thickness have 12.73 and 13.00% growth, respectively, in the LV between the groups of young and senescent . This discrepancy may be originated from the simplification of the geometry we chose or the prediction of the myocardium stiffness. Future experiments will build on these limitations to improve our results by incorporating more realistic geometry and taking more advanced techniques to determine LV stiffness.
We have established the first mathematical model to study age-related temporal-spatial LV remodeling by adopting thick-wall theory and stretch-induced tissue growth theory. IInputs of the mathematical model were real experimental data including temporal profiles of LV mass, collagen content change, and pressure across LV, which were obtained by over 140 mice. The established model captured the major property of LV remodeling with age and yielded predicted results of LV geometry progressions in mice comparable to experimental results.
The authors acknowledge grant and contract support from NSF CAREER award #0644646 (to HCH), NHLBI HHSN268201000036C (N01-HV-00244), NIH R01 HL75360, Veteran's Administration Merit Award, and the Max and Minnie Tomerlin Voelcker Fund (to MLL), and NIH 1R03EB009496, and NIH SC2HL101430 (to YFJ.).
This article has been published as part of BMC Systems Biology Volume 6 Supplement 3, 2012: Proceedings of The International Conference on Intelligent Biology and Medicine (ICIBM) - Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/6/S3.
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.