Tutorials › Quantitative Systems Pharmacology (QSP) › Positive Feedback and Self-Amplifying Inflammation

Inflammation and Autoimmunity · Tutorial 305 of 799

Positive Feedback and Self-Amplifying Inflammation

Learn to distinguish stable amplification from threshold-dependent self-sustaining activation using simple cytokine feedback models and their equilibrium conditions.

Advanced 11 min read

What You'll Learn

  • Formulate cytokine production and clearance as a dynamical feedback model
  • Calculate local loop gain and the resulting amplification of a cytokine response
  • Identify stability conditions for a linearized positive-feedback loop
  • Derive the threshold and bistability conditions for a saturating cooperative loop
  • Assess how cytokine clearance can shift a model between inflammatory regimes
  • Recognize data patterns and modeling assumptions that can mislead feedback analysis

Why Positive Feedback Can Change an Inflammatory Response

A cytokine can do more than act on target cells: it can also increase the production of itself, directly or through intervening cells. That arrangement forms a positive-feedback loop. A small increase in cytokine concentration may then produce a larger increase in its production, potentially amplifying an initiating signal or, under some conditions, sustaining activation after that signal has passed.

The key distinction is between amplification and self-sustaining activation. A loop can increase the size of a response while the response still decays when its initiating input is removed. Sustained activation requires stronger conditions: feedback production must be sufficient relative to cytokine removal over the relevant range of concentrations. Building on cytokine source–sink balance and soluble cytokine antagonism, this tutorial makes that comparison explicit in dynamical models.

Definition: Positive cytokine feedback occurs when an increase in a cytokine state increases a process that further increases that state. In a one-state model, let \(C(t)\) be cytokine concentration, \(q_b\) basal production, \(q_{\mathrm{ext}}(t)\) an external or initiating input, \(q_{\mathrm{fb}}(C)\) feedback-dependent production, and \(kC\) first-order removal. Then $$ \frac{dC}{dt}=q_b+q_{\mathrm{ext}}(t)+q_{\mathrm{fb}}(C)-kC. $$ Here production terms have units of concentration per time, \(k\) has units of inverse time, and \(C\) has units of concentration.

This model assumes a well-mixed compartment, a single cytokine state, constant first-order removal, and a specified production response to \(C\). It treats the feedback pathway as acting without an explicitly modeled delay. Those assumptions make the loop interpretable, but they are not claims that real inflammatory tissue is homogeneous or that its cellular intermediates are instantaneous.

Local Amplification and Loop Gain

Near a steady state \(C^*\), small changes can be analyzed by approximating the feedback production curve by its local slope. At equilibrium, net production equals removal. Define the dimensionless local loop gain \(L\) as the slope of feedback production divided by the removal rate:

$$ L=\frac{q_{\mathrm{fb}}'(C^*)}{k}. $$

If \(L\) is positive but less than one, a small perturbation is amplified yet remains locally stable: deviations from the steady state decay. If \(L\) reaches or exceeds one, the local linear approximation no longer predicts a stable return to that equilibrium. The resulting behavior depends on the full nonlinear production curve, not on loop gain alone.

Result: For a stable equilibrium with \(L<1\), a small constant increase \(\delta q\) in production changes the equilibrium concentration by approximately \(\delta C=\delta q/(k-q_{\mathrm{fb}}'(C^*))\). Relative to the change \(\delta q/k\) without feedback, the amplification factor is \(1/(1-L)\).

This result is local: it describes small perturbations near the stated equilibrium. It does not guarantee that a large perturbation returns to that equilibrium, nor does it establish that the feedback function remains linear over a large concentration range.

Worked Example: Stable Amplification in a Linearized Loop

Suppose a hypothetical tissue has \(q_b=10\) nM/day, first-order removal \(k=1\) per day, and feedback production \(q_{\mathrm{fb}}(C)=gC\), where \(g=0.60\) per day. Assume no time-varying external input. The model is $$ \frac{dC}{dt}=10+0.60C-1.00C. $$

At steady state, set the derivative to zero and solve: $$ 0=10-(1.00-0.60)C^*,\qquad C^*=\frac{10\ \mathrm{nM/day}}{0.40\ \mathrm{day}^{-1}}=25\ \mathrm{nM}. $$ Without feedback, the equilibrium would be \(10/1.00=10\) nM. The positive loop therefore raises the steady concentration by a factor of \(25/10=2.50\).

The local gain is \(L=g/k=0.60/1.00=0.60\), so the equilibrium is stable because \(L<1\). A small additional production input of \(1\) nM/day would change the concentration by approximately $$ \delta C=\frac{1\ \mathrm{nM/day}}{1.00-0.60\ \mathrm{day}^{-1}}=2.50\ \mathrm{nM}, $$ rather than \(1.00\) nM without feedback. This is amplification, not autonomous persistence: if basal production were also removed, this purely linear model with \(g<k\) would decay toward zero.

Nonlinear Feedback and Thresholds

Cytokine-induced production may be nonlinear. A common phenomenological choice is a Hill function, which rises with concentration and approaches a finite maximum:

$$ q_{\mathrm{fb}}(C)=q_{\max}\frac{C^n}{K^n+C^n}. $$

Here \(q_{\max}\) is the maximum feedback-dependent production rate, \(K\) is the concentration at which that production is half-maximal, and \(n\) controls the steepness of the response. The Hill function is a modeling choice; using it does not establish a particular molecular binding mechanism.

To see how a threshold can arise, consider a hypothetical loop with no basal production after the initiating input is removed, and with \(n=2\). Define \(x=C/K\), \(a=q_{\max}/(kK)\), and dimensionless time \(\tau=kt\). The model becomes:

$$ \frac{dx}{d\tau}=a\frac{x^2}{1+x^2}-x. $$

At zero concentration, the state remains zero. For a positive equilibrium, divide the equilibrium equation by \(x\), giving \(a x/(1+x^2)=1\). Rearranging yields \(x^2-ax+1=0\). Its roots are real and distinct when \(a>2\). For this case, the zero state is stable, the smaller positive root is an unstable threshold, and the larger positive root is stable. Initial conditions below the threshold return to zero; those above it move toward the high-concentration state. Thus a sufficiently large transient trigger can, in this idealized model, leave the system in sustained activation.

Derivation: For positive \(x\), equilibrium requires \(a x/(1+x^2)=1\). The function \(x/(1+x^2)\) has derivative \((1-x^2)/(1+x^2)^2\), so its maximum for \(x>0\) is \(1/2\), at \(x=1\). Therefore no positive equilibrium exists when \(a<2\); at \(a=2\) the two positive roots meet at \(x=1\); and when \(a>2\), two positive equilibria are possible. This threshold condition follows from feedback strength relative to removal in this specific model.

Worked Example: Finding the Activation Threshold

Use a hypothetical Hill feedback loop with \(K=10\) pM, \(k=0.50\) per day, \(q_{\max}=12.5\) pM/day, \(n=2\), and zero basal production after the initiating input ends. The dimensionless feedback strength is $$ a=\frac{q_{\max}}{kK} =\frac{12.5\ \mathrm{pM/day}}{(0.50\ \mathrm{day}^{-1})(10\ \mathrm{pM})} =2.50. $$

The positive equilibria are $$ x_{\pm}=\frac{a\pm\sqrt{a^2-4}}{2} =\frac{2.50\pm\sqrt{2.25}}{2} =\frac{2.50\pm1.50}{2}, $$ so \(x_-=0.50\) and \(x_+=2.00\). Since \(C=Kx\), the threshold is \(5.00\) pM and the high equilibrium is \(20.0\) pM.

The direction of change verifies the threshold interpretation. At \(x=0.40\), $$ \frac{dx}{d\tau}=2.50\frac{0.40^2}{1+0.40^2}-0.40 =0.3448-0.40=-0.0552, $$ so the concentration decreases. At \(x=0.60\), $$ \frac{dx}{d\tau}=2.50\frac{0.60^2}{1+0.60^2}-0.60 =0.6618-0.60=0.0618, $$ so it increases. A transient input that leaves the system at \(6.00\) pM after the input ends places it above the \(5.00\) pM threshold; the model predicts movement toward \(20.0\) pM. A smaller residual concentration, such as \(4.00\) pM, is below threshold and declines. These are model predictions for the stated hypothetical parameters, not universal inflammatory concentrations.

Worked Example: How Clearance Can Remove a Sustained State

Keep the hypothetical \(q_{\max}=12.5\) pM/day and \(K=10\) pM, but increase the removal rate from \(0.50\) to \(0.70\) per day. Then $$ a=\frac{12.5}{(0.70)(10)}=1.786. $$

Because \(a<2\), the positive-equilibrium equation has no real roots. More directly, for every \(x>0\), \(x/(1+x^2)\leq 1/2\), so $$ a\frac{x}{1+x^2}\leq\frac{1.786}{2}=0.893<1. $$ Consequently \(a x^2/(1+x^2)-x<0\) at every positive \(x\): concentration declines after the external input is removed. In this model, the higher clearance shifts the loop from a regime with a threshold and high stable state to one without a positive sustained state.

Assumptions, Diagnostics, and Alternative Models

A feedback interpretation is only as reliable as its model structure and evidence. A time series in which cytokine concentration rises alongside inflammatory activity does not by itself show that the cytokine caused its own continued production. Both could be responding to a third process. Evidence for a feedback mechanism is stronger when the proposed direction is supported by mechanistic knowledge and by perturbations that alter a loop component in a way the model predicts.

Useful checks include whether the model reproduces the observed time scale and concentration range, whether estimated production and removal rates have plausible units and signs, and whether the data include measurements around initiation and decay. A steep Hill curve can make a fitted model appear to have a threshold even when the data do not distinguish it from a smoother response. If observations cover only a narrow concentration range, a linear approximation may be more defensible than claiming a threshold or bistability.

The linearized model is useful for local response size and stability near one operating point. A saturating Hill model can represent a bounded response and, when sufficiently cooperative, threshold-like behavior. More detailed multi-state models can represent distinct cytokine-producing and responding cell populations, transport between compartments, or explicit delays. Those added states require additional data and assumptions; they do not automatically make a model more informative.

  • Check the sign of the proposed loop. The modeled feedback contribution should increase with the cytokine state over the range where positive feedback is claimed.
  • Check the source–sink balance. Strong feedback does not imply persistence if removal dominates throughout the relevant concentration range.
  • Check equilibrium stability, not just equilibrium existence. A mathematical steady state may repel nearby trajectories rather than describe a sustained outcome.
  • Check whether the initiating input is truly absent. Continued external production can maintain cytokine concentration without self-sustaining feedback.
Common pitfall: “Cytokine concentration rose, so positive feedback must have made the response self-sustaining.” This conclusion does not follow from a rise alone. The increase could reflect a continuing external input, reduced removal, or another upstream driver. To claim self-sustaining activation, specify a feedback model, remove or account for the external input, and show that the model has a stable activated state reachable under the stated conditions.

These simple models have important limitations. A single concentration can conceal multiple cell types, compartments, and mediator interactions; constant clearance may be inaccurate; and unmodeled delays can alter the time course. Parameters inferred from concentration data alone may also be difficult to separate: increased production and decreased removal can both raise observed concentrations. The equations organize hypotheses and yield testable predictions, but they do not establish the biological loop without appropriate evidence.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. In the linear model, what is the local loop gain when \(g=0.30\) per day and \(k=0.60\) per day? Is the equilibrium locally stable?
  2. For the quadratic Hill model with \(a=1.6\), can a positive sustained equilibrium exist after the external input is removed? Justify your answer.
  3. For \(a=3\), calculate the two positive equilibria \(x_-\) and \(x_+\). Which is the threshold?
  4. Why does observing a high cytokine concentration after a stimulus not, by itself, demonstrate self-sustaining feedback?

Solutions. 1. \(L=g/k=0.30/0.60=0.50\), so the linear equilibrium is locally stable because \(L<1\). 2. No. For the specified quadratic model, positive equilibria require \(a\geq2\); here \(a=1.6<2\), so concentration declines at every positive state after the input ends. 3. \(x_{\pm}=(3\pm\sqrt{9-4})/2=(3\pm\sqrt{5})/2\), giving \(x_-\approx0.382\) and \(x_+\approx2.62\). The smaller root is the unstable threshold and the larger root is the stable activated state. 4. A continuing external source, reduced removal, or another driver could maintain the concentration. Demonstrating self-sustaining feedback requires evidence for the loop and a model or experiment showing persistence after the initiating input is removed.

Clinical Biostats