Computational and Mathematical Methods in Medicine

Volume 2012 (2012), Article ID 736394, 16 pages

doi:10.1155/2012/736394

Research Article

Computational Modeling of Microabscess Formation

Alexandre Bittencourt Pigozzo,¹ Gilson Costa Macedo,² Rodrigo Weber dos Santos,¹ and Marcelo Lobosco¹

¹Graduate Program in Computational Modeling, UFJF, Rua José Lourenço Kelmer s/n, Campus Universitário, Bairro São Pedro, 36036-900 Juiz de Fora, MG, Brazil

²Graduate Program in Biological Sciences, UFJF, Rua José Lourenço Kelmer s/n, Campus Universitário, Bairro São Pedro, 36036-900 Juiz de Fora, MG, Brazil

Received 29 June 2012; Revised 28 August 2012; Accepted 5 September 2012

Academic Editor: Vladimir Brusic

Copyright © 2012 Alexandre Bittencourt Pigozzo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Bacterial infections can be of two types: acute or chronic. The chronic bacterial infections are characterized by being a large bacterial infection and/or an infection where the bacteria grows rapidly. In these cases, the immune response is not capable of completely eliminating the infection which may lead to the formation of a pattern known as microabscess (or abscess). The microabscess is characterized by an area comprising fluids, bacteria, immune cells (mainly neutrophils), and many types of dead cells. This distinct pattern of formation can only be numerically reproduced and studied by models that capture the spatiotemporal dynamics of the human immune system (HIS). In this context, our work aims to develop and implement an initial computational model to study the process of microabscess formation during a bacterial infection.

1. Introduction

The immune system is one of the most important and complex system of our organism. Despite great advances in recent years that shed light on its understanding and unravel the underlying key mechanisms behind its functions, there are still many functions of the human immune system (HIS) that are not well understood. Computational models of HIS dynamics can contribute to a better understanding of the relationship between cells and molecules of the HIS.

In this study, we developed a mathematical model of some cells and molecules of the HIS to reproduce the spatiotemporal dynamics of the initial formation of microabscesses during an immune response to a bacteria.

To reproduce these dynamics, we introduce a mathematical model composed of a system of partial differential equations (PDEs) that extends our previous models [1, 2] and defines the dynamics of representative cells and molecules of the HIS during the immune response to a bacteria. The model presented is descriptive, mechanistic, and deterministic; therefore, it enables the understanding of how different complex phenomena, structures, and elements interact during an immune response. In addition, the model’s parameters reflect the physiological features of the system, making the model appropriate for general use.

The remainder of this paper is organized as follows. First, the necessary biological background is presented. Next, related works are briefly discussed. This exposition is followed by a description of the mathematical model proposed in this work and the numerical scheme used to implement it. Then simulation results obtained from the proposed model are discussed, and, finally, our conclusions and plans for future work are presented.

2. Biological Background

The initial response of the host to a diverse array of biological stressors including bacterial infection, burns, trauma, and invasive surgery is an inflammatory response. Despite the growing understanding of the cellular and molecular mechanisms of inflammation, the complexity of the inflammatory response has challenged therapeutic development [3, 4]. A key reason for this conundrum has been speculated to be the difficulty of predicting the impact of manipulating individual components of the highly complex, nonlinear, and redundant inflammatory response [5]. Thus, progress would require a greater understanding of how components are organized. This makes systems biology based approaches appealing [6].

Most inflammatory reactions begin as a suppurative or purulent exudation process when the first line of cellular defense, the neutrophils, accumulate in the area. A suppurative process is characterized by the presence of pus (neutrophils mixed with cellular debris). Classically, there are three requisites of suppuration:

(i)presence of neutrophils that release proteolytic enzymes,(ii)necrosis of some types,(iii)liquefaction.

This suppurative process may lead to the formation of microabscesses. A microabscess is a localized collection of dead cells, body fluids, microbes, and other cells of the HIS. The process of formation of a microabscess begins when a cell of the HIS encounters bacteria and warn other cells that there is a stranger in the host. Its “warn” is in the form of a class of biochemicals called cytokines, which beckon other HIS cells to come to the point of infection and surround the enemy.

Most often all this goes unnoticed because the first few immune system cells phagocytize (engulf and digest) the invaders and the battle is finished. But every so often, an invader has a trick to escape the immune response and cannot be killed by the phagocytes. Those bacteria continue to grow and to spew out whatever they do. More and more immune system cells, mainly neutrophils, congregate at the infection site trapping the pathogens in the center. If this microabscess is close enough to the body surface, it can be seen as a blob of pus under the skin. When a microabscess is well developed, it has a wall or capsule of fibrous connective tissue separating it from the surrounding tissue, helping to prevent any microbes present in the microabscess from spreading to other areas of the body. Thus, microabscesses can be considered as a natural strategy used to fight against infection.

Microabscesses are found in many different diseases, for instance, the authers in [7–10] present animal studies detailing the formation of liver microabscess and microabscess by different types of infections. Epidermal microabscess formation by neutrophils was also evaluated in [11–14]. Infection of the heart by bacteria (bacterial myocarditis [15]) or by viruses (viral myocarditis [16]) is also correlated with microabscess formation by neutrophils. The interaction between tumor cells and inflammatory cells plays an important role in cancer initiation and progression and was investigated in [17] for the case of tumor-infiltrating neutrophils in pancreatic neoplasia, where the pattern of microabscess formation by neutrophils was reported once again.

3. Related Work

This section presents and discusses other mathematical and computational models of the immune response. Essentially, two distinct approaches are used: agent-based models and ordinary differential equations (ODEs). These models have some features in common with our model. All models include representative cells and molecules of the innate immune system. For example, neutrophils, macrophages, and proinflammatory cytokines are modeled in the majority of models. Some models as our model consider the important interactions between endothelial cells, tissue cells, and cytokines. Despite some similarities with our model, none of the works focus on modeling microabscesses.

4. Models Based on Agents

In [18, 19], it was developed an agent-based model of the dynamics of some cells, such as polymorphonuclear leukocytes (PMNs) and mononuclear cells and molecules, such as TNF-α and IL-1, during the initial inflammatory response in the interface endothelium/blood at the capillary level. Some characteristics of the model are as follows:

(i)all the cells are represented as agents whose behavior is close to the real;(ii)it considers the interactions between endothelial cells and circulating inflammatory cells at the blood/blood vessel-lining interface;(iii)the initial injury number (IIN) defines the number of tissue cells that are dead initially;(iv)the oxygen concentration is one important variable in the model;(v)the total tissue damage is represented as a deficit in the oxygen variable;(vi)the injury state of an endothelial cell depends on the available oxygen concentration; (vii)proinflammatory mediators and endothelial cells surface adhesion molecules are modeled by state variables;(viii)it considers a generic pathogen that causes the infection.

The work aims to reproduce the time course of the early inflammatory response associated with the Systemic Inflammatory Syndrome Response (SIRS)/Multiple Organ Failure (MOF) from massive trauma or large exposure to endotoxin. The objective of the simulations were to compare the results with the soluble TNF-receptor experiment [20] where the soluble TNF-receptor is tested as a therapeutic treatment for the sepsis. The author claims that his results generally replicate the results of several large-scale clinical trials of cytokine-directed antimediator agents.

In [6], an agent-based modeling (ABM) framework is proposed to study the nonlinear dynamics of acute inflammatory responses to LPS. Their work uses an agent-based approach to elucidate molecular interactions involved in the NF- signaling pathway, coupled with the spatial orientation of various inflammation specific molecules and cell populations such as macrophages and T-helper cells. In their model, the propagation of LPS signaling across the system is studied by considering the coupling between extracellular signals and transcriptional response. Some relevant characteristics considered in the model are as follows:

(i)molecular interactions, (ii)cellular heterogeneity,(iii)LPS/TLR4 signal transduction pathway,(iv)transcriptional response.

The proposed in silico model is evaluated through its ability to successfully reproduce a self-limited inflammatory response as well as a series of scenarios: a persistent (non)infectious response or innate immune tolerance and potentiation effects followed by perturbations in intracellular signaling molecules and cascades.

5. Models Based on ODEs

The model of [21] studies immunomodulatory strategies for treating cases of severe sepsis. They introduced and evaluated the concept of conducting a randomized clinical trial in silico based on simulated patients generated from a mechanistic mathematical model of bacterial infection, acute inflammatory response, global tissue dysfunction, and a therapeutic intervention. Trial populations are constructed to reflect heterogeneity in bacterial load and virulence as well as propensity to mount and modulate an inflammatory response. They constructed a cohort of trial patients submitted to therapy with one of three different doses of a neutralizing antibody directed against tumor necrosis factor (anti-TNF) for , , or hrs. Their focus was to assess the feasibility of using differential equation models to improve the design of clinical trials. This paper replicates in silico the general findings from actual clinical trials—that it is very difficult to design a treatment strategy that is effective over a broad range of sepsis patients.

6. Hybrid Models

In [22], a hybrid model that coupled an agent-based model [18, 19] and a system dynamics/differential equation model [21] was created using the System Dynamics tool within Netlogo [23]. They developed an interface between the agent-based and system dynamics models. The area of initial infection and subsequent interactions at this point were simulated with the agent-based model and certain aspects primarily related to the production and life cycle of circulating inflammatory cells was modeled with the system dynamics model. They implemented a simple system dynamics model that focused on systemic polymorphoneutrophils (PMN) production, maturation, sequestration, and release. Some of the flows inside the system dynamics model were influenced by the conditions within the agent-based model. The primary role of the system dynamics model was to manifest a delay between the elevation of the cytokines in the tissue and the increase in PMNs in the circulating blood. The objective of the work was to reproduce some of the results of the agent-based model [18, 19] using the hybrid model.

7. Mathematical Model

Our main objective is to develop a parameterized mathematical model of the human innate immune system that simulates the immune response occurring in a generic tissue. To achieve this goal, we first build a model of the immune response to LPS [1, 2]. In this work, we extend this model to reproduce the spatiotemporal dynamics of a bacterial infection and the process of microabscess formation.

The mathematical model simulates the temporal and spatial behavior of bacteria (), dead bacteria (), macrophages, neutrophils (), apoptotic neutrophils (), proinflammatory cytokines (), healthy tissue cells (), and dead tissue cells (). Macrophages are present in two states of readiness: resting () and hyperactivated (). We must stress that the equations modeling proinflammatory cytokines are generic in the sense that they model the role of distinct cytokines taking part in the inflammatory process. Equation parameters can be adjusted to model the role of a specific proinflammatory cytokine.

The relationships among all of the model’s components are described next. Neutrophils, resting macrophages, and active macrophages phagocytose the bacteria. The neutrophils then undergo apoptosis, which may or may not be induced by the phagocytosis process. In this different state, apoptotic neutrophils cannot perform phagocytosis or produce proinflammatory cytokines; as a result, apoptotic neutrophils are eliminated from the body after being phagocytosed by active macrophages. Apoptotic neutrophils will die after a period of time, releasing cytotoxic granules and degradation enzymes in the medium that cause tissue damage destroying healthy tissue cells. Active neutrophils and bacteria also cause tissue damage by producing toxic products that are not modeled here. The infection site is “cleaned” by resting and active macrophages that do the phagocytosis of dead tissue cells. Healthy tissue cells in contact with bacteria, neutrophils, and active macrophages produce proinflammatory cytokines. The proinflammatory cytokines increase the permeability of the blood vessels; consequently, more neutrophils and monocytes are recruited to the infected tissue. In addition, the proinflammatory cytokines act as a chemoattractant substance to the resting macrophages, active macrophages, and neutrophils.

