Skip to main content
  • Research article
  • Open access
  • Published:

Qualitative modelling of the interplay of inflammatory status and butyrate in the human gut: a hypotheses about robust bi-stability



Gut microbiota interacts with the human gut in multiple ways. Microbiota composition is altered in inflamed gut conditions. Likewise, certain microbial fermentation products as well as the lipopolysaccharides of the outer membrane are examples of microbial products with opposing influences on gut epithelium inflammation status. This system of intricate interactions is known to play a core role in human gut inflammatory diseases. Here, we present and analyse a simplified model of bidirectional interaction between the microbiota and the host: in focus is butyrate as an example for a bacterial fermentation product with anti-inflammatory properties.


We build a dynamical model based on an existing model of inflammatory regulation in gut epithelial cells. Our model introduces both butyrate as a bacterial product which counteracts inflammation, as well as bacterial LPS as a pro-inflammatory bacterial product. Moreover, we propose an extension of this model that also includes a feedback interaction towards bacterial composition. The analysis of these dynamical models shows robust bi-stability driven by butyrate concentrations in the gut. The extended model hints towards a further possible enforcement of the observed bi-stability via alteration of gut bacterial composition. A theoretical perspective on the stability of the described switch-like character is discussed.


Interpreting the results of this qualitative model allows formulating hypotheses about the switch-like character of inflammatory regulation in the gut epithelium, involving bacterial products as constitutive parts of the system. We also speculate about possible explanations for observed bimodal distributions in bacterial compositions in the human gut. The switch-like behaviour of the system proved to be mostly independent of parameter choices. Further implications of the qualitative character of our modeling approach for the robustness of the proposed hypotheses are discussed, as well as the pronounced role of butyrate compared to other inflammatory regulators, especially LPS, NF- κB and cytokines.


The human gut harbours a plethora of microbiota able to digest and process a vast amount of dietary compounds. These are needed for the human metabolism and have regulatory implications for the immune system [15, 25]. Dysbiosis is defined as a drastic change in the composition of the gut microbiota, and is often associated with disease: major examples are Irritable Bowel Syndrome (IBS) and Inflammatory Bowel Disease (IBD) [27, 28], metabolic syndrome and type-2 diabetes [31, 54], obesity and gut problems in the elderly [7, 36]. Generally speaking, dysbiosis has been hypothesized to involve abnormal inflammatory status or at least long or even chronic, systemic low-grade inflammation [6].

Bacterial products are known to play decisive roles in triggering or down-regulating host inflammatory responses. For example, lipopolysaccharides (LPS) which are part of certain bacteria’s cell walls, are potent pro-inflammatory stimuli. Short-chain fatty acids (SCFA), as another example, are key metabolic products from fibre processing bacteria in the gut lumen [43, 53]. They have multiple effects for human health [5, 26]. Butyrate, in particular, plays a major role in the regulation of inflammation as an anti-inflammatory signal for the epithelial cells of the large intestine [8]. It is derived from microbial fermentation of dietary fibres in the colon. Anti-carcinogenic, anti-inflammatory, and barrier-protective activity in the distal gut are among many of its health-promoting effects [19, 40]. Butyrate also accounts for ≥ 70% of the energy used by healthy colonocytes [32], and it also functions as a histone deacetylase inhibitor inside the nucleus to epigenetically regulate gene expression and cell fate [10, 19]. Given its impact on a variety of mechanisms, there is a growing interest in butyrate and the microbes that produce this compound. Indeed, butyrate-producing bacteria are important for a healthy colon and, if reduced, contribute to emerging diseases such as IBD [55] and IBS [41]. For example, butyrate has been shown to down-regulate inflammation responses in Crohn’s disease through inhibition of NF- κB activation in immune cells [47]. The NF- κB-signalling system is implicated in the regulation of a variety of genes during immune and inflammatory responses, including those encoding pro-inflammatory cytokines such as TNF, IL-1 and IL-6. Low concentrations of butyrate have also been detected in IBD patients [52]. Treatment with butyrate enemas has been shown to reduce inflammation in this patient group [45]. Inflammatory signalling in healthy volunteers has also been shown to be affected by butyrate concentrations [20]. In addition, butyrate is believed to play a role in maintaining intestinal barrier function, since a deficit in butyrate causes tight junction lesions [40], and it decreases the permeability in intestinal cell models [38]. IBS is an example of a disease characterized by a systemic low-grade inflammation. IBS patients have also been shown to suffer from increased intestinal permeability, suggesting a deteriorated intestinal barrier function, also known as leaky gut [39, 58]. Increased intestinal permeability would allow translocation of endotoxic compounds, such as LPS, that trigger a sustained immune response. These processes are proposed to subsequently result in a low-grade systemic inflammatory response [27].

Inflammation is proposed to mediate a decreased butyrate uptake into gut epithelial cells. The specific carrier-mediated transport systems involved in the transport of butyrate from the colonic lumen into colonic epithelial cells are mono-carboxylate transporter 1 (MCT1) and sodium-coupled mono-carboxylate transporter 1 (SMCT1) [17, 44]. These transporters are down-regulated in patients with IBD [42, 49]. The reduced expression of MCT1 and SMCT1 in IBD has been suggested to be a consequence of intestinal inflammation, since treatment of human intestinal epithelial cells with pro-inflammatory cytokines (e.g., IFN- γ and TNF- α) down-regulates MCT1 expression, leading to butyrate deficiency [50].

Dysbiosis of the gut microbiota has been frequently related to reduced bacterial diversity [24]. In addition, bimodal distributions for certain phylogenetic subgroups of gut bacteria have been previously described [30], yielding the hypothesis that these species are decisive indicators of malfunction in the gut microbiota’s ecosystem. An additional level of complexity relates to the mucin layer in the large intestine, where the microbiota is mainly involved in fibre processing [53]. Mucus layer status has been proposed to determine distinct gut microbiota ecosystems called enterotypes [1, 13].

Here, we introduce and analyse a drastically simplified, qualitative model for the intricate interplay between bacteria, their products and gut inflammatory status. The key components of this dynamical model are the interactions between bacterial SCFA production, where butyrate is a key player, LPS as pro-inflammatory signal, intracellular uptake of butyrate, barrier function and inflammatory status of the gut epithelial cell layer, represented by NF- κB and its regulators. We implement two models, a core model and an extended version, as systems of ordinary differential equations which are then investigated by fixed-point analyses. Our modelling approach aims at discovering general systems characteristics.


Simplified inflammation model by Yde et al. [57]

We base our investigation on the simplistic inflammation model by [57]. It is a system with four interacting elements which will be summarized by the following symbols and their corresponding classes of agents in the inflammatory regulatory system:

  1. N

    indicating pro-inflammatory elements such as NF- κB or upstream in the regulatory system

  2. R

    representing repressors, such as IkB, IkBe, A20, Cesanne

  3. T

    representing pro-inflammatory cytokines, such as TNF, IL-1, IL-6

  4. L

    equals S in the original model formulation by [57], which is an activator of T. In their model this corresponds to the bacterial load, which in turn is represented by LPS, a pro-inflammatory outer membrane compound.

Using these species (N,R,T,L), our model is based on the approach proposed by [57] with some adapted parameters, as defined in the following system of differential equations:

$$ \dot{N} = k_{a} \ \frac{T^{3}}{T^{3}+k_{s}^{3}} \ (1-N) - k_{ai} \ R \ \frac{N}{N+k_{r}} $$
$$ \dot{R} = k_{b} \ N - k_{br} \ R $$
$$ \dot{T} = k_{p} \ \frac{N^{2}}{N^{2}+k_{n}^{2}} - k_{t} \ T + k_{tl} \ L $$
$$ \dot{L} = k_{bl} - k_{l} \ L $$

Equations (1-4) represent the model by [57] with some simplifications: Diffusion is not considered, as we neglect the spatial aspect. Self-limitation is added for all elements to prevent unrealistic negative values. As pointed out by [57], both [48] and [12] already analysed the NF- κB inflammatory system with respect to limit boundaries and activation threshold points: These are adequately modelled in this model, based on a Hill-function. The model in its non-spatial version is bistable, with one fixed point in the origin and another fixed point with all components (N,T,R,L) strictly positive [57].

Core model implementation and analysis

Implementation of the core model

The simplified inflammation model by [57] is now further developed, and integrated with the basic interactions of butyrate, to constitute our core model. The first step is to assemble a list of interactions between elements of the inflammatory system and butyrate, in gut epithelial cells. The precise implementation of the interactions between certain elements of the model, such as butyrate, LPS, or parts of the inflammatory system, is to be understood as a coarse-grained, simplified approach to the underlying regulatory processes. Based on the biological knowledge as reviewed in the “Background” section, the relevant interactions for our proposed core model are defined: Butyrate leaky barrier Butyrate is supposed to strengthen the intestinal barrier and thus in our model it blocks LPS influx from the gut lumen into epithelial cells. Butyrate NF- κB Butyrate blocks the process of NF- κB entering the nucleus during ongoing inflammation. Inflammation (cytokines) butyrate transport In the inflamed state, butyrate transporter proteins (as e.g. MCT1) are down-regulated.

From these interactions, together with the non-spatial part of the simplified model of inflammation by [57], the core model is now built: Butyrate B interacts with and reduces the influx of LPS L, by strengthening the barrier function of the gut epithelial cells. Therefore, it reduces the activation of the cytokines T, as influx of LPS L is reduced. Butyrate also blocks the entrance of NF- κB into the nucleus. Cytokines T on the other hand reduce the expression of the transporter for butyrate, hence diminishing its transport into the epithelial cells from the lumen. All interactions of the model by [57] together with the new interactions via butyrate can be collected in an interaction graph, see Fig. 1. The core model is now constructed as a five-dimensional ODE system:

Fig. 1
figure 1

Interaction graph of the core model. LPS (L) in the intestinal lumen environment enters the epithelial cell through leaky barrier, butyrate (B) via transporter channel. Butyrate enhances barrier function of epithelial cells and blocks LPS to enter. Transport of lumen butyrate into the cell is blocked by cytokines (T). L, T, NF- κB (N) and repressors (R) are representing the inflammatory model variables in the model by [57] which interact with butyrate (B), as described above and by B blocking NF- κB’s transport into the nucleus (modelled as blocking NF- κB). Black arrows represent transport, green arrows activation, red arrows inhibition

LPS (kbl) and butyrate in the lumen (Bout=kBo) concentrations are modelled as linear diffusion parameters since they are assumed constant. Also NF- κB (N) is recruited from the cytosol by the term (N0N) which we normalize as in [57] to (1−N) in Eq. 5. Butyrate (B) acts as an inhibitor for the NF- κB transport into the nucleus, by a linear term kbn (see Eq. 5). NF- κB is activated through cytokines and itself activating cytokines in turn. Butyrate in the epithelial cell can be provided by the reservoir (kBo) in the gut lumen, modelled as a linear diffusion process with parameter kd, or by a transporter (first term on r.h.s. of Eq. 9), with a reduction caused by L via T, with the parameters ktl (Eq. 7, last term r.h.s.) and kbt (Eq. 9, first term r.h.s.). The influx of L will be reduced by the blocking effect of butyrate on the leaky barrier with parameter klb (Eq. 8, first term r.h.s.). kB is a linear degradation rate for butyrate in the lumen, which we assume to be constant over time. For the inhibition of T on B (inflammatory negative effect on butyrate transport) a Hill function of second order is used (see Eq. 9, first term r.h.s.), as this allows to consider possible cooperative effects. The same consideration holds for Eq. 1.

This is now our core model:

$$ {}\dot{N} = k_{a} \ \frac{T^{3}}{T^{3}+k_{s}^{3}} \ (1-N) - (k_{ai} \ R + k_{bn} \ B) \ \frac{N}{N+k_{r}} $$
$$ \dot{R} = k_{b} \ N - k_{br} \ R $$
$$ \dot{T} = k_{p} \ \frac{N^{2}}{N^{2}+k_{n}^{2}} - k_{t} \ T + k_{tl} \ L $$
$$ \dot{L} = k_{bl} \ \frac{k_{lb}^{2}}{B^{2}+ k_{lb}^{2}} - k_{l} \ L $$
$$ \dot{B} = k_{Bo} \ \frac{k_{bt}^{2}}{T^{2}+k_{bt}^{2}} + k_{d} \ k_{Bo} - k_{B} \ B $$

In Eqs. 5 to 9, the ODE system for the respective components in Fig. 1 is defined. While Eqs. 5 to 8 reflect essentially the model by [57], Eq. 9 represents its extension with butyrate-related interactions. In Eq. 5, a limitation factor in the last two terms was added, to avoid negative values of N. According to [33], in Eq. 9, a linear dependency upon the butyrate influx to the epithelial cell was assumed, as mediated by the transporter (e.g. MCT1; first term r.h.s.) and also by the diffusion through the epithelial cell (second term r.h.s.). This is due to enhanced expression of transporter proteins, as a regulatory effect by butyrate itself [33]. The model parameters are partly already reduced in number to avoid unnecessary redundancy. Scaling of variables is not considered here, since the model allows every absolute value for the variables in general, except for N.

Fixed point analysis of the core model

A fixed point analysis of our core model, as defined in Eqs. 5 to 9, was conducted. From Fig. 2a, kBo appears as a suitable parameter for the switch between inflamed state (B low) and healthy state (B high). Around the average value (here kBo≈130) a bistable system exists. There are three fixed points of which the upper and lower fixed points are stable. The middle fixed point separates both the upper and lower layer of the hysteresis, and hence is unstable. Below a threshold (kBo≈ 55) the system changes to a low concentration of butyrate within the epithelial cell (B low; inflamed state), and above a threshold (kBo≈ 200) the healthy within-cell concentration can be recovered (non-inflamed state).

Fig. 2
figure 2

