Skip to main content

Mathematical modeling of left ventricular dimensional changes in mice during aging

Abstract

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.

Introduction

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 [68]. With physiological aging, the LV undergoes monotonic structural changes that include increased wall thickness, chamber diameter, and mass [68].

We have shown diastolic dysfunction at the organ level in mice during cardiac aging [9], and diastolic dysfunction is caused by increased myocardial stiffness at the tissue level [10]. 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 [12]. 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 [11]. 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.

Methods

With physiological aging, the LV undergoes tonic structural and functional changes including increased mean LV wall thickness, chamber diameter, mass, concentric remodeling, and decline in LV diastolic function in humans [68]. We assumed mice age-related remodeling follow similar rules. Our model was established based on the following 2 assumptions.

1) The LV was modeled as a concentric spherical shell with interior and exterior radii R i and R o , respectively, at the no-load state. Due to the transmural pressure difference P on the inner and outer surfaces, the deformed inner and outer radii were r i and r o , 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.

(2) 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.

Thick-wall LV model

Based on assumption 1), deformation of the LV wall obeys Hook's law, which can be expressed in tensor form: i j = 1 + ν E σ i j - ν E δ i j σ k k , where is strain tensor, σ is stress tensor, v is Poisson's ratio, E is material's Young's modulus, δ ij is Kronecker delta, i, j, k {R, θ, φ} are the spatial indices for radial, meridional, and circumferential directions. The strain in the radial direction is:

R R = r R - 1 = 1 E [ σ R R - ν ( σ θ θ + σ ϕ ϕ ) ] .
(1)

The stresses for elastic, isotropic, homogeneous, and spherical thick-wall contains the following analytical form[17]:

σ R R ( R ) = P R i 3 R o 3 - R i 3 ( 1 - R o 3 R 3 ) ,
(2)
σ θ θ ( R ) = σ ϕ ϕ ( R ) = P R i 3 R o 3 - R i 3 ( 1 + 1 2 R o 3 R 3 ) ,
(3)

in which, P denotes pressure difference applied to the inner and exterior walls, and R [R i , R o ]. Applying equations (2) and (3) to equation (1), we had

r R = 1 + P E 1 ( R o / R i ) 3 - 1 [ ( 1 - 2 ν ) - ( 1 + ν ) ( R o / R ) 3 ] .
(4)

Integrate both sides of Eq. (4) from R i to R o , we obtained

r o - r i = { 1 - P E 1 ( R o / R i ) 3 - 1 [ ( 1 + ν ) ( R o / R i ) 2 + R o / R i 2 - ( 1 - 2 ν ) ] } ( R o - R i ) .
(5)

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 [18]. 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].

The Young's modulus E increases with aging due to changes in the LV constituents. We previously reported that the fractions of the two constituents, i.e. myocytes and collagen changes with aging.[9, 11] In this study, we adopted a simplified version of the linear mixture theory of composite material, Young's modulus obeys the following equation [20, 21]:

E ( t ) = E c v c ( t ) + E m [ 1 - v c ( t ) ] ,
(6)

in which E c and E m are two constants representing Young's moduli of collagen and muscle, respectively. Function v c (t) is the volume fraction of collagen, which can be expressed as:

υ c ( t ) = V c ( t ) V c ( t ) + V m ( t ) = M c ( t ) / ρ c M c ( t ) / ρ c + M m ( t ) / ρ m .
(7)

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 [22]. 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].

Due to the fact that both collagen and muscle mass evolve with aging, the volume of the wall changed accordingly. The mass evolution enforces the following equations:

M c ( t ) = ρ c υ c ( t ) R i ( t ) R o ( t ) 4 π R 2 d R ,
(8)
M m ( t ) = ρ m [ 1 - υ c ( t ) ] R i ( t ) R o ( t ) 4 π R 2 d R .
(9)

Take the derivative of equations (8) and (9) with respect to time, we obtained