Below, we provide the equations derived from the model. Equation (1) provides the bacteria differential equation:

In this equation, denotes the reproduction term of the bacteria, where is the rate of reproduction and is a function of the total density of cells in a discretized area of the two-dimensional space at a specific time step. The variable is defined as

denotes the decay of bacteria, where is the rate of decay. denotes the phagocytosis of bacteria by neutrophils, where is the rate of this phagocytosis. denotes the phagocytosis of bacteria by resting macrophages, where is the rate of this phagocytosis. denotes the phagocytosis of bacteria by active macrophages, where is the rate of this phagocytosis. denotes bacteria diffusion, where represents the diffusion coefficient and is calculated in the following way:

The function models the probability of a cell being pushed from a site due to the pressure exerted by neighboring cells [24, 25]. This population pressure is modelled by a Hill equation [26]. It increases with the total density of cells occupying the same position in space and has a saturation in a high density of cells.

The function is defined as

and are constant values.

The function returns the percentage of free space in a discretized area of the two-dimensional space and its use is motivated by some important biological concepts such as quorum sensing/volume sensing [24, 25, 27–29]. The idea is that cells have a set of cell density sensing mechanisms and changes its behavior in crowded regions. In the context of our model it is used to limit the density of cells that occupy a discretized area of our two-dimensional domain. The function is defined as

The variable total represents the maximum density of cells that fits in a discretized area of the tissue.

The differential equation corresponding to dead bacteria (BD) is given as follows:

Here, note that , , and were defined previously. denotes the phagocytosis of dead bacteria by active macrophages, where is the rate of phagocytosis. denotes the phagocytosis of dead bacteria by resting macrophages, where is the rate of phagocytosis. denotes dead bacteria diffusion, where represents the diffusion coefficient and the function was defined previously.

The differential equation corresponding to the resting macrophage (RM) is given as follows:

denote the increase in endothelium permeability and its effects on monocyte extravasation. The permeability of blood vessel endothelium is modeled by a Hill equation [26], which also has been used to model drug dose-response relationships [30]. The idea is to model the increase in the permeability of the endothelium in accordance with the number of proinflammatory cytokines deposited on the endothelium.

The calculation of involves the following parameters: (a) , the maximum endothelium permeability induced by the proinflammatory cytokine; (b) , the minimum endothelium permeability induced by the proinflammatory cytokine; (c) keqCH, the number of proinflammatory cytokines that exert of the maximum effect on permeability. denotes the source term of macrophages, which is related to the number of monocytes that will enter into the tissue from the blood vessels. This number depends on the endothelium permeability and on the number of monocytes appearing in the blood ().

denotes resting macrophage apoptosis, where is the apoptosis rate. denotes the activation of resting macrophages, where is the rate of activation. denotes resting macrophage diffusion, where represents the diffusion coefficient and the function was defined previously. denotes resting macrophage chemotaxis, where is the chemotaxis rate and is calculated in the following way:

The differential equation corresponding to the active macrophage (AM) is given as follows:

Here, note that was defined previously. Above, , , and denote the active macrophage apoptosis, diffusion, and chemotaxis, respectively, whereas , , and are the apoptosis rate, diffusion coefficient, and chemotaxis rate, respectively.

The neutrophil differential equation () is given as follows:

In this equation, denotes the increase in endothelium permeability and its effects on neutrophil extravasation. In the top equation, is the maximum endothelium permeability induced by proinflammatory cytokines, is the minimum endothelium permeability induced by proinflammatory cytokines, and keqCH is the number of proinflammatory cytokines that exert of the maximum effect on endothelium permeability.

Here, denotes neutrophil apoptosis, where is the rate of apoptosis. denotes the neutrophil apoptosis induced by phagocytosis, where represents the rate of this induced apoptosis. represents the source term of neutrophil, that is, the number of neutrophils entering the tissue from the blood vessels. This number depends on the endothelium permeability () and on the number of neutrophils in the blood ().

denotes neutrophil diffusion, where represents the diffusion coefficient and the function was defined previously. denotes neutrophil chemotaxis, where is the chemotaxis rate and was defined previously.

The differential equation corresponding to the apoptotic neutrophil (ND) is given as follows:

Here, note that and were defined previously, whereas denotes the phagocytosis of the apoptotic neutrophil carried out by active macrophages, and is the rate of this phagocytosis. denotes the neutrophil necrosis, where is the rate of necrosis. denotes apoptotic neutrophil diffusion, where represents the diffusion coefficient and the function was defined previously.

The differential equation for the proinflammatory cytokine (CH) is given in as follows:

In this equation, denotes the proinflammatory cytokine decay, where is the decay rate. denotes the proinflammatory cytokine production by the neutrophils, where is the production rate. denotes the proinflammatory cytokine production by active macrophages, where is the production rate. denotes the proinflammatory cytokine production by healthy tissue cells in contact with bacteria, where is the production rate. denotes the proinflammatory cytokine diffusion, where represents the diffusion coefficient and the function was defined previously.

The differential equation corresponding to the healthy tissue (HT) is given as follows:

denotes the tissue damage caused by the release of toxic products from necrotic neutrophils. denotes the tissue damage caused by bacteria, where is the rate of damage.

The differential equation corresponding to the dead tissue (TD) is given as follows:

and were defined previously. denotes the phagocytosis of dead tissue cells by resting macrophages, where is the rate of phagocytosis. denotes the phagocytosis of dead tissue cells by active macrophages, where is the rate of phagocytosis.

