Volume 6 Supplement 3

## The International Conference on Intelligent Biology and Medicine (ICIBM): Systems Biology

# Mathematical modeling of left ventricular dimensional changes in mice during aging

- Tianyi Yang
^{1, 2}, - Ying Ann Chiao
^{1, 3, 4}, - Yunji Wang
^{2}, - Andrew Voorhees
^{1, 5}, - Hai-Chao Han
^{1, 5}, - Merry L Lindsey
^{1, 3}and - Yu-Fang Jin
^{1, 2}Email author

**6(Suppl 3)**:S10

**DOI: **10.1186/1752-0509-6-S3-S10

© Yang et al.; licensee BioMed Central Ltd. 2012

**Published: **17 December 2012

## 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 [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 [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 [6–8]. 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

*∈*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*

_{ i },

*R*

_{ o }]. Applying equations (2) and (3) to equation (1), we had

*R*

_{ i }to

*R*

_{ o }, we obtained

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

*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*

_{ 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:

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

*π*, the following equation about free radii was obtained:

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

*r*

_{ o }>

*R*

_{ o }and

*r*

_{ i }

*> R*

_{ i }. Additionally, the ratio between deformed and free exterior radii $\lambda \left(t\right):=\frac{{r}_{o}\left(t\right)}{{R}_{o}\left(t\right)}$ 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,

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 ${\lambda}^{D}\left(t\right)-1=\frac{r{o}^{D}\left(t\right)-R{o}^{D}\left(t\right)}{R{o}^{D}\left(t\right)}$ 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}\left(t\right)=1.21+1.50\left\{1-exp[-\left(t-7.50\right)\mathsf{\text{/}}{\tau}_{c}\right\}$, with ${\tau}_{c}=7.23$ months and other units are the same to relevant ones in $M\left(t\right)$. 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}\left(t\right)=M\left(t\right)-{M}_{c}\left(t\right)$. 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.

Constant parameters and initial values adopted in numerical calculation and the source of the data.

Parameters (unit) | Values | References |
---|---|---|

Initial value | 2.0 | Our experiments on over 148 mice |

Muscle density | 1.06 | |

Collagen density | 1.70 | [22] |

Muscle elastic modulus | 5.0×10 | [28] |

Collagen elastic modulus | 3.2×10 | |

Poisson' ratio | 0.47 | [31] |

Time step | 3 |

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

*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:

The radii ${r}_{i}\left(t\right)$ and ${r}_{o}\left(t\right)$ were real positive numbers subject to the condition *d*(*t*) > *c*^{3}(*t*).

We studied the variation of free and deformed radii in the time domain [*t*_{0}, *t*_{0} + *T*] and evenly divided this region T by a total of N small constant intervals of *h*. Equivalently, the domain could be rewritten as [*t*_{0}, *t*_{0} + *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 $\lambda \left({t}_{0}\right)=\frac{{r}_{o}\left({t}_{0}\right)}{{R}_{o}\left({t}_{0}\right)}$, note that *α*(*t* _{0}) = 1. Iterate steps 2-5 until the final step *N*.

*R*

_{ i }at the next time step.

*c*(

*n*) and

*d*(

*n*) using their definitions by substituting the newly derived

*R*

_{ i }and

*R*

_{ o }.

*r*

_{ i },

*r*

_{ o }and λ

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

*σ*

_{ RR }and mean

*σ*

_{ θθ }(or

*σ*

_{ φφ }) across the wall were shown in Figures 4E and 4F calculating by

As is seen in Eq. (2), ${\sigma}_{RR}\le 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

### Analysis of constraints of remodeling equation

*R*

_{ o }with

*α*and using equation for $\dot{\alpha}$, the inequality became $f\left(t\right)\le \frac{b\left(t\right)}{{R}_{o}\left({t}_{0}\right)\left[{R}_{\mathsf{\text{o}}}^{2}\left(t\right)-{R}_{i}^{2}\left(t\right)\right]}$, which led to

*b*(

*t*) was a universal input function in our system and did not depend on the selection of parameters.

*R*

_{ o }(

*t*

_{0}), i.e. initial value of

*R*

_{ o }, was a constant. However, ${R}_{o}^{2}\left(t\right)$ and ${R}_{o}^{2}\left(t\right)-{R}_{i}^{2}\left(t\right)$ 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, $\frac{b\left(t\right)}{{R}_{o}^{3}\left({t}_{0}\right)}$ is a natural high upper bound of

*f*(

*t*). By defining $S\left(t\right)={R}_{o}^{2}\left(t\right)-{R}_{i}^{2}\left(t\right)$, 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 }(

*t*

_{0}). 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, $\frac{b\left(t\right)}{{R}_{o}\left({t}_{0}\right)\left[{R}_{o}^{2}\left(t=30\u1e41ons\right)-{R}_{i}^{2}\left(t=30\phantom{\rule{0.3em}{0ex}}mons\right)\right]}$ 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,

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.

*f*(

*t*) calculated in equations (30) for the set of parameters listed above.

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.

## Declarations

### 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.

## Authors’ Affiliations

## References

- Roger VL, et al: Heart disease and stroke statistics--2011 update. Circulation. 2011, 123 (4): e18-e209. 10.1161/CIR.0b013e3182009701.PubMed CentralView ArticlePubMed
- Howden LM, Meyer JA: Age and sex composition: 2010. 2010,
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- Lakatta EG: Cardiovascular aging in health. Clin Geriatr Med. 2000, 16 (3): 419-44. 10.1016/S0749-0690(05)70021-5.View ArticlePubMed
- 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.View ArticlePubMed
- 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
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- Zhong L, Ghista DN, Tan RS: Left ventricular wall stress compendium. Comput Methods Biomech Biomed Engin. 2011, 1-
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- Greaves GN, et al: Poisson's ratio and modern materials. Nat Mater. 2011, 10 (11): 823-37. 10.1038/nmat3134.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View Article
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- Mendez J, Keys A: Density and composition of mammalian muscle. Metabolism. 1960, 9: 184-188.
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.PubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.PubMed

## Copyright

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.