c ( t ) / ρ c = 4 π υ c ( t ) [ R o 2 ( t ) o ( t ) - R i 2 ( t ) i ( t ) ] + 4 π 3 [ R o 3 ( t ) - R i 3 ( t ) ] υ ˙ c ( t ) ,
(10)
m ( t ) / ρ m = 4 π [ 1 - υ c ( t ) ] [ R o 2 ( t ) o ( t ) - R i 2 ( t ) i ( t ) ] - 4 π 3 [ R o 3 ( t ) - R i 3 ( t ) ] υ ˙ c ( t ) .
(11)

Add up Eqs. (10-11) and divide both sides by 4π, the following equation about free radii was obtained:

R o 2 ( t ) o ( t ) - R i 2 ( t ) i ( t ) = c ( t ) + ρ c ρ m m ( t ) 4 π ρ c .
(12)

Wall incompressibility yields:

r o 3 ( t ) - r i 3 ( t ) = R o 3 ( t ) - R i 3 ( t ) .
(13)

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.

Stretch-induced tissue growth model

The dynamic equation that governs the temporal behavior of the exterior radius in the free state was defined as α ( t ) := R o ( t ) R o ( t 0 ) , which is the ratio of free exterior radius to its initial value. Further, since the deformed wall thickness is smaller than that in the free state, volume conservation between free state and deformed state of LV wall guaranties r o > R o and r i > R i . Additionally, the ratio between deformed and free exterior radii λ ( t ) := r o ( t ) R o ( t ) is always larger than 1 (λ(t) > 1) since r o > R o . We proposed that the rate of change of free exterior radius had a power law relation kernel with the ratio λ(t) and the curve followed an overall exponential decay,

α ˙ ( t ) = f ( t ) = τ R - 1 [ λ D ( t ) - 1 ] e x p [ - ( t - t 0 ) / τ f ] ,
(14)

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 λ D ( t ) -1= r o D ( t ) - R o D ( t ) R o D ( t ) 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.

Input functions and parameters

