Tutorials › Quantitative Systems Pharmacology (QSP) › Cholesterol Synthesis, LDL Receptor Clearance, and Bile-Acid Cycling

Metabolic Disease · Tutorial 347 of 799

Cholesterol Synthesis, LDL Receptor Clearance, and Bile-Acid Cycling

Learn to write and check a cholesterol balance model that separates LDL uptake, cholesterol disposal, bile-acid cycling, and net fecal loss.

Advanced 13 min read

What You'll Learn

  • Write hepatic cholesterol and plasma LDL balance equations with consistent flux units.
  • Represent LDL-receptor-mediated clearance separately from non-receptor removal.
  • Distinguish biliary cholesterol secretion from its net loss after intestinal return.
  • Model bile-acid secretion, reabsorption, and fecal loss as distinct fluxes.
  • Derive how enterohepatic recycling changes the relationship between pool size and net bile-acid loss.
  • Diagnose common balance errors and compare reduced and compartmental model formulations.

Why a Cholesterol Balance Model Needs Several Compartments

A measured LDL cholesterol concentration is the outcome of several flows, not a direct readout of one process. Hepatic synthesis adds cholesterol to the liver pool; LDL-receptor-mediated uptake transfers cholesterol from plasma LDL to the liver; and biliary routes remove cholesterol from the body, with some secreted material returning from the intestine. A model that combines these processes into one “clearance” term can conceal which flux changed.

The model below tracks amounts rather than concentrations. This makes the balance explicit: each state changes by inflows minus outflows. It is a reduced model, not a complete account of lipoprotein transport or sterol chemistry. In particular, the LDL input term aggregates the production and remodeling processes described in the earlier material on lipoprotein remodeling.

Definition: Let \(H\) be hepatic cholesterol amount and \(L\) be plasma LDL cholesterol amount, both measured in mmol of cholesterol. Let \(B\) be the body bile-acid pool and \(G\) the intestinal bile-acid amount, measured in mmol of bile-acid molecules. Since each bile-acid molecule is derived from one cholesterol molecule, cholesterol-to-bile-acid conversion is represented in molar equivalents. Fluxes are in mmol/day, and rate constants are in day\(^{-1}\).

The model distinguishes gross biliary cholesterol secretion from net cholesterol loss. Let \(J_C\) be the gross flux of cholesterol into bile, and let \(f_C\) be the fraction of that flux that returns to the hepatic cholesterol pool after intestinal absorption. Thus, \(0\leq f_C\leq1\), and the net loss through this route is \((1-f_C)J_C\). This reduced representation treats return as an immediate fraction; a gut compartment could be added to represent transit delay.

Balance Equations and Assumptions

Let \(S(H)\) be hepatic cholesterol synthesis, \(A\) the absorbed dietary cholesterol input, \(J_B\) the flux of hepatic cholesterol converted to bile acids, \(J_V\) hepatic cholesterol export in very-low-density lipoprotein and other apoB-containing particles, and \(J_O\) other hepatic cholesterol use or loss. Let \(J_{\mathrm{in}}\) denote LDL cholesterol appearance in plasma, including production and remodeling into LDL. The LDL-receptor uptake rate constant is \(k_R R\), where \(R\) is a dimensionless relative receptor-activity state; \(k_N\) is the non-receptor LDL removal rate constant.

$$ \frac{dH}{dt} = S(H)+k_R R L+A-J_B-(1-f_C)J_C-J_V-J_O $$
$$ \frac{dL}{dt} = J_{\mathrm{in}}-(k_R R+k_N)L $$

The hepatic uptake term \(k_RRL\) appears as an outflow from plasma LDL and an equal inflow to the hepatic cholesterol pool. This transfer is internal to the combined liver-plus-plasma cholesterol system: it changes the location of cholesterol, not its total amount. By contrast, synthesis and absorbed dietary cholesterol are external inputs to the modeled system, while net biliary loss, export, and other uses are outflows.

For enterohepatic bile-acid cycling, let \(k_s\) be the rate constant for secretion of bile acids from the body pool into the gut, \(k_a\) the intestinal reabsorption rate constant, and \(k_f\) the fecal loss rate constant. A two-compartment representation is:

$$ \frac{dB}{dt}=J_B+k_aG-k_sB, \qquad \frac{dG}{dt}=k_sB-(k_a+k_f)G $$

In this formulation, secretion and reabsorption are internal transfers. Only fecal loss, \(k_fG\), removes bile acids from the body pool. The equations assume well-mixed pools, first-order transport, constant rate parameters over the modeled period, and no explicit transit delay. They also assume that \(J_B\) is expressed in molar equivalents consistent with the bile-acid pools.

Physiological feedback can be represented through the fluxes and receptor state. For example, synthesis may decrease as hepatic cholesterol rises, while receptor activity may increase when hepatic cholesterol is low. One possible model structure is:

$$ S(H)=\frac{S_{\max}}{1+(H/K_S)^n}, \qquad \frac{dR}{dt} = \frac{k_{\mathrm{syn},R}}{1+(H/K_R)^m} -k_{\mathrm{deg},R}R $$

Here \(K_S\) and \(K_R\) are positive pool scales, \(n\) and \(m\) are positive exponents, and the receptor synthesis and degradation parameters are positive. These functions encode a qualitative negative relationship between hepatic cholesterol and synthesis or receptor expression; they are candidate model structures, not universal parameterizations. Data and the scientific question determine whether feedback needs to be modeled dynamically or can be approximated as a fixed rate over a short interval.

Steady-State Consequences of Recycling

At a steady state, \(dG/dt=0\), so intestinal bile-acid amount satisfies \(G^*=k_sB^*/(k_a+k_f)\). Substituting this into the bile-acid pool balance gives:

$$ 0=J_B+k_a\frac{k_sB^*}{k_a+k_f}-k_sB^* =J_B-\frac{k_sk_f}{k_a+k_f}B^* $$

The last equality follows because \(k_sB^* - k_ak_sB^*/(k_a+k_f)=k_sk_fB^*/(k_a+k_f)\). Therefore, the effective net loss rate from the body bile-acid pool is \(k_sk_f/(k_a+k_f)\), not the secretion rate constant \(k_s\). The fraction of bile acids reabsorbed on a gut passage is \(r=k_a/(k_a+k_f)\); the fraction lost is \(1-r=k_f/(k_a+k_f)\). Consequently, steady-state synthesis balances fecal loss: \(J_B=(1-r)k_sB^*\). Repeated recycling can make the total secretion flux large while net loss remains comparatively small.

For LDL, if \(J_{\mathrm{in}}\), \(R\), \(k_R\), and \(k_N\) are constant and \(k_RR+k_N>0\), the steady-state amount is \(L^*=J_{\mathrm{in}}/(k_RR+k_N)\). This relation does not imply that LDL concentration is determined by receptor activity alone: appearance and non-receptor removal also matter. If concentration rather than amount is modeled, distribution volume must be specified to convert between them.

Worked Examples

Worked Example: Checking a Hepatic Cholesterol Balance

Suppose a hypothetical steady-state model uses mmol/day and has synthesis \(S=1.8\), LDL-receptor uptake \(k_RRL=1.2\), dietary absorption \(A=0.3\), bile-acid synthesis \(J_B=0.6\), gross biliary cholesterol secretion \(J_C=0.6\), return fraction \(f_C=0.5\), apoB-particle export \(J_V=1.5\), and other use \(J_O=0.9\). Check whether the hepatic pool is at steady state.

First calculate net biliary cholesterol loss: \( (1-f_C)J_C=(1-0.5)(0.6)=0.30\) mmol/day. Total hepatic inputs are \(1.8+1.2+0.3=3.3\) mmol/day. Total outflows are \(0.6+0.30+1.5+0.9=3.3\) mmol/day. Thus:

$$ \frac{dH}{dt}=3.3-3.3=0\ \text{mmol/day} $$

The specified fluxes are consistent with a steady hepatic cholesterol amount. This calculation checks mass balance only; it does not establish that the assumed fluxes are physiologically identifiable or uniquely correct.

Worked Example: Separating LDL-Receptor and Non-Receptor Clearance

In a hypothetical plasma LDL model, LDL appearance is \(J_{\mathrm{in}}=2.4\) mmol/day, \(k_RR=0.080\) day\(^{-1}\), and \(k_N=0.040\) day\(^{-1}\). Calculate steady-state LDL amount and the receptor-mediated share of total removal.

The total removal rate constant is \(0.080+0.040=0.120\) day\(^{-1}\). Hence:

$$ L^*=\frac{2.4\ \text{mmol/day}}{0.120\ \text{day}^{-1}} =20\ \text{mmol} $$

Receptor-mediated removal is \(0.080\times20=1.6\) mmol/day; non-receptor removal is \(0.040\times20=0.80\) mmol/day. Their sum, \(2.4\) mmol/day, matches appearance. The receptor-mediated share is \(1.6/2.4=0.667\), or about 66.7%.

If receptor-mediated clearance falls by 25% while the other quantities remain fixed, its rate constant becomes \(0.060\) day\(^{-1}\). The new total rate constant is \(0.060+0.040=0.100\) day\(^{-1}\), giving \(L^*=2.4/0.100=24\) mmol. The model predicts a 20% higher steady-state LDL amount, not a 25% increase, because the non-receptor route remains active.

Worked Example: Calculating Bile-Acid Recycling and Fecal Loss

Consider a hypothetical bile-acid pool of \(B=12\) mmol with \(k_s=1.0\) day\(^{-1}\), \(k_a=4.0\) day\(^{-1}\), and \(k_f=1.0\) day\(^{-1}\). Calculate the gut amount at steady state, reabsorbed flux, fecal loss, and the synthesis required to maintain the body pool.

The gut balance at steady state gives \(G=k_sB/(k_a+k_f)\), so:

$$ G=\frac{(1.0\ \text{day}^{-1})(12\ \text{mmol})} {4.0\ \text{day}^{-1}+1.0\ \text{day}^{-1}} =2.4\ \text{mmol} $$

The secreted flux is \(k_sB=12\) mmol/day. Reabsorption is \(k_aG=(4.0)(2.4)=9.6\) mmol/day, and fecal loss is \(k_fG=(1.0)(2.4)=2.4\) mmol/day. Their difference, \(12-9.6=2.4\) mmol/day, is the net loss. Therefore \(J_B=2.4\) mmol/day is required at steady state. The reabsorbed fraction is \(4/(4+1)=0.80\), so 80% of each modeled gut passage returns and 20% is lost.

Diagnostics, Alternatives, and Limitations

A balance model should be checked before its parameters are interpreted. At a proposed steady state, substitute the states and parameters into every differential equation and verify that each net rate is zero within rounding. Also check units: a rate constant in day\(^{-1}\) multiplied by an amount in mmol gives a flux in mmol/day. An amount cannot be added directly to a flux.

  • Check internal transfers. LDL-receptor uptake must be a loss from \(L\) and an equal gain to \(H\). Bile-acid secretion and reabsorption likewise cancel when the \(B\) and \(G\) equations are summed.
  • Check net versus gross loss. Gross biliary cholesterol secretion is not necessarily net body loss. If a return fraction is used, the hepatic equation must include \((1-f_C)J_C\), not \(J_C\) alone.
  • Check the model’s time behavior. With constant LDL input and positive total clearance, \(L\) approaches its steady state. If simulated LDL grows without bound under those conditions, inspect signs, units, and whether clearance has been omitted.
  • Check physical constraints. Amounts and fluxes should not become negative under the modeled conditions. A fitted feedback function that predicts negative synthesis or receptor activity signals an unsuitable functional form or parameter range.

The fractional-return approximation for biliary cholesterol is compact when the main question concerns net balance. An explicit intestinal compartment is preferable when timing, absorption kinetics, or transient perturbations matter. Likewise, the two-compartment bile-acid model distinguishes secretion, return, and loss; a single effective loss term can be adequate at steady state but cannot show the individual cycle fluxes. Adding compartments creates more parameters, so it is useful only when the available data or question can support that detail.

The model is also limited by its aggregation. It does not distinguish cholesterol among intracellular pools, lipoprotein subclasses, or individual bile-acid species. A steady-state fit can conceal compensating changes in synthesis and clearance: the same pool may be consistent with different combinations of input and output fluxes. A pool measurement alone therefore does not identify every flux. Independent flux or perturbation information is needed to discriminate among mechanisms.

Common Mistake

A frequent error is to treat every molecule secreted into bile as permanently eliminated. For example, one might subtract the entire \(J_C\) from hepatic cholesterol and also count intestinal return as a separate input without defining which fraction returns. The correction is to specify the bookkeeping once: either use net loss \((1-f_C)J_C\) in the reduced hepatic balance, or write explicit gut transport equations and include the returning flux. Do not do both, and do not count an internal transfer as a net loss from the whole system.

Key takeaway: Cholesterol balance depends on both pool sizes and fluxes. LDL-receptor uptake redistributes cholesterol from plasma to liver; biliary secretion is not synonymous with net loss; and bile-acid recycling reduces fecal loss relative to the gross secretion flux.

Check Your Understanding

Try each problem, then compare with the solutions below.

  1. A model has \(S=1.0\), LDL uptake \(=0.7\), and dietary input \(A=0.2\) mmol/day. Its outflows are \(J_B=0.3\), net biliary loss \(=0.1\), \(J_V=0.8\), and \(J_O=0.7\) mmol/day. Is the hepatic pool at steady state?
  2. If \(J_{\mathrm{in}}=1.5\) mmol/day, \(k_RR=0.05\) day\(^{-1}\), and \(k_N=0.10\) day\(^{-1}\), what is \(L^*\), and what fraction of LDL removal is receptor-mediated?
  3. For \(k_a=3\) day\(^{-1}\) and \(k_f=1\) day\(^{-1}\), what fraction of bile acids entering the gut is reabsorbed? If gross secretion is 8 mmol/day, what are reabsorption and fecal-loss fluxes?
  4. Why is an explicit gut compartment more informative than a fractional-return term when modeling delayed changes in biliary secretion?

Solutions. 1. Inputs total \(1.0+0.7+0.2=1.9\) mmol/day; outflows total \(0.3+0.1+0.8+0.7=1.9\) mmol/day, so \(dH/dt=0\). 2. \(L^*=1.5/(0.05+0.10)=10\) mmol. Receptor-mediated removal is \(0.05(10)=0.5\) mmol/day, so its share is \(0.5/1.5=0.333\), or 33.3%. 3. The reabsorbed fraction is \(3/(3+1)=0.75\). Reabsorption is \(0.75(8)=6\) mmol/day and fecal loss is \(0.25(8)=2\) mmol/day. 4. An explicit gut compartment represents material in transit and its time-dependent return or loss; a fixed fraction collapses those events into an immediate net effect.

Clinical Biostats