Hysteresis curves as evoked by parameters butyrate in the lumen or LPS in the lumen. a Hysteresis driven by lumen butyrate, as observed for the core model. Abscissa: butyrate in the gut lumen kBo, ordinate: butyrate in the epithelial cell B. The hysteresis curve is a plot of all fixed points, i.e. no further change in concentrations: \(\dot {N}=\dot {R}=\dot {T}=\dot {L}=\dot {B}=0\), of Eqs. 5 to 9, for each value of kBo. Parameters: ka=12,ks=1.0,kai=5,kr=0.5,kbn=4.7,kb=5,kbr=0.5,kn=0.2,kp=7,kt=0.2,ktl=2.8,kbl=1.02,klb=3.4,kl=0.7,kbt=2.1,kd=0.125,kB=2.9. The limit points on the hysteresis of the core model have normal form coefficients a left/right: a=−0.1079160/−0.01114689. b Hysteresis driven by lumen LPS kbl, as observed for the core model. Abscissa: LPS in the gut lumen kbl, ordinate: butyrate in the epithelial cell B. Parameters as for (a), but with fixed kBo=130. c Schematic view of the transformation of the butyrate distribution in the lumen over the hysteresis to the bimodal distribution of within-cell butyrate, as depicted on the y-axis. d Hysteresis observed for the extended model. This model includes a positive feedback loop between butyrate in the epithelial mucin production SCFA producing microbiota. The area of bi-stability is broadened for this model

Figure 3 shows the model dynamics as driven by pulsed changes of the butyrate concentration in the intestinal lumen (kBo). Here, driving means that a change in the investigated parameter can result in the switching of the bistable system from one of its stable fixed points to the other. Such a change can be observed from one fixed point of the steady state system to the other, mediated by kBo-pulses in both directions (brown line). Bout=kBo can trigger the switch between both fixed points around Bout=130. Above a threshold (kBo=200 between time =50-170) the non-inflamed state is triggered, leading to a significant reduction of the variables N, R, T, and L. Below a threshold (kBo=50 with time = 200 - 270) the system switches back to an inflamed state, where N, R, T and L are established on higher values again. The following special fixed points were found: Inflamed state fixed point at N =0.178369, R = 1.783690, T =20.365623, L = 0.347144, B =6.079745; non-inflamed state fixed point at N =0.00002003, R = 0.00020025, T = 0.09248014, L = 0.00660570, B =50.382992.

Fig. 3
figure 3

Bistable dynamics of the core model as driven by lumen butyrate. Pulsing of the butyrate concentration in the lumen kBo around an average value and its influence on the switch from the inflamed state to the healthy state. The brown line represents the concentration of butyrate in the lumen kBo and N (green), R (red), T (blue), L (yellow), B (orange) lines, respectively

A constant average level of butyrate supply (by bacterial fermentation processes) is assumed available in the lumen, for the state space region of bi-stability of the core model. In such a situation, butyrate concentration in the intestinal lumen (kBo) is a proper parameter for a switch between both observed fixed points around an intermediate “normal” operating mode.

We also investigated if lumen LPS, kbl, could act as driver of the core model’s bistability. This however is only partly the case, as shown in Fig. 2b. When starting with high levels of within-cell butyrate, increasing levels of lumen LPS lead to a switch of the system to an inflamed state, with small within-cell butyrate values. However, there is no possibility to use decreasing lumen LPS values to re-establish higher within-cell butyrate values and a non-inflamed state again (as the necessary negative values for this parameter are not biologically defined). For these analyses, lumen butyrate was kept stable at kBo=130.

The qualitative bistable behaviour of the core model (5 to 9) is stable under variations of parameters in general, as shown in more detail below. The boundary of the N,R,T,L,B-system is repulsive due to self-limitation so that the hysteresis is generally attracting. Infinity is repulsive and the origin is repulsive such that at least one finite equilibrium always exists. To obtain more realistic relations between inflamed state (NF- κB concentration high) and non-inflamed state (NF- κB concentration low), the experimental results from [3] were used, with a ratio of inflamed/non-inflamed NF- κB levels in the nucleus of ≈100:1. To obtain comparable values in the core model, parameter choices were adapted for the model part based on the work by [57].

Additional file 1: Figures S1 A and S1 B show how NF- κB (N) or the cytokines (T) are effected by lumen butyrate driven hysteresis respectively.

Transformation of an average Gaussian distribution with the hysteresis results in bi-modality

In Fig. 2c it is exemplified how the observed hysteresis in our model could be used to transform an average Gaussian distribution of lumen butyrate levels into a bi-modal distribution of within-cell butyrate levels: Starting with the hysteresis from model 5 to 9, with the three layers (upper (blue)- middle (red) and lower (green)), note that the middle layer is unstable whereas the upper and lower layers are stable. We assume a Gaussian distribution of the butyrate distribution (p(x), pink curve) in the average gut lumen. The transformation of p(x) over the upper and lower layer of the hysteresis y=h(x) to the new distribution of butyrate in the epithelial cell, q(y), with Eq. 10, gives for each layer of the hysteresis a single peak in the corresponding colour (lightblue for the upper layer, lightgreen for the lower layer). Both transformed distributions (lightblue and lightgreen) constitute the overall bi-modal distribution of the butyrate concentration in the epithelial cell. Grey horizontal and vertical lines show corresponding sections for the transformation process. It is assumed that upper and lower layers occur with the same probability. A proper fitting of the curve would have to rely on fitting the probability for the upper or lower stable layer using experimental data.

Transforming a probability distribution p(x) by a function h(x)=y we get a new distribution q(y) by:

$$ q(y) = \frac{p(x)}{h'(x)} \, \ \ \ where \ \ y=h(x). $$

Here we can set h(x) as the hysteresis, p(x) as the Gaussian distribution. Since the hysteresis is not everywhere unique we have to split the transformation for the upper layer and lower layer and weight each of them by an appropriate factor if necessary. Furthermore, the hysteresis has two points x1,x2 where h(xi) (i{1,2}) tends to infinity. These points are biologically not stable since a jump towards the stable layer is more likely to happen than a fixation at this point. These points can therefore be neglected for our computations.

Consequences for the predicted bi-modal within-cell butyrate distribution in a population become first evident when considering our proposed extended model now in the following section.

Extended model implementation and analysis

An extension of the core model is proposed, with feedback relations describing interactions between mucin and butyrate-producing bacteria, which colonise the mucus layer of the gut epithelium. The fermentation products of SCFA-producers, among them butyrate, are indispensable for functioning of gut epithelial cells. The gut epithelium produces a mucin layer. This mucin layer is the place where many SCFA-producing bacteria colonise [14, 22]. Here are the three additional interactions: butyrate → mucin layer High levels of butyrate in the epithelial cells promote mucin layer production [2, 14]. mucin layer → SCFA producers A more pronounced mucin layer promotes colonisation by SCFA producers. SCFA producers → butyrate SCFA-producing bacteria produce lumen butyrate.

These three additional interactions constitute a positive feedback-loop for butyrate inside the epithelial cells. As there is a two-step relation (butyrate → mucin → mucin-adhered bacteria), an extra quadratic term with necessary self-limitation is used. The quadratic term is justified by the double feedback (mucin & microbiota) and the self-limitation due to crowding limitation. For implementation of these additional interactions, a modified version of Eq. 9 is proposed as follows:

$$\begin{array}{@{}rcl@{}} \dot{B} = & k_{Bo} \ \frac{k_{bt}^{2}}{T^{2}+k_{bt}^{2}} \cdot \frac{B^{2}}{\left(B^{2}+k_{2B}^{2}\right)} + k_{d} \ k_{Bo} - k_{B} \ B \end{array} $$

Here, k2B represents the threshold for the Hill curve. This extra term also models the fact that at zero level of butyrate there is no recovering of the butyrate system, which is neglected here, since a completely zero level is not biologically reasonable, and the origin is repulsive.

Fixed-point analysis results for the proposed extended model turned out similar to those of the core model. However, Fig. 2d demonstrates an increased separation of the bistable layers due to the additional feedback.

Considering the transformation of a gut lumen Gaussian distribution of butyrate towards a bimodal within-cell distribution within a population, these findings can now be interpreted with regard to the proposed extended model: a negative correlation between inflammation and the thickness of the mucin layer is hypothesized. The higher the inflammation markers, the thinner the mucin layer and vice versa. This corresponds to the initially defined interaction – a positive correlation – between within-cell butyrate and the mucin layer [2, 14]. Consequently, a bi-modal distribution of within-cell butyrate can be hypothesized to result in a bi-modal distribution of mucin layer thicknesses, and, thus, could result in a bi-modal distribution of SCFA-producing bacteria, within a population.

In depth mathematical analyses of the core model

Existence of interior fixed points and stability of the hysteresis

In this section, we elucidate the more theoretical background of our statement about robust hysteresis, and thereby bistability, in the investigated system (here: the core model). We show that at least one fixed point exists in the positive orthant N,L,R,T,B≥0, in Fig. 2a, Eqs. 5 - 9. We denote the maximum and minimum of all constants k from Eqs. 5 to 9 by:

$$\begin{array}{@{}rcl@{}} u: \ &=&\ \max{k_{\ast} \ >0,} \end{array} $$
$$\begin{array}{@{}rcl@{}} d: \ &=& \ \min{k_{\ast}}\ >0. \end{array} $$

For lumen butyrate concentration, we assume that kBo>0. Then, the following limit relations for N,L,R,T,B>0 can be observed:

  1. (1)


    1. (a)

      case N→1: \(\Rightarrow \dot {N} = - \left (k_{ai} \ R + k_{bn} \ B\right) < 0\)

    2. (b)

      case N→0: \(\Rightarrow \dot {N} = k_{a} \ \frac {T^{3}}{T^{3}+k_{s}^{3}} > 0\)

  2. (2)


    1. (a)

      case R: \(\Rightarrow \dot {R} \sim - k_{br} \ R < 0\)

    2. (b)

      case R→0: \(\Rightarrow \dot {R} = k_{b} \ N > 0\)

  3. (3)


    1. (a)

      case T: \(\Rightarrow \dot {T} \sim - k_{t} \ T < 0\)

    2. (b)

      case T→0: \(\Rightarrow \dot {T} = k_{p} \ \frac {N^{2}}{N^{2}+k_{n}^{2}} + k_{tl} \ L > 0\)

  4. (4)


    1. (a)

      case \(L > \frac {k_{bl}}{k_{l}} \)\(\Rightarrow \dot {L} < 0\)

    2. (b)

      case L→0: \(\Rightarrow \dot {L} = k_{bl} \frac {k_{lb}^{2}}{B^{2}+ k_{lb}^{2}} > 0\)

  5. (5)


    1. (a)

      case B: \(\Rightarrow \dot {B} \sim - k_{B} \ B < 0\)

    2. (b)

      case B→0: \(\Rightarrow \dot {B} = k_{Bo} \ \frac {k_{bt}^{2}}{T^{2}+k_{bt}^{2}} + k_{d} \ k_{Bo} > 0\)

We construct a box P, which is mapped into itself under the ODE of Eqs. 5 to 9. Hence by Brouwer’s fixed point theorem [4] there exists at least one fixed point in P°. We show that there is no fixed point on the boundary of the box P. First observe that for

$$ N,L,R,T,B \ \leq \ \epsilon\ > \ 0 \ $$

with ε sufficiently small, both \(\dot {L}, \dot {B} > 0\) since d>0. Hence, there is no fixed point in the ball Bε=(N,R,T,L,B)ε around the origin {N,R,T,L,B}=0. We set \(L_{0}=\frac { k_{bl}}{k_{l}}+\delta \), where δ>0. Let K>L0 be large enough such that for LL0 case (a) from limit relations 2, 3 and 5 are satisfied for R,T,BK. Based on these assumptions, we analyse the case, where B,Lε and the different cases with N,R,T→0:

$$\begin{array}{@{}rcl@{}} T=0, & N,R >0 & \Rightarrow \dot{T} >0 \end{array} $$
$$\begin{array}{@{}rcl@{}} N=0, & R,T >0 & \Rightarrow \dot{N} >0 \end{array} $$
$$\begin{array}{@{}rcl@{}} R=0, & N,T >0 & \Rightarrow \dot{R} >0 \end{array} $$
$$\begin{array}{@{}rcl@{}} R,N=0, & T >0 & \Rightarrow \dot{N} >0, \dot{R} =0 \end{array} $$
$$\begin{array}{@{}rcl@{}} R,T=0, & N >0 & \Rightarrow \dot{T} >0, \dot{R} >0 \end{array} $$
$$\begin{array}{@{}rcl@{}} N,T=0, & R >0 & \Rightarrow \dot{N} =0, \dot{T} >0 \end{array} $$
$$\begin{array}{@{}rcl@{}} N,T,R=0, & & \Rightarrow \dot{N} =0, \dot{T} >0, \dot{R} =0 \end{array} $$

From relations 15 to 21 it can be deduced that the positive orthant N,L,R,T,B≥0 is invariant and that there is no fixed point on its border. Furthermore, we can now determine the box P by:

$$\begin{array}{@{}rcl@{}} N, R, T, L, B &\geq& 0 \end{array} $$
$$\begin{array}{@{}rcl@{}} L &\leq& L_{0} \end{array} $$
$$\begin{array}{@{}rcl@{}} N &\leq& 1 \end{array} $$
$$\begin{array}{@{}rcl@{}} R, T, B &\leq& K \ . \end{array} $$

From Eqs. 23 to 25 and cases (a) from limit relations 1 to 5, we can see that the upper border of P is repulsive since the vectorfield defined by Eqs. 5 to 9 points into Po. Since δ was not bound from above, the box P can be arbitrarily large in the positive orthant except for N which is bound by definition. Hence, the previous numerically obtained bifurcation curve from Fig. 2a has at least one fixed point. In order to investigate the special form of the hysteresis we have to use a continuation program and evaluate its stability. For the parameter set of Fig. 2a, we can give the sign of the normal form coefficient a of the normal form approximation at the left and right limit points with active parameter kBo (values for a see Fig. 2a, as computed by matcont [29]). Both are negative (a<0), indicating a stable top and bottom layer of the hysteresis. There are no other bifurcation points for the hysteresis for the used parameters. For other choices of parameters, we could, in principle, also get stability changes within the hysteresis layers, indicating bifurcation points with emerging limit cycles. Brouwer’s fixed point theorem only provides the existence for one or more fixed points in the positive orthant for our model in general. From a biological perspective, these conclusions can be interpreted as statements about the robustness of bistability in the investigated system and their necessary preconditions. These preconditions are motivated by experimental evidence.

Limit behaviour for various parameters

As in Fig. 2a we can see that for large kBo the upper layer of the hysteresis converges to a line through the origin after being bowed away from it by the presence of inflammation for low kKo values. A short limit analysis for kBo is available through:

$$\begin{array}{*{20}l} {}\dot{B}= 0 & \ \Rightarrow \ B= \frac{k_{Bo}}{k_{B}}\left(\frac{k_{bt}^{2}}{T^{2}+k_{bt}^{2}} + k_{d}\right),\\ &\qquad\qquad for \ \ \ k_{Bo} \rightarrow \infty \ \Rightarrow \ B \rightarrow \infty \\ \end{array} $$
$$\begin{array}{*{20}l} {}\dot{L}=0 & \ \Rightarrow \ L= \frac{k_{bl} \cdot k_{lb}^{2}}{\left(B^{2}+ k_{lb}^{2}\right)\cdot k_{l}}, & for \ \ \ B \rightarrow \infty \ \ \Rightarrow L \rightarrow 0 \\ \end{array} $$
$$\begin{array}{*{20}l} \dot{N}=0 & \Rightarrow \dot{N} \leq k_{a} - k_{bn} B \frac{N}{N+k_{r}},& \ for \ B \rightarrow \infty \! \Rightarrow\! N \!\rightarrow \! 0 \\ \end{array} $$
$$\begin{array}{*{20}l} \dot{R}=0 & \ \ \Rightarrow{R} =N \frac{k_{b}}{k_{br}}, & \ \ for \ B \rightarrow \infty \ \& \ (1.8) \ \ \Rightarrow \ \ R \rightarrow 0 \\ \end{array} $$
$$\begin{array}{*{20}l} {}\dot{T}=0, & \ \ \ B \rightarrow \infty \ \& \ (1.7)\ \&\ (1.8) & \Rightarrow \ \dot{T} \leq 0 \ \ \Rightarrow \ \ T \rightarrow 0 \\ \end{array} $$
$$\begin{array}{*{20}l} {}\dot{B} = 0, &\quad k_{Bo} \rightarrow \infty \ \& \ (1.6) - (1.10)&\Rightarrow B \approx k_{Bo} \ \frac{\left(1 + k_{d}\right) }{k_{B}} \ \ \end{array} $$

For large values of lumen butyrate, kBo, the upper layer of the hysteresis follows a negative linear trend, with slope \( \frac {(1 + k_{d}) }{k_{B}} \), which is larger than that at the origin with ongoing inflammation. This results in the bi-modal distribution, as explained in Fig. 2c, with different height and width as predicted by Eq. 10. Biologically, this translates again as a statement about robust bistability in our investigated system, as being rather independent of parameter choices in our modelling approach.


Robust bistability from the interplay of two antagonists: inflammation and butyrate

The cross-talk between butyrate and the inflammatory state of human gut epithelium cells comprises complex interactions, including feedback between different elements, and therefore constitutes a problem eligible for a systems biological analysis. We built a core ODE dynamical model for this set of interactions based on previous knowledge and additional terms. Choice of parameters for our ODE model was based both on previous modelling work [57] and further prior knowledge, but no dynamical data, rendering all results essentially qualitative in character. Our analyses showed a clear bistable behaviour for the system, characterised by a hysteresis-like fixed point curve (Figs. 2 a, b and d). The found bistability is largely independent of the specific choices of dynamic parameters. Our model predicts that butyrate concentration in the gut (kBo) is a possible driver for the switch between the two stable, non-zero fixed points in the system, a non-inflamed state with high inner butyrate levels as opposed to the inflamed state with low inner butyrate concentrations. However, other model-specific parameters can also control this switch-like behaviour, for example, lumen LPS. Essentially, the system, even in its extended form, consists of two antagonistic parts, inflammation and butyrate. Each of these parts suppresses the other, resulting in the observed robust bistability. Both lumen butyrate and lumen LPS can act as drivers of this system, determining which of the two switch-like states (the two stable non-zero fixed points) would dominate. However, only butyrate can drive the system both into the inflamed state, by decreasing lumen butyrate levels, and back into the non-inflamed state, by sufficiently increasing lumen butyrate levels, at least during a pulse.

Inflammation is modelled qualitatively, as a simplistic interplay between three species, NF- κB (N), the cytokines (T), and a repressor (R) with only one non-zero stable fixed point. In our modelling approach these species are of comparatively re-active character, responding to, rather than driving, the switch between inflamed or non-inflamed status. However, the modelled system is analysed as if cut-out from the larger regulatory network of the organ, here, the human gut. In reality, there might be additional regulatory interactions, and certain situations in which these get strong and decisive. Other than NF- κB, which is an internal cellular signal, cytokines of different kinds also act as inter-cellular signals. Such pro-inflammatory or anti-inflammatory cytokines are known to originate from other places in the body, and, by regulating the NF- κB inflammatory system, also drive the system towards an inflamed (pro-inflammatory cytokines) or an uninflamed (anti-inflammtory cytokines) state. A virus or bacterial infection of the gut, for example, could lead to additional signals resulting in high NF- κB levels, and hence drive the system towards inflammation. Biologics treatment of inflammation via anti-TNF antibodies, an upstream regulator of NF- κB, is another example for an alternative driver in the analysed system [35]. In summary, many elements of the system have the potential to induce a bi-stable behaviour. However, our simplified model shows that lumen butyrate on its own, as well as lumen LPS, can drive a bi-stable behaviour. More data is needed to elucidate the behaviour of the system under a wide range of physiological conditions.

The proposed core model constitutes a very simplified account of the interaction between inflammation and fermentation products from gut bacteria (butyrate in particular). It cannot therefore serve to generate quantitative conclusions about concentration thresholds or time-scales. However, from our point of view, it still can provide useful insights into regulatory mechanisms of inflammation in the gut and its potential consequences. As opposed to the model in [57] that shows one non-zero fixed point, our core model has two non-zero stable fixed points (and an unstable third one). The observed bi-stability would be lost if the system was parametrised in a way such that interactions in the inflammation-regulatory subsystem (N, T, R) would be significantly damped or lost. Experimental evidence suggests that the modelled interactions are not marginal in a real scenario and, therefore, our model provides a simplistic but also realistic account of this system.

Our model also links to possible ways to interfere with the two observed antagonistic states. We refer to treatment of inflammation with butyrate enemas as a therapy (compare e.g. [51]) or appropriate dietary interventions, e.g. butyrate-promoting dietary fibres (compare e.g. [9]), especially for situations of low-grade inflammation.

Model-based insights on the bimodality of the gut microbiota

Assuming an average Gaussian distribution of gut lumen butyrate concentrations in a collection of systems (cohort, population), our core model predicts a bimodal transformation for the epithelial cell inner butyrate levels. Reasoning about possible consequences of this characteristic of our modelled system leads over to a more speculative part of our work, discussed in the following. The proposed extended model builds on additional interactions that are hypothesized to operate in the gut-microbiota system. Experimental evidence indicates that an inflamed gut epithelium produces less of a mucus layer, and also results in a damaged mucus layer [21]. As butyrate counteracts an inflamed state in epithelial cells, it might indirectly promote a more pronounced/intact mucus layer. Butyrate is also an important energy source for gut epithelial cells [23] and might therefore indirectly increase mucin production. Important SCFA-producing bacteria are dependent on the mucus layer as their habitat [11]. Hence, these proposed interactions result in a positive feedback loop from gut lumen butyrate, via butyrate within the cell, and back, mediated by the mucus layer. This positive feedback loop would also constitute an interaction to within-cell butyrate levels, as their control of inflammation allows mucus production, and thereby influence indirectly microbiota composition in the gut lumen. Simulations in [46] showed a positive feedback between mucin producing epithelial cells and mucin feeding bacteria. This line of argument leads us to propose a possible explanation for the bimodality of gut microbiome species compositions as described in a recent work by [30]: A possible effect of the bi-stability in the intestinal inflammatory system could have consequences also for the composition of the gut microbiome. In [30] five taxonomic groups were classified as bimodal with regard to their abundance distribution. Also [18] described an essential part of mucin-related microbiota as bimodally distributed. Some of these are SCFA- and specifically butyrate-producers [11, 37, 53, 56]. Moreover, even the key species of the proposed three enterotypes by [1] correspond to the bimodal types as observed in [30].

Conclusions and outlook

Our model shows robust bistability as a result of the feedback between two antagonists, microbial-derived butyrate and gut epithelial inflammation. These observations warrant further mathematical and experimental studies on the emergence of bistability. Additional experimental data is required to parametrise our modeling approach, and would need to involve a dynamical screening of the concentration of key components or derived variables thereof. In vitro experiments with gut epithelial cell-derived cell cultures are the basis to further understand key pathways and signals involved in the butyrate-inflammation interplay. Organoid studies on the same topic are started, as well as ex vivo studies involving Ussing chamber experiments [16], where barrier function can directly be related to inflammatory signals. Finally, in the clinical setting, carefully managed butyrate bolus experiments, as well as tightly controlled fibre-rich diets, are already under way. All these levels of experimental evidence are necessarily dynamic in nature as to enable systems biological analysis of mechanistic hypotheses.


Aim, design and software for analysis

The aim of the study is to manifest the fundamental antagonistic behaviour between inflammation of the gut epithelial cell layer and specific bacterial products in the large intestine, represented here by butyrate. We analyse the ODE-model based on [57] extended by suitable interaction terms with butyrate, as described in detail in the “Results” section, model development. First, we observe by simulations with Matlab [34] a bistable behaviour, visible as hysteresis between antagonistic variables in our mathematical model. A switch-like behaviour is modelled also by Simulink (available as add-on for Matlab) showing the operation of the hysteresis. A continuation analysis is implemented using the matcont (package for Matlab) [29], giving bifurcation points for determining the stability of the hysteresis. A plausible stability analysis for the variables is given, indicating the global stability of the top and bottom layer of the hysteresis. A probability transformation of an assumed Gaussian distribution of butyrate in the lumen, mediated by the hysteresis, results in a qualitative similarity to observed bimodal distributions within the gut intestinal microbiome as observed in [30].

Availability of data and materials

Simulation and analysis scripts are made available as tar-archive (Additional file 2) within the section Additional files.


  1. Arumugam M, Raes J, Pelletier E, Paslier DL, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto J-M, Bertalan M, Borrue N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S, Torrents D, Ugarte E, Zoetendal EG, Wang J, Guarner F, Pedersen O, de Vos WM, Brunak S, Dore J, Consortium M, Weissenbach J, Ehrlich SD, Bork P. Enterotypes of the human gut microbiome. Nature. 2011; 473:174–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Barcelo A, Claustre J, Moro F, Chayvialle JA, Cuber JC, Plaisancie P. Mucin secretion is modulated by luminal factors in the isolated vascularly perfused rat colon. Gut. 2000; 46:218–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Blaecke A, Delneste Y, Herbault N, Jeannin P, Bonnefoy J-Y, Beck A, Aubry J-P. Measurement of Nuclear factor-kappa B translocation on lipopolysaccharide-activated human dendritic cells by confocal microscopy and flow cytometry. Cytometry. 2002; 48:71–9.

    Article  CAS  PubMed  Google Scholar 

  4. Brouwer LEJ. Über Abbildung von Mannigfaltigkeiten. Math Ann. 1912; 71(4):598.

    Google Scholar 

  5. Canfora EE, Jocken JW, Blaak EE. Short-chain fatty acids in control of body weight and insulin sensitivity. Nat Rev Endocr. 2015; 11:577–91.

    Article  CAS  Google Scholar 

  6. Chassaing B, Gewirtz AT. Gut microbiota, low-grade inflammation, and metabolic syndrome. Toxicol Pathol. 2014; 42(1):49–53.

    Article  PubMed  Google Scholar 

  7. Claesson MJ, Jeffery IB, Conde S, Power SE, O’Connor EM, Cusack S, Harris HM, Coakley M, Lakshminarayanan B, O’Sullivan O, Fitzgerald GF, Deane J, O’Connor M, Harnedy N, O’Connor K, O’Mahony D, van Sinderen D, Wallace M, Brennan L, Stanton C, Marchesi JR, Fitzgerald AP, Shanahan F, Hill C, Ross RP, O’Toole PW. Gut microbiota composition correlates with diet and health in the elderly. Nature. 2012; 488(7410):178–84.

    Article  CAS  PubMed  Google Scholar 

  8. Corrêa R-O, Fachi JL, Vieira A, Sato FT, Vinolo MAR. Regulation of immune cell function by short-chain fatty acids. Clin Transl Immunol. 2016; 5:73.

    Article  Google Scholar 

  9. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, Ling AV, Devlin AS, Varma Y, Fischbach MA, Biddinger SB, Dutton RJ, Turnbaugh PJ. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014; 505(7484):559–63.

    Article  CAS  Google Scholar 

  10. Davie JR. Inhibition of histone deacetylase activity by butyrate. J Nutr. 2003; 133(7):2485–93.

    Article  Google Scholar 

  11. den Abbeele PV, Belzer C, Goossens M, Kleerebezem M, Vos WMD, Thas O, Weirdt RD, Kerckhof1 F-M, de Wiele TV. Butyrate-producing clostridium cluster XIVa species specifically colonize mucins in an in vitro gut model. ISME J. 2013; 7:949–61.

  12. Enesa K, Zakkar M, Chaudhury H, Luong LA, Rawlinson L, Mason JC, Haskard DO, Dean JLE, Evans PC. Nf-kappaB suppression by the deubiquitinating enzyme Cezanne: a novel negative feedback loop in pro-inflammatory signaling. J Biol Chem. 2008; 283:7036–45.

    Article  CAS  PubMed  Google Scholar 

  13. Filippo CD, Cavalieri D, Paola MD, Ramazzotti M, Poullet JB, Massart S, Collini S, Pieraccini G, Lionetti P. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci USA. 2010; 107(33):14691–6.

    Article  PubMed  Google Scholar 

  14. Finnie IA, Dwarakanath AD, Taylor BA, Rhodes JM. Colonic mucin synthesis is increased by sodium butyrate. Gut. 1995; 36:93–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Flint HJ, Scott KP, Duncan SH, Louis P, Forano E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes. 2012; 3(4):289–306.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ganda-Mall JP, Löfvendahl L, Lindqvist CM, Brummer RJ, Keita A, Schoultz I. Differential effects of dietary fibres on colonic barrier function in elderly individuals with gastrointestinal symptoms. Sci Rep. 2018; 8(1):13404.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Goncales P, Martel F. Butyrate and colorectal cancer: The role of butyrate transport. Curr Drug Metab. 2013; 14(9):994–1008.

    Article  Google Scholar 

  18. Hakansson A, Molin G. Gut microbiota and inflammation. Nutrients. 2011; 3:637–82.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer R-J. Review article: the role of butyrate on colonic function. Aliment Pharmacol Ther. 2008; 15(2):104–19.

    Article  Google Scholar 

  20. Hamer HM, Jonkers DMAE, Vanhoutvin SALW, Troost FJ, Rijkers G, de Brune A, Bast A, Venema K, Brummer RJ. Effect of butyrate enemas on inflammation and antioxidant status in the colonic mucosa of patients with ulcerative colitis in remission. Clin Nutr. 2010; 29(6):738–44.

    Article  CAS  PubMed  Google Scholar 

  21. Hansson GC. Role of mucus layers in gut infection and inflammation. Curr Opin Microbiol. 2012; 15(1):57–62.

    Article  CAS  PubMed  Google Scholar 

  22. Hatayama H, Iwashita J, Kuwajima A, Abe T. The short chain fatty acid, butyrate, stimulates MUC2 mucin production in the human colon cancer cell line, LS174T. Biochem Biophys Res Commun. 2007; 356:599–603.

    Article  CAS  PubMed  Google Scholar 

  23. Jung T-H, Park JH, Jeon W-M, Han K-S. Butyrate modulates bacterial adherence on LS174T human colorectal cells by stimulating mucin secretion and MAPK signaling pathway. Nutr Res Pract. 2015; 9(4):343–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Karkman A, Lehtimäki J, Ruokolainen L. The ecology of human microbiota: dynamics and diversity in health and disease. Ann N Y Acad Sci. 2017; 1399(1):78–92.

    Article  PubMed  Google Scholar 

  25. Kau AL, Ahern PP, Griffin NW, Goodman AL, Gordon JI. Human nutrition, the gut microbiome, and immune system: envisioning the future. Nature. 2012; 474(7351):327–36.

    Article  Google Scholar 

  26. Koh A, De Vadder F, Kovatcheva-Datchary P, Bäckhed F. From dietary fiber to host physiology: Short-chain fatty acids as key bacterial metabolites. Cell. 2016; 165:1332–45.

    Article  CAS  Google Scholar 

  27. König J, Brummer RJ. Alteration of the intestinal microbiota as a cause of and a potential therapeutic option in irritable bowel syndrome. Benef Microbes. 2014; 5(3):247–61.

    Article  PubMed  Google Scholar 

  28. Kostic AD, Xavier RJ, Gevers D. The microbiome in inflammatory bowel diseases: Current status and the future ahead. Gastroenterol. 2014; 146(6):1489–99.

    Article  CAS  Google Scholar 

  29. Kuznetsov YA. Elements of Applied Bifurcation Theory. New York: Springer; 2004.

    Book  Google Scholar 

  30. Lahti L, Salojärvi J, Salonen A, Scheffer M, de Vos WM. Tipping elements in the human intestinal ecosystem. Nat Commun. 2014; 5(4344):1–10.

    Google Scholar 

  31. Larsen N, Vogensen FK, van den Berg FWJ, Nielsen DS, Andreasen AS, Pedersen BK, Al-Soud WA, Sørensen SJ, Hansen LH, Jakobsen M. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PloS ONE. 2010; 5:9085.

    Article  Google Scholar 

  32. Macfarlane S, Macfarlane GT. Regulation of short-chain fatty acid production. Proc Nutr Soc. 2003; 62(1):67–72.

    Article  CAS  PubMed  Google Scholar 

  33. Mathewson ND, Jenq R, Mathew AV, Koenigsknecht M, Hanash A, Toubai T, Oravecz-Wilson K, Wu S-R, Sun Y, Rossi C, Fujiwara H, Byun J, Shono Y, Lindemans C, Calafiore M, Schmidt TC, Honda K, Young VB, Pennathur S, van den Brink M, Reddy P. Gut microbiome-derived metabolites modulate intestinal epithelial cell damage and mitigate graft-versus-host disease. Nature. 2016; 17(5):505–15.

    CAS  Google Scholar 

  34. MathWorks T. MATLAB and Statistics Toolbox Release 2012b. Natick: The MathWorks, Inc.; 2012.

  35. Monaco C, Nanchahal J, Taylor P, Feldmann M. Anti-TNF therapy: past, present and future. Int Immunol. 2015; 27:55–62.

    Article  CAS  PubMed  Google Scholar 

  36. Moreno-Indias I, Cardona F, Tinahones FJ, Queipo-Ortuño MI. Impact of the gut microbiota on the development of obesity and type 2 diabetes mellitus. Front Microbiol. 2014; 5:190.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Patrick S, Blakely GW, Houston S, Moore J, Abratt VR, Bertalan M, Cerdeño-Tárraga AM, Quail MA, Corton N, Corton C, Bignell A, Barron A, Clark L, Bentley SD, Parkhill J. Twenty-eight divergent polysaccharide loci specifying within- and amongst-strain capsule diversity in three strains of Bacteroides fragilis. Microbiol. 2010; 156(Pt 11):3255–69.

    Article  CAS  Google Scholar 

  38. Peng L, Li Z-R, Green RS, Holzman IR, Lin J. Butyrate enhances the intestinal barrier by facilitating tight junction assembly via activation of AMP-activated protein kinase in Caco-2 cell monolayers. J Nutr. 2009; 139(9):1619–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Piche T, Barbara G, Aubert P, des Varannes SB, Dainese R, Nano JL, Cremon C, Stanghellini V, Giorgio RD, Galmiche JP, Neunlist M. Impaired intestinal barrier integrity in the colon of patients with irritable bowel syndrome: involvement of soluble mediators. Gut. 2009; 58(2):196–201.

    Article  CAS  PubMed  Google Scholar 

  40. Plöger S, Stumpff F, Penner GB, Schulzke J-D, Gäbel G, Martens H, Shen Z, Günzel D, Aschenbach JR. Microbial butyrate and its role for barrier function in the gastrointestinal tract. Ann N Y Acad Sci. 2012; 1258:52–9.

    Article  PubMed  Google Scholar 

  41. Pozuelo M, Panda S, Santiago A, Mendez S, Accarino A, Santos J, Guarner F, Azpiroz F, Manichanh C. Reduction of butyrate- and methane-producing microorganisms in patients with Irritable Bowel Syndrome. Sci Rep. 2015; 5:12693.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Preter VD, Arijs I, Windey K, Vanhove W, Vermeire S, Schuit F, Rutgeerts P, Verbeke K. Impaired butyrate oxidation in ulcerative colitis is due to decreased butyrate uptake and a defect in the oxidation pathway. Inflamm Bowel Dis. 2012; 18(6):1127–36.

    Article  PubMed  Google Scholar 

  43. Raos-Covian D, Ruas-Madiedo P, Margolles A, Gueimonde M, de los Reyes-Gavilan CG, Salazar N. Intestinal short chain fatty acids and their link with diet and human health. Front Microbiol. 2016; 10(3389):2016–00185.

    Google Scholar 

  44. Ritzhaupt A, Wood IS, Ellis A, Hosie KB, Shirazi-Beechey S. Identification and characterization of a monocarboxylate transporter (MCT1) in pig and human colon: its potential to transport lactate as well as butyrate. J Physiol. 1998; 513(3):719–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Scheppach W, Sommer H, Kirchner T, Paganelli GM, Bartram P, Christi S, Richter F, Dusel G, Kasper H. Effect of butyrate enemas on the colonic mucosa in distal ulcerative colitis. Gastroenterol. 1992; 103(1):51–6.

    Article  CAS  Google Scholar 

  46. Schluter J, Foster KR. The evolution of mutualism in gut microbiota via host epithelial selection. PLoS Biol. 2012; 10(11):1001424.

    Article  Google Scholar 

  47. Segain JP, Blétiere DRDL, Bourreille A, Leray V, Gervois N, Rosales C, Ferrier L, Bonnet C, Blottiere HM, Galmiche JP. Butyrate inhibits inflammatory responses through NFkB inhibition: implications for Crohn’s disease. Gut. 2000; 47(3):397–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Tay S, Hughey JJ, Lee TK, Lipniacki T, Quake SR, Covert MW. Single-cell NF-kappaB dynamics reveal digital activation and analogue information processing. Nature. 2010; 466:267–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Thibault R, Blachier F, Darcy-Vrillon B, de Coppet P, Bourreille A, Segain JP. Butyrate utilization by the colonic mucosa in inflammatory bowel diseases: a transport deficiency. Inflamm Bowel Dis. 2010; 16(4):684–95.

    Article  PubMed  Google Scholar 

  50. Thibault R, Coppet PD, Daly K, Bourreille A, Cuff M, Bonnet C, Mosnier J-F, Galmiche J-P, Shirazi-Beechey S, Segain J-P. Down-regulation of the monocarboxylate transporter 1 is involved in butyrate deficiency during intestinal inflammation. Gastroenterol. 2007; 133(6):1916–27.

    Article  CAS  Google Scholar 

  51. Vanhoutvin SALW, Troost FJ, Hamer HM, Lindsey PJ, Koek GH, Jonkers DMAE, Kodde A, Venema K, Brummer RJM. Butyrate-induced transcriptional changes in human colonic mucosa. PLoS ONE. 2009; 4(8):6759.

    Article  Google Scholar 

  52. Vernia P, Gnaedinger A, Hauck W, Breuer RI. Organic anions and the diarrhea of inflammatory bowel disease. Dig Dis Sci. 1988; 33(11):1353–8.

    Article  CAS  PubMed  Google Scholar 

  53. Vital M, Howe AC, Tiedjea JM. Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data. mBio. 2014; 5(2):00889–14.

    Article  Google Scholar 

  54. Vrieze A, Van Nood E, Holleman F, Salojärvi J, Kootte RS, Bartelsman JFWM, Dallinga-Thie GM, Ackermans MT, Serlie MJ, Oozeer R, Derrien M, Druesne A, Van Hylckama Vlieg JET, Bloks VW, Groen AK, Heilig HGHJ, Zoetendal EG, Stroes ES, de Vos WM, Hoekstra JBL, Nieuwdorp M. Transfer of intestinal microbiota from lean donors increases insulin sensitivity in individuals with metabolic syndrome. Gastroenterol. 2012; 143:913–9167.

    Article  CAS  Google Scholar 

  55. Wang W, Chen L, Zhou R, Wang X, Song L, Huang S, Wang G, Xia B. Increased proportions of bifidobacterium and the lactobacillus group and loss of butyrate-producing bacteria in inflammatory bowel disease. J Clin Microbiol. 2014; 52(2):398–406.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Weir TL, Manter DK, Sheflin AM, Barnett BA, Heuberger AL, Ryan EP. Stool microbiome and metabolome differences between colorectal cancer patients and healthy adults. PLoS ONE. 2013; 8(8):70803.

    Article  Google Scholar 

  57. Yde P, Mengel B, Jensen MH, Krishna S, Trusina A. Modeling the NF-kB mediated inflammatory response predicts cytokine waves in tissue. BMC Syst Biol. 2011; 5:115.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Zhou Q, Zhang B, Verne GN. Intestinal membrane permeability and hypersensitivity in the Irritable Bowel Syndrome. Pain. 2009; 146(1-2):41–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


Thanks to Olaf Wolkenhauer for his immediate hints on modelling software and Leo Lathi for helpful clarifications. We acknowledge the Nutrition-Gut-Brain Interactions Research Centre (NGBI) at Örebro University, Sweden, for valuable scientific discussion and research infrastructure offered for this project. Two reviewers and the editor assisted with detailed valuable comments and concrete proposals to improve discussion of our results and the manuscript. Stephen Rush kindly reviewed our manuscript with regard to language.


This work has been funded by the Swedish Knowledge Foundation, grant number #20110225.

Availability of data and materials

Part of the dynamic model used in this work has been published by others, as indicated. All data have been simulated using the models developed in this work. Scripts for these simulations are made available in the Additional Files section.

Author information

Authors and Affiliations



DR conceived the study and supervised model development, simulations and analysis. GN devised the project, the main conceptual ideas and proof outline and run all simulations. RW, IR and TM contributed to model development and discussion of model analysis results. All authors worked on finalising the manuscript and have approved the final manuscript.

Corresponding author

Correspondence to Dirk Repsilber.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1

Figure S1. Hysteresis effects driven by lumen butyrate on further variables in the core model. A: for NF- κB (N). B: for cytokines (T). Parameters as for Fig. 2a. (JPG 40 kb)

Additional file 2

MATLAB scripts. tar archive with MATLAB scripts for simulating the hysteresis and for transforming a Gaussian distribution via the simulated hysteresis. (TAR 966 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Neumann, G., Wall, R., Rangel, I. et al. Qualitative modelling of the interplay of inflammatory status and butyrate in the human gut: a hypotheses about robust bi-stability. BMC Syst Biol 12, 144 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: