Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Genomics**|**v.13(Suppl 6); 2012**|**PMC3481436

Formats

Article sections

- Abstract
- Background
- Results
- Discussion
- Conclusions
- Methods
- Competing interests
- Authors' contributions
- Supplementary Material
- References

Authors

Related links

BMC Genomics. 2012; 13(Suppl 6): S21.

Published online 2012 October 26. doi: 10.1186/1471-2164-13-S6-S21

PMCID: PMC3481436

Yunji Wang,^{1,}^{3} Tianyi Yang,^{1,}^{3} Yonggang Ma,^{2,}^{3,}^{4} Ganesh V Halade,^{2,}^{3,}^{4} Jianqiu Zhang,^{1} Merry L Lindsey,^{2,}^{3,}^{4} and Yu-Fang Jin^{}^{1,}^{3}

Yu-Fang Jin: ude.astu@nij.gnafuy

Selected articles from the IEEE International Workshop on Genomic Signal Processing and Statistics (GENSIPS) 2011

Ranadip Pal, Yufei Huang and Yidong Chen

This supplement has not been supported by sponsorship or other external funding.

IEEE International Workshop on Genomic Signal Processing and Statistics (GENSIPS) 2011

4-6 December 2011

San Antonio, TX, USA

Copyright ©2012 Wang et al.; licensee 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.

This article has been cited by other articles in PMC.

About 6 million Americans suffer from heart failure and 70% of heart failure cases are caused by myocardial infarction (MI). Following myocardial infarction, increased cytokines induce two major types of macrophages: classically activated macrophages which contribute to extracellular matrix destruction and alternatively activated macrophages which contribute to extracellular matrix construction. Though experimental results have shown the transitions between these two types of macrophages, little is known about the dynamic progression of macrophages activation. Therefore, the objective of this study is to analyze macrophage activation patterns post-MI.

We have collected experimental data from adult C57 mice and built a framework to represent the regulatory relationships among cytokines and macrophages. A set of differential equations were established to characterize the regulatory relationships for macrophage activation in the left ventricle post-MI based on the physical chemistry laws. We further validated the mathematical model by comparing our computational results with experimental results reported in the literature. By applying Lyaponuv stability analysis, the established mathematical model demonstrated global stability in homeostasis situation and bounded response to myocardial infarction.

We have established and validated a mathematical model for macrophage activation post-MI. The stability analysis provided a possible strategy to intervene the balance of classically and alternatively activated macrophages in this study. The results will lay a strong foundation to understand the mechanisms of left ventricular remodelling post-MI.

Myocardial infarction is defined by pathology as myocytes necrosis and apoptosis due to prolonged ischemia. Since myocytes cannot divide and replace themselves, myocytes in the infarct area deprived of oxygen die and are replaced by a collagen scar. There is a series of cellular and molecular activities respond to MI in the myocardium. Myocytes apoptosis appears in the first 6 to 8 hours post-MI, and necrosis occurs in 12 hrs to 4 days post-MI [1]. Necrosis of myocytes results in significantly elevated interlukin-1 (IL-1), tumour necrotic factor-*β *(TNF-*β*), IL-10, and monocyte chemotactic protein-1 (MCP-1) levels. MCP-1 is a strong chemoattractants that recruit and confine monocytes to the injury site. It's been reported that over 95% of monocytes differentiate to macrophages [2]. There are two major types of macrophages post-MI: classically activated macrophages (M1) and alternatively activated macrophages (M2). Porcheray has reported a switch between M1 and M2 macrophages with *in vitro *stimuli including IL-4, IL-10, and TNF-*β *[3]. In addition, biomarkers of M1 and M2 macrophages show a temporal *in vivo *transition [3]. Since M1 and M2 macrophages are responsible for extracellular matrix (ECM) destruction and construction, respectively, the transition and dynamic balance between two macrophage phenotypes might lead to the balance between ECM destruction and construction, and thus determine the ECM remodeling post-MI.[4] Therefore, characterizing macrophage activation pattern is essential to better understand the ECM remodeling post-MI.

A large amount of experimental research has been conducted to elucidate the underlying mechanisms of macrophage activation, and an abundant accumulation of experimental results define on macrophage responses to different stimuli. There is a need, however, to systemically analyze the accumulated data and integrate the results into a framework that will allow a more complete understanding. To address this need, several mathematical models have been established to characterize the effects of macrophages on wound healing, inflammatory responses, and collagen synthesis post-MI [5-9]. However, most models do not consider the effect of macrophage activation patterns and ignore the differences between macrophage phenotypes. Therefore, the aim of this study was to establish and validate a set of ordinary differential equations to characterize macrophage activation patterns post-MI. Since our mathematical model was established based on *in vivo *and *in vitro *experimental results, all parameters in the model were determined by the averages of the experimental data.