M c (t), M m (t), and P(t) are temporal functions established by fitting experimental data. We measured the LV mass of 148 wild type mice at different ages by necropsy. We adopted a smooth monotonic function that saturates with time to fit the experimental data for LV mass. The total LV mass function used in our simulation is M(t) = 76.88 + 8.92 {1 - exp[- (t - 8.18)/τ LV } (solid line in Figure 1), the mass unit is milligram, the unit of time t is months, and τ LV = 6.08 months. One of our recent studies showed that the collagen content (volume fraction) in wild type C57/BL6J mice is doubled in senescent mice (~1%) compared to young mice (~0.5%). Our previous studies also demonstrated the same trend of collagen content increase [11]. In our model, we chose 1% collagen content for young at 7.5 months, and approaching 2% for old mice after 30 months. By adopting the converted collagen content, LV mass function M(t) above, and corresponding densities (see Table 1), the temporal curve for collagen mass was M c t = 1 . 21 + 1 . 50 1 - e x p [ - t - 7 . 50 / τ c , with τ c = 7 . 23 months and other units are the same to relevant ones in M t . Since our model assumed that the LV consisted of two constituents (assumption 4), muscle mass were calculated by subtracting collagen mass from the total LV mass as M m t = M t - M c t . Figures 2 showed the temporal profiles of the mass of the total LV, collagen, and myocytes, as well as collagen volume fraction with age. Based on our computational prediction, total LV mass underwent a 12.8% increase from 7.5 to 30 months. The mass growth for collagen was about 121% in parallel with a twofold increase in volume fraction as shown in Figure 2D due to a higher density in relative to muscle. The myocytes mass grew about 11.1% in senescent mice compared to young group.

Figure 1
figure 1

Fitted temporal profile of the LV mass with age. Red rectangles denoted the average values of LV mass at the ages of 8, 15, 20, 24, and 29.5 months collected over 140 wild type C57BL6J mice. Mass was measured by necropsy. The solid curve was the fitted temporal profile based on our experimental data.

Table 1 Constant parameters and initial values adopted in numerical calculation and the source of the data.
Figure 2
figure 2

Computational results for mass functions for LV, collagen, and muscle, as well as volume fraction for collagen. These functions were used as inputs to our mathematical model. (A) The total LV mass fitted with data generated in autopsy experiments; (B) Collagen temporal profile fitted with known trend from experiments; (C) Temporal profile of muscle mass calculated as the difference between total and collagen mass; (D) Volume fraction of collagen, the values for young (6-9 months) and senescent (>25 months) were obtained from published results of mouse.

Time-dependent pressure difference between the interior and exterior surfaces of the LV was another input function. We collected the pressure measurement in the literature and derived its temporal progression as P t = 4 . 5 + 3 . 0 1 - e x p - t - 2 . 5 / τ p with unit mmHg (Figure 3A), where τ p = 14 . 0 months [11, 19]. The pressure ranged from 4.5 mmHg at 2.5 months to 7.1 mmHg at 30 months. We also showed the calculated Young's modulus of LV wall by assuming linear and homogeneously mixed composite made of muscle and collagen (Figure 3B). The modulus increased by 82% in senescent mice compared to young mice. Though there is only an increase from 1% to 2% increase of collagen volume ratio, it gives rise to significant enlargement in the modulus of the composite because the elastic modulus of collagen is much larger compared that of the muscle.

Figure 3
figure 3

Temporal profile of pressure difference of LV wall and predicted Young's modulus of the myocardium. (A) Time-dependent pressure difference across the LV. The curved is derived by fitting experimental data; (B) Elastic modulus of LV, which was calculated as isotropic and homogeneous composite made up of muscle and collagen. The modulus was calculated by taking account of temporal profiles of muscle and collagen in LV.

Further, from equations (2-3), the pressure difference P 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.

Discretization of equations

For convenience in the writing of equations, we defined b ( t ) = c ( t ) + ρ c ρ m m ( t ) 4 π ρ c , c ( t ) = { 1 - P ( t ) E ( t ) 1 R i o ( t ) 3 - 1 [ ( 1 + ν ) R i o ( t ) 2 + R i o ( t ) 2 - ( 1 - 2 ν ) ] } [ R o ( t ) - R i ( t ) ] , in which R i o ( t ) = R o ( t ) / R i ( t ) , and d ( t ) = R o 3 ( t ) - R i 3 ( t ) . Now Equations (5), (12), and (13) can be rewritten as:

R o 2 ( t ) o ( t ) - R i 2 ( t ) i ( t ) = b ( t ) ,
(15)
r o ( t ) - r i ( t ) = c ( t ) ,
(16)
r o 3 ( t ) - r i 3 ( t ) = d ( t ) .
(17)

Given R i (t), R o (t), P(t) and, E(t) variables c(t) and d(t) could be determined and r i (t) and r o (t) can be solved as:

r i ( t ) = - c ( t ) + [ 4 d ( t ) / c ( t ) - c 2 ( t ) ] / 3 2 ,
(18)
r o ( t ) = c ( t ) + [ 4 d ( t ) / c ( t ) - c 2 ( t ) ] / 3 2 .
(19)

The radii r i t and r o t were real positive numbers subject to the condition d(t) > c3(t).

We studied the variation of free and deformed radii in the time domain [t0, t0 + T] and evenly divided this region T by a total of N small constant intervals of h. Equivalently, the domain could be rewritten as [t0, t0 + Nh]. The following is the procedure to calculate time-dependent radii.

1. Read initial values R i (t 0), R o (t 0), calculate r i (t 0) and r o (t 0) and λ ( t 0 ) = r o ( t 0 ) R o ( t 0 ) , note that α(t 0) = 1. Iterate steps 2-5 until the final step N.

2. According to discretized equation (14) using Euler method, perform the following compute:

α ( n + 1 ) = α ( n ) + λ D ( n ) - 1 τ R e x p ( - n h / τ f ) h ,
(20)
R o ( n + 1 ) = α ( n + 1 ) R o ( t 0 ) .
(21)

3. Utilize discretized version of Eq. (15) taking account of Eq. (14) to calculate the R i at the next time step.

R i ( n + 1 ) = R i ( n ) + [ R o ( n ) R i ( n ) ] 2 [ λ D ( n ) - 1 ] e x p ( - n h / τ f ) h τ R R o ( t 0 ) - b ( n ) R i 2 ( n ) h .
(22)

4. Update values for c(n) and d(n) using their definitions by substituting the newly derived R i and R o .

c ( n + 1 ) = { 1 - P ( n + 1 ) E ( n + 1 ) 1 R i o ( n + 1 ) 3 - 1 [ ( 1 + ν ) R i o ( n + 1 ) 2 + R i o ( n + 1 ) 2 ( 1 - 2 ν ) ] } [ R o ( n + 1 ) - R i ( n + 1 ) ] ,
(23)
d ( n + 1 ) = R o 3 ( n + 1 ) - R i 3 ( n + 1 ) .
(24)

5. Compute new values for r i , r o and λ

r i ( n + 1 ) = - c ( n + 1 ) + [ 4 d ( n + 1 ) / c ( n + 1 ) - c 2 ( n + 1 ) ] / 3 2 ,
(25)
r o ( n + 1 ) = c ( n + 1 ) + [ 4 d ( n + 1 ) / c ( n + 1 ) - c 2 ( n + 1 ) ] / 3 2 ,
(26)
λ ( n + 1 ) = r o ( n + 1 ) R o ( n + 1 ) .
(27)

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.

Results

Changes of EDD, wall thickness and stress

The numerical simulation on LV dimensional changes including EDD and wall thickness were shown in Figures 4. The interior radius increased from 2.00 mm to about 2.05 mm (see Figure 4B), which is consistent with the 2.60% increase in our mice experiments for EDD from young to senescence. None of the percentage difference between strainless and deformed radii exceeds 0.1% for all time of concern, indicating a proper adoption of small deformation Hook's law (see Figure 4C). The wall thickness started from 0.92 mm, ended with 0.98 mm and had an increase of 0.053 mm, or a 5.72% ascent (Figure 4D), which is of the same magnitude as the 13.20% increase in our experiments.

Figure 4
figure 4

Temporal curves of LV dimensions and stresses derived by solving the simultaneous equations. The parameters in Eq. (14) adopted are τ R = 0 . 22 mons, D = 2 . 5 and τ f = 5 . 6 mons. (A) Temporal profile of the strainless radius on the external surface of LV; (B) The internal radius change with age; (C) Percentage change between deformed and free states. They are almost identical; (D) Wall thickness that has a 5.72% thickening from young to old mice; (E) Absolute value of negative spatial mean strain in radial direction; (F) Circumferential (or meridional) stress averaged across the wall.

The temporal profiles of the mean σ RR and mean σ θθ (or σ φφ ) across the wall were shown in Figures 4E and 4F calculating by

< σ R R > = 1 R o - R i R i R o P R i 3 R o 3 - R i 3 ( 1 - R o 3 R 3 ) d R = - R i ( R o + 2 R i ) 2 ( R o 2 + R o R i + R i 2 ) P
(28)
< σ θ θ > = < σ ϕ ϕ > = 1 R o - R i R i R o P R i 3 R o 3 - R i 3 ( 1 + R o 3 2 R 3 ) d R = R i ( R o 2 + R o R i + 4 R i 2 ) 4 ( R o 3 - R i 3 ) P
(29)

As is seen in Eq. (2), σ R R 0 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).

Pressure-volume relationship with aging

End diastolic portion of pressure-volume loop in a cardiac cycle provides important information of passive filling for LV and therefore describes passive properties of myocardium. For example, the slope of end diastolic pressure-volume relationship (EDPVR) gives the measure of LV stiffness. We plotted end diastolic pressure-volume curve with age for C57/BL6J mice (see Figure 5), in which EDV= ( 4 / 3 ) π r i 3 is the volume of the cavity. The arrow indicates the aging direction. The ratio of changes between pressure and end diastolic volume goes up with age, together with increased pressure and volume with age. It is worth noting that the end diastolic P-V relation shown in Figure 5 is a temporal progression with age, which is different from the typical one cardiac cycle P-V loop.

Figure 5
figure 5

End diastolic pressure-volume relationship with age. The arrow implies aging direction. Both pressure and end diastolic volume go up with age.

Analysis of constraints of remodeling equation

We also analyzed the candidates of the tissue growth function f(t) based on assumption 2). From Eq. (12), we had: R i 2 ( t ) i ( t ) = R o 2 ( t ) o ( t ) -b ( t ) . So i ( t ) 0 is equivalent as o ( t ) b ( t ) R o - 2 ( t ) . Apply dynamic equation, i.e. equation (14), o ( t ) = R o ( t 0 ) α ˙ ( t ) , the inequality became: f ( t ) b ( t ) R o ( t 0 ) R o 2 ( t ) . On the other hand, wall thickness growth also gave o ( t ) i ( t ) . Since i ( t ) = R o 2 ( t ) o ( t ) - b ( t ) R i 2 ( t ) , o ( t ) R o 2 ( t ) o ( t ) - b ( t ) R i 2 ( t ) , i.e. o ( t ) b ( t ) R o 2 - R i 2 , substituting R o with α and using equation for α ˙ , the inequality became f ( t ) b ( t ) R o ( t 0 ) [ R o 2 ( t ) - R i 2 ( t ) ] , which led to

b ( t ) R o ( t 0 ) R o 2 ( t ) f ( t ) b ( t ) R o ( t 0 ) [ R o 2 ( t ) - R i 2 ( t ) ] ,
(30)

where b(t) was a universal input function in our system and did not depend on the selection of parameters. R o (t0), i.e. initial value of R o , was a constant. However, R o 2 ( t ) and R o 2 ( t ) - R i 2 ( t ) could only be determined after the simultaneous equations being solved. Therefore, these two terms were subject to the parameters adopted in the mathematical model. However, we could find universal upper and lower bounds for the growth rate function, f(t). As we know, R o (t) increases with aging. Hence, b ( t ) R o 3 ( t 0 ) is a natural high upper bound of f(t). By defining S ( t ) = R o 2 ( t ) - R i 2 ( t ) , it is easy to see S(t) went up with time because both inner, outer radii and wall thickness enlarge. So S(t) was maximized at 30 months which was the terminal simulation time span. To achieve the maximum surface area difference for two concentric spheres kept at a constant volume, the smaller the inner/outer radius, the bigger the value S is. Since both interior and exterior radii increase, the lowest possible radius at the internal surface was R i (t = 30 mons) = R i (t0). Use the mass of each component at the final time and corresponding density, the final volume was determined and Parameter R o (t = 30 mon) was calculated, which is the exterior radius in month 30 assuming the interior radius is the same as that in month 7.5, i.e. its initial value. Hence, b ( t ) R o ( t 0 ) [ R o 2 ( t = 30 o n s ) - R i 2 ( t = 30 m o n s ) ] offers a lower upper bound. Therefore, if the following inequalities are satisfied, the selection of the radius growth rate function f(t) should give us proper remodeling trend,

b ( t ) R o 3 ( t 0 ) f ( t ) b ( t ) R o ( t 0 ) [ R o 2 ( t = 30 m o n s ) - R i 2 ( t = 30 m o n s ) ] .
(31)

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.

As is shown in Figure 6, the green line with circle markers and red line with hexagram markers are the universal boundaries determined by equation (31). The blue curve with star markers and the black one with square markers are true upper and lower bounds respectively for f(t) calculated in equations (30) for the set of parameters listed above.

Figure 6
figure 6

Boundaries for the LV growth rate function f t to generate growing trends for internal radius and wall thickness. The dotted cyan line in the middle is our f t for τ R = 0 . 22 mons, D = 2 . 5 and τ f = 5 . 6 mons. The star-marked blue line and square-marked black line form boundaries for our selection of parameters. Between the green curve with circle markers and red one with hexagram markers is the region satisfying conditions determined by equations (30-31).

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.

Discussion

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 [25]. 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 [26], 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 [27]. 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.

Conclusions

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.

References

  1. Roger VL, et al: Heart disease and stroke statistics--2011 update. Circulation. 2011, 123 (4): e18-e209. 10.1161/CIR.0b013e3182009701.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Howden LM, Meyer JA: Age and sex composition: 2010. 2010,

    Google Scholar 

  3. Blackman DJ, et al: Revascularisation for acute coronary syndromes in older people. Age Ageing. 2003, 32 (2): 129-135. 10.1093/ageing/32.2.129.

    Article  PubMed  Google Scholar 

  4. Ruan Q, Nagueh SF: Effect of age on left ventricular systolic function in humans: a study of systolic isovolumic acceleration rate. Exp Physiol. 2005, 90 (4): 527-34. 10.1113/expphysiol.2005.030007.

    Article  PubMed  Google Scholar 

  5. Watanabe S, et al: Influence of aging on cardiac function examined by echocardiography. Tohoku J Exp Med. 2005, 207 (1): 13-9. 10.1620/tjem.207.13.

    Article  PubMed  Google Scholar 

  6. Gates PE, et al: Left ventricular structure and diastolic function with human ageing. Relation to habitual exercise and arterial stiffness. Eur Heart J. 2003, 24 (24): 2213-20. 10.1016/j.ehj.2003.09.026.

    Article  PubMed  Google Scholar 

  7. Lakatta EG: Cardiovascular aging in health. Clin Geriatr Med. 2000, 16 (3): 419-44. 10.1016/S0749-0690(05)70021-5.

    Article  CAS  PubMed  Google Scholar 

  8. Lakatta EG: Age-associated cardiovascular changes in health: impact on cardiovascular disease in older persons. Heart Fail Rev. 2002, 7 (1): 29-49. 10.1023/A:1013797722156.

    Article  PubMed  Google Scholar 

  9. Chiao YA, et al: Multi-Analyte Profiling Reveals Matrix Metalloproteinase-9 and Monocyte Chemotactic Protein-1 as Plasma Biomarkers of Cardiac Aging/Clinical Perspective. Circulation: Cardiovascular Genetics. 2011, 4 (4): 455-462. 10.1161/CIRCGENETICS.111.959981.

    PubMed Central  CAS  Google Scholar 

  10. Conrad CH, et al: Myocardial Fibrosis and Stiffness With Hypertrophy and Heart Failure in the Spontaneously Hypertensive Rat. Circulation. 1995, 91 (1): 161-170. 10.1161/01.CIR.91.1.161.

    Article  CAS  PubMed  Google Scholar 

  11. Lin J, et al: Age-related cardiac muscle sarcopenia: Combining experimental and mathematical modeling to identify mechanisms. Exp Gerontol. 2008, 43: 296-306. 10.1016/j.exger.2007.12.005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Swartz M, et al: Mechanical stress is communicated between different cell types to elicit matrix remodeling. PNAS. 2001, 98 (11): 6180-6185. 10.1073/pnas.111133298.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Zhong L, Ghista DN, Tan RS: Left ventricular wall stress compendium. Comput Methods Biomech Biomed Engin. 2011, 1-

    Google Scholar 

  14. Han HC, et al: A left ventricle model to predict post-revascularization ejection fraction based on cine magnetic resonance images. J Biomech Eng. 2002, 124 (1): 52-5. 10.1115/1.1428555.

    Article  PubMed  Google Scholar 

  15. Rachev A, Stergiopulos N, Meister JJ: A model for geometric and mechanical adaptation of arteries to sustained hypertension. J Biomech Eng. 1998, 120 (1): 9-17. 10.1115/1.2834313.

    Article  CAS  PubMed  Google Scholar 

  16. Rachev A: Theoretical study of the effect of stress-dependent remodeling on arterial geometry under hypertensive conditions. J Biomech. 1997, 30 (8): 819-27. 10.1016/S0021-9290(97)00032-8.

    Article  CAS  PubMed  Google Scholar 

  17. Wong AY, Rautaharju PM: Stress distribution within the left ventricular wall approximated as a thick ellipsoidal shell. Am Heart J. 1968, 75 (5): 649-62. 10.1016/0002-8703(68)90325-6.

    Article  CAS  PubMed  Google Scholar 

  18. Greaves GN, et al: Poisson's ratio and modern materials. Nat Mater. 2011, 10 (11): 823-37. 10.1038/nmat3134.

    Article  CAS  PubMed  Google Scholar 

  19. Batkai S, et al: Decreased age-related cardiac dysfunction, myocardial nitrative stress, inflammatory gene expression, and apoptosis in mice lacking fatty acid amide hydrolase. Am J Physiol Heart Circ Physiol. 2007, 293 (2): H909-18. 10.1152/ajpheart.00373.2007.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Hori M, Nemat-Nasser S: On two micromechanics theories for determining micro-macro relations in heterogeneous solids. Mechanics of Materials. 1999, 31 (10): 667-682. 10.1016/S0167-6636(99)00020-4.

    Article  Google Scholar 

  21. Jin YF, et al: Combining experimental and mathematical modeling to reveal mechanisms of macrophage-dependent left ventricular remodeling. BMC Syst Biol. 2011, 5: 60-10.1186/1752-0509-5-60.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Martin RB, Boardman DL: The effects of collagen fiber orientation, porosity, density, and mineralization on bovine cortical bone bending properties. J Biomech. 1993, 26 (9): 1047-54. 10.1016/S0021-9290(05)80004-1.

    Article  CAS  PubMed  Google Scholar 

  23. Mendez J, Keys A: Density and composition of mammalian muscle. Metabolism. 1960, 9: 184-188.

    CAS  Google Scholar 

  24. Vinnakota KC, Bassingthwaighte JB: Myocardial density and composition: a basis for calculating intracellular metabolite concentrations. American Journal of Physiology - Heart and Circulatory Physiology. 2004, 286 (5): H1742-H1749. 10.1152/ajpheart.00478.2003.

    Article  CAS  PubMed  Google Scholar 

  25. Pacher P, et al: Left ventricular pressure-volume relationship in a rat model of advanced aging-associated heart failure. Am J Physiol Heart Circ Physiol. 2004, 287 (5): H2132-7. 10.1152/ajpheart.00405.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Yang B, Larson DF, Watson R: Age-related left ventricular function in the mouse: analysis based on in vivo pressure-volume relationships. Am J Physiol. 1999, 277 (5 Pt 2): H1906-13.

    CAS  PubMed  Google Scholar 

  27. Lindsey ML, et al: Age-dependent changes in myocardial matrix metalloproteinase/tissue inhibitor of metalloproteinase profiles and fibroblast function. Cardiovasc Res. 2005, 66 (2): 410-9. 10.1016/j.cardiores.2004.11.029.

    Article  CAS  PubMed  Google Scholar 

  28. Ghista DN, Advani SH, Rao BN: In vivo elastic modulus of the left ventricle: its determination by means of a left ventricular vibrational model and its physiological significance and clinical utility. Med Biol Eng. 1975, 13 (2): 162-70. 10.1007/BF02477723.

    Article  CAS  PubMed  Google Scholar 

  29. Graham JS, et al: Structural changes in human type I collagen fibrils investigated by force spectroscopy. Exp Cell Res. 2004, 299 (2): 335-42. 10.1016/j.yexcr.2004.05.022.

    Article  CAS  PubMed  Google Scholar 

  30. Yang L, et al: Mechanical properties of native and cross-linked type I collagen fibrils. Biophys J. 2008, 94 (6): 2204-11. 10.1529/biophysj.107.111013.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Vannah WM, Childress DS: Modelling the mechanics of narrowly contained soft tissues: the effects of specification of Poisson's ratio. J Rehabil Res Dev. 1993, 30 (2): 205-9.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yu-Fang Jin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Y.F.J and M.L.L designed the research; Y.A.C, A.V, H.C.H, and M.L.L performed the experiment, Y.F.J, T.Y, and Y.W performed the computational analysis and simulation. T.Y, H.C.H, M.L.L, and Y.F.J analyzed the results and wrote the manuscript.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Yang, T., Chiao, Y.A., Wang, Y. et al. Mathematical modeling of left ventricular dimensional changes in mice during aging. BMC Syst Biol 6 (Suppl 3), S10 (2012). https://doi.org/10.1186/1752-0509-6-S3-S10

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1752-0509-6-S3-S10

Keywords