Tutorials › Quantitative Systems Pharmacology (QSP) › Setting Compartment Boundaries at Exchange Interfaces

Defining Model Scope and Boundaries · Tutorial 20 of 799

Setting Compartment Boundaries at Exchange Interfaces

Learn to define compartments around exchange interfaces and represent the flows, boundary conditions, and feedback needed to preserve the biology relevant to a QSP decision.

Intermediate 12 min read

What You'll Learn

  • Define an exchange interface by the flows and causal links crossing it, not only by anatomical labels.
  • Write compartment balances that distinguish internal exchange from external input and loss.
  • Choose between a coupled compartment, a prescribed boundary flux, and a fixed-concentration boundary.
  • Check boundary choices with units, conservation, and plausible perturbations.
  • Identify how circulating biomarkers can carry tissue signals and feed back across an interface.

Why Exchange Interfaces Matter

A compartment boundary is a modeling choice about what is inside the system and what crosses its edge. A tissue may exchange drug with plasma, release a biomarker into blood, and respond to circulating signals. If a boundary is drawn only around the tissue, some of those processes become boundary flows; if it includes tissue and plasma, the same processes become internal exchanges. The boundary should therefore be chosen to represent the flows and feedback needed for the decision, not simply to match an anatomical outline.

For example, a model intended to predict tissue drug exposure needs a representation of drug entry into and exit from that tissue. A model intended to interpret a circulating biomarker may also need the tissue source and the transport from tissue to blood. The central test is whether the chosen boundary lets the model account for relevant material crossing the interface and for any causal influence that crosses in the opposite direction.

Definition: An exchange interface is a modeled or conceptual surface across which material or information passes between compartments or between the modeled system and its surroundings. A compartment boundary specifies which states and processes are represented inside the model and how each relevant crossing is represented. A crossing may be an internal exchange, an externally prescribed input or output, or a condition imposed at the boundary.

Representing a Boundary With Flows

For a compartment containing an amount \(A_i\) of a substance, a balance states that its rate of change equals the rates entering minus the rates leaving. An internal exchange appears with opposite signs in the balances for the two compartments it connects. An external input or loss appears in only the balance for the modeled system. This distinction is essential: an internal transfer changes where material is located, but does not by itself change the total amount in the combined system.

A useful illustrative law for drug exchange between plasma and tissue is a permeability-limited net flux. Let \(C_p\) and \(C_t\) denote plasma and tissue concentrations, \(PS\) a permeability-surface-area product with units of volume per time, and \(K_{t:p}\) the tissue-to-plasma partition ratio at equilibrium. Define \(J_{p\to t}\) as positive when drug moves from plasma to tissue:

$$ J_{p\to t}=PS\left(C_p-\frac{C_t}{K_{t:p}}\right). $$

The flux has units of amount per time. It is zero when \(C_t=K_{t:p}C_p\), positive when the plasma-side term is larger, and negative when net movement is in the reverse direction. This is one possible exchange representation, not a universal law: a model may instead use measured transport, a flow-limited description, or a prescribed flux if the evidence and decision support that choice.

With amounts \(A_p=V_pC_p\) and \(A_t=V_tC_t\), where \(V_p\) and \(V_t\) are the corresponding volumes, example balances are:

$$ \frac{dA_p}{dt}=R_{\mathrm{in}}-CL\,C_p-J_{p\to t}, \qquad \frac{dA_t}{dt}=J_{p\to t}-L_t. $$

Here \(R_{\mathrm{in}}\) is an external input rate, \(CL\,C_p\) is an external elimination rate when \(CL\) is clearance, and \(L_t\) is any tissue loss rate. The notation makes the boundary accounting visible: exchange is subtracted from plasma and added to tissue, while input and elimination cross the boundary of the combined plasma-tissue system.

Worked Example: Checking a Drug Exchange Boundary

Suppose a hypothetical model has \(V_p=3\) L, \(V_t=12\) L, \(C_p=2\) mg/L, \(C_t=0.25\) mg/L, \(PS=0.6\) L/h, and \(K_{t:p}=0.5\). At this instant, let \(R_{\mathrm{in}}=0.4\) mg/h, \(CL=0.2\) L/h, and \(L_t=0\).

The net plasma-to-tissue flux is

$$ J_{p\to t}=0.6\left(2-\frac{0.25}{0.5}\right) =0.6(1.5)=0.9\ \mathrm{mg/h}. $$

The elimination rate is \(CL\,C_p=(0.2\ \mathrm{L/h})(2\ \mathrm{mg/L})=0.4\) mg/h. Therefore,

$$ \frac{dA_p}{dt}=0.4-0.4-0.9=-0.9\ \mathrm{mg/h}, \qquad \frac{dA_t}{dt}=0.9-0=0.9\ \mathrm{mg/h}. $$

The plasma amount is decreasing at this instant because more drug leaves plasma for tissue than enters it, after accounting for input and elimination. Tissue amount is increasing at the same rate. The combined amount changes at \(-0.9+0.9=0\) mg/h from these two compartment rates, but the model also has external input and elimination, which cancel here: \(0.4-0.4=0\) mg/h. The opposing exchange terms are consistent with conservation.

Deriving the Conservation Check

A simple derivation provides a strong diagnostic for compartment boundaries. Add the two balances above. Since the same internal flux appears once with a minus sign and once with a plus sign, it cancels:

$$ \frac{d(A_p+A_t)}{dt} = \left(R_{\mathrm{in}}-CL\,C_p-J_{p\to t}\right) + \left(J_{p\to t}-L_t\right) = R_{\mathrm{in}}-CL\,C_p-L_t. $$

This identity holds because \(J_{p\to t}\) is an internal transfer and the amounts use the same substance accounting on both sides of the interface. If the flux does not cancel, check whether its sign was reversed in one equation, whether the two terms represent the same quantity, or whether an unmodeled transformation or loss has been folded into one of them. If the model represents only one side of an interface, the flux need not cancel within that model: it is then a boundary input or output for the modeled system.

This check is not proof that every biological process is represented correctly. It establishes that the equations are internally consistent with the stated accounting. It cannot determine whether the chosen boundary includes all decision-relevant sources, sinks, or feedback.

Choosing the Boundary Condition

A boundary condition describes what the model assumes about a crossing. The useful choice depends on whether the other side of the interface must respond to the modeled state and intervention.

RepresentationWhat is specifiedAppropriate interpretation
Coupled compartmentsA flux depends on states on both sides of the interface.Both sides are modeled because their changing states or feedback matter.
Prescribed fluxA time-varying amount-per-time input or output is supplied.The crossing is relevant, but its value is treated as external to the modeled states.
Fixed concentrationThe concentration on one side is held at a specified value.The unmodeled side is approximated as a sufficiently stable reservoir for the intended use.

These are modeling assumptions, not interchangeable mathematical conveniences. A fixed concentration does not imply that the corresponding amount is fixed if the volume changes; it asserts a concentration condition at the interface. A prescribed flux gives an amount-transfer rate and may be inappropriate if that rate should change in response to a modeled concentration. A coupled representation requires states and parameters for both sides, but can represent state-dependent exchange directly.

Worked Example: A Tissue-Derived Circulating Biomarker

Consider a hypothetical tissue that produces a biomarker and releases it into plasma. Let \(A_t\) and \(A_p\) be tissue and plasma amounts in ng, with concentrations \(C_t=A_t/V_t\) and \(C_p=A_p/V_p\). Suppose the tissue-to-plasma flux is represented by \(J_{t\to p}=G(C_t-K_{t:p}C_p)\), where \(G=0.8\) L/h and \(K_{t:p}=2\). At one instant, \(C_t=10\) ng/L and \(C_p=2\) ng/L. Let tissue production be 10 ng/h, local tissue loss be 2 ng/h, and plasma clearance be \(k_pA_p\), with \(k_p=0.5\) h\(^{-1}\) and \(V_p=3\) L.

The net tissue-to-plasma flux is

$$ J_{t\to p}=0.8(10-2\times2)=0.8(6)=4.8\ \mathrm{ng/h}. $$

The plasma amount is \(A_p=V_pC_p=(3\ \mathrm{L})(2\ \mathrm{ng/L})=6\) ng, so plasma clearance is \(k_pA_p=(0.5\ \mathrm{h}^{-1})(6\ \mathrm{ng})=3\) ng/h. The instantaneous rates are therefore \(dA_t/dt=10-2-4.8=3.2\) ng/h and \(dA_p/dt=4.8-3=1.8\) ng/h. Their sum is 5 ng/h, equal to tissue production minus local loss and plasma clearance: \(10-2-3=5\) ng/h. The circulating concentration is not, by itself, the tissue production rate; it reflects production, exchange, and removal in the modeled system.

An interface can also carry feedback in the reverse direction. Suppose circulating biomarker concentration suppresses tissue production according to the illustrative relation \(q(C_p)=12/[1+C_p/(4\ \mathrm{ng/L})]\) ng/h. At \(C_p=2\) ng/L, production is \(12/(1+2/4)=8\) ng/h. If \(C_p\) rises to 4 ng/L, production becomes \(12/(1+4/4)=6\) ng/h. A model that replaces the tissue source with a constant 8 ng/h would miss this two ng/h change in production under that perturbation. If the feedback is relevant to the decision, plasma concentration must be available to the tissue production rule, whether or not the full plasma dynamics are modeled. In this example, a coupled plasma state is one way to represent that dependency; a justified approximation might instead supply the relevant plasma signal as an external input.

Worked Example: Choosing Between a Fixed Boundary and a Coupled Tissue

A hypothetical decision asks whether a drug intervention changes tissue exposure enough to alter a tissue response. In the baseline setting, plasma concentration is 1 mg/L and a tissue measurement is 0.4 mg/L. A team proposes to impose \(C_p=1\) mg/L as a fixed boundary and model tissue exchange from that value.

The fixed boundary can be adequate for a question limited to tissue behavior under a stable plasma concentration of 1 mg/L. It is not adequate, without further justification, for comparing interventions expected to change plasma concentration. For example, if the intervention raises \(C_p\) to 2 mg/L, a fixed value of 1 mg/L will not pass that change into the tissue exchange calculation. The boundary condition would suppress an intended effect of the intervention by construction.

For the intervention comparison, a coupled plasma-tissue representation is the direct choice if plasma concentration changes dynamically and influences tissue entry. If a separate model supplies the intervention-specific plasma concentration profile, the tissue model may instead treat that profile as a prescribed boundary input. In either case, the boundary representation must transmit the relevant change rather than silently hold the baseline value fixed. The choice is about which response is modeled inside and which is supplied from outside, not about whether one formulation is universally superior.

Diagnostics, Limitations, and Common Mistakes

Several checks help reveal a poorly specified exchange boundary. First, verify the units of every flux: concentration multiplied by volume per time gives amount per time. Second, confirm the sign convention by testing a simple case: a positive net flux into a compartment must increase that compartment’s amount and decrease the donor compartment’s amount. Third, add the balances for all compartments in a closed exchange group and check that internal transfers cancel. Finally, ask whether a planned perturbation changes a boundary condition that has been fixed or prescribed. If it should, the model may be excluding a relevant response.

A boundary can also be too broad or too narrow for the question. Including both sides of an interface can require additional states and assumptions; treating one side as external can be more economical, but may hide state dependence. Measurements at only one side may not identify the exchange law uniquely. A good fit to a measured concentration does not establish that the flux mechanism or boundary condition is correct. Such choices should be assessed against the intended prediction and available evidence, rather than inferred from fit alone.

  • Pitfall: treating exchange as disappearance. If drug moves from plasma into tissue, subtracting the flux from plasma but not adding it to tissue makes the amount vanish from the combined system. The correction is to represent the same internal flux with opposite signs, or explicitly define the crossing as an external loss if the destination is outside the model.
  • Pitfall: fixing a boundary that should respond. A constant plasma concentration may be convenient, but it cannot transmit a predicted intervention-induced plasma change. Use a coupled state or a prescribed intervention-specific concentration profile when the change matters.
  • Pitfall: assuming a biomarker measures its source directly. A blood concentration can depend on tissue production, transport, distribution, and clearance. Interpret it as a direct proxy for tissue production only under assumptions that justify ignoring or controlling those other processes.

The practical comparison is therefore not “more compartments are always better” versus “fewer compartments are always better.” A coupled compartment is useful when both sides’ states or cross-interface feedback are needed. A prescribed boundary is useful when the crossing is important but can be supplied independently. A fixed condition is useful when the unmodeled side behaves sufficiently like a stable reservoir for the intended question. Each can fail when its implied assumptions do not hold.

Key takeaway: Place a compartment boundary where the model can account for every relevant crossing. Represent internal exchange with consistent, opposite-signed fluxes; represent crossings from outside with explicit boundary conditions; and retain cross-interface feedback whenever the decision depends on it.

Check Your Understanding

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

  1. A drug flux is 3 mg/h from plasma into tissue. What signs should this flux have in the plasma and tissue amount balances, and what is its contribution to the combined amount rate?
  2. A tissue releases a biomarker into plasma, but the model includes only plasma. Give one condition under which a prescribed tissue-to-plasma flux could be reasonable, and one reason it might fail.
  3. For an exchange law \(J=PS(C_p-C_t/K)\), with \(PS=0.5\) L/h, \(C_p=1.2\) mg/L, \(C_t=0.4\) mg/L, and \(K=0.5\), calculate the flux and interpret its direction.
  4. A model fixes plasma concentration while comparing interventions that change plasma drug exposure. Explain the boundary problem and give two possible representations that could address it.

Solutions. 1. The flux contributes \(-3\) mg/h to the plasma balance and \(+3\) mg/h to the tissue balance, so its contribution to the combined amount rate is zero. 2. A prescribed flux could be reasonable if tissue release is known or supplied independently and is not expected to respond to the modeled plasma or intervention states. It may fail if release changes with tissue state or circulating feedback, because a fixed external flux cannot represent that dependence. 3. \(J=0.5[1.2-(0.4/0.5)]=0.5(0.4)=0.2\) mg/h. The positive sign means net movement from plasma to tissue under the stated convention. 4. The fixed concentration prevents the intervention-induced plasma change from affecting tissue exchange. One option is to model plasma and tissue as coupled compartments; another is to use an intervention-specific plasma concentration profile as a prescribed boundary input to the tissue model.

Clinical Biostats