The mathematical model presented here introduced some modifications to our previous model [2] with the aim to reproduce the microabscess formation. We included equations for the dynamics of the tissue to take into account some effects of infection such as tissue damage and production of cytokines by tissue cells. We also replaced the LPS equation in our previous model [2] by the bacteria equation with a term for reproduction of bacteria. Besides we modified the calculus of the diffusion and chemotaxis terms [24] (a) to limit the number of cells that are allowed to stay at the same time in the same area of the domain and (b) to reduce the efficiency of the diffusion and chemotaxis processes in overcrowded regions. More specifically, the method we implemented incorporates the following general mechanisms which may lead to dispersal of the population [24].

(i)Population pressure: we assume that a high cell density results in increased probability of a cell being pushed from a site, for example, due to the pressure exerted by neighboring cells. This is achieved phenomenologically with the function of our mathematical model and the changes in diffusion and chemotaxis calculus. (ii)Limited space: here we assume that no more cells can enter a site above a total cell density. In our model, this is achieved with the incorporation of the function in the diffusion and chemotaxis calculations. (iii)Gradient detection: cells may detect and respond to a local gradient in the cell density and as a consequence cells can move to higher concentrations of the attractant substance.

8. Implementation

The numerical method used to solve the mathematical model was the Finite Difference Method [31], a method commonly used in the numeric discretization of PDEs.

A complex part of the resolution of the PDEs is the resolution of the convective term, the chemotaxis term. The development of numerical methods to approximate convective terms (in most cases not linear) have been subject of intense researches [32–35].

Different numerical approaches have been proposed for the discretization of the chemotaxis term [36, 37]. Our implementation is based on the finite difference method for the spatial discretization and the explicit Euler method for the time evolution. The discretization of the chemotaxis term () uses the First-Order Upwind scheme [38]. Therefore, the precision of our numerical implementation is first-order in time (explicit Euler) and first-order in space (upwind scheme). The upwind scheme discretizes the hyperbolic PDEs through the use of differences with bias in the direction given by the signal of the characteristics’ speeds. The upwind scheme uses an adaptive or solution-sensitive stencil to numerically simulate more precisely the direction of information propagation.

In two-dimension, the upwind scheme approximates the chemotaxis term as the sum of the flux in the direction () with the flux in the direction (). is the sum of the at the point with the at the point and is the sum of the at the point with the at the point in Algorithm 1.

float chemotaxis (float** , float **ch, int , int , float **)

{

    flux_left = 0;

    flux_right = 0;

    if ( > 0)

    {

             if (()

             {

                    flux_left = − (ch – ch[])* [] *()* )/delta;

            }

            else

            {

                    flux_left = − (ch[] – ch[])* *()*()/delta;

            }

    }

    if ( < (size − 1))

    {

            if ((ch[ – ch[]) > 0)

            {

                    flux_right = (ch[] [] – ch[])*[] []*()*()/delta;

            }

            else

            {

                    flux_right = (ch[] [] – ch[])*[]*(])*()/delta;

            }

    }

    res = (flux_left + flux_right)/delta;

    flux_up = 0;

    flux_down = 0;

    if ( > 0)

    {

            if ((ch[] [] – ch[] []) > 0)

            {

                    flux_up = − (ch[] [] – ch[][])* [][]*(][])*()/delta;

            }

            else

            {

                    flux_up = − (ch[] [] – ch[][])* [ []* ()*()/delta;

            }

    }

    if ( < (size − 1))

    {

            if ((ch[] [] – ch[] []) > 0)

            {

                    flux_down = (ch[] [] – ch[] [])*] []*)*] [])/delta;

            }

            else

            {

                    flux_down = (ch[] [] – ch[] [])*[] []*(] [])*/delta;

            }

    }

    res = (flux_up + flux_down)/delta;

    return (res + res);

}

Algorithm 1

In this code, represents the discretization of the proinflammatory cytokine, represents the discretization of neutrophils, is the total density of cells in a position of the space, and are the positions in space, and and are the spatial discretizations in and directions, respectively. The test made is to define what is the signal of the characteristic speed, where the speed of the movement of is given by the term . This value is then used to choose between two schemes of finite differences: forward or backward.

We decided to implement our own numerical method to solve the systems of PDEs because (a) we have the possibility to parallelize the code and (b) most of the numerical libraries offer few functions that are suitable to our problem. The sequential code was implemented in C.

9. Numerical Experiments

We performed several simulations in order to verify that the model’s results are in agreement with what is described in the literature. Our objective was to reproduce some characteristics of the microabscess such as an accumulation of dead cells and bacteria in the infection site.

The model’s initial conditions and parameters are given in Tables 1 and 2, respectively.

Table 1: Initial conditions.

Parameter
Value
Unit

 cells/mm³

 cells/mm³

 cells/mm³

 cells/mm³

N₀

 cells/mm³

cells/mm³

 cells/mm³

cells/mm³

 cells/mm³

Table 2: Parameters.

Parameter
Value
Unit
Reference

0.05
Adimensional
Estimated*

35
Cells/mm³
Estimated*

Total
70
Cells/mm³
Estimated*

15000
Cells/mm³
[44]**

250000
Cells/mm³
[44]**

1
1/day
[45]**

1
1/day
Estimated*

0.001
1/day
Estimated*

0.01
1/day
Estimated*

keqCH
5
Cells/mm³
Estimated*

4
1/day
[46]

0.05
1/(cells/mm³)·day
[47]

0.6
1/(cells/mm³)·day
[47]**

0.8
1/(cells/mm³)·day
[44]

0.6
1/(cells/mm³)·day
[47]**

0.6
1/(cells/mm³)·day
[47]**

