Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Modeling Combined Renal and Hepatic Impairment

Renal and Hepatic Impairment · Tutorial 631 of 799

Modeling Combined Renal and Hepatic Impairment

Learn to combine renal and hepatic disease effects mechanistically, keep blood- and plasma-based clearances consistent, and diagnose double counting in a joint PBPK model.

Advanced 11 min read

What You'll Learn

  • Represent renal and hepatic impairment with separate, physiologically meaningful covariates.
  • Convert hepatic blood clearance to plasma clearance using the blood-to-plasma concentration ratio.
  • Incorporate altered protein binding once in each affected clearance pathway.
  • Calculate joint clearance and exposure under a stated set of linear model assumptions.
  • Recognize when multiplying impairment ratios double counts a shared disease effect.
  • Check joint predictions against pathway contributions, units, and model assumptions.

Why Combined Impairment Requires a Joint Model

A patient may have both renal and hepatic disease, but the effects on drug disposition are not necessarily independent. Renal dysfunction can alter filtration and secretion; hepatic dysfunction can change blood flow, functional liver mass, and enzyme capacity. Either condition may also change binding, and binding can affect more than one clearance pathway. A combined prediction therefore needs to represent the relevant physiological changes together, rather than simply applying two overall impairment multipliers.

The key accounting rule is to add pathway clearances only after expressing them on the same concentration basis. A hepatic clearance calculated relative to blood concentration and a renal clearance calculated relative to plasma concentration cannot be added directly in general. Binding must also enter through the appropriate pathway equations, not as an additional multiplier applied after those equations have already included it.

Definition: A joint renal-hepatic impairment model predicts drug disposition for a patient with both renal and hepatic disease by specifying the disease-dependent physiological or drug-specific parameters for each relevant elimination pathway. For a linear model, total clearance is the sum of pathway clearances expressed on one common concentration basis.

Specify the Joint State and the Clearance Basis

Represent the patient's disease state by \(z=(r,h)\), where \(r\) denotes the renal impairment state and \(h\) the hepatic impairment state. These labels may index measured or estimated quantities such as glomerular filtration rate, renal secretion capacity, hepatic blood flow, functional liver mass, enzyme capacity, plasma unbound fraction, or blood-to-plasma ratio. They are not themselves universal severity multipliers: each model must define how the state changes its parameters.

For clarity, use plasma concentration as the common basis for total clearance. Let \(C_p\) and \(C_b\) be plasma and blood concentrations, respectively, and define \(R_b(z)=C_b/C_p\). Let \(f_{u,p}(z)\) and \(f_{u,b}(z)\) be the unbound fractions in plasma and blood. They are distinct quantities; do not assume they are equal unless justified for the drug and model.

In the well-stirred liver model, hepatic clearance is first calculated on a blood basis. Here \(Q_H(z)\) is hepatic blood flow and \(CL_{\mathrm{int},H}(z)\) is hepatic intrinsic clearance, with units consistent with blood concentration.

$$ CL_{H,b}(z)= \frac{Q_H(z)f_{u,b}(z)CL_{\mathrm{int},H}(z)} {Q_H(z)+f_{u,b}(z)CL_{\mathrm{int},H}(z)} $$

Because the hepatic elimination rate is \(CL_{H,b}C_b\), and \(C_b=R_bC_p\), the equivalent hepatic clearance on a plasma basis is \(CL_{H,p}=R_bCL_{H,b}\). For a simplified renal model with filtration and active secretion, the plasma-based renal clearance is \(CL_{R,p}=f_{u,p}[GFR+CL_{\mathrm{int},R\mathrm{sec}}]\). The secretion term is defined here as an intrinsic clearance referenced to unbound plasma concentration. More detailed models may represent transport, reabsorption, or saturation explicitly instead of using this simplification.

$$ CL_{\mathrm{total},p}(z)= R_b(z)CL_{H,b}(z) + f_{u,p}(z)\left[GFR(z)+CL_{\mathrm{int},R\mathrm{sec}}(z)\right] + CL_{\mathrm{other},p}(z) $$

The final term represents any additional elimination pathway, also expressed on a plasma basis. If no other pathway is modeled, set it to zero. A total clearance defined on a blood basis could be used instead, but then every pathway—including the renal pathway—must be converted to that basis. The choice of basis does not change the predicted elimination rate when conversions are done consistently.

Why Pathway Clearances Add

For linear elimination at a given state, the rate of drug removal through each pathway is its clearance times the concentration on which that clearance is defined. The total removal rate is the sum of the pathway rates. This provides both the justification for adding pathway clearances and a direct check on the concentration basis.

1
Write each pathway rate.
Hepatic removal is \(CL_{H,b}C_b\); renal removal in this simplified model is \(CL_{R,p}C_p\).
2
Express both rates using plasma concentration.
Substitute \(C_b=R_bC_p\), giving hepatic removal \(R_bCL_{H,b}C_p\), while renal removal is already plasma based.
3
Factor out the common concentration.
The sum of rates is \((R_bCL_{H,b}+CL_{R,p})C_p\). Dividing by \(C_p\) gives the total plasma-based clearance.

This derivation assumes the pathway rates are additive and linear in concentration over the range considered. It also assumes the pathway equations describe distinct removal routes. If a fitted “nonrenal clearance” already includes a pathway represented separately, adding both terms counts that route twice.

Worked Example: Calculating Joint Clearance on a Plasma Basis

Suppose a hypothetical patient with both impairments has \(Q_H=90\) L/h, \(f_{u,b}=0.060\), \(CL_{\mathrm{int},H}=300\) L/h, \(R_b=1.50\), \(f_{u,p}=0.100\), \(GFR=4.0\) L/h, and \(CL_{\mathrm{int},R\mathrm{sec}}=6.0\) L/h. Assume there are no other modeled pathways. Calculate total clearance.

First calculate hepatic blood clearance:

$$ CL_{H,b}= \frac{90(0.060)(300)}{90+(0.060)(300)} = \frac{1620}{108} = 15.0\ \text{L/h} $$

Convert this to a plasma basis: \(CL_{H,p}=1.50(15.0)=22.5\) L/h. The renal plasma-based clearance is \(CL_{R,p}=0.100(4.0+6.0)=1.00\) L/h. Therefore,

$$ CL_{\mathrm{total},p}=22.5+1.00=23.5\ \text{L/h} $$

Adding \(15.0\) L/h directly to \(1.00\) L/h would mix blood- and plasma-based clearances. That calculation would be valid only if \(R_b=1\), or if the renal term were converted to the blood basis instead.

Enter Disease Effects Once, at the Mechanistic Level

For a joint prediction, assign each disease-related change to the parameter or pathway it plausibly affects. Renal impairment may reduce \(GFR\) or secretion capacity. Hepatic impairment may affect \(Q_H\), \(CL_{\mathrm{int},H}\), or both. A binding change can affect hepatic and renal clearance through their respective unbound fractions, while a change in \(R_b\) changes the conversion between blood- and plasma-based hepatic clearance. The resulting model may therefore show a joint effect that is not the product of two simple impairment ratios.

A central safeguard is to check how each input was obtained. If a disease-dependent hepatic clearance estimate already incorporates altered binding, do not additionally multiply it by a binding correction. Likewise, if an empirically estimated total nonrenal clearance already includes hepatic elimination, do not add a separate modeled hepatic pathway on top of it. A parameter's meaning and concentration basis must be explicit before it is combined with other terms.

Worked Example: Applying Shared Binding Changes Only Once

Consider a hypothetical drug with no other modeled clearance. In a reference state, let \(Q_H=90\) L/h, \(f_{u,b}=0.040\), \(CL_{\mathrm{int},H}=300\) L/h, and \(R_b=1.20\). Let \(f_{u,p}=0.080\), \(GFR=5.0\) L/h, and \(CL_{\mathrm{int},R\mathrm{sec}}=5.0\) L/h. In a joint impairment state, suppose \(Q_H=60\) L/h, \(f_{u,b}=0.080\), \(CL_{\mathrm{int},H}=180\) L/h, \(R_b=1.20\), \(f_{u,p}=0.160\), \(GFR=3.0\) L/h, and \(CL_{\mathrm{int},R\mathrm{sec}}=4.0\) L/h. These values are illustrative model inputs, not clinical estimates.

In the reference state, \(CL_{H,b}=90(0.040)(300)/[90+(0.040)(300)]=10.588\) L/h, so \(CL_{H,p}=1.20(10.588)=12.706\) L/h. Renal clearance is \(0.080(5.0+5.0)=0.800\) L/h, giving total plasma clearance \(13.506\) L/h.

In the joint state, \(CL_{H,b}=60(0.080)(180)/[60+(0.080)(180)]=864/74.4=11.613\) L/h. Thus \(CL_{H,p}=1.20(11.613)=13.935\) L/h. Renal clearance is \(0.160(3.0+4.0)=1.120\) L/h, and total plasma clearance is \(13.935+1.120=15.055\) L/h.

The higher unbound fractions are already included in the hepatic and renal equations. Multiplying either calculated pathway clearance by the same binding change again would double count it. The example also shows why the joint result cannot be inferred from labels such as “renal impairment” and “hepatic impairment” alone: the net change depends on the specified inputs and their interactions in the equations.

Exposure, Diagnostics, and Model Alternatives

Under linear, time-invariant conditions after an intravenous dose, with plasma-based clearance defined consistently, plasma AUC is dose divided by total plasma clearance. In the first worked example, a hypothetical 100 mg intravenous dose gives \(AUC=100/23.5=4.255\) mg·h/L. This relationship would need modification if clearance were nonlinear, changed over time, or if the dose were not intravenous and bioavailability were not accounted for.

Several checks can reveal a faulty joint model. Confirm that every clearance has the same concentration basis and units. Reconstruct the total elimination rate from pathway rates at a common concentration. Inspect whether binding or a disease covariate appears both inside a mechanistic pathway equation and as an extra correction factor. Check that empirically estimated “total” or “nonrenal” terms do not overlap with explicitly modeled pathways. Finally, ask whether each disease-dependent parameter has a defensible physiological or drug-specific interpretation.

A simple alternative is to multiply renal and hepatic impairment ratios derived separately. This may be a useful approximation when the ratios are defined on a common basis, are based on comparable conditions, and their effects are reasonably independent. It is not a general mechanistic rule: a shared binding change, altered blood-to-plasma partition, or nonlinear well-stirred hepatic response can make the joint effect depart from that product. A pathway-based PBPK model makes these dependencies visible, but requires more inputs and can be poorly identified when the data are sparse.

The simplified equations also have limits. They do not automatically represent dialysis, saturable transport, biliary elimination, time-varying disease, or changes in distribution and bioavailability. Such processes need explicit representation when relevant. A mathematically consistent clearance sum is necessary, but it does not establish that the physiological inputs or model structure are adequate for a particular patient.

Common Mistakes

  • Mixing concentration bases: adding hepatic blood clearance directly to renal plasma clearance is generally invalid. Convert one pathway so both use the same concentration basis.
  • Applying binding twice: if \(f_{u,b}\) or \(f_{u,p}\) is already in the pathway equation, do not add another binding multiplier unless it represents a distinct, justified effect.
  • Multiplying impairment ratios automatically: separate ratios need not combine multiplicatively when their underlying physiology interacts.
  • Double counting a pathway: an empirical nonrenal term and an explicit hepatic pathway cannot both represent the same hepatic elimination.
Key takeaway: Build the combined prediction from disease-dependent pathway inputs, include each effect once, and express all pathway clearances on one stated concentration basis before summing them.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A hepatic blood clearance is 8.0 L/h and \(R_b=1.25\). What is its plasma-based value?
  2. Why can hepatic blood clearance not generally be added directly to renal plasma clearance?
  3. In the simplified model, let \(f_{u,p}=0.12\), \(GFR=2.0\) L/h, and \(CL_{\mathrm{int},R\mathrm{sec}}=3.0\) L/h. Calculate renal clearance.
  4. A model uses an empirical nonrenal clearance that already includes hepatic elimination, then adds a mechanistic hepatic clearance. What is the concern?
  5. State one reason that multiplying independently estimated renal and hepatic impairment ratios may fail.

Solutions. (1) \(CL_{H,p}=R_bCL_{H,b}=1.25(8.0)=10.0\) L/h. (2) The clearances are referenced to different concentrations; convert the hepatic term using \(R_b\), or convert the renal term to blood basis, before adding. (3) \(CL_{R,p}=0.12(2.0+3.0)=0.600\) L/h. (4) The hepatic pathway may be counted twice, so the meaning of the empirical term must be clarified or the overlapping term removed. (5) The disease effects may interact through shared binding, blood-to-plasma partition, or another parameter, so the joint effect need not equal the product of separate ratios.