Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Reporting Model Structure and Assumptions

PBPK in Regulatory Submissions · Tutorial 789 of 799

Reporting Model Structure and Assumptions

Learn to describe a PBPK model’s structure and assumptions precisely enough for a reviewer to understand what was simulated and how.

Intermediate 12 min read

What You'll Learn

  • Distinguish model structure from parameter values and evidence about them.
  • Describe compartments, connections, and blood-flow pathways in a model diagram.
  • State governing equations and explain how they represent mass balance.
  • Specify the simulated population, dosing regimen, initial conditions, and boundaries.
  • Make simplifying assumptions explicit and identify their implications.
  • Recognize reporting gaps that prevent a reviewer from reconstructing the analysis.

Why Structure and Assumptions Need to Be Reported

A PBPK prediction is interpretable only when a reviewer can tell what system was represented and what rules governed its simulation. A statement such as “a PBPK model predicted exposure in the target population” leaves crucial questions unanswered: which tissues were represented, how drug moved between them, how the dose entered the system, and what was assumed about the population and physiology? A model diagram alone is not enough if its compartments and arrows are not defined; equations alone are not enough if the simulated regimen and initial conditions are unclear.

This tutorial focuses on describing those elements. It does not replace an executive summary or a discussion of the evidence supporting a model. The reporting task here is to make the model’s implemented structure and assumptions legible: a reviewer should be able to trace the path from dose, through the represented system, to the predicted quantity.

Definition: Model structure is the set of represented compartments, connections, state variables, and governing relationships that determine how simulated quantities change. A model assumption is a simplifying condition imposed on that representation, such as instantaneous tissue partitioning or a fixed physiological flow. A clear report distinguishes these structural choices from parameter values, which specify quantities within the chosen structure.

Describe the Model as a Connected System

A useful model diagram identifies each compartment, the substance or quantity represented in it, and the direction of each relevant transfer. It should show the route of administration and the route by which drug reaches systemic circulation, as well as elimination or metabolism pathways represented in the model. Labeling a box “liver” is not sufficient if the equations treat it as a single well-mixed compartment; the report should say so. Likewise, distinguish compartments that represent blood from those representing tissue, and state whether a diagram depicts anatomical connections, mathematical connections, or both.

The diagram and its accompanying description should agree with the implemented model. If the model includes a gut absorption compartment, a portal pathway, and a liver compartment, the reported structure should make clear how material moves among them. If a tissue is omitted, combined with another compartment, or represented by a lumped residual compartment, say so. These descriptions help a reviewer determine whether a reported prediction corresponds to the system implied by the model rather than to a different, more detailed system suggested by a generic diagram.

Worked Example: Turning a Diagram into a Structural Description

Consider a hypothetical oral-drug model with a gut depot, liver, systemic arterial blood, venous blood, and a perfusion-limited muscle compartment. The gut depot empties into portal blood; portal blood reaches the liver; liver outflow joins systemic venous return; arterial blood supplies muscle; and venous blood returns to the central circulation. The liver is the only compartment with drug metabolism.

A reviewable description states the compartments and the paths explicitly: “An oral dose enters a gut depot, transfers to portal blood, and passes through the liver before reaching systemic circulation. Arterial blood supplies the muscle compartment, whose venous outflow returns to the central circulation. Drug is metabolized in the liver; no muscle metabolism is represented.” This conveys more than a list of organ names because it identifies both the route and the implemented elimination pathway.

The description should also state the model’s level of detail. For example, if the liver is represented as one well-mixed compartment, that is a structural choice, not a claim that the anatomical liver is physically homogeneous. The reviewer can then interpret the diagram and equations as a simplified representation rather than mistake the boxes for a complete anatomical map.

State Equations and the Meaning of Simplifications

For each dynamic compartment, report the state variable and its governing equation, or provide an equivalent specification that makes the implemented relationships unambiguous. A common starting point is mass balance: the rate of change of drug amount equals the rate entering minus the rate leaving, including any metabolism or other modeled loss. For a perfusion-limited tissue compartment \(i\), one possible equation is:

$$ V_i\frac{dC_i}{dt} = Q_i\left(C_a-\frac{C_i}{K_{p,i}}\right) - CL_{\mathrm{int},i}C_i $$

Here \(C_i\) is the tissue concentration, \(V_i\) is tissue volume, \(Q_i\) is blood flow to the tissue, \(C_a\) is arterial blood concentration, and \(K_{p,i}\) is the tissue-to-blood partition coefficient used in the model. Under the stated equilibrium assumption, the venous concentration leaving the tissue is \(C_i/K_{p,i}\). The term \(CL_{\mathrm{int},i}C_i\) represents a modeled metabolic loss when \(CL_{\mathrm{int},i}\) is expressed as a clearance acting on the tissue concentration. If metabolism is not included in that tissue, the loss term is absent. The report must define its own state variables and concentration conventions; symbols alone do not establish what the implementation means.

The equation follows directly from mass balance. The amount in the tissue is \(V_iC_i\), so its rate of change is \(V_i\,dC_i/dt\) when \(V_i\) is fixed. Blood enters at rate \(Q_iC_a\) and leaves at rate \(Q_iC_i/K_{p,i}\). Subtracting the leaving rate and the metabolic loss gives the equation above. This derivation also exposes assumptions: constant tissue volume, the stated relationship between tissue and venous concentrations, and the specified flow and metabolism terms.

Worked Example: Interpreting a Tissue Equation Numerically

Suppose a hypothetical muscle compartment has \(V_i=10\) L, \(Q_i=2\) L/h, \(C_a=3\) mg/L, \(C_i=6\) mg/L, and \(K_{p,i}=5\). Let the modeled metabolic clearance be \(CL_{\mathrm{int},i}=0.10\) L/h. Assume the compartment is perfusion-limited, its volume is fixed, and the equation above applies.

First calculate the venous concentration under the model’s partition assumption: \(C_i/K_{p,i}=6/5=1.2\) mg/L. The net blood-flow contribution is therefore:

$$ Q_i(C_a-C_i/K_{p,i}) = 2\ \mathrm{L/h}\,(3-1.2)\ \mathrm{mg/L} = 3.6\ \mathrm{mg/h}. $$

The modeled metabolic loss is \(CL_{\mathrm{int},i}C_i=(0.10\ \mathrm{L/h})(6\ \mathrm{mg/L})=0.60\ \mathrm{mg/h}\). Thus the net rate of increase in tissue amount is \(3.6-0.60=3.0\) mg/h, and:

$$ \frac{dC_i}{dt} = \frac{3.0\ \mathrm{mg/h}}{10\ \mathrm{L}} = 0.300\ \mathrm{mg\,L^{-1}\,h^{-1}}. $$

At this instant, the modeled tissue concentration is increasing at 0.300 mg/L per hour. This is a local rate, not a prediction that the concentration will continue to increase linearly: \(C_a\), \(C_i\), and the resulting flows may change as the system evolves. The arithmetic also checks units: each term on the right side of the amount-balance equation is mass per time.

Define the Population, Regimen, and Simulation Conditions

A report should identify the population represented by the simulation in clinically meaningful terms, including relevant physiological characteristics that define the modeled group. It should say whether the model represents, for example, adults with normal organ function or a specified patient group, rather than relying on a broad label such as “the clinical population.” The report should not imply that a simulation represents people or circumstances that were not specified in its setup.

For each simulated regimen, give the administered amount, route, timing, duration or number of administrations, and any relevant formulation or input assumptions. State whether the listed amount is the administered dose or the amount treated as available to the model. Also identify the simulated output and its units, such as a plasma concentration-time profile or an exposure summary. These details are needed to interpret what the model’s result describes.

Initial conditions and boundaries are equally important. Specify the initial amount or concentration in each dynamic compartment, any pre-dose or steady-state condition, and the simulation time interval. State how doses are introduced, what concentration or flux conditions apply at the system boundary, and how modeled material leaves the system. If a compartment begins empty, say so rather than leaving that condition implicit.

Worked Example: Reporting an Oral Dosing Schedule and Its Input

Suppose a hypothetical adult simulation uses a 100 mg oral dose at 0, 24, and 48 hours. The model represents each dose in a gut depot, uses first-order transfer with \(k_a=0.70\ \mathrm{h^{-1}}\), and assumes that 80% of the administered dose is available for absorption. The gut-depot amount is initially zero before the first dose, and the simulation ends at 72 hours.

A clear description states that doses are added to the depot at the three specified times and that the modeled systemic input rate from the depot is \(F_{\mathrm{abs}}k_aA_{\mathrm{dep}}\), where \(F_{\mathrm{abs}}=0.80\) and \(A_{\mathrm{dep}}\) is the depot amount. Between doses, the depot follows:

$$ \frac{dA_{\mathrm{dep}}}{dt}=-k_aA_{\mathrm{dep}}. $$

For the first dose, the initial depot amount immediately after administration is 100 mg. The initial modeled systemic input rate is therefore \((0.80)(0.70\ \mathrm{h^{-1}})(100\ \mathrm{mg})=56\ \mathrm{mg/h}\). The dose schedule must also specify how doses at 24 and 48 hours are applied: each adds 100 mg to the amount then remaining in the depot. The total amount made available for absorption across the three doses is \(3(0.80)(100\ \mathrm{mg})=240\) mg, under this assumed absorption fraction.

This description makes clear that 100 mg is the administered amount, not the initial systemic amount; it identifies the depot initial condition, dose events, absorption rule, and simulation horizon. It also avoids implying that all 300 mg enters systemic circulation.

Make Assumptions and Boundaries Reviewable

List assumptions that materially define the simulated system, rather than relying on a generic statement that “standard PBPK assumptions” were used. Examples include instantaneous tissue-blood equilibrium, well-mixed compartments, fixed or time-varying physiology, linear versus saturable processes, and whether a pathway or process is omitted. For each, briefly state where it applies and what it excludes. A concise, specific statement is more useful than an unqualified claim that the model is physiologically realistic.

The boundary conditions and assumptions should be internally consistent with the equations. For example, a model that assumes fixed tissue volume should not silently use a time-varying volume in its mass balance. A model that represents venous concentration as tissue concentration divided by a partition coefficient should explain that relationship rather than leaving the reviewer to infer it from the diagram. If external input, physiological conditions, or dosing change during a simulation, state when and how those changes are imposed.

Reporting principle: Describe each simplifying assumption in terms of what the model does, where the assumption applies, and what is not represented. This makes the scope of the prediction visible without treating an assumption as evidence that the simplification is harmless.

Check for Completeness and Common Failure Modes

A practical review asks whether a reader can reconstruct the model’s conceptual operation from the report. The diagram, equations, population description, regimen, and initial and boundary conditions should be mutually consistent. Their purpose is not to reproduce every implementation detail in prose, but to ensure that key choices are identifiable and that the predicted quantity has a clear meaning.

  • Unlabeled arrows. A diagram shows connections but not whether an arrow represents blood flow, absorption, or another transfer. Label the process and its direction.
  • Compartments without definitions. A box named for an organ may conceal whether it represents tissue, blood, or a combined space. State the represented quantity and mixing assumption.
  • Implicit dose handling. A nominal dose is reported without clarifying its route, input compartment, timing, or fraction treated as available to the model. Specify each item.
  • Unstated starting conditions. A concentration-time profile is shown, but the reader cannot tell whether it begins from zero, after a prior dose, or at steady state. Report the initial state.
  • Equations and diagrams that disagree. A pathway is shown but has no corresponding transfer in the equations, or an equation contains a process absent from the diagram. Reconcile the two descriptions.

A common mistake is to assume that naming a recognized PBPK structure is enough: “a whole-body model with standard compartments.” That phrase does not specify which compartments were included, how they connect, or which processes were simplified or omitted. The correction is to describe the implemented structure and its conditions explicitly. This tutorial’s focus is structural reporting; a detailed accounting of where parameter values came from belongs in the separate reporting of parameter provenance and supporting evidence.

The main alternatives are not simply “more detail” and “less detail.” A model may represent many organs separately or use a reduced structure in which some tissues are grouped. Either can be described adequately if its compartments, connections, equations, and intended simulated population are clear. A more detailed diagram does not by itself demonstrate that the model is appropriate for its use, and a concise model description does not excuse leaving assumptions unstated. Structural transparency lets a reviewer understand what was simulated; it does not, on its own, establish the adequacy of the choices.

Check Your Understanding

Answer the questions, then compare your reasoning with the solutions below.

  1. In a tissue equation, \(V_i\,dC_i/dt\) has units of mass per time. What units must \(Q_i(C_a-C_i/K_{p,i})\) have?
  2. A report says “three oral doses were simulated” but gives no dose times or initial depot condition. Name two missing details that affect interpretation.
  3. A model diagram shows a tissue metabolism arrow, but the tissue equation contains no loss term. What reporting issue should be resolved?
  4. In the oral-input example, what is the initial input rate following a 50 mg dose under the same \(F_{\mathrm{abs}}\) and \(k_a\)?

Solutions. (1) \(Q_i\) has units L/h and the concentration difference has units mg/L, so their product has units mg/h, consistent with the amount-balance equation. (2) The report should give the timing or interval of the doses and the initial depot amount or condition; it should also clarify the dose route and input representation. (3) The diagram and equation are inconsistent as reported. Either include and define the modeled metabolism term or correct the diagram to show that tissue metabolism is not represented. (4) The rate is \(F_{\mathrm{abs}}k_aA_{\mathrm{dep}}=(0.80)(0.70\ \mathrm{h^{-1}})(50\ \mathrm{mg})=28\ \mathrm{mg/h}\) immediately after the dose enters an otherwise empty depot.