0.6
1/(cells/mm³)·day
[47]**

0.55
1/(cells/mm³)·day
[44]

0.24
1/(cells/mm³)·day
[47]

0.25
1/(cells/mm³)·day
[44]

0.8
1/(cells/mm³)·day
[44]

0.01
1/day
[44]

0.67
1/day
[44]

0.05
1/day
[44]

0.0033
1/day
[44]

0.07
1/day
[44]

12
1/day
[44]

0.05
mm²/day
[44]

5
mm²/day
[44]

5
mm²/day
[44]

10
mm²/day
[44]

0.001
mm²/day
[44]

6
mm²/day
[44]

10
mm²/day
[44]

5
mm²/day
[44]

7
mm²/day
[44]

1
1/(cells/mm³)·day
[48]*

1
1/(cells/mm³)·day
[48]*

0.2
1/(cells/mm³)·day
[48]*

RMact
0.4
1/(cells/mm³)·day
[44]

In Table 2, parameters marked with * were adjusted to qualitatively reproduce the results obtained in several studies of the immune response to LPS. In the case of the bacteria (results not shown here), we adjust the equation parameters in order to obtain an exponential decrease in time as shown in [39]. The results of the concentration of proinflammatory cytokines in time (results not shown here) are qualitatively similar to those obtained in some experimental works [40–42]. The parameters marked with ** were based on the values given in the references but were adjusted due to the use of distinct units (e.g., from liter to ) or to fit in a tissue. In this paper, we obtained parameter values for humans whenever they were available. The variables and in (5) and (4) were defined based on the concentrations of neutrophils and macrophages per liter given in [43] and were adjusted (a) due to the use of distinct units (e.g., from liter to ) and (b) to fit in a tissue. The variable total represents the maximum density of cells that fits in a discretized area of the tissue.

In the next sections, we will show the results of the simulation performed with the parameters given in Table 2. In this simulation, we considered a two-dimensional domain representing a tissue with 25  of area and a simulation time of 5 days. In our model, the exchange between the vascular system (arterioles and vessels) and tissue was assumed to occur only at the points , , and , In this point, immune cells (neutrophils and macrophages) that are in the blood stream can enter into the tissue. The communication between blood vessels and tissue is modeled by permeabilities that vary in time and may depend on the concentration of different cells and molecules (in our model, the endothelium permeability of neutrophils and macrophages depends on the concentration of the proinflammatory cytokine).

10. Bacteria

In the case of bacteria (Figure 1), we observe that initially the bacteria diffuses through the tissue causing tissue damage without its presence to be noticed.

Figure 1: Temporal evolution and spatial distribution of bacteria.

As soon as resting macrophages residents in the tissue recognize the bacteria they start to produce proinflammatory cytokines that will diffuse through the tissue reaching the blood vessel. Once proinflammatory cytokines interact with the endothelial cells an increase in the endothelium permeability occurs allowing neutrophils and monocytes to migrate to the tissue.

The bacteria starts to die a lot due to the presence of huge numbers of neutrophils. However, the immune response can not completely eliminate bacteria due to the formation of the microabscess pattern. In the microabscess, there are bacteria and a huge concentration of dead cells around it (Figures 4, 5, and 6). In this context, the cleaning process realized by macrophages is very important to allow neutrophils to reach bacteria and eliminate them. Macrophages are responsible for phagocyte dead cells that accumulated in the microabscess. The pattern of microabscess could have lasted longer if we had considered the formation of fibrous tissue around the microabscess.

11. Neutrophil

Neutrophils are initially attracted to the tissue by proinflammatory cytokines produced by activated resident macrophages (Figure 2). Once a neutrophil encounters bacteria, it phagocytizes bacteria and starts to produce proinflammatory cytokines that will attract more neutrophils and macrophages. The cytokine gradient will guide the movement of neutrophils and macrophages in the direction of the highest bacteria concentration.

Figure 2: Temporal evolution and spatial distribution of neutrophil.

After the microabscess formation, the immune system cells lose contact with a high number of bacteria. As these cells tend to move following the cytokine gradient, we can observe an accumulation of them around the microabscess.

After a significant number of macrophages phagocyte dead cells, the neutrophils can encounter the bacteria and phagocyte them from the border to the center of the microabscess area. As a consequence, a reduction in the microabscess area is observed (Figure 1), which indicates that the immune response is succeeding in controlling the infection.

12. Cytokine

The cytokines in Figure 3 are produced primarily by resident macrophages that are the first to recognize the bacteria presence. The cytokines will increase the endothelium permeability allowing neutrophils to migrate to the tissue. The arrived neutrophils will produce even more cytokines that will guide the movement of neutrophils and macrophage cells in the direction of high concentrations of bacteria.

Figure 3: Temporal evolution and spatial distribution of cytokine.

Figure 4: Temporal evolution and spatial distribution of apoptotic neutrophil.

Figure 5: Temporal evolution and spatial distribution of dead bacteria.

Figure 6: Temporal evolution and spatial distribution of damaged tissue cells.

During the formation of the microabscess and after it, the production of cytokines is higher in the regions where neutrophils and macrophages have contact with the bacteria that is surrounding the microabscess.

13. Apoptotic Neutrophils

In Figure 4, it can be observed that initially the neutrophils that came from the blood vessel closer to the site of infection died in large number than the neutrophils that came from other sites. Then, after microabscess formation, a lot of neutrophils start to die around the entire microabscess. This phenomenon continues until the microabscess disappears.

14. Dead Bacteria

Figure 5 shows that initially more bacteria died near blood vessels. After the formation of the microabscess, the bacteria starts to die around the entire microabscess since this is the area where the immune response is acting.

15. Tissue Cells

In Figure 6, it can be observed that a lot of tissue damage by bacteria during the period the immune system took to mount an effective immune response. The number of dead tissue cells then reduces, because of the phagocytosis realized by macrophages.

Figure 7 shows the evolution of the healthy tissue area destroyed by the bacteria.

Figure 7: Temporal evolution and spatial distribution of healthy tissue cells.

16. Microabscess Area

Figure 8 shows the microabscess area using a set of level curves. We defined the microabscess as an area where the concentration of bacteria plus concentration of dead bacteria plus the concentrations of damaged tissue and apoptotic neutrophil is higher. These results show that our model was capable to reproduce the formation of the microabscess in agreement with the observed characteristics of a microabscess [7, 14, 49–54].

Figure 8: Level curves highlighting the microabscess area at day 5 of the immune response.

17. Conclusions and Future Works

In this work, we presented a computational model for the dynamics of representative types of cells and molecules of the HIS during an immune response to a bacteria. Despite the simplifications and limitations of the model, our results showed that we were able to reproduce an initial microabscess formation. The spatial results show a collection of dead tissue cells, dead bacteria, and apoptotic neutrophil in the microabscess region. This distinct pattern of formation can only be reproduced by spatiotemporal models, such as PDEs models.

As future work, we plan to perform a detailed sensitivity analysis of our microabscess formation model. A previous work [55] has given us some hints about the most sensitivity parameter of the model. We also plan to validate our model against experimental data.

We also plan to modify many aspects of the model to make it more real. For example, we plan to consider a more adequate model to represent the structure of the tissue and its constituents. The tissue can be better characterized as a multiphasic porous medium subjected to stress and deformation variations mainly during the inflammatory process. This porous medium would comprise fluids, extracellular matrix, cells, and molecules. We also plan to model the mechanical behaviour of each of these phases and the mechanical interactions between them.

We have interest in developing models for processes such as vasodilation, coagulation, and others and analyse its effects on the mechanical behaviour of immune system cells and the consequences for the immune response.

With the aim to investigate better the formation of the microabscess, we plan to add another features that contribute to this formation such as the effects of extracellular pH on immune response [56]. Acidic pH predominates at inflammatory loci and other sites of immune activity. Investigations on neutrophils demonstrate mainly inhibition of chemotaxis, respiratory activity, and bactericidal capacity at reduced pH. Besides diminished extracellular pH may play a role in suppressing cytokine production and cytotoxic activities by pulmonary macrophages [56].

Besides we plan to add to the model the process of fibrous tissue formation around the microabscess. We plan to investigate what factors determine if the fibrous tissue will be produced or not. The production of fibrous tissue as well as the coagulation process are ways of the immune system to prevent the bacteria to spread throughout the body doing damage with possible serious consequences, for example, SIRS/MOF. In particular, we are interested in modeling the participation of macrophages, fibroblasts, tissue cells, endothelial cells, and many mediators in the process of fibrous tissue formation. An important step in this process is the production of collagen by fibroblasts induced by the cytokine TGF- produced by macrophages [57, 58]. The macrophages has many roles in the processes of wound healing and tissue repair. For example, during the coagulation process, macrophages and endothelial cells are responsible for the production of diverse growth factors and chemotcatic substances that attracts and stimulates the proliferation of tissue cells initiating tissue repair [59].

Acknowledgments

The authors would like to thank FAPEMIG, CNPq, CAPES, and UFJF for supporting this study.

References

A. B. Pigozzo, G. C. Macedo, R. W. Dos Santos, and M. Lobosco, “Implementation of a computational model of the innate immune system,” in 10th International Conference on Artificial Immune Systems (ICARIS '11), vol. 6825, pp. 95–107, July 2011.

A. B. Pigozzo, G. C. Macedo, R. W. dos Santos, and M. Lobosco, “On the computational modelling of the innate immune system,” BMC Bioinformatics. In press

E. J. Kerschen, J. A. Fernandez, B. C. Cooley et al., “Endotoxemia and sepsis mortality reduction by non-anticoagulant-activated protein C,” Journal of Experimental Medicine, vol. 204, no. 10, pp. 2439–2448, 2007.

V. K. Klaitman and Y. Almog, “Corticosteroids in sepsis: a new concept for an old drug,” Israel Medical Association Journal, vol. 5, no. 1, pp. 51–55, 2003.

Y. Vodovotz, G. Constantine, J. Rubin, M. Csete, E. O. Voit, and G. An, “Mechanistic simulations of inflammation: current state and future prospects,” Mathematical Biosciences, vol. 217, no. 1, pp. 1–10, 2009.

X. Dong, P. T. Foteinou, S. E. Calvano, S. F. Lowry, and I. P. Androulakis, “Agent-based modeling of endotoxin-induced acute inflammatory response in human blood leukocytes,” PLoS ONE, vol. 5, no. 2, Article ID e9249, 2010.

K. B. Seydel, T. Zhang, and S. L. Stanley, “Neutrophils play a critical role in early resistance to amebic liver abscesses in severe combined immunodeficient mice,” Infection and Immunity, vol. 65, no. 9, pp. 3951–3953, 1997.

M. C. Rigothier, H. Khun, P. Tavares, A. Cardona, M. Huerre, and N. Guillén, “Fate of Entamoeba histolytica during establishment of amoebic liver abscess analyzed by quantitative radioimaging and histology,” Infection and Immunity, vol. 70, no. 6, pp. 3208–3215, 2002.

J. C. Lin, F. Y. Chang, C. P. Fung et al., “Do neutrophils play a role in establishing liver abscesses and distant metastases caused by Klebsiella pneumoniae?” PLoS ONE, vol. 5, no. 11, Article ID e15005, 2010.

Y. Ebe, G. Hasegawa, H. Takatsuka et al., “The role of Kupffer cells and regulation of neutrophil migration into the liver by macrophage inflammatory protein-2 in primary listeriosis in mice,” Pathology International, vol. 49, no. 6, pp. 519–532, 1999.

M. Schön, D. Denzer, R. C. Kubitza, T. Ruzicka, and M. P. Schön, “Critical role of neutrophils for the generation of psoriasiform skin lesions in flaky skin mice,” Journal of Investigative Dermatology, vol. 114, no. 5, pp. 976–983, 2000.

A. Wetzel, T. Wetzig, U. F. Haustein et al., “Increased neutrophil adherence in psoriasis: role of the human endothelial cell receptor Thy-1 (CD90),” Journal of Investigative Dermatology, vol. 126, no. 2, pp. 441–452, 2006.

W. Grayson, “The HIV-positive skin biopsy,” Journal of Clinical Pathology, vol. 61, no. 7, pp. 802–817, 2008.

M. I. Cybulsky, I. J. Cybulsky, and H. Z. Movat, “Neutropenic responses to intradermal injections of Escherichia coli. effects on the kinetics of polymorphonuclear leukocyte emigration,” American Journal of Pathology, vol. 124, no. 1, pp. 1–9, 1986.

J. Bhavsar, D. Gordon, and M. Shea, “Listeria myocarditis with an atypical intracavitary vegetation/thrombus,” Heart, vol. 96, no. 20, p. 1637, 2010.

K. Oka, K. Oohira, Y. Yatabe et al., “Fulminant myocarditis demonstrating uncommon morphology—a report of two autopsy cases,” Virchows Archiv, vol. 446, no. 3, pp. 259–264, 2005.

M. D. Reid, O. Basturk, D. Thirabanjasak et al., “Tumor-infiltrating neutrophils in pancreatic neoplasia,” Modern Pathology, vol. 24, no. 12, pp. 1612–1619, 2011.

G. An, “Agent-based computer simulation and sirs: building a bridge between basic science and clinical trials,” Shock, vol. 16, no. 4, pp. 266–273, 2001.

G. An, “In silico experiments of existing and hypothetical cytokine-directed clinical trials using agent-based modeling,” Critical Care Medicine, vol. 32, no. 10, pp. 2050–2060, 2004.

C. J. Fisher, J. M. Agosti, S. M. Opal et al., “Treatment of septic shock with the tumor necrosis factor receptor:Fc fusion protein,” The New England Journal of Medicine, vol. 334, no. 26, pp. 1697–1702, 1996.

G. Clermont, J. Bartels, R. Kumar, G. Constantine, Y. Vodovotz, and C. Chow, “In silico design of clinical trials: a method coming of age,” Critical Care Medicine, vol. 32, no. 10, pp. 2061–2070, 2004.

W. Wakeland, L. Macovsky, and G. An, “A hybrid simulation model for studying acute inammatory response,” in Proceedings of the Spring Simulation Multiconference (SpringSim '07), vol. 2, pp. 39–46, Society for Computer Simulation International, San Diego, Calif, USA, 2007.

NetLogo, http://ccl.northwestern.edu/netlogo/

K. J. Painter and J. A. Sherratt, “Modelling the movement of interacting cell populations,” Journal of Theoretical Biology, vol. 225, no. 3, pp. 327–339, 2003.

K. J. Painter, “Continuous models for cell migration in tissues and applications to cell sorting via differential chemotaxis,” Bulletin of Mathematical Biology, vol. 71, no. 5, pp. 1117–1147, 2009.

S. Goutelle, M. Maurin, F. Rougier et al., “The Hill equation: a review of its capabilities in pharmacological modelling,” Fundamental and Clinical Pharmacology, vol. 22, no. 6, pp. 633–648, 2008.

K. J. Painter and T. Hillen, “Volume-filling and quorum-sensing in models for chemosensitive movement,” Canadian Applied Mathematics Quarterly, vol. 10, no. 4, pp. 501–543, 2002.

H. M. Byrne and M. R. Owen, “A new interpretation of the Keller-Segel model based on multiphase modelling,” Journal of Mathematical Biology, vol. 49, no. 6, pp. 604–626, 2004.

Z. A. Wang, “On chemotaxis models with cell population interactions,” Mathematical Modelling of Natural Phenomena, vol. 5, no. 3, pp. 173–190, 2010.

J. G. Wagner, “Kinetics of pharmacologic response I. Proposed relationships between response and drug concentration in the intact animal and man,” Journal of Theoretical Biology, vol. 20, no. 2, pp. 173–201, 1968.

R. J. LeVeque, “Finite Difference Methods for Ordinary and Partial Differential Equations Steady State and Time Dependent Problems,” Society for Industrial and Applied Mathematics. In press

A. Harten, “High resolution schemes for hyperbolic conservation laws,” Journal of Computational Physics, vol. 135, no. 2, pp. 260–278, 1997.

B. P. Leonard, “Simple high-accuracy resolution program for convective modelling of discontinuities,” International Journal for Numerical Methods in Fluids, vol. 8, no. 10, pp. 1291–1318, 1988.

C. W. Shu and S. Oshert, “Efficient implementation of essentially non-oscillatory shock-capturing schemes, II,” Journal of Computational Physics, vol. 83, no. 1, pp. 32–78, 1989.

G. A. Sod, “A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws,” Journal of Computational Physics, vol. 27, no. 1, pp. 1–31, 1978.

A. Marrocco, “Numerical simulation of chemotactic bacteria aggregation via mixed finite elements,” Mathematical Modelling and Numerical Analysis, vol. 37, no. 4, pp. 617–630, 2003.

F. Filbet, “A finite volume scheme for the Patlak-Keller-Segel chemotaxis model,” Numerische Mathematik, vol. 104, no. 4, pp. 457–488, 2006.

M. M. Hafez and J. J. Chattot, Innovative Methods For Numerical Solution of Partial Differential Equations, World Scientific Publishing Company, 2002.

Y. Li, A. Karlin, J. D. Loike, and S. C. Silverstein, “A critical concentration of neutrophils is required for effective bacterial killing in suspension,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 12, pp. 8289–8294, 2002.

R. De Waal Malefyt, J. Abrams, B. Bennett, C. G. Figdor, and J. E. De Vries, “Interleukin 10(IL-10) inhibits cytokine synthesis by human monocytes: an autoregulatory role of IL-10 produced by monocytes,” Journal of Experimental Medicine, vol. 174, no. 5, pp. 1209–1220, 1991.

I. P. Oswald, T. A. Wynn, A. Sher, and S. L. James, “Interleukin 10 inhibits macrophage microbicidal activity by blocking the endogenous production of tumor necrosis factor α required as a costimulatory factor for interferon γ-induced activation,” Proceedings of the National Academy of Sciences of the United States of America, vol. 89, no. 18, pp. 8676–8680, 1992.

G. D. Martich, R. L. Danner, M. Ceska, and A. F. Suffredini, “Detection of interleukin 8 and tumor necrosis factor in normal humans after intravenous endotoxin: the effect of antiinflammatory agents,” Journal of Experimental Medicine, vol. 173, no. 4, pp. 1021–1024, 1991.

M. H. Beers, R. S. Porter, and T. V. Jones, The Merck Manual, Merck & Co., 18th edition, 2006.

B. Su, W. Zhou, K. S. Dorman, and D. E. Jones, “Mathematical modelling of immune response in tissues,” Computational and Mathematical Methods in Medicine, vol. 10, pp. 1748–6718, 2009.

T. H. Price, H. D. Ochs, R. Gershoni-Baruch, J. M. Harlan, and A. Etzioni, “In vivo neutrophil and lymphocyte function studies in a patient with leukocyte adhesion deficiency type II,” Blood, vol. 84, no. 5, pp. 1635–1639, 1994.

R. Kumar, G. Clermont, Y. Vodovotz, and C. C. Chow, “The dynamics of acute inflammation,” Journal of Theoretical Biology, vol. 230, no. 2, pp. 145–155, 2004.

A. Reynolds, J. Rubin, G. Clermont, J. Day, Y. Vodovotz, and G. Bard Ermentrout, “A reduced mathematical model of the acute inflammatory response: I. Derivation of model and analysis of anti-inflammation,” Journal of Theoretical Biology, vol. 242, no. 1, pp. 220–236, 2006.

A. Andoh, H. Takaya, T. Saotome et al., “Cytokine regulation of chemokine (IL-8, MCP-1, and RANTES) gene expression in human pancreatic periacinar myofibroblasts,” Gastroenterology, vol. 119, no. 1, pp. 211–219, 2000.

C. W. Hanke, H. R. Higley, D. M. Jolivette, N. A. Swanson, and S. J. Stegman, “Abscess formation and local necrosis after treatment with Zyderm or Zyplast Collagen Implant,” Journal of the American Academy of Dermatology, vol. 25, no. 2, part 1, pp. 319–326, 1991.

G. A. MacDonald, J. K. Greenson, E. A. DelBuono et al., “Mini-microabscess syndrome in liver transplant recipients,” Hepatology, vol. 26, no. 1, pp. 192–197, 1997.

L. W. Lamps, C. W. Pinson, D. S. Raiford, Y. Shyr, M. A. Scott, and M. K. Washington, “The significance of microabscesses in liver transplant biopsies: a clinicopathological study,” Hepatology, vol. 28, no. 6, pp. 1532–1537, 1998.

A. D. Gable, D. K. Marsee, D. A. Milner, and S. R. Granter, “Suppurative inflammation with microabscess and pseudocyst formation is a characteristic histologic manifestation of cutaneous infections with rapid-growing Mycobacterium species,” American Journal of Clinical Pathology, vol. 130, no. 4, pp. 514–517, 2008.

J. C. Lin, F. Y. Chang, C. P. Fung et al., “Do neutrophils play a role in establishing liver abscesses and distant metastases caused by Klebsiella pneumoniae?” PLoS ONE, vol. 5, no. 11, Article ID e15005, 2010.

H. S. Alsaif, S. K. Venkatesh, D. S. G. Chan, and S. Archuleta, “CT appearance of pyogenic liver abscesses caused by Klebsiella pneumoniae,” Radiology, vol. 260, no. 1, pp. 129–138, 2011.

P. Rocha, A. Pigozzo, B. Quintela, G. Macedo, R. Santos, and M. Lobosco, “Modelling the innate immune system,” in Bio-Inspired Computational Algorithms and their Applications, InTech, 2012.

A. Lardner, “The effects of extracellular pH on immune function,” Journal of Leukocyte Biology, vol. 69, no. 4, pp. 522–530, 2001.

E. P. Amento, N. Ehsani, H. Palmer, and P. Libby, “Cytokines and growth factors positively and negatively regulate interstitial collagen gene expression in human vascular smooth muscle cells,” Arteriosclerosis and Thrombosis, vol. 11, no. 5, pp. 1223–1230, 1991.

M. R. Duncan, K. S. Frazier, S. Abramson et al., “Connective tissue growth factor mediates transforming growth factor β- induced collagen synthesis: downregulation by cAMP,” The FASEB Journal, vol. 13, no. 13, pp. 1774–1786, 1999.

A. C. Newby and A. B. Zaltsman, “Fibrous cap formation or destruction—the critical importance of vascular smooth muscle cell proliferation, migration and matrix formation,” Cardiovascular Research, vol. 41, no. 2, pp. 345–360, 1999.