We have collected experimental data from adult C57 mice and built a framework to represent the regulatory relationship among cytokines and macrophages. Based on this framework, we established a set of nonlinear differential equations to characterize the regulatory relationship for macrophage activation in the left ventricle post-myocardial infarction using physical chemistry laws. Our framework and the mathematical model were established based on the following three assumptions.

1) All monocytes that migrate to the infarct region are differentiated to unactivated macrophages [10].

2) All activated macrophages are differentiated from unactivated macrophage since previous studies have shown that <5% of macrophages undergo mitotic division [11,12].

3) All parameters and coefficients in this model are constant.

In this framework, myocytes and monocytes were considered as inputs to the system. Cellular densities of M1 and M2 macrophages were considered as the outputs of the system. We chose IL-1, IL-10, and TNF-*α *as three molecules which regulate macrophage activations in this mathematical model since they were well-recognized as stimuli for macrophage activation [13-15]. M1 macrophages and myocytes secrete IL-1 and TNF-*α*. M2 macrophages secrete IL-10 [3]. Further, TNF-*α *and IL-1 promote M1 activation and IL_10 promotes M2 activation [16]. IL-10 inhibits TNF-*α*, IL-1, and itself [17]. The input-output and regulating relationship were shown in Figure Figure11.

Temporal profiles of monocytes and myocytes densities were used as inputs to our mathematical framework (Figure (Figure1).1). Myocytes density in healthy adult mice was 6 × 10^{9 }cells/ml as an initial value. Myocytes numbers monotonically decreases post-MI and is directly associated with LV wall thickness. We have measured the LV wall thickness at days 0, 1, 3, 5, and 7 post-MI. The temporal profile of myocytes was determined by combining the initial value and the monotone progression trend (the crosses) as shown in Figure 2(A).

Macrophages density in the left ventricle of healthy adult mice is 2000 cells/ml, which will be used as initial values of unactivated macrophage density in this study [18,19]. Yang et al has measured temporal profiles of macrophages post-MI in mice at days 1, 2, 4, 7, 14, 21, and 28 [20]. The temporal profile of macrophages was obtained by fitting the experimental results as a continuous function as shown in Figure 2(B).

Based on our assumptions 1 and 2, all macrophages were differentiated from monocytes and emigrated from infarct area to the lymph node system. The estimation of unactivated macrophage based on the experimental results [18] is shown as follows,

$${M}_{un}\left(t\right)={M}_{un}\left(0\right){e}^{-\mu t}+{\int}_{0}^{t}{e}^{-\mu \left(t-s\right)}Mds,$$

(1)

where *M *denotes the differentiation rates of monocytes, *M _{un }*denotes the unactivated macrophage density, and

The mathematical model of macrophage activations is a set of nonlinear differential equations represented by cellular densities (cell number/ml) of *M _{un}*, M1 and M2, and concentrations (pg/ml) of chemical factors such as IL-1, IL-10, and TNF-

$$\frac{d{M}_{un}}{dt}=\stackrel{differentiation}{\overbrace{M}}\stackrel{activated\mathsf{\text{}}to\mathsf{\text{}}M1}{\overbrace{-{k}_{2}{M}_{un}\frac{I{L}_{1}}{I{L}_{1}+{c}_{IL1}}-{k}_{3}{M}_{un}\frac{{T}_{\alpha}}{{T}_{\alpha}+{c}_{{T}_{\alpha}}}}}-\stackrel{activated\mathsf{\text{}}to\mathsf{\text{}}M2}{\overbrace{{k}_{4}{M}_{un}\frac{I{L}_{10}}{I{L}_{10}+{c}_{IL10}}}}-\stackrel{emigration}{\overbrace{\mu {M}_{un}}},$$

(2)

