Embedding Saturable Metabolism in a PBPK Tissue

From an Enzyme Rate to a Tissue Mass Balance

A Michaelis–Menten rate law describes how quickly an enzyme consumes substrate at a specified concentration. A PBPK tissue model must do more: it must account for substrate delivered by blood, substrate leaving in venous blood, and substrate metabolized within the tissue. These terms are rates of parent-drug amount, so they must be combined in a single mass balance.

The central distinction is between tissue concentration, which drives the metabolic rate in this model, and blood concentration, which determines delivery and removal by perfusion. They are not interchangeable. A tissue-to-blood partition coefficient links them under the perfusion-limited assumption used here.

With these definitions, the tissue amount balance is

Substituting \(A_T=V_TC_T\) and \(C_v=C_T/K_{T:b}\), and treating \(V_T\) as constant, gives the concentration equation

The first term is net perfusion: \(QC_a\) enters and \(QC_v\) leaves. The second is metabolism within the tissue. The earlier material on the Michaelis–Menten elimination rate law establishes the saturation form; here its role is specifically as a local sink in a tissue balance.

Assumptions and Unit Checks

This equation is appropriate when the tissue is represented as a single, well-mixed compartment; blood flow is the limiting exchange process; and tissue and venous blood are in instantaneous partition equilibrium. The metabolic rate is assumed to depend on \(C_T\), with fixed \(V_{\max,T}\) and \(K_{m,T}\) over the interval modeled. The equation does not assert that every biological tissue satisfies these assumptions.

Concentrations used in the partition relation and in the rate law must be on compatible bases. Here \(C_a\) and \(C_v\) are blood concentrations, while \(C_T\) and \(K_{m,T}\) are tissue concentrations. If an enzyme rate is parameterized by an unbound intracellular concentration instead, the model needs an appropriate relationship between that concentration and \(C_T\); it is not valid simply to relabel one concentration as the other.

Check units before calculating. If flow is in L/min and concentration in mg/L, each perfusion term is in mg/min. The metabolic term is also in mg/min when \(V_{\max,T}\) is in mg/min and \(C_T/(K_{m,T}+C_T)\) is dimensionless. Therefore the right side is amount per time, matching \(V_T\,dC_T/dt\). Dividing by \(V_T\), in L, gives concentration change per time, such as mg/L/min.

Sometimes enzyme capacity is provided per tissue volume rather than for the whole tissue. If \(v_{\max,T}\) is in mg/L/min, the total tissue capacity is \(V_{\max,T}=v_{\max,T}V_T\). Use either total capacity or the appropriately volume-scaled rate, not both; multiplying a total \(V_{\max,T}\) by tissue volume again overcounts metabolism.

Why the Concentrations Must Stay Distinct

The partition relation connects the concentration that determines venous outflow to the tissue state. It does not say that tissue concentration equals arterial concentration. The difference matters both for transport and for metabolism. Using \(C_a\) as the metabolic concentration would make the metabolic sink respond to incoming blood rather than to the modeled tissue concentration.

Mass conservation is explicit in this calculation. The amount entering is \(QC_a\); the amount leaving in blood is \(QC_v\); and metabolism removes \(v_{\mathrm{met}}\) of parent drug per unit time. No exchange term is an additional elimination process: it transfers parent drug between blood and tissue. Across a whole-body model, internal transfers cancel when compartment balances are added, while metabolic parent-drug loss remains a sink.

Steady State and a Short Derivation

A steady tissue concentration satisfies \(dC_T/dt=0\). Therefore, net blood delivery must exactly balance the metabolic rate:

The left side decreases as \(C_T\) increases, because the venous concentration rises. The right side increases with \(C_T\), approaching \(V_{\max,T}\). This gives a useful existence and uniqueness check for positive parameters and \(C_a>0\). Define the difference between net delivery and metabolism as \(g(C_T)\). At zero concentration, \(g(0)=QC_a>0\). As \(C_T\) becomes large, the perfusion term decreases without bound while metabolism approaches a finite maximum, so \(g(C_T)\) becomes negative. Also,

Consequently, \(g\) crosses zero once: there is a unique positive steady concentration under these assumptions. The negative slope also means a small increase in \(C_T\) above that steady value makes the balance favor depletion, while a small decrease favors accumulation.

Diagnostics, Alternatives, and Limitations

A useful first diagnostic is a term-by-term unit check: both perfusion and metabolism must be amount per time before they are subtracted. Next, check the sign of the balance. If \(QC_a\) is less than \(QC_v+v_{\mathrm{met}}\), the tissue concentration should decrease; if it is greater, the tissue concentration should increase. A computed negative concentration or a steady state that fails to balance the two sides signals an algebra, coding, or unit error.

Also inspect the concentration mapping. If \(K_{T:b}=1\), the model implies equal tissue and venous blood concentrations; if it differs from one, they differ by that factor. Confusing the partition coefficient’s direction is a common source of error: with the definition here, \(C_v=C_T/K_{T:b}\), not \(K_{T:b}C_T\). Finally, verify that the \(K_{m,T}\) concentration basis matches the concentration used in the rate law.

This perfusion-limited formulation is one modeling choice, not a universal tissue description. A permeability-limited formulation represents exchange across a barrier explicitly and may be preferable when perfusion alone cannot describe equilibration. A linear metabolic sink can be a useful approximation at concentrations well below \(K_{m,T}\), since \(C_T/(K_{m,T}+C_T)\) is then approximately \(C_T/K_{m,T}\); it does not preserve the capacity-limited behavior at higher concentrations. These alternatives change assumptions or model structure, so they should not be substituted without checking whether their conditions suit the tissue and data.

Common mistake: “The tissue is supplied by blood, so use arterial blood concentration in the Michaelis–Menten term.” That reasoning confuses delivery with the local concentration driving metabolism. In this model the rate law uses \(C_T\); arterial and venous blood concentrations belong in the perfusion term, with venous concentration linked to tissue concentration through \(K_{T:b}\). A different model can use a different enzyme-site concentration, but it must define and justify that mapping explicitly.