Tutorials › Quantitative Systems Pharmacology (QSP) › Progression from Insulin Resistance to Hyperglycemia

Metabolic Disease · Tutorial 350 of 799

Progression from Insulin Resistance to Hyperglycemia

Learn to construct and analyze a coupled progression model in which insulin resistance and beta-cell capacity shape glucose while glucose toxicity feeds back on both.

Advanced 12 min read

What You'll Learn

  • Define glucose regulation as a coupled model of glucose, insulin sensitivity, beta-cell capacity, and glucose toxicity.
  • Derive the glucose steady state and the critical sensitivity-capacity product for a chosen glucose threshold.
  • Distinguish positive toxicity production from a net increase in toxicity.
  • Calculate a self-consistent steady state when toxicity feeds back on sensitivity and capacity.
  • Identify model assumptions, diagnostics, and limitations when interpreting progression.

From Compensation to Progressive Hyperglycemia

Insulin resistance does not imply that glucose must immediately rise. If beta cells increase insulin secretion enough to compensate, glucose can remain near its regulated range even as insulin sensitivity falls. Hyperglycemia emerges when the combined ability of sensitivity and beta-cell capacity to support glucose control becomes inadequate. A further feedback can then arise: sustained high glucose may contribute to glucose toxicity, which in turn impairs sensitivity and beta-cell capacity.

The model here represents that progression with four state variables: glucose concentration \(G\), insulin sensitivity \(S\), beta-cell capacity \(b\), and a dimensionless glucose-toxicity burden \(T\). It is a reduced-order model, not a detailed account of every tissue flux or secretion pathway. Its purpose is to make the joint feedback explicit and to show how the model’s assumptions determine its predictions.

Definition: \(S\) is insulin sensitivity relative to a reference value, and \(b\) is beta-cell capacity relative to a reference value; both are positive and dimensionless. \(T\geq 0\) represents a slowly changing burden of glucose-associated impairment. \(G\) is glucose concentration in mmol/L, and time \(t\) is measured in days. In this model, the product \(Sb\) summarizes how much insulin-mediated compensation is available relative to the reference state.

A Reduced-Order Progression Model

Assume glucose adjusts more rapidly than sensitivity, capacity, and toxicity. Represent its adjustment toward a quasi-steady level \(G_{\mathrm{eq}}\) by a first-order relaxation equation. Let sensitivity and capacity recover toward reference values in the absence of toxicity, while toxicity increases their loss. Let glucose above a specified threshold \(G_{\mathrm{tox}}\) produce toxicity, and let toxicity also decay.

$$ \frac{dG}{dt}=\frac{G_{\mathrm{eq}}(S,b)-G}{\tau_G}, \qquad G_{\mathrm{eq}}(S,b)=G_{\mathrm{base}}+\frac{A}{Sb} $$
$$ \frac{dS}{dt}=\frac{S_0-S}{\tau_S}-k_{ST}TS, \qquad \frac{db}{dt}=\frac{b_0-b}{\tau_b}-k_{bT}Tb $$
$$ \frac{dT}{dt}=k_T\max(G-G_{\mathrm{tox}},0)-\frac{T}{\tau_T} $$

Here \(G_{\mathrm{base}}\) is a glucose baseline in mmol/L, \(A\) is a positive glucose increment in mmol/L, and \(\tau_G,\tau_S,\tau_b,\tau_T\) are positive time constants in days. \(S_0\) and \(b_0\) are positive reference values. The coefficients \(k_{ST}\) and \(k_{bT}\) have units of inverse days per unit of \(T\); \(k_T\) has units of \(T\) per day per mmol/L. The maximum function means toxicity production is zero at or below \(G_{\mathrm{tox}}\), and proportional to the glucose excess above it. The term \(-T/\tau_T\) describes decay at every glucose level.

The inverse dependence of \(G_{\mathrm{eq}}\) on \(Sb\) is a compact phenomenological assumption: less sensitivity or less capacity raises the glucose level required by the model’s balance. It is not a universal physiological law or a substitute for fitting a more detailed flux model. The reference values, parameters, and thresholds must be chosen for the purpose and data at hand. In examples below, \(7.0\) mmol/L is an illustrative hyperglycemia cutoff, not a claim that this model establishes a clinical diagnosis.

What the Equations Imply

If \(S\) and \(b\) are approximately constant, the glucose equation approaches \(G_{\mathrm{eq}}\). Solving the equilibrium relation for the product \(Sb\) gives a direct measure of how much combined compensation is required to keep glucose below a chosen level. For an illustrative cutoff \(G_H>G_{\mathrm{base}}\), the model predicts \(G_{\mathrm{eq}}\geq G_H\) when \(Sb\leq A/(G_H-G_{\mathrm{base}})\).

$$ G_{\mathrm{eq}}\geq G_H \quad\Longleftrightarrow\quad G_{\mathrm{base}}+\frac{A}{Sb}\geq G_H \quad\Longleftrightarrow\quad Sb\leq\frac{A}{G_H-G_{\mathrm{base}}} $$

This result shows why considering sensitivity or beta-cell capacity alone can be misleading. A decline in \(S\) can be offset by a rise in \(b\), but only while the product remains sufficient. The model’s toxicity feedback changes that product over time, so a static threshold calculation is a starting point rather than a full progression analysis.

Worked Example: Finding the Compensation Threshold

Suppose \(G_{\mathrm{base}}=4.8\) mmol/L, \(A=1.2\) mmol/L, and \(G_H=7.0\) mmol/L. Treat \(S\) and \(b\) as fixed over the period of interest. Calculate the product threshold and glucose when \(S=0.65\) and \(b=0.70\).

The critical product is \(A/(G_H-G_{\mathrm{base}})=1.2/(7.0-4.8)=1.2/2.2=0.5455\). The specified product is \(Sb=0.65(0.70)=0.4550\), below that threshold.

$$ G_{\mathrm{eq}}=4.8+\frac{1.2}{0.4550} =4.8+2.637=7.437\ \mathrm{mmol/L} $$

Under these assumptions, the equilibrium glucose exceeds the illustrative cutoff. This calculation does not show that either parameter alone caused the result, and it does not establish that the glucose state is stable over longer periods: it holds \(S\) and \(b\) fixed and omits evolving toxicity.

For a fixed glucose concentration above \(G_{\mathrm{tox}}\), toxicity production is positive, but the net rate of change depends on the current toxicity burden. Setting \(dT/dt=0\) gives the toxicity equilibrium at that fixed glucose. Below this equilibrium, \(T\) increases; above it, \(T\) decreases. At or below \(G_{\mathrm{tox}}\), production is zero and positive \(T\) decays.

$$ T^*=\tau_T k_T\max(G-G_{\mathrm{tox}},0) $$

For fixed \(G\), the equation is linear in \(T\), so its solution is \(T(t)=T^*+[T(0)-T^*]e^{-t/\tau_T}\). This justifies the direction-of-change interpretation: the difference \(T(t)-T^*\) retains its initial sign while shrinking exponentially.

Worked Example: Production Is Positive but Toxicity Falls

Consider the illustrative fixed-glucose case \(G=8\) mmol/L, \(G_{\mathrm{tox}}=6\) mmol/L, \(k_T=0.05\) \(T\) per day per mmol/L, \(\tau_T=10\) days, and current \(T=2\). The production rate is \(0.05(8-6)=0.10\) \(T\) per day. The decay rate is \(2/10=0.20\) \(T\) per day.

$$ \frac{dT}{dt}=0.10-0.20=-0.10\ T\text{ per day}, \qquad T^*=10(0.05)(8-6)=1 $$

Thus toxicity production is positive, yet the net burden decreases because the current \(T=2\) exceeds the fixed-glucose equilibrium \(T^*=1\). If glucose stayed fixed for 10 days, the model gives \(T(10)=1+(2-1)e^{-10/10}=1.368\). The burden moves toward 1 rather than increasing simply because glucose remains above threshold.

Feedback Equilibrium and Progression

At a steady state with fixed parameters, the sensitivity and capacity equations can also be solved explicitly as functions of \(T\). Setting each derivative to zero yields \(S=S_0/(1+\tau_S k_{ST}T)\) and \(b=b_0/(1+\tau_b k_{bT}T)\). These expressions make the feedback clear: at a positive toxicity burden, both states are below their reference values. Together with the glucose equation and \(T^*\), they define a coupled equilibrium, not three independent calculations.

$$ S^*=\frac{S_0}{1+\tau_S k_{ST}T^*}, \qquad b^*=\frac{b_0}{1+\tau_b k_{bT}T^*}, \qquad G^*=G_{\mathrm{base}}+\frac{A}{S^*b^*} $$

The derivation follows by rearranging, for example, \(0=(S_0-S)/\tau_S-k_{ST}TS\): multiplying by \(\tau_S\) gives \(S_0=S(1+\tau_S k_{ST}T)\), and division gives the stated expression. A positive-toxicity equilibrium must also satisfy \(T^*=\tau_T k_T(G^*-G_{\mathrm{tox}})\) and \(G^*>G_{\mathrm{tox}}\). If the resulting candidate violates that inequality, it is not an equilibrium on the positive-production branch; the zero-production branch must instead be checked.

Worked Example: Solving a Coupled Positive-Toxicity Equilibrium

Let \(G_{\mathrm{base}}=4.8\) mmol/L, \(A=2.3\) mmol/L, \(S_0=b_0=1\), \(G_{\mathrm{tox}}=6.5\) mmol/L, \(\tau_S=\tau_b=20\) days, \(k_{ST}=0.002\) and \(k_{bT}=0.003\) per day per unit \(T\), \(\tau_T=10\) days, and \(k_T=0.02\) \(T\) per day per mmol/L. Seek an equilibrium on the branch \(G^*>6.5\).

The sensitivity and capacity equilibria are \(S^*=1/(1+0.04T^*)\) and \(b^*=1/(1+0.06T^*)\). Since \(\tau_Tk_T=0.20\) \(T\) per mmol/L, the toxicity balance implies \(T^*=0.20(G^*-6.5)\), or \(G^*=6.5+5T^*\). Substitution into the glucose equilibrium gives

$$ G^*=4.8+2.3(1+0.04T^*)(1+0.06T^*) =7.1+0.230T^*+0.00552(T^*)^2 $$

Equating this expression to \(6.5+5T^*\) gives \(0.00552(T^*)^2-4.770T^*+0.600=0\). The smaller root is approximately \(T^*=0.126\); the larger root is approximately \(T^*=864.0\). Using the smaller root gives \(G^*=6.5+5(0.126)=7.130\) mmol/L, which satisfies the assumed positive-production condition. It also gives \(S^*=1/(1+0.04(0.126))=0.9950\) and \(b^*=1/(1+0.06(0.126))=0.9925\), to four significant figures.

This algebra identifies a self-consistent equilibrium candidate, not a proof that it is dynamically stable or that the large-root candidate is physiologically meaningful. Stability must be assessed from the coupled dynamics, and plausible parameter ranges must be justified independently.

Assumptions, Diagnostics, and Alternatives

The model’s interpretation depends on several assumptions: glucose relaxes toward the stated equilibrium on the timescale \(\tau_G\); \(S\), \(b\), and \(T\) are adequate summaries of slower progression; glucose toxicity is driven by excess above a fixed threshold; and the loss terms are proportional to both \(T\) and the state being impaired. The model also assumes the parameters remain constant during the modeled period. These assumptions are simplifications, not consequences of the equations.

When fitting or applying this structure, check whether it can reproduce both transient and steady-state behavior. Useful warning signs include systematic residual patterns in glucose over time, implausible fitted values of \(S\), \(b\), or \(T\), and parameter estimates that change greatly when the observation window or initial conditions change. Because glucose depends on the product \(Sb\), glucose observations alone may not distinguish low sensitivity from low capacity; independent measurements or informative experimental challenges may be needed to constrain them. A threshold-driven toxicity term can also make predictions sensitive to small changes near \(G_{\mathrm{tox}}\). Examine that sensitivity rather than treating the threshold as precisely known.

A static calculation using \(Sb\) can be useful for locating a compensation boundary, but it cannot describe delayed feedback or recovery. A more mechanistic flux model can represent secretion, insulin action, and glucose production separately, but requires additional states, parameters, and data to identify them. This reduced model is useful for organizing a progression hypothesis and testing its qualitative implications; it cannot, by itself, identify a biological cause, establish an individual diagnosis, or replace validation against relevant data.

Common pitfall: “Glucose is above \(G_{\mathrm{tox}}\), so \(T\) must be increasing.” That confuses positive production with positive net change. The correct comparison is between current \(T\) and \(T^*=\tau_Tk_T(G-G_{\mathrm{tox}})\) at fixed glucose: \(T\) increases below that value and decreases above it.

Check Your Understanding

Use the model definitions and state any assumptions needed for your calculation.

  1. If \(G_{\mathrm{base}}=5.0\) mmol/L, \(A=1.0\) mmol/L, \(S=0.50\), and \(b=0.80\), what is \(G_{\mathrm{eq}}\)? Does it exceed an illustrative cutoff of \(7.0\) mmol/L?
  2. At fixed \(G=9\) mmol/L, \(G_{\mathrm{tox}}=7\) mmol/L, \(k_T=0.03\) \(T\) per day per mmol/L, and \(\tau_T=5\) days, calculate \(T^*\). Is \(T\) increasing or decreasing when its current value is \(0.4\)?
  3. Explain why an observed rise in glucose does not, by itself, identify whether sensitivity or beta-cell capacity declined.
  4. Which model assumption is challenged if glucose remains far from \(G_{\mathrm{eq}}\) over a period much longer than \(\tau_G\)?

Solutions. 1. \(Sb=0.50(0.80)=0.40\), so \(G_{\mathrm{eq}}=5.0+1.0/0.40=7.50\) mmol/L, above the illustrative cutoff. 2. \(T^*=5(0.03)(9-7)=0.300\). Because current \(T=0.4>T^*\), decay exceeds production and \(T\) decreases. 3. The glucose equilibrium depends on the product \(Sb\); different combinations of \(S\) and \(b\) can give the same product and thus the same predicted glucose. 4. That observation challenges the assumed rapid first-order relaxation toward \(G_{\mathrm{eq}}\), or suggests that the model’s target or timescale is misspecified.

Clinical Biostats