$$\frac{d{M}_{1}}{dt}=\stackrel{activation\mathsf{\text{}}effects\mathsf{\text{}}by\mathsf{\text{}}I{L}_{1}\mathsf{\text{}}and\mathsf{\text{}}TN{F}_{\alpha}}{\overbrace{{k}_{2}{M}_{un}\frac{I{L}_{1}}{I{L}_{1}+{c}_{IL1}}+{k}_{3}{M}_{un}\frac{{T}_{\alpha}}{{T}_{\alpha}+{c}_{{T}_{\alpha}}}}}+\stackrel{transation\mathsf{\text{}}between\mathsf{\text{}}M1\mathsf{\text{}}and\mathsf{\text{}}M2}{\overbrace{{{k}^{\prime}}_{1}{M}_{2}-{k}_{1}{M}_{1}}}-\stackrel{emigration}{\overbrace{\mu {M}_{1}}},$$

(3)

$$\frac{d{M}_{2}}{dt}=\stackrel{activation\mathsf{\text{}}effects\mathsf{\text{}}by\mathsf{\text{}}I{L}_{10}}{\overbrace{{k}_{4}{M}_{un}\frac{I{L}_{10}}{I{L}_{10}+{c}_{IL10}}}}+\stackrel{transation\mathsf{\text{}}between\mathsf{\text{}}M1\mathsf{\text{}}and\mathsf{\text{}}M2}{\overbrace{{k}_{1}{M}_{1}-{{k}^{\prime}}_{1}{M}_{2}}}-\stackrel{emigration}{\overbrace{\mu {M}_{2}}},$$

(4)

$$\frac{dI{L}_{10}}{dt}=\stackrel{secretion\mathsf{\text{}}by\mathsf{\text{}}M2}{\overbrace{{k}_{5}{M}_{2}\frac{{c}_{1}}{{c}_{1}+I{L}_{10}}}}-\stackrel{degradation}{\overbrace{{d}_{IL10}I\phantom{\rule{0.3em}{0ex}}{L}_{10}}},$$

(5)

$$\frac{d{T}_{\alpha}}{dt}=\stackrel{secretion\mathsf{\text{}}by\mathsf{\text{}}M1\mathsf{\text{}}and\mathsf{\text{}}myocytes}{\overbrace{\left({k}_{6}{M}_{1}+\lambda Mc\right)\frac{c}{c+I{L}_{10}}}}-\stackrel{degradation}{\overbrace{{d}_{T\alpha}{T}_{\alpha}}},$$

(6)

$$\frac{dI{L}_{1}}{dt}=\stackrel{secretion\mathsf{\text{}}by\mathsf{\text{}}M1\mathsf{\text{}}and\mathsf{\text{}}myocytes}{\overbrace{\left({k}_{7}{M}_{1}+\lambda Mc\right)\frac{c}{c+I{L}_{10}}}}-\stackrel{degradation}{\overbrace{{d}_{IL1}I{L}_{1}}},$$

(7)

where *M _{un}*,

Equation 2 determines the density of unactivated macrophages in the infarct area. For the construction part, the unactivated macrophages are differentiated from monocytes as shown in Figure 2(C). For the destruction part, the unactivated macrophages are activated to *M*_{1 }or *M*_{2}. Additionally, inactivated macrophages do not die locally in the scar tissue but die out in the lymph node system [27].

Equation 3 determines the activation rate of M1 macrophages. For the construction part, IL-1 and TNF-*α *promote M1 activation. Parameters *k*_{2 }and *k*_{3 }denote the activation rates of M1 macrophages by IL-1 and TNF-*α*. Hill equations are used to represent the promotion effects of IL-1 and TNF-*α *and parameters *c*_{IL1 }and ${c}_{{T}_{\alpha}}$ are the effectiveness of IL-1 and TNF-*α *promotion on M1 calculated based on the experimental results [3,29]. Steinmuller has shown the transition between M1 and M2 phenotypes *in vivo *[21]. Correspondingly, we use parameter *k*_{1 }to denote the transition from M1 to M2 and parameter ${k}_{1}^{\prime}$ for the transition from M2 to M1 [21]. The destruction part includes emigration of macrophage (*μ*)and transition from M1 to M2 macrophages (*k*_{1}) [27].

Equation 4 determines the activation rate of M2 macrophages. The construction part is denoted by activation of M2 macrophages promoted by IL-10, and transition from M1 to M2. IL-10 promotes M2 activation and this activation rate has been approximated by parameters *k*_{4 }based on the experimental results [3]. The transition rate from M1 to M2 is denoted as *k*_{1}. The destruction part includes emigration of M2 macrophages (*μ*) and transition from M2 to M1 macrophages (${k}_{1}^{\prime}$), similarly as described in equation 3.

Equation 5 determines the secretion rate of IL-10. For construction part, IL-10 is secreted by M2 macrophages, and parameter *k*_{5 }denotes the secretion rate of IL-10 by M2 macrophages [22,30]. The destruction part includes the self-inhibition and degradation of IL-10. A Hill equation is employed to represent this self-inhibition effect and parameter *c*_{1 }denotes the self-inhibition effect of IL-10 post-MI [6,17]. Parameter ${d}_{I{L}_{10}}$ denotes the decay rate of IL-10 determined by its half-life time [13].

Equation 6 determines the deposition rate of TNF-*α*, which is secreted by both M1 macrophages and myocytes [6,23,28,31]. We used *in vitro *results from Meng to determine the secretion rate of TNF-*α *by M1 and secretion rate of TNF-*α *by myocytes (*λ*) is determined by Horio's experimental results [23,31]. The inhibition of IL-10 is presented by a Hill equation where parameter c represents the effectiveness of IL-10 inhibiting TNF-*α *[25]. The destruction part is denoted by the degradation of TNF-*α*. Parameter ${d}_{{T}_{\alpha}}$ is the decay rate of TNF-*α *determined by its half-life time [15,26].

Equation 7 determines the deposition rate of IL-1. IL-1 is secreted by both M1 macrophage and myocytes. Parameter *k*_{7 }denotes the secretion rate of IL-1 in cultured rat cardiac myocytes [31]. The inhibition of IL-10 is presented by a Hill equation similarly as in equation 6 [25]. In the destruction part, parameter *d*_{IL1 }represents the decay rate of IL-1 determined by its half-life time [14].

If there is no myocardial infarction, monocytes differentiation and myocytes apoptosis should be at a very low level, and the studied macrophage activation pathway should maintain homeostasis. We have calculated the equilibrium point of the system without any input and performed Lyapunov stability analysis. Our analysis showed that without any monocytes differentiation and myocytes secretion, the system would stay at the origin when *t *→ ∞.

In the case of myocardial infarction, myocytes apoptosis and necrosis triggered inflammatory responses and significant monocytes differentiation, which will drive the system to a new equilibrium point. Correspondingly, the cell densities of M1 and M2 increase post-MI. We have obtained a steady state as *E *= [20, 1200, 3500, 0.73, 1.1, 5.9] from our computational simulations. The steady states match with the experimental measurements collected from healthy adults without myocardial infarction [32]. In addition, the stable region of the established mathematical model depends on the strength of the input.

Computational simulations of macrophage activation were carried out by solving the nonlinear differential equations with MATLAB. The initial conditions of unactivated, M1 and M2 macrophage densities were chosen as *M _{un }*= 2000 cells/ml and

Our computational results demonstrated that from days 0 to day 3 post-MI, cellular densities of the M1 phenotype increased at a faster rate than the M2 phenotype. At day 10, the M2 phenotype dominates over the M1 phenotype. This prediction agrees with the results reported by Troidl [4]. Additionally, temporal profiles of IL-1 and TNF-*α *significantly increased from days 0 to 1 post-MI in our computational simulations, which match the experimental results reported by Sumitra [33]. Comparison between the computational and experimental results demonstrates a similar trend for the temporal profiles, suggesting the effectiveness of our model.

This study established a mathematical model for macrophage activation in the left ventricle post-MI by combining experimental and computational approaches. This is the first mathematical model focusing on the dynamic interactions among cytokines, myocytes, monocytes, and macrophages. Computational predictions based on this mathematical model match with experimental measurements, suggesting effectiveness of the model. In addition, our stability analysis provided insight for the activation pattern of macrophages post-myocardial infarction. In our mathematical model, there are two inputs, myocytes apoptosis and monocytes differentiation. In this study, we have predicted a stable equilibrium for homeostasis, which means without myocardial infarction or following small injury stimuli, macrophage densities and concentrations of IL-1, IL-10, and TNF-α should stay at the equilibrium. After MI, the monocytes differentiated into macrophage and apoptotic myocytes secreted significant amounts of cytokines to activate the macrophages. The strength of the monocytes differentiation and myocytes apoptosis (inputs of the system) drive the system to different states while all states will be bounded due to the bounded strength of the inputs.

However, there exist some differences between computational predictions and real experimental results. To address this issue and the variation in different experiments, a stochastic parameter distribution will need to be introduced to replace the constant parameters. In addition, more detailed measurements on monocytes and myocytes and concentration temporal profiles of IL-1, IL-10 and TNF-*α *from C57 mice will be carried out in our future research, which will help to solidify our equations.

Our study has established framework for macrophage activation and used ordinary differential equations to model the cellular interactions between macrophage activation types post-MI. The results on stability analysis can be used as a useful tool to predict the behaviour of biological systems.

To incorporate the experimental data, curve fitting algorithm was applied to obtain temporal continuous density profiles of myocytes and monocytes based on discrete experimental data.

To analyse the stability of the proposed mathematical model, we have calculated the equilibrium point of the system and performed the Lyapunov stability analysis of the system.

In our mathematical model, equations (2-7) are six first-order equations with input M and *M _{c}*. Now we use

$$V\left(x\right)=\frac{1}{2}{\sum}_{i=1}^{6}{x}_{i}^{2}$$

(8)

and obtained its derivative as the following equation

$$\begin{array}{llll}\hfill \dot{V}\left(x\right)& ={\sum}_{i=1}^{6}{x}_{i}{\u1e8b}_{i}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & =-\mu {x}_{1}^{2}-{x}_{1}\left[{k}_{2}{x}_{1}\frac{{x}_{6}}{{x}_{6}+{c}_{IL1}}+{k}_{3}{x}_{1}\frac{{x}_{5}}{{x}_{5}+{c}_{{T}_{\alpha}}}+{k}_{4}{x}_{1}\frac{{x}_{4}}{{x}_{4}+{c}_{IL10}}\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & -\left(\mu +{k}_{1}\right){x}_{2}^{2}+{k}_{2}{x}_{2}{x}_{1}\frac{{x}_{6}}{{x}_{6}+{c}_{IL1}}+{k}_{3}{x}_{1}{x}_{2}\frac{{x}_{5}}{{x}_{5}+{c}_{{T}_{\alpha}}}+{{k}^{\prime}}_{1}{x}_{3}{x}_{2}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & -\left({{k}^{\prime}}_{1}+\mu \right){x}_{3}^{2}+{k}_{1}{x}_{2}{x}_{3}+{k}_{4}{e}_{1}{x}_{3}\frac{{x}_{4}}{{x}_{4}+{c}_{IL10}}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & +{k}_{5}{x}_{3}{x}_{4}\frac{{c}_{1}}{{c}_{1}+{x}_{4}}-{d}_{{T}_{\alpha}}{x}_{5}^{2}+{k}_{6}{x}_{2}{x}_{5}\frac{c}{c+{x}_{4}}-{d}_{IL1}{x}_{6}^{2}+{k}_{7}{x}_{2}{x}_{6}\frac{c}{c+{x}_{4}}.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{2em}{0ex}}& \hfill \end{array}$$

(9)

By applying the boundary of the Hill equations, we got

$$\begin{array}{llll}\hfill \dot{V}\left(x\right)& \le -\left(\mu +{k}_{2}+{k}_{3}+{k}_{4}\right){x}_{1}^{2}-\left(\mu +{k}_{1}\right){x}_{2}^{2}+{k}_{2}{x}_{2}{x}_{1}+{k}_{3}{x}_{1}{x}_{2}+{{k}^{\prime}}_{1}{x}_{3}{x}_{2}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & -\left({{k}^{\prime}}_{1}+\mu \right){x}_{3}^{2}+{k}_{1}{x}_{2}{x}_{3}+{k}_{4}{x}_{1}{x}_{3}-{d}_{IL10}{x}_{4}^{2}+{k}_{5}{x}_{3}{x}_{4}-{d}_{{T}_{\alpha}}{x}_{5}^{2}+{k}_{6}{x}_{2}{x}_{5}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & -{d}_{IL1}{x}_{6}^{2}+{k}_{7}{x}_{2}{x}_{6}.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{2em}{0ex}}& \hfill \end{array}$$

(10)

Applying the parameters in Table 1 to equation (10), we got

$$\begin{array}{llll}\hfill \dot{V}\left(x\right)& \le -1.6{x}_{1}^{2}-0.275{x}_{2}^{2}-0.25{x}_{3}^{2}-2.5{x}_{4}^{2}-55{x}_{5}^{2}-10.5{x}_{6}^{2}+1.1{x}_{1}{x}_{2}+0.3{x}_{1}{x}_{3}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}0.125{x}_{2}{x}_{3}+0.0005{x}_{3}{x}_{4}+0.0007{x}_{2}{x}_{5}+0.0006{x}_{2}{x}_{6}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \le -0.02{x}_{1}^{2}-0.02{x}_{2}^{2}-0.05{x}_{3}^{2}-0.5{x}_{4}^{2}-5{x}_{5}^{2}-0.5{x}_{6}^{2}-1.28{\left({x}_{1}-0.429{x}_{2}\right)}^{2}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{1em}{0ex}}-0.3{\left({x}_{1}-0.5{x}_{3}\right)}^{2}-0.25{\left({x}_{3}-0.25{x}_{2}\right)}^{2}-2{\left({x}_{4}-0.000125{x}_{3}\right)}^{2}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{1em}{0ex}}-50{\left({x}_{5}-0.000007{x}_{2}\right)}^{2}-10{\left({x}_{6}-0.00003{x}_{2}\right)}^{2}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & <0.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{2em}{0ex}}& \hfill \end{array}$$

(11)

Since the derivative of Lyapunov function is negative semi-definite and the semi-definite is satisfied with all states equal to zero, the system is globally asymptotically stable without any input. Given any bounded differentiate rate of monocytes and myocytes density as input, the system will have bounded states.

The authors declare that they have no competing interests.

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

**Table 1**. Pre-determined parameters from literature search.

Click here for file^{(41K, pdf)}

Based on “Mathematical modeling of macrophage activation in left ventricular remodeling post-myocardial infarction”, by Yunji Wang, Yu-Fang Jin, Yonggang Ma, Ganesh V Halade and Merry L Linsey which appeared in *Genomic Signal Processing and Statistics (GENSIPS), 2011 IEEE International Workshop on. *© 2011 IEEE [38].

The authors acknowledge grant and contract support from NHLBI HHSN268201000036C (N01-HV-00244), NIH R01 HL75360, Veteran's Administration Merit Award, and the Max and Minnie Tomerlin Voelcker Fund (to M.L.L.), and NIH 1R03EB009496, and NIH SC2HL101430 (to Y.F.J.).

This article has been published as part of *BMC Genomics *Volume 13 Supplement 6, 2012: Selected articles from the IEEE International Workshop on Genomic Signal Processing and Statistics (GENSIPS) 2011. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/13/S6.

- Haunstetter A, Izumo S. Apoptosis: Basic Mechanisms and Implications for Cardiovascular Disease. Circ Res. 1998;82(11):1111–1129. doi: 10.1161/01.RES.82.11.1111. [PubMed] [Cross Ref]
- Li Q, Park PW, Wilson CL, Parks WC. Matrilysin Shedding of Syndecan-1 Regulates Chemokine Mobilization and Transepithelial Efflux of Neutrophils in Acute Lung Injury. Cell. 2002;111(5):635–646. doi: 10.1016/S0092-8674(02)01079-6. [PubMed] [Cross Ref]
- Porcheray F, Viaud S, Rimaniol AC, Léone C, Samah B, Dereuddre-Bosquet N, Dormont D, Gras G. Macrophage activation switching: an asset for the resolution of inflammation. Clinical & Experimental Immunology. 2005;142(3):481–489. [PubMed]
- Troidl C, Mollmann H, Nef H, Masseli F, Voss S, Szardien S, Willmer M, Rolf A, Rixe J, Troidl K. Classically and alternatively activated macrophages contribute to tissue remodelling after myocardial infarction. Journal of Cellular and Molecular Medicine. 2009;13:3485–3496. doi: 10.1111/j.1582-4934.2009.00707.x. [PubMed] [Cross Ref]
- Waugh Helen V, S JA. Modeling the effects of treating diabetic wounds with engineered skin substitutes. Wound Repair and Regeneration. 2007;15(4):556–565. doi: 10.1111/j.1524-475X.2007.00270.x. [PubMed] [Cross Ref]
- Day J, Friedman A, Schlesinger L. Modeling the immune rheostat of macrophages in the lung in response to infection. Proc Natl Acad Sci U S A. 2009;106(27):11246–11251. doi: 10.1073/pnas.0904846106. [PubMed] [Cross Ref]
- Vodovotz Y, Clermont G, Chow C, An G. Mathematical models of the acute inflammatory response. Curr Opin Crit Care. 2004;10:383–390. doi: 10.1097/01.ccx.0000139360.30327.69. [PubMed] [Cross Ref]
- Munoz-Alicea R-M, Marcano-Velazquez M. A mathematical model for macrophage, T-cell and mycobacterium tuberculosis interactions. 1999.
- Owen MR, Sherratt JA. Modelling the macrophage invasion of tumours: Effects on growth and composition. IMA J Math Appl Med Biol. 1998;15(2):165–185. doi: 10.1093/imammb/15.2.165. [PubMed] [Cross Ref]
- Krause S, Rehli M, Kreutz M, Schwarzfischer L, Paulauskis J, Andreesen R. Differential screening identifies genetic markers of monocyte to macrophage maturation. J Leukoc Biol. 1996;60:540–545. [PubMed]
- Mantovani A. Macrophage diversity and polarization: in vivo veritas. Blood. 2006;108(2):408–409. doi: 10.1182/blood-2006-05-019430. [Cross Ref]
- Burke B, L C. The Macrophage. 2. Oxford: Oxford University Press; 2002.
- Jin Y, Lindsey M. Multi-Scale Modeling and Analysis of Left Ventricular Remodeling Post Myocardial Infarction: Integration of Experimental and Computational Approaches. Application of Machine Learning. 2009. pp. 267–280.
- Kudo S, Mizuno K, Hirai Y, Shimizu T. Clearance and Tissue Distribution of Recombinant Human Interleukin 1α in Rats. Cancer Research. 1990;50(18):5751–5755. [PubMed]
- Tsutsumi Y, Kihira T, Tsunoda S, Kamada H, Nakagawa S, Kaneda Y, Kanamori T, Mayumi T. Molecular design of hybrid tumor necrosis factor-alpha III: polyethylene glycol-modified tumor necrosis factor-alpha has markedly enhanced antitumor potency due to longer plasma half-life and higher tumor accumulation. Journal of Pharmacology and Experimental Therapeutics. 1996;278(3):1006. [PubMed]
- Mocellin S, Panelli MC, Wang E, Nagorsen D, Marincola FM. The dual role of IL-10. Trends in Immunology. 2003;24(1):36–43. doi: 10.1016/S1471-4906(02)00009-1. [PubMed] [Cross Ref]
- Moore KW, de Waal Malefyt R, Coffman RL, O'Garra A. INTERLEUKIN-10 AND THE INTERLEUKIN-10 RECEPTOR. Annual Review of Immunology. 2001;19(1):683–765. doi: 10.1146/annurev.immunol.19.1.683. [PubMed] [Cross Ref]
- Jin Y-F, Han H-C, Berger J, Dai Q, Lindsey M. Combining experimental and mathematical modeling to reveal mechanisms of macrophage-dependent left ventricular remodeling. BMC Systems Biology. 2011;5(1):60. doi: 10.1186/1752-0509-5-60. [PMC free article] [PubMed] [Cross Ref]
- Gao X, Xu Q, Kiriazis H, Dart A, Du X. Mouse model of post-infarct ventricular rupture: time course, strain- and gender-dependency, tensile strength, and histopathology. Cardiovascular Research. 2005;65:469–477. doi: 10.1016/j.cardiores.2004.10.014. [PubMed] [Cross Ref]
- Yang F, Liu YH, Yang XP, Xu J, Kapke A, Carretero OA. Myocardial infarction and cardiac remodelling in mice. Exp Physiol. 2002;87(5):547–555. doi: 10.1113/eph8702385. [PubMed] [Cross Ref]
- Steinmuller C, Franke-Ullmann G, Lohmann-Matthes M-L, Emmendorffer A. Local Activation of Nonspecific Defense against a Respiratory Model Infection by Application of Interferon-gamma. Comparison between Rat Alveolar and Interstitial Lung Macrophages. Am J Respir Cell Mol Biol. 2000;22(4):481–490. [PubMed]
- Edwards JP, Zhang X, Frauwirth KA, Mosser DM. Biochemical and functional characterization of three activated macrophage populations. Journal of Leukocyte Biology. 2006;80:1298–1307. doi: 10.1189/jlb.0406249. [PMC free article] [PubMed] [Cross Ref]
- Meng F, Lowell CA. Lipopolysaccharide (LPS)-induced Macrophage Activation and Signal Transduction in the Absence of Src-Family Kinases Hck, Fgr, and Lyn. The Journal of Experimental Medicine. 1997;185(9):1661–1670. doi: 10.1084/jem.185.9.1661. [PMC free article] [PubMed] [Cross Ref]
- Mantovani A, Sica A, Sozzani S, Allavena P, Vecchi A, Locati M. The chemokine system in diverse forms of macrophage activation and polarization. Trends in Immunology. 2004;25(12):677–686. doi: 10.1016/j.it.2004.09.015. [PubMed] [Cross Ref]
- Marino S, Hogue IB, Ray CJ, Kirschner DE. A methodology for performing global uncertainty and sensitivity analysis in systems biology. Journal of theoretical biology. 2008;254(1):178–196. doi: 10.1016/j.jtbi.2008.04.011. [PMC free article] [PubMed] [Cross Ref]
- Oliver J, Bland L, Oettinger C, Arduino M, McAllister S, Aguero S, Favero M. Cytokine kinetics in an in vitro whole blood model following an endotoxin challenge. Lymphokine and cytokine research. 1993;12(2):115. [PubMed]
- Rohde L, Ducharme A, Arroyo L, Aikawa M, Sukhova G, Lopez-Anaya A, McClure K, Mitchell P, Libby P, Lee R. Matrix metalloproteinase inhibition attenuates early left ventricular enlargement after experimental myocardial infarction in mice. Circulation. 1999;99:3063–3070. doi: 10.1161/01.CIR.99.23.3063. [PubMed] [Cross Ref]
- Maass DL, Hybki DP, White J, Horton JW. The Time Course of Cardiac NF-[kappa]B Activation and TNF-[alpha] Secretion by Cardiac Myocytes After Burn Injury: Contribution to Burn-Related Cardiac Contractile Dysfunction. Shock. 2002;17(4):293–299. doi: 10.1097/00024382-200204000-00009. [PubMed] [Cross Ref]
- Onozaki K, Urawa H, Tamatani T, Iwamura Y, Hashimoto T, Baba T, Suzuki H, Yamada M, Yamamoto S, Oppenheim J. Synergistic interactions of interleukin 1, interferon-beta, and tumor necrosis factor in terminally differentiating a mouse myeloid leukemic cell line (M1). Evidence that interferon-beta is an autocrine differentiating factor. The Journal of Immunology. 1988;140(1):112–119. [PubMed]
- Edwards JP, Zhang X, Frauwirth KA, Mosser DM. Biochemical and functional characterization of three activated macrophage populations. Journal of Leukocyte Biology. 2006;80(6):1298–1307. doi: 10.1189/jlb.0406249. [PMC free article] [PubMed] [Cross Ref]
- Horio T, Nishikimi T, Yoshihara F, Nagaya N, Matsuo H, Takishita S, Kangawa K. Production and Secretion of Adrenomedullin in Cultured Rat Cardiac Myocytes and Nonmyocytes: Stimulation by Interleukin-1α and Tumor Necrosis Factor-α Endocrinology. 1998;139(11):4576–4580. doi: 10.1210/en.139.11.4576. [PubMed] [Cross Ref]
- Steensberg A, Fischer CP, Keller C, Møller K, Pedersen BK. IL-6 enhances plasma IL-1ra, IL-10, and cortisol in humans. American Journal of Physiology - Endocrinology And Metabolism. 2003;285(2):E433–E437. [PubMed]
- Sumitra M, Manikandan P, Nayeem M, Manohar B, 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. doi: 10.1007/s11010-005-3856-8. [PubMed] [Cross Ref]
- Frangogiannis NG, Mendoza LH, Lindsey ML, Ballantyne CM, Michael LH, Smith CW, Entman ML. IL-10 Is Induced in the Reperfused Myocardium and May Modulate the Reaction to Injury. The Journal of Immunology. 2000;165(5):2798–2808. [PubMed]
- Yang F, Liu Y, Yang X, Xu J, Kapke A, Carretero O. Myocardial infarction and cardiac remodelling in mice. Experimental Physiology. 2002;87(5):547–555. doi: 10.1113/eph8702385. [PubMed] [Cross Ref]
- Patti G, D'Ambrosio A, Mega S, Giorgi G, Zardi EM, Zardi DM, Dicuonzo G, Dobrina A, Di Sciascio G. Early interleukin-1 receptor antagonist elevation in patients with acute myocardial infarction. Journal of the American College of Cardiology. 2004;43(1):35–38. doi: 10.1016/j.jacc.2003.07.032. [PubMed] [Cross Ref]
- Drimal J, Knezl V, Paulovicova E, Drimal D. Enhanced early after-myocardial infarction concentration of TNF-alpha subsequently increased circulating and myocardial adrenomedullin in spontaneously hypertensive rats. General physiology and biophysics. 2008;27(1):12. [PubMed]
- Wang Y, Jin Y, Ma Y, Halade GV, Linsey ML. Mathematical modeling of macrophage activation in left ventricular remodeling post-myocardial infarction. Genomic Signal Processing and Statistics (GENSIPS), 2011 IEEE International Workshop on: 4-6 December 2011. 2011. pp. 202–205. [Cross Ref]

Articles from BMC Genomics are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |