# Combining experimental and mathematical modeling to reveal mechanisms of macrophage-dependent left ventricular remodeling

- Yu-Fang Jin
^{1}Email author, - Hai-Chao Han
^{2}, - Jamie Berger
^{3}, - Qiuxia Dai
^{3}and - Merry L Lindsey
^{3}

**5**:60

**DOI: **10.1186/1752-0509-5-60

© Jin et al; licensee BioMed Central Ltd. 2011

**Received: **19 March 2010

**Accepted: **5 May 2011

**Published: **5 May 2011

## Abstract

### Background

Progressive remodeling of the left ventricle (LV) following myocardial infarction (MI) can lead to congestive heart failure, but the underlying initiation factors remain poorly defined. The objective of this study, accordingly, was to determine the key factors and elucidate the regulatory mechanisms of LV remodeling using integrated computational and experimental approaches.

### Results

By examining the extracellular matrix (ECM) gene expression and plasma analyte levels in C57/BL6J mice LV post-MI and ECM gene responses to transforming growth factor (TGF-β_{1}) in cultured cardiac fibroblasts, we found that key factors in LV remodeling included macrophages, fibroblasts, transforming growth factor-β_{1}, matrix metalloproteinase-9 (MMP-9), and specific collagen subtypes. We established a mathematical model to study LV remodeling post-MI by quantifying the dynamic balance between ECM construction and destruction. The mathematical model incorporated the key factors and demonstrated that TGF-β_{1} stimuli and MMP-9 interventions with different strengths and intervention times lead to different LV remodeling outcomes. The predictions of the mathematical model fell within the range of experimental measurements for these interventions, providing validation for the model.

### Conclusions

In conclusion, our results demonstrated that the balance between ECM synthesis and degradation, controlled by interactions of specific key factors, determines the LV remodeling outcomes. Our mathematical model, based on the balance between ECM construction and destruction, provides a useful tool for studying the regulatory mechanisms and for predicting LV remodeling outcomes.

## Background

Myocardial infarction (MI) is a leading cause of congestive heart failure (CHF) [1, 2]. In response to the MI stimulus, the left ventricle (LV) undergoes structural and functional adaptations that collectively have been termed as LV remodeling [3]. Adverse LV remodeling progresses to CHF in about 25% of post-MI patients, but the mechanisms that drive this progression remain poorly understood. During LV remodeling, both extracellular matrix (ECM) degradation and synthesis increase [4]. When ECM degradation dominates over synthesis, LV rupture can occur. When ECM synthesis dominates over degradation rates, fibrosis can occur. Fibrosis increases myocardial stiffness and further depresses LV function to culminate in CHF [5, 6]. Therefore, understanding what regulates the balance between ECM degradation and synthesis post-MI is critical to understand the mechanisms of LV remodeling and may allow us to target specific early diagnostic indicators to better guide treatment protocols.

Previous studies have shown that matrix metalloproteinases (MMPs) regulate ECM degradation and fibroblasts regulate ECM synthesis [7–9]. MMP-9, transforming growth factor-β_{1} (TGF-β_{1}), tissue inhibitor of metalloproteinase-1 (TIMP-1), and collagen I levels are significantly elevated from day 1 to day 7 post-MI [4, 10, 11]. These increases are concomitant with increased infiltration of macrophages and activation of fibroblasts [12]. LV remodeling is a complex process that involves the spatiotemporal interactions among many biological components that remains poorly understood, in part due to the lack of complete sets of experimental data and computational models. Therefore, the objectives of this study were to 1) identify candidate biomarkers of LV remodeling post-MI from ECM gene expression and plasma analyte analyses, and 2) establish a mathematical model that incorporates experimental results to predict LV remodeling outcomes following different interventions. This model would provide a tool to elucidate LV regulatory mechanisms, estimate un-measurable variables, and predict outcomes following multiple therapeutic scenarios.

## Results

### Identifying Key Factors

ECM Gene Array (Data are Mean ± SD normalized levels.)

Gene | Control | Remote | Infarct |
---|---|---|---|

Cdh3 | 1.0E-05 ± 9.8E-06 | 7.4E-05 ± 5.5E-05 | 2.8E-04 ± 1.3E-04 |

Col1a1 | 1.5E-02 ± 2.6E-03 | 1.9E-01 ± 1.0E-01 | 8.2E-01 ± 3.4E-01 |

Col2a1 | 3.1E-06 ± 6.3E-07 | 1.1E-05 ± 1.1E-05 | 2.3E-04 ± 1.8E-04 |

Col3a1 | 5.0E-03 ± 1.5E-03 | 8.0E-02 ± 5.5E-02 | 2.3E-01 ± 9.7E-02 |

Col5a1 | 8.0E-3 ± 2.0E-3 | 2.3E-2 ± 1.1E-2 | 4.1E-2 ± 1.6E-2 |

Ctgf | 2.5E-02 ± 4.3E-03 | 7.6E-02 ± 5.0E-02 | 1.7E-01 ± 9.1E-02 |

Fn1 | 9.0E-03 ± 2.5E-03 | 4.3E-02 ± 2.7E-02 | 1.6E-01 ± 6.9E-02 |

Mmp2 | 3.5E-03 ± 1.2E-03 | 1.2E-02 ± 7.9E-03 | 1.8E-02 ± 7.5E-03 |

Mmp14 | 1.9E-03 ± 4.6E-04 | 6.8E-03 ± 3.2E-03 | 1.5E-02 ± 6.8E-03 |

Ncam1 | 6.2E-04 ± 1.3E-04 | 1.8E-03 ± 1.3E-03 | 4.1E-03 ± 1.7E-03 |

Postn | 5.3E-03 ± 1.6E-03 | 1.5E-01 ± 9.2E-02 | 2.9E-01 ± 1.2E-01 |

Sparc | 1.3E-02 ± 3.7E-03 | 7.0E-02 ± 3.7E-02 | 1.1E-01 ± 4.5E-02 |

Spp1 | 1.3E-05 ± 5.7E-06 | 3.1E-04 ± 4.5E-04 | 2.6E-03 ± 2.2E-03 |

Tgfbi | 2.0E-03 ± 3.3E-04 | 3.3E-03 ± 1.2E-03 | 7.1E-03 ± 3.6E-03 |

Thbs1 | 7.2E-03 ± 6.7E-03 | 1.2E-02 ± 7.6E-03 | 4.0E-02 ± 1.6E-02 |

Thbs2 | 9.0E-03 ± 2.3E-03 | 1.1E-02 ± 3.6E-03 | 3.9E-02 ± 1.4E-02 |

Timp1 | 2.0E-04 ± 9.4E-05 | 4.1E-03 ± 3.1E-03 | 6.1E-03 ± 3.3E-03 |

Multi-analyte Profiling of Control and 7 day Post-MI Plasma

Control n = 6 | 7 d MI n = 7 | P value | |
---|---|---|---|

Clusterin (μg/mL) | 330 ± 40 | 510 ± 200 | 0.046 |

Cystatin-C (ng/mL) | 360 ± 20 | 530 ± 100 | 0.008 |

Eotaxin (pg/mL) | 1500 ± 200 | 1900 ± 400 | 0.049 |

Fibrinogen (mg/mL) | 12 ± 2 | 19 ± 5 | 0.012 |

Haptoglobin (μg/mL) | 83 ± 30 | 180 ± 50 | 0.002 |

Macrophage Inflammatory Protein-1 gamma (ng/mL) | 24 ± 4 | 33 ± 9 | 0.042 |

Matrix Metalloproteinase-9 (ng/mL) | 71 ± 9 | 96 ± 20 | 0.009 |

Myeloperoxidase (ng/(mL) | 58 ± 10 | 100 ± 30 | 0.006 |

Osteopontin (ng/mL) | 250 ± 60 | 390 ± 100 | 0.041 |

Serum Amyloid Protein (μg/mL) | 21 ± 3 | 33 ± 9 | 0.013 |

TIMP-1 (ng/mL) | 0.8 ± 0.1 | 3.7 ± 2.0 | 0.002 |

_{1}. Fibroblast ECM array analysis showed that TGF-β

_{1}stimulation of cardiac fibroblasts up regulated 5 genes and down regulated 7 genes, which are shown in Table 3. Collagen I expression and SPARC expression levels were doubled, and TIMP-1 expression increased 5-fold in response to TGF-β

_{1}stimulation at the concentration of 10 ng/mL. These experimental data indicated a primary regulatory effect of TGF-β

_{1}on fibroblast ECM production. Interestingly, the 5 up regulated genes are among the 17 ECM genes that were significantly expressed in the LV, indicating that the cardiac fibroblast is likely the major tissue source for these ECM genes (Table 1).

Fibroblast ECM array in serum free control and 10 ng/ml TGF-β stimulated fibroblasts

Serum Free | TGF-β | p Value | |
---|---|---|---|

ECM/Growth Factors | |||

Col1a1 | 4.842 ± 1.399 | 9.614 ± 3.324 | 0.028 |

Col5a1 | 0.301 ± 0.085 | 0.605 ± 0.230 | 0.038 |

Fbln1 | 0.012 ± 0.005 | 0.008 ± 0.004 | 0.032 |

Sparc | 2.709 ± 0.204 | 4.932 ± 0.379 | <0.001 |

Tgfbi | 0.022 ± 0.016 | 0.013 ± 0.015 | 0.009 |

Cell Adhesion Molecules | |||

Ncam1 | 0.043 ± 0.016 | 0.169 ± 0.054 | 0.011 |

Pecam1 | 0.000060 ± 0.000012 | 0.000022 ± 0.000011 | 0.015 |

Sgce | 0.187 ± 0.035 | 0.146 ± 0.018 | 0.030 |

Vcam1 | 0.301 ± 0.136 | 0.161 ± 0.088 | 0.031 |

MMPs/TIMPs | |||

Mmp7 | 0.0000075 ± 0.0000011 | 0.0000056 ± 0.0000004 | 0.026 |

Timp1 | 0.052 ± 0.022 | 0.272 ± 0.136 | 0.035 |

Timp2 | 0.326 ± 0.109 | 0.188 ± 0.050 | 0.027 |

We also examined the correlations between LV wall thickness with 6 genes that were over expressed post-MI. The R^{2} values were 0.76 for collagen 1α1, 0.64 for collagen 2α1, 0.75 for collagen 5α1, 0.60 for periostin, 0.61 for osteopontin, and 0.63 for TGF-β_{1}.

_{1}, MMP-9, and collagen. Based on these experimental results, we developed a framework of the interaction loops among the identified key factors (Figure 1).

### Linking Experimental Results to Mathematical Modeling Framework

In the ECM construction pathway, collagen is secreted by fibroblasts. Growth and secretion of fibroblasts are stimulated by TGF-β_{1} (Table 3). In ECM destruction pathway, MMP-9 is a key factor by cleaving collagen. The major source of MMP-9 is the macrophage, which infiltrates into the infarct region post-MI. The majority of macrophages in the LV post-MI are differentiated from peripheral blood monocytes stimulated by chemoattractants including TGF-β_{1}[17]. The major source of TGF-β_{1} is the macrophage. Meanwhile, there are interactions between the ECM construction and destruction pathways: a) MMP-9 regulates ECM construction by activating TGF-β_{1} which stimulates collagen synthesis; b) TGF-β_{1} induces TIMP-1, which inhibits collagen degradation by blocking MMP-9 activity. Linking these key factors with their sources and effects allows us to develop the mathematical model for quantitative analysis.

### Mathematical Modeling

We established a set of nonlinear differential equations to model the temporal interactions among the key factors identified by our experimental results. The model incorporated the following variables: macrophage cell density (M_{Φ}, cells/mm^{3}), fibroblast cell density (F, cells/mm^{3}), collagen concentration (C, μg/μL), TGF-β_{1} concentration (T_{β}, pg/μL), and activated MMP-9 concentration (M_{9A}, pg/μL). Rates of cell number change were determined by the summation of constructive effects (migration rate or proliferation rate) and destructive effects (death rate or removal rate). Rates of chemical factors (TGF-β_{1}, MMP-9, collagen, etc) change were determined by the net difference between the synthesis rate and degradation rate.

Four assumptions were used based on experimental results: 1) All monocytes that migrate to the infarct region are differentiated to macrophages [18]; all activated macrophages are differentiated from peripheral blood monocytes since previous studies have shown that <5% of macrophages undergo mitotic division [19]; 2) The major source of fibroblasts is the resident cell and the contribution of circulating fibroblasts is ignored [20]; 3) The majority (80%) of TGF-β_{1} secreted at the injured site becomes activated [21, 22]; 4) ECM proteins in the infarct region are secreted by local cells. Our experimental data showed that ECM gene expression was higher in the infarct region than the expression in the remote region, and ECM gene levels in the remote region were higher than that in the LV of controls, indicating a local source of ECM production.

_{β}),

*F*

_{ g }(

*T*

_{ β }), and Fc(T

_{β}), were established based on

*in vivo*, and

*in vitro*experiments [17, 23–26]. We employed the function Fc(C) in a form given by Waugh and colleagues [27, 28]. The forms of these functions we took in this study were described in equation 8.

Pre-determined parameters in the mathematical model

Symbol | Biological meaning | Value | Units | Ref |
---|---|---|---|---|

d | Macrophage removal rate**(eqn 1) | 0.6 | day | [29] |

ρ | maximal macrophage density (eqn 2) | 2500 | cells/mm | [28] |

ρ | maximal fibroblast density(eqn 2) | 1250 | cells/mm | [28] |

ρ | maximal collagen density (eqn 2) | 3300 | μg/mm | [28] |

k | Fibroblast growth rate* (eqn 3) | 0.924 | day | [62] |

d | Fibroblast apoptosis rate (eqn 3) | 0.12 | day | [33] |

k | Macrophage TGF-β production rate (eqn 4) | 0.07 | pg/cell/day | [34] |

k | Fibroblast TGF-β production rate (eqn 4) | 0.004 | pg/cell/day | [35] |

d | TGF-β degradation rate | 15 | day | [36] |

K | Macrophage secretion MMP9 rate (eqn 5) | 3 | pg/cell/day | estimated |

d | MMP-9 degradation rate(eqn 5) | 0.875 | day | |

k | Kinetic reaction speed (eqn 5) | 3 × 10 | 1/(μg/mm | |

k | Kinetic reaction speed (eqn 5) | 4 × 10 | s | |

k | Kinetic reaction speed (eqn 6) | 0.004 | s | |

k | Fibroblast collagen production rate (eqn 6) | 20 | μg/cell/day | [41] |

Equation 1 determines the rate of macrophage (M_{Φ}) infiltration. The function M(T_{β}) describes the migration rate of macrophages to the scar tissue [17]. Since we assumed that all monocytes differentiated into macrophages (assumption 1), M(T_{β}) also represents the migration of monocytes stimulated by TGF-β_{1}. Parameter d_{MΦ} denotes the emigration rate of macrophages [29], since macrophages do not die locally in the scar tissue but emigrate to the lymph node system for disposal.

_{mc}= 0.05 represents the rate of myocyte cell death since ischemic myocytes undergo necrosis in the infarct region post-MI. The crowding effect of macrophages, fibroblasts, and collagen was considered by calculating the normalized density with respect to their maximum density in scar tissue [28]. Temporal profiles of total crowding effects, crowding effects contributed by macrophages, fibroblasts, and collagen, and the parameter

*k*

_{ mem }(

*t*) were shown in Figure 3.

Equation 3 determines the rate of fibroblast (F) density changes based on the assumption that majority of fibroblasts come from the proliferation of resident cells (assumption 2). The function *F*_{g}(*T*_{
β
}) denotes the stimulating effects of TGF-β_{1} on the growth rate of fibroblasts [17]. Parameter d_{F} represents the apoptosis rate of the fibroblast [33].

Equation 4 determines the rate of TGF-β_{1} concentration change, wherein k_{FT} denotes the TGF-β_{1} secretion rate of fibroblasts [34] and k_{MΦT} denotes the TGF-β_{1} secretion rate of macrophages [35] since a major source of TGF-β_{1} in the scar tissue is the activated macrophage. Parameter d_{Tβ} represents the degradation rate of TGF-β_{1}, which can be calculated from the half life data [36].

TGF-β_{1} gene levels demonstrated temporal progression at the early stage post-MI. Gene expression profile of TGF-β_{1} increased post-MI, peaked at day 2, and returned to normal levels after day 7 in mice post-MI [37]. Since majority of TGF-β_{1} secreted in the infarct is activated (Assumption 3), gene expression profile can be used as an activation pattern of TGF-β_{1}. The function, *ut* denotes temporal profile of TGF-β_{1} activation post-MI [38–40] and works as the input in the simulation. It's worth to mention that TGF-β1 activation levels were different with respect to different scar size. A temporal profile of *ut* representing small scar was shown in Figure 3.

Equation 5 determines the rate of activated MMP-9 concentration change. Proteolytic collagen degradation with activated MMP-9 is described in equation 5a, where M_{9A}, C, CM_{9}, and CID denote activated MMP-9, collagen, binding of MMP-9 and collagen, and degraded collagen peptide concentration, respectively. MMP-9 is inhibited primarily by TIMP-1, and TIMP-1 is induced by TGF-β_{1}. Thus, we established an inhibition function h(T) = 1/(1+T_{β}/T_{βN}) with T_{βN} = 6.0 pg/μL to represent the inhibition effect.

Equation 6 determines the rate of collagen concentration changes. Collagen secretion rate by fibroblasts was denoted by parameter*K*_{
FC
}. Meanwhile, the function, *F*_{
c
}(*T*_{
β
}), characterizes effects of TGF-β_{1} on collagen secretion rate by fibroblasts [23, 41]. Function, *F*_{
c
}(*c*), denotes the effect of collagen density on fibroblast secretion rate [27, 28].

Equation 7 determines the concentration change of CM_{9}, based on the theoretical model for collagen degradation by MMPs proposed by Popel's group [42, 43].

### Computational simulations

Computational simulations of scar formation (collagen deposition) were carried out by solving the nonlinear differential equations with MATLAB. Initial conditions of the fibroblast and macrophage densities were chosen as F(0) = 20 cells/mm^{3} and M_{Φ}(0) = 5 cells/mm^{3}. Accordingly, T_{β}(0) = 0.21 pg/μL, M_{9A}(0) = 7.1 pg/μL, C(0) = 839.5 μg/μL, CM_{9}(0) = 447.6 μg/μL were calculated by the equilibriums of equations 1-7 for normal LV. All the simulations shown in this study used the same initial conditions. The initial conditions were chosen based on measurements in the normal LV for both the control and MI groups (MI induced at day 0). The simulations covered the LV remodeling process from day 0 to day 30 post-MI.

### Model validation

_{1}stimuli corresponding to reduced, normal, and elevated post-MI activation strength. Others have reported an early increase of MMP-9 levels of 78 ± 19 pg/μL for small infarcts and 195 ± 63 pg/μL for large infarcts [45, 46]. These experimental results agreed with the MMP-9 predictions by our mathematical model, which peak at 80 pg/μL for small infarcts, and 220 pg/μL for large infarcts, respectively. MMP-9 profiles were shown in Figure 5 and details of the simulation setup for Figure 5 were explained as follows.

#### Effects of TGF-β_{1} levels

We employed activation of TGF-β_{1}, *u*_{
T
}in equation 4, at reduced, normal, and elevated post-MI strength. The activation function peaked at 15, 30, and 60 pg/μL/day, according to active expression levels observed in small, median, and very large infarcts, respectively. Temporal profiles of cell densities of macrophages, fibroblasts, concentrations of MMP-9, collagen, and TGF-β_{1} responses to the stimuli *u*_{
T
}are shown in Figure 5. For the TGF-β_{1} stimulus level setting at 30 pg/μL/day (a normal post-MI level in mice), TGF-β_{1} peaked after day 2 post-MI, macrophage density peaked at day 3, MMP-9 concentration peaked at day 4, and all were returned to normal levels at 30 days post-MI. In contrast, fibroblast density and collagen density continued to increase beginning at day 4, reached a stable value after day 20 post-MI, and then remained at a higher equilibrium level at day 30 post-MI (blue solid line in Figure 5). This computational predictions agreed with the experimental observations on the peak time of macrophages [9, 12], and progression trend of MMP-9 [10] and collagen for stable LV remodeling [47].

In the case of reduced TGF-β_{1} levels (with an activation peak at 15 pg/μL/day in Figure 3), less monocytes were attracted to the infarct. Thereby, less macrophages and MMP-9 appeared at the early stage. The reduced MMP-9 levels slowed down collagen degradation, leading to a higher collagen concentration (black dotted lines in Figure 5. This simulation agreed with the temporal profiles of macrophage and fibroblast density [9, 12], and relative low expression levels of MMP-9 for moderate size of infarcts compared to large infarcts [46, 48].

When the TGF-β_{1} stimulus strength was elevated to a level 2-fold higher than normally seen post-MI in mice, more macrophages infiltrated to the infracted region in the early days (Figure 5 red dash lines). Elevated macrophage infiltration led to high levels of crowding effect and high concentrations of MMP-9 for collagen degradation. Therefore, fibroblast growth was inhibited earlier due to the crowding effect, which led to less collagen secretion and a net negative collagen deposition. These results indicated an increased early susceptibility to LV rupture between day 2 to 6 post-MI and prolonged LV remodeling for those mice that did not rupture. Though we know that LV ruptures frequently occurred at day 2-4 post-MI in a normal remodeling process [49], this simulation has not been examined experimentally and raises some interesting hypotheses.

#### Effects of MMP-9 interventions on ECM destruction

## Discussion

This study is the first investigation to integrate mathematical modeling with ECM and fibroblast gene array data and plasma analytes to predict ECM remodeling post-MI. We have integrated *in vivo, in vitro*, and *in silico* approaches to dissect the complicated interactions among multiple regulatory factors in LV remodeling. As such, this study provides a promising modeling approach for elucidating the complicated LV remodeling process. The most important findings of this study are highlighted as follows. 1) The balance between ECM construction and destruction kinetics is likely the key determinant of scar formation, and interventions to either the construction or destruction side lead to different remodeling outcomes. 2) Dynamic interactions among key factors in LV remodeling determine LV remodeling outcomes post-MI. Altering initial conditions, intervention strengths, or intervention times have significant effects on LV remodeling outcomes, and these effects could be simulated accurately by our model. 3) Collagen, MMP-9, TGF-β_{1}, and TIMP-1 are critical biomarker candidates of LV remodeling outcomes.

Our experimental results on microarray and plasma data provided the foundations to build our computational framework. We examined 84 ECM genes and chose the 17 genes that were most highly over-expressed in the infarct region compared to both control and non-infarcted groups (>2.5-fold over-expression). The expression levels of several of these factors were further verified by our plasma data at protein level. One interesting finding was that SPP1 (osteopontin) gene expression levels increased 206-fold in the infarct region compared to the control group at day 7 post-MI, suggesting strong macrophage activation. The plasma protein levels of osteopontin increased from 250 ± 60 ng/mL in controls to 390 ± 100 ng/mL in post-MI samples, which adds support for the critical role of macrophages in our mathematical model. Therefore, the primary selection of the most highly changed genes allowed us to focus on the most significant factors at gene level and predict the possible interactions at protein level and cellular level.

An interesting observation was that MMP-9 mRNA levels did not increase in our gene array analysis, but MMP-9 protein levels increased in the plasma data. It is well known that MMP-9 protein levels and activation were increased at day 7 post-MI. This conflicting phenomena is caused by pre-formed MMP-9 proteins stored in leukocytes, which do not rely on increased gene expression to induce MMP-9 secretion and activation [50]. In addition, the shift from a normal ventricle composed of cardiac myocytes, endothelial cells, and smooth muscle cells to an infarcted ventricle comprised primarily of cardiac fibroblasts and inflammatory cells could result in a quantitative no gain in gene expression level but qualitative increase in MMP-9 function due to increased activation and increased substrate availability. Karl Weber's group showed that MMP-9 protein levels increased early post-MI, but the mRNA levels were not elevated [51]. We have also shown that MMP-9 activation increases upon reperfusion in a dog model of cardiac ischemia/reperfusion [7]. These results support our observations in this study and highlight the fact that MMP-9 is primarily regulated at the post-translational level. Furthermore, the fact that MMP-9 gene levels do not change at day 7 post-MI also suggest that MMP-9 should be evaluated at the protein level.

There is limited data available to construct suitable functions of *M*(*T*_{
β
}), *F*_{
g
}(*T*_{
β
}), *F*_{
c
}(*T*_{
β
}), and *F*_{
c
}(*c*). Wahl et al presented a set of *in vitro* data on monocyte chemotaxis induced by TGF-β with concentration ranges from 0.05 pg/mL to 10 pg/mL [17]. Their experiments also provided production of fibroblast growth activity induced with TGF-β concentration range from 10 pg/mL to 10 ng/mL. These data were used in this study to construct the functions of *M* (*T*_{
β
}) and *F*_{
g
}(*T*_{
β
}). Roberts et al have presented a set of data on *in vitro* collagen formation stimulated by TGF-β with the concentration range from 0.0825 pg/mL to 2.5 pg/mL [23]. Ignotz and Massague have examined the effects of TGF-β on collagenases in chick embryo fibroblasts stimulated by TGF-β at a concentration range of 1.25 ng/mL to 25 ng/mL [41]. We have done experiments to quantify the ECM production of cardiac fibroblast stimulated with TGF-β (10 ng/mL). These experimental data were used to build and justify the function, *F*_{
c
}(*T*_{
β
}). Loftis and colleagues have studied effects of collagen density on cardiac fibroblast behavior and showed elevated fibroblast activities stimulated with higher collagen concentrations (750 μg/mL - 1250 μg/mL) [52]. Wu and colleagues showed that propeptides at lower levels caused about 80% decrease in collagen synthesis compared to control [53]. In this study, we used an *in silico* function *F*_{
c
}(*c*) provided by Waugh and colleagues [27, 28]. These experiments were done by different groups with different conditions and the *in vitro* experiments might not reflect the *in vivo* interactions. In addition, we established these functions using polynomial interpolation. It's likely that these functions are not the optimal forms to describe the *in vivo* interactions; however, they provide us a baseline for further study.

Our computational results demonstrated that altering the strength of TGF-β_{1} altered LV remodeling outcomes. Elevated TGF-β_{1} levels at the early stage (day 3 post-MI) led to elevated macrophage density and MMP-9 levels, decreased fibroblast secretion of collagen and collagen deposition, and thereby, prolonged the progression of remodeling. Wetzler and colleagues have shown that the prolonged persistence of macrophages at the late phase (after day 7 post-MI) impairs the wound healing process [54]. It has also been shown that elevated TGF-β_{1} levels delays wound healing post-MI [55, 56]. Our simulation results are in agreement with these previously published studies.

Simulations of different TGF-β_{1} strengths also shed light on the regulation scheme of ECM construction and destruction. For ECM construction regulation, active TGF-β_{1} stimulated fibroblast proliferation and collagen secretion, which increased the crowding effect. The increased crowding coefficient ramped down fibroblast proliferation and TGF-β_{1} secretion (negative feedback), which slowed down the stimulus for monocytes to migrate into the infarct region. The decrease of monocytes number led to less macrophage infiltration, which then reduced the crowding coefficient. Meanwhile, reduced macrophages lead to less TGF-β_{1} secretion by macrophage, which further slowed down collagen synthesis. For collagen destruction regulation, TGF-β_{1} induced macrophage infiltration, which lead to elevated MMP-9 secretion, elevated collagen degradation, and thereby reduced crowding coefficients. Smaller crowding coefficients lead to elevated fibroblast proliferation and collagen secretion. Notably, there are two types of negative feedback schemes in the mathematical model: degradation (apoptosis or emigration) rates associated with proteins (cells) and the crowding effects. Degradation rates are constants and determine how fast the proteins (cells) can respond to stimuli. Crowding effects are time varying impacts imposed by the environment. Through these regulation schemes, a dynamic balance of collagen construction and destruction can be maintained to generate a stable scar. Furthermore, profiles of crowding effects elucidated the transition from the normal LV to scar tissue with respect to cell types and collagen concentrations (Figure 3).

It is worth mentioning that there exist biological negative feedbacks in our mathematical model. For macrophage density regulation, there was a positive feedback loop containing macrophage and TGF-β_{1}: TGF-β_{1} stimulated monocytes migration, leading to macrophage infiltration; macrophages secreted TGF-β_{1} which might attract more macrophages to the infarct site. We observed elevated macrophage density and MMP-9 concentrations corresponding to increased TGF-β_{1} levels (Figure 5). This positive feedback loop was inhibited by emigration of macrophages and degradation of TGF-β_{1}. In addition, Wahl et al pointed out that the monocyte chemotactic activity increased in response to low concentrations of TGF-β_{1} stimuli, while the chemotactic activity decreased in response to higher concentrations of TGF-β_{1}[17]. Therefore, as TGF-β_{1} levels continuously increased, infiltration speed of macrophages decreased as shown in figure 2 (M(T_{β})), suggesting a secondary biological inhibition scheme of the TGF-β_{1} -- macrophage positive feedback loop.

There are a few limitations of the mathematical model that resulted to a large degree from the model assumptions. More research is needed to address these limitations and further enhance the models. First, our model calls for accurate determination of MMP-9 activation and inhibition functions. Complete time-course measurements of TIMP-1 and the other three TIMPs (for MMP-9 inhibition) and MMP-3 (for MMP-9 activation) would provide additional details on the regulators of MMP-9 function. Second, the large differences between simulation results and experimental measurement of MMP-9 concentrations before day 4 post-MI in Figure 4 is likely due to the effect of neutrophils, an early source of MMP-9, indicating the need to investigate the role of neutrophils in the early stage (days 1-5) post-MI [7]. This inconsistency of the simulated MMP-9 concentrations and measurements from experimental results is one of the approaches for us to employ more biomarkers and key factors in model development. Third, Interleukin IL-1, IL-6, IL-10, and tumor necrosis factor-α profiles need to be included to better study macrophage activation phenotypes [57]. Fourth, the effects of other chemoattractants on monocyte migration on the LV remodeling need further investigation. The predicted fibroblast density at days 14-30 post-MI was lower than experimental measurements. This might be caused by assumption 2, that the increase of fibroblast numbers post-MI was contributed solely by proliferation of local cells. We will investigate the migration rates of fibroblasts in the future to refine our mathematical model. Regardless of the above limitations, we were able to use our model to compute the progression of macrophages, fibroblasts, and collagen density, MMP-9 and TGF-β_{1} concentrations post-MI.

This systems biology study for LV remodeling can be expanded to include proteomics and cardiac functions in future studies. We employed plasma data in this study since plasma proteins reflect the process of LV remodeling and plasma data are more directly translatable to the clinic. However, we are well aware that measuring tissue protein levels will provide a more direct evaluation of LV remodeling. Further investigation on the ECM proteomics in cardiac samples has been planned in our future research to establish a more complete mathematical model for LV remodeling. Though it is beyond the scope of the current paper, we have previously reported some data on cardiac function [58–60] and our future study will integrate cardiac function into the mathematical model.

The computational model for post-MI LV remodeling developed here illustrated the dynamic interactions among critical factors in LV remodeling. This is the first mathematical model focusing on the protein and cellular interactions post-MI. Thus, this model provides a strong foundation for future studies to build a more comprehensive model that takes into account a more complete set of parameters. The model also provides a tool to guide experimental designs by identifying candidate factors to intervene, the proper intervention time, and doses for effective interventions to achieve the most beneficial outcomes. As an example, we have shown the effects of MMP-9 intervention time and doses on LV remodeling outcomes in this study. Though this model was established based on *in vitro* data and *in vivo* data from mice, the modeling approach can be used to develop models for other scenarios such as the LV remodeling of human MI and LV remodeling under reperfusion conditions.

## Conclusions

In conclusion, we developed a set of differential equations to quantitatively model the dynamic interactions and temporal changes of the key components identified from our experimental results. Predictions of the mathematical model fell well within experimental measurements, particularly with regard to macrophage infiltration and matrix remodeling. This mathematical model provides a powerful tool to better understand how the dynamic balance between ECM construction and ECM destruction influences LV remodeling outcomes.

## Methods

### Mice

All animal procedures were conducted in accordance with the Guide for the Care and Use of Laboratory Animals (National Research Council, 1996) and were approved by the Institutional Animal Care and Use Committee at the University of Texas Health Science Center at San Antonio. Male C57BL/6J wild type adult mice (n = 13) at age 8.0 ± 0.5 months were used. One group (n = 6) served as unoperated controls, while the other group (n = 7) underwent coronary artery ligation for 7 days as described previously [9] and were sacrificed at day 7 post-myocardial infarction.

### In Vivo Procedures

Blood was collected from the jugular vein and placed in a heparinized tube for plasma collection. Tissue was collected for the gene array analysis as described previously[9].

### Microarray and Plasma Analysis

Total RNA was isolated using the TRIzol plus Total RNA purification kit (Invitrogen). The RT^{2} qPCR Primer Array for Extracellular Matrix and Adhesion Molecules (SuperArray catalog APMM-013A) was used for the gene array. Results were analyzed based on the ΔΔCt method with normalization of raw data to the GAPDH housekeeper gene. Data are presented as average 2^{-ΔCT} levels.

For the fibroblast ECM microarray analysis, cardiac fibroblasts were isolated from adult C57BL/6J mice and stimulated with or without 10 ng/ml TGF-β_{1} for 24 hours[61]. Plasma samples (80 μL) were analyzed for 67 antigens using the quantitative immunoassay panel for mice provided by Rules Based Medicine in Austin, TX (http://www.rulesbasedmedicine.com).

### Statistical Analysis

Control, remote, and infarct LV groups were analyzed by ANOVA, with the Bonferroni post hoc test. Control and MI plasma was analyzed by Students t-test. Unstimulated and TGF-β_{1} stimulated fibroblast groups were analyzed by paired t-test. A p < 0.05 was considered statistically significant.

## Declarations

### Acknowledgements

The authors acknowledge grants support from NIH 1R03EB009496, NIH SC2HL101430, NSF 0649172, and AT&T foundation (to YFJ), from NSF 0644646, and NSF 0602834 (to HCH), and NHLBI HHSN268201000036C (N01-HV-00244), NIH R01 HL75360, AHA Grant-in-Aid 0855119F, and the Morrison Fund (to MLL).

## Authors’ Affiliations

## References

- Hellermann JP, Jacobsen SJ, Redfield MM, Reeder GS, Weston SA, Roger VL: Heart failure after myocardial infarction: clinical presentation and survival. Eur J Heart Fail. 2005, 7: 119-125. 10.1016/j.ejheart.2004.04.011View ArticlePubMedGoogle Scholar
- Anavekar NS, McMurray JJV, Velazquez EJ, Solomon SD, Kober L, Rouleau JL, White HD, Nordlander R, Maggioni A, Dickstein K, et al.: Relation between Renal Dysfunction and Cardiovascular Outcomes after Myocardial Infarction. N Engl J Med. 2004, 351: 1285-1295. 10.1056/NEJMoa041365View ArticlePubMedGoogle Scholar
- Cohn JN, Ferrari R, Sharpe N: Cardiac remodeling--concepts and clinical implications: a consensus paper from an international forum on cardiac remodeling. Behalf of an International Forum on Cardiac Remodeling. J Am Coll Cardiol. 2000, 35: 569-582. 10.1016/S0735-1097(99)00630-0View ArticlePubMedGoogle Scholar
- Whittaker P: Collagen and ventricular remodeling after acute myocardial infarction: concepts and hypotheses. Basic Res Cardiol. 1997, 92: 79-81.View ArticlePubMedGoogle Scholar
- Zhang M, Shah MA: Role of reactive oxygen species in myocardial remodeling. Current Heart Failure Reports. 2007, 4: 26-30. 10.1007/s11897-007-0022-5View ArticlePubMedGoogle Scholar
- Rumberger JA: Ventricular Dilatation and Remodeling After Myocardial Infarction. Mayo Clin Proc. 1994, 69: 664-674.View ArticlePubMedGoogle Scholar
- Lindsey M, Wedin K, Brown MD, Keller C, Evans AJ, Smolen J, Burns AR, Rossen RD, Michael L, Entman M: Matrix-Dependent Mechanism of Neutrophil-Mediated Release and Activation of Matrix Metalloproteinase 9 in Myocardial Ischemia/Reperfusion. Circulation. 2001, 103: 2181-2187.View ArticlePubMedGoogle Scholar
- Rohde LE, Ducharme A, Arroyo LH, Aikawa M, Sukhova GH, Lopez-Anaya A, McClure KF, Mitchell PG, Libby P, Lee RT: Matrix metalloproteinase inhibition attenuates early left ventricular enlargement after experimental myocardial infarction in mice. Circulation. 1999, 15: 3063-3070.View ArticleGoogle Scholar
- Lindsey ML, Escobar GP, Dobrucki LW, Goshorn DK, Bouges S, Mingoia JT, McClister DM, Su H, Gannon J, MacGillivray C, et al.: Matrix metalloproteinase-9 gene deletion facilitates angiogenesis after myocardial infarction. Am J Physiol Heart Circ Physiol. 2006, 290: H232-239.View ArticlePubMedGoogle Scholar
- Webb CS, Bonnema DD, Ahmed SH, Leonardi AH, McClure CD, Clark LL, Stroud RE, Corn WC, Finklea L, Zile MR, Spinale FG: Specific Temporal Profile of Matrix Metalloproteinase Release Occurs in Patients After Myocardial Infarction: Relation to Left Ventricular Remodeling. Circulation. 2006, 114: 1020-1027. 10.1161/CIRCULATIONAHA.105.600353View ArticlePubMedGoogle Scholar
- Vanhoutte D, Schellings M, Pinto Y, Heymans S: Relevance of matrix metalloproteinases and their inhibitors after myocardial infarction: A temporal and spatial window. Cardiovasc Res. 2006, 69: 604-613. 10.1016/j.cardiores.2005.10.002View ArticlePubMedGoogle Scholar
- Yang F, Liu YH, Yang XP, Xu J, Kapke A, Carretero OA: Myocardial infarction and cardiac remodelling in mice. Exp Physiol. 2002, 87: 547-555. 10.1113/eph8702385View ArticlePubMedGoogle Scholar
- N Ashizawa KG, Do YS, Nunohiro T, Giachelli CM, Meehan WP, Tuan TL, Hsueh WA: Osteopontin is produced by rat cardiac fibroblasts and mediates A(II)-induced DNA synthesis and collagen gel contraction. J Clin Invest. 1996, 98: 2218-2227. 10.1172/JCI119031View ArticleGoogle Scholar
- Lenga Y, Koh A, Perera AS, McCulloch CA, Sodek J, Zohar R: Osteopontin Expression Is Required for Myofibroblast Differentiation. Circ Res. 2008, 102: 319-327. 10.1161/CIRCRESAHA.107.160408View ArticlePubMedGoogle Scholar
- Zahradka P: Novel Role for Osteopontin in Cardiac Fibrosis. Circ Res. 2008, 102: 270-272. 10.1161/CIRCRESAHA.107.170555View ArticlePubMedGoogle Scholar
- Borg TK, Markwald R: Periostin: More Than Just an Adhesion Molecule. Circ Res. 2007, 101: 230-231. 10.1161/CIRCRESAHA.107.159103View ArticlePubMedGoogle Scholar
- Wahl SM, Hunt DA, Wakefield LM, McCartney-Francis N, Wahl LM, Roberts AB, Sporn MB: Transforming growth factor type beta induces monocyte chemotaxis and growth factor production. Proceedings of the National Academy of Sciences of the United States of America. 1987, 84: 5788-5792. 10.1073/pnas.84.16.5788PubMed CentralView ArticlePubMedGoogle Scholar
- Krause SW, Rehli M, Kreutz M, Schwarzfischer L, Paulauskis JD, Andreesen R: Differential screening identifies genetic markers of monocyte to macrophage maturation. J Leuko Biol. 1996, 60: 510-545.Google Scholar
- Burke B LC: The Macrophage. 2002, Oxford: Oxford University Press, 2,Google Scholar
- Quan TE, Cowper S, Wu SP, Bockenstedt LK, Bucala R: Circulating fibrocytes: collagen-secreting cells of the peripheral blood. The International Journal of Biochemistry & Cell Biology. 2004, 36: 598-606. 10.1016/j.biocel.2003.10.005View ArticleGoogle Scholar
- Yehualaeshet T, O'Connor R, Green-Johnson J, Mai S, Silverstein Roy, Murphy-Ullrich JE, Khalil N: Activation of Rat Alveolar Macrophage-Derived Latent Transforming Growth Factor β-1 by Plasmin Requires Interaction with Thrombospondin-1 and its Cell Surface Receptor, CD36. Am J Pathol. 1999, 155: 841-851. 10.1016/S0002-9440(10)65183-8PubMed CentralView ArticlePubMedGoogle Scholar
- Khalil N, Whitman C, Zuo L, Danielpour D, Greenberg A: Regulation of alveolar macrophage transforming growth factor-beta secretion by corticosteroids in bleomycin-induced pulmonary inflammation in the rat. JClin Invest. 1993, 92: 1812-1818. 10.1172/JCI116771.View ArticleGoogle Scholar
- Roberts AB, Sporn MB, Assoian RK, Smith JM, Roche NS, Wakefield LM, Heine UI, Liotta LA, Falanga V, Kehrl JH: Transforming growth factor type beta: rapid induction of fibrosis and angiogenesis in vivo and stimulation of collagen formation in vitro. Proceedings of the National Academy of Sciences of the United States of America. 1986, 83: 4167-4171. 10.1073/pnas.83.12.4167PubMed CentralView ArticlePubMedGoogle Scholar
- Ignotz RA, Massagué J: Transforming growth factor-beta stimulates the expression of fibronectin and collagen and their incorporation into the extracellular matrix. Journal of Biological Chemistry. 1986, 261: 4337-4345.PubMedGoogle Scholar
- Leroy EC: Increased Collagen Synthesis by Scleroderma Skin Fibroblasts In Vitro: a possibledefect in the regulation or activation of the scleroderma fibroblast. J Clin Invest. 1974, 54: 880-889. 10.1172/JCI107827PubMed CentralView ArticlePubMedGoogle Scholar
- Naugle JE OE, Zhang X, Mase SE, Pilati CF, Maron MB, Folkesson HG, Horne WI, Doane KJ, Meszaros JG: Type VI collagen induces cardiac myofibroblast differentiation: implications for postinfarction remodeling. Am J Physiol Heart Circ Physiol. 2006, 290: H323-330.View ArticlePubMedGoogle Scholar
- Waugh H, Sherratt J: Macrophage Dynamics in Diabetic Wound Dealing. Bulletin of Mathematical Biology. 2006, 68: 197-207. 10.1007/s11538-005-9022-3View ArticlePubMedGoogle Scholar
- Helen V, Waugh JAS: Modeling the effects of treating diabetic wounds with engineered skin substitutes. Wound Repair and Regeneration. 2007, 15: 556-565. 10.1111/j.1524-475X.2007.00270.xView ArticleGoogle Scholar
- Bellingan GJ, Caldwell H, Howie SE, Dransfield I, Haslett C: In vivo fate of the inflammatory macrophage during the resolution of inflammation: inflammatory macrophages do not die locally, but emigrate to the draining lymph nodes. J Immunol. 1996, 157: 2577-2585.PubMedGoogle Scholar
- Olsen L, Sherratt JA, Maini PK: A mathematical model for fibro-proliferative wound healing disorders. Bull Math Biol. 1996, 58: 787-808. 10.1007/BF02459482View ArticlePubMedGoogle Scholar
- Munoz-Alicea RN-MP, Marcano-Velazquez M: A mathematical model for macrophage, T-cell and mycobacterium tuberculosis interactions. Book A mathematical model for macrophage, T-cell and mycobacterium tuberculosis interactions. 1999, Editor ed.^eds, City: University of Puerto Rico,Google Scholar
- Banerjee I, Fuseler JW, Price RL, Borg TK, Baudino TA: Determination of cell types and numbers during cardiac development in the neonatal and adult rat and mouse. Am J Physiol Heart Circ Physiol. 2007, 293: H1883-1891. 10.1152/ajpheart.00514.2007View ArticlePubMedGoogle Scholar
- Darby IA, Bisucci T, Hewitson TD, MacLellan DG: Apoptosis is increased in a model of diabetes-impaired wound healing in genetically diabetic mice. The International Journal of Biochemistry & Cell Biology. 1997, 29: 191-200. 10.1016/S1357-2725(96)00131-8View ArticleGoogle Scholar
- Huang M, Sharma S, Zhu L, Keane M, Luo J, Zhang L, Burdick M, Lin Y, Dohadwala M, Gardner B, et al.: IL-7 inhibits fibroblast TGF-beta production and signaling in pulmonary fibrosis. J Clin Invest. 2002, 109: 931-937.PubMed CentralView ArticlePubMedGoogle Scholar
- Cobbold CA, Sherratt JA: Mathematical Modelling of Nitric Oxide Activity in Wound Healing can explain Keloid and Hypertrophic Scarring. Journal of Theoretical Biology. 2000, 204: 257-288. 10.1006/jtbi.2000.2012View ArticlePubMedGoogle Scholar
- Zhang H, Ahmad M, Gronowicz G: Effects of transforming growth factor-beta 1 (TGF-[beta]1) on in vitro mineralization of human osteoblasts on implant materials. Biomaterials. 2003, 24: 2013-2020. 10.1016/S0142-9612(02)00616-6View ArticlePubMedGoogle Scholar
- Vandervelde S, van Luyn MJA, Rozenbaum MH, Petersen AH, Tio RA, Harmsen MC: Stem cell-related cardiac gene expression early after murine myocardial infarction. Cardiovasc Res. 2007, 73: 783-793. 10.1016/j.cardiores.2006.11.030View ArticlePubMedGoogle Scholar
- Sun Y, Weber KT: Cardiac remodeling by fibrous tissue: role of local factors and circulating hormones. Ann Med. 1998, 30: 3-8.PubMedGoogle Scholar
- Czarkowska-P¸czek B, Przybylski J, Marciniak A, Pawłowska M, Juskowa J, Foroncewicz B, Mucha K: Proteolytic Enzymes Activities in Patients After Myocardial Infarction Correlate with Serum Concentration of TGF-β. Inflammation. 2004, 28: 279-284. 10.1007/s10753-004-6051-2View ArticlePubMedGoogle Scholar
- Schaan B, Quadros A, Sarmento-Leite R, De Lucca G, Bender A, Bertoluci M: 'Correction:' Serum transforming growth factor beta-1 (TGF-beta-1) levels in diabetic patients are not associated with pre-existent coronary artery disease. Cardiovascular Diabetology. 2007, 6: 19- 10.1186/1475-2840-6-19PubMed CentralView ArticlePubMedGoogle Scholar
- Ignotz R, Massague J: Transforming growth factor-beta stimulates the expression of fibronectin and collagen and their incorporation into the extracellular matrix. J Biol Chem. 1986, 261: 4337-4345.PubMedGoogle Scholar
- Karagiannis ED, Popel AS: A Theoretical Model of Type I Collagen Proteolysis by Matrix Metalloproteinase (MMP) 2 and Membrane Type 1 MMP in the Presence of Tissue Inhibitor of Metalloproteinase 2. J Biol Chem. 2004, 279: 39105-39114. 10.1074/jbc.M403627200View ArticlePubMedGoogle Scholar
- Vempati P, Karagiannis ED, Popel AS: A Biochemical Model of Matrix Metalloproteinase 9 Activation and Inhibition. J Biol Chem. 2007, 282: 37585-37596. 10.1074/jbc.M611500200View ArticlePubMedGoogle Scholar
- Fang L, Gao XM, Moore XL, Kiriazis H, Su Y, Ming Z, Lim YL, Dart AM, Du XJ: Differences in inflammation, MMP activation and collagen damage account for gender difference in murine cardiac rupture following myocardial infarction. Journal of Molecular and Cellular Cardiology. 2007, 43: 535-544. 10.1016/j.yjmcc.2007.06.011View ArticlePubMedGoogle Scholar
- Kelly D, Cockerill G, Ng LL, Thompson M, Khan S, Samani NJ, Squire IB: Plasma matrix metalloproteinase-9 and left ventricular remodelling after acute myocardial infarction in man: a prospective cohort study. Eur Heart J. 2007, 28: 711-718. 10.1093/eurheartj/ehm003PubMed CentralView ArticlePubMedGoogle Scholar
- Kaden JJ, Dempfle CE, Sueselbeck T, Brueckmann M, Poerner TC, Haghi D, Haase KK, Borggrefe M: Time-Dependent Changes in the Plasma Concentration of Matrix Metalloproteinase 9 after Acute Myocardial Infarction. Cardiology. 2003, 99: 140-144. 10.1159/000070670View ArticlePubMedGoogle Scholar
- Sumitra M, Manikandan P, Nayeem M, Manohar BM, Lokanadam B, Vairamuthu S, Subramaniam S, Puvanakrishnan R: Time course studies on the initiation of complement activation in acute myocardial infarction induced by coronary artery ligation in rats. Molecular and cellular biochemistry. 2005, 268: 149-158. 10.1007/s11010-005-3856-8View ArticlePubMedGoogle Scholar
- Xiao-Ming G, Ziqiu M, Yidan S, Lu F, Helen K, Qi X, Anthony MD, Xiao-Jun D: Infarct size and post-infarct inflammation determine the risk of cardiac rupture in mice. International journal of cardiology. 2010, 143: 20-28. 10.1016/j.ijcard.2009.01.019View ArticleGoogle Scholar
- Gao XM, Xu Q, Kiriazis H, Dart AM, Du XJ: Mouse model of post-infarct ventricular rupture: time course, strain- and gender-dependency, tensile strength, and histopathology. Cardiovascular Research. 2005, 65: 469-477. 10.1016/j.cardiores.2004.10.014View ArticlePubMedGoogle Scholar
- Kjeldsen L, Sengelov H, Lollike K, Nielsen MH, Borregaard N: Isolation and Characterization of Gelatinase Granules from Human Neutrophils. Blood. 1994, 83: 1640-1649.PubMedGoogle Scholar
- Cleutjens JP, Kandala JC, Guarda E, Guntaka RV, Weber KT: Regulation of collagen degradation in the rat myocardium after infarction. Journal of molecular and cellular cardiology. 1995, 27: 1281-1292. 10.1016/S0022-2828(05)82390-9View ArticlePubMedGoogle Scholar
- Loftis MJ, Sexton D, Carver W: Effects of collagen density on cardiac fibroblast behavior and gene expression. Journal of Cellular Physiology. 2003, 196: 504-511. 10.1002/jcp.10330View ArticlePubMedGoogle Scholar
- Wu CH, Donovan CB, Wu GY: Evidence for pretranslational regulation of collagen synthesis by procollagen propeptides. Journal of Biological Chemistry. 1986, 261: 10482-10484.PubMedGoogle Scholar
- Wetzler C, Kampfer H, Stallmeyer B, Pfeilschifter J, Frank S: Large and Sustained Induction of Chemokines during Impaired Wound Healing in the Genetically Diabetic Mouse: Prolonged Persistence of Neutrophils and Macrophages during the Late Phase of Repair. 2000, 115: 245-253.Google Scholar
- Ertl G, Frantz S: Healing after myocardial infarction. Cardiovasc Res. 2005, 66: 22-32. 10.1016/j.cardiores.2005.01.011View ArticlePubMedGoogle Scholar
- Nakajima H, Nakajima HO, Salcher O, Dittie AS, Dembowsky K, Jing S, Field LJ: Atrial but Not Ventricular Fibrosis in Mice Expressing a Mutant Transforming Growth Factor-{beta}1 Transgene in the Heart. Circ Res. 2000, 86: 571-579.View ArticlePubMedGoogle Scholar
- Day J, Friedman A, Schlesinger LS: Modeling the immune rheostat of macrophages in the lung in response to infection. Proceedings of the National Academy of Sciences. 2009, 106 (27): 11246-51. 10.1073/pnas.0904846106.View ArticleGoogle Scholar
- Lin J, Lopez EF, Jin Y, Van Remmen H, Bauch T, Han HC, Lindsey ML: 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.005PubMed CentralView ArticlePubMedGoogle Scholar
- Lindsey ML, Goshorn DK, Squires CE, Escobar GP, Hendrick JW, Mingoia JT, Sweterlitsch SE, Spinale FG: Age-dependent changes in myocardial matrix metalloproteinase/tissue inhibitor of metalloproteinase profiles and fibroblast function. Cardiovasc Res. 2005, 66: 410-419. 10.1016/j.cardiores.2004.11.029View ArticlePubMedGoogle Scholar
- Zamilpa R, Lopez EF, Chiao YA, Dai Q, Escobar GP, Hakala K, Weintraub ST, Lindsey ML: Proteomic analysis identifies in vivo candidate matrix metalloproteinase-9 substrates in the left ventricle post-myocardial infarction. PROTEOMICS. 2010, 10: 2214-2223. 10.1002/pmic.200900587PubMed CentralView ArticlePubMedGoogle Scholar
- Squires CE, Escobar GP, Payne JF, Leonardi RA, Goshorn DK, Sheats NJ, Mains IM, Mingoia JT, Flack EC, Lindsey ML: Altered fibroblast function following myocardial infarction. Journal of Molecular and Cellular Cardiology. 2005, 39: 699-707. 10.1016/j.yjmcc.2005.07.008View ArticlePubMedGoogle Scholar
- Johan D, Heilborn KBAH: Inhibited proliferation of fibroblasts derived from chronic diabetic wounds and normal dermal fibroblasts treated with high glucose is associated with increased formation of L-lactate. Wound Repair and Regeneration. 1998, 6: 135-141. 10.1046/j.1524-475X.1998.60207.xView ArticleGoogle Scholar
- , : Loss of BRCA2 promotes prostate cancer cell invasion through up-regulation of matrix metalloproteinase-9. Cancer Science. 2008, 99: 553-563. 10.1111/j.1349-7006.2007.00719.xView ArticleGoogle Scholar
- Eberhardt W, Akool ELS, Rebhan J, Frank S, Beck KF, Franzen R, Hamada FMA, Pfeilschifter J: Inhibition of cytokine-induced MMP-9 expression by PPARalpha agpnists is indirect and is due to a no-mediated reduction of mRNA stability. J Biol Chem. 2002, M202008200,Google Scholar
- Eberhardt W, Akool ELS, Rebhan J, Frank S, Beck KF, Franzen R, Hamada FMA, Pfeilschifter J: Inhibition of Cytokine-induced Matrix Metalloproteinase 9 Expression by Peroxisome Proliferator-activated Receptor alpha Agonists Is Indirect and Due to a NO-mediated Reduction of mRNA Stability. J Biol Chem. 2002, 277: 33518-33528. 10.1074/jbc.M202008200View ArticlePubMedGoogle Scholar
- Monaco S, Sparano V, Gioia M, Sbardella D, Pierro DD, Marini S, Coletta M: Enzymatic processing of collagen IV by MMP-2 (gelatinase A) affects neutrophil migration and it is modulated by extracatalytic domains. Protein Sci. 2006, 15: 2805-2815. 10.1110/ps.062430706PubMed CentralView ArticlePubMedGoogle Scholar
- Iris Metzmacher PR, Abel Martin, Friess Wolfgang: In vitro binding of matrix metalloproteinase-2 (MMP-2), MMP-9, and bacterial collagenase on collagenous wound dressings. Wound Repair and Regeneration. 2007, 15: 549-555. 10.1111/j.1524-475X.2007.00263.xView ArticlePubMedGoogle Scholar
- Atkins P: Physical Chemistry for the Life Sciences. 2006, Oxford, Oxford University Press,Google Scholar

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