Tutorials › Quantitative Systems Pharmacology (QSP) › Connecting Biological Scales Without Duplicating States

Defining Model Scope and Boundaries · Tutorial 22 of 799

Connecting Biological Scales Without Duplicating States

Learn to assign ownership to model states and choose scale-appropriate mappings that preserve meaning, units, and biological totals.

Intermediate 11 min read

What You'll Learn

  • Distinguish model states from quantities derived or passed across scales.
  • Choose interface quantities with clear units and biological meaning.
  • Convert per-cell molecular quantities into population totals without double counting.
  • Map cellular production into a tissue-level source term.
  • Aggregate concentrations using amounts and compartment volumes.
  • Diagnose inconsistent units, overlapping populations, and duplicate states.

Why Connecting Scales Requires Care

A QSP model may represent a molecular interaction, the response of a cell, a tissue-level biomarker, and a whole-body clinical measurement. Connecting these levels is not simply a matter of placing more variables in one system. Each connection must specify what quantity crosses the interface, how it is converted, and which model component owns it. Otherwise, a model can count the same biological material or process twice, or pass a quantity whose meaning changes unnoticed between scales.

The earlier material on integrating evidence across biological scales motivates using different types of evidence together. Here the focus is model architecture: defining explicit mappings between representations. The earlier material on compartment boundaries provides the idea of an exchange interface; a scale interface also needs to connect quantities such as per-molecule occupancy, per-cell activity, and tissue-level amount or concentration.

Definition: A scale interface is a specified mapping that passes a quantity, or a function of quantities, between model components representing different biological scales. A valid interface identifies the source quantity, its units and biological meaning, the mapping rule, and the receiving quantity. A duplicate state occurs when two model components independently represent the same biological quantity or pool as if each were a distinct state.

A state is a quantity the model tracks as changing over time, such as the amount of a biomarker in a compartment. A mapping may instead produce a derived quantity, such as the concentration calculated from an amount and volume, or a per-cell molecular signal passed to a cell-response function. A parameter, such as the number of cells in a fixed tissue sample, specifies a model property rather than a changing biological pool. These distinctions help determine whether a quantity should be tracked, calculated, or supplied from outside the model.

Specify Ownership, Units, and Direction

For each biological quantity, assign one primary owner: the component in which that quantity is represented as a state. Other components may receive a mapped value, but should not independently track the same pool unless the model explicitly defines how the representations relate. For example, an amount and the concentration calculated from that amount are two representations of one pool, not two independent pools.

A useful interface record states: (1) the source and receiving scales; (2) the quantity being passed; (3) the source and receiving units; (4) the conversion or aggregation rule; and (5) whether the receiving quantity is an input, a derived quantity, or a new state. If material physically transfers between modeled compartments, the interface should also make clear which component loses it and which gains it. A signal that influences a response, by contrast, need not represent material transfer.

Write the mapping explicitly, even when it is simple. In general, if the source quantities are collected in \(x\), an interface can be written as \(z=g(x)\), where \(g\) is the defined mapping and \(z\) is the quantity received. The mapping must preserve the meaning needed by the receiving component. It need not preserve the numerical value: a number of molecules per cell and a total number of molecules in a tissue are different quantities linked by a cell count.

Worked Example: Passing Molecular Occupancy to a Cell Response

Suppose a hypothetical tissue contains \(1.00\times10^6\) cells. Each cell has \(2.00\times10^5\) copies of a receptor, and the molecular component estimates that 30.0% of receptors are occupied. Assume the cells have the same receptor count and occupancy. The cell component uses occupancy as a signal for a response score, defined for this example as 2.50 response units per unit occupancy.

Define the interface. The molecular component supplies the dimensionless occupied fraction, \(0.300\). The cell component uses this fraction as an input to its response calculation; it does not create a second time-varying state for the same receptor occupancy.

The occupied receptors per cell are

$$ 2.00\times10^5\ \frac{\text{receptors}}{\text{cell}}\times 0.300 =6.00\times10^4\ \frac{\text{occupied receptors}}{\text{cell}}. $$

Across the stated population, the implied total is

$$ 6.00\times10^4\ \frac{\text{receptors}}{\text{cell}} \times1.00\times10^6\ \text{cells} =6.00\times10^{10}\ \text{occupied receptors}. $$

The response score is \(2.50\times0.300=0.750\) response units. It is a distinct downstream output, not another representation of the number of occupied receptors. The population total is valid only under the stated uniformity assumption; if cells differ, the total must be calculated from the cell-specific occupancies or an appropriately weighted average.

From Cells to Tissue: Map Rates Before Concentrations

When cells produce or consume a substance, the scale interface often maps a per-cell rate into a total tissue rate. The units provide a direct check: multiplying a rate per cell by a cell count should produce an amount per time. If only a fraction reaches the receiving pool, that fraction must be included and its interpretation stated. A tissue concentration then requires a further conversion using the volume of the pool in which the amount is represented.

This separation prevents a frequent ambiguity: a production rate is not a concentration, and a concentration is not an amount. The total rate can be passed as a source term to a tissue amount balance. If concentration is the tissue state instead, the corresponding source term must be divided by the represented volume. These are alternative representations of the same source contribution, not two sources to add together.

Worked Example: Converting Cellular Secretion Into a Tissue Source

In a hypothetical tissue, cell density is \(2.00\times10^8\) cells/L and tissue volume is 0.200 L. Each cell secretes a biomarker at \(4.00\times10^{-18}\) mol/cell/h. Assume the density and rate are uniform and that all secreted biomarker enters the modeled tissue pool.

The represented cell count is

$$ 2.00\times10^8\ \frac{\text{cells}}{\text{L}}\times0.200\ \text{L} =4.00\times10^7\ \text{cells}. $$

The tissue-level source rate is therefore

$$ 4.00\times10^7\ \text{cells} \times4.00\times10^{-18}\ \frac{\text{mol}}{\text{cell}\cdot\text{h}} =1.60\times10^{-10}\ \text{mol/h}. $$

If the tissue model tracks concentration in mol/L, this source contributes

$$ \frac{1.60\times10^{-10}\ \text{mol/h}}{0.200\ \text{L}} =8.00\times10^{-10}\ \text{mol/L/h}. $$

The total source rate and the concentration source term describe the same secretion at different levels of representation. The model should use the form consistent with its tissue state, not add both as independent production terms. If only half the secreted biomarker entered the modeled pool, the mapped total rate would instead be \(0.50\times1.60\times10^{-10}=8.00\times10^{-11}\) mol/h.

Aggregation Across Tissue Regions

A whole-body or pooled measurement may summarize several tissue regions. For concentrations, a simple arithmetic average is generally inappropriate when the regions have different volumes. The conserved quantity is amount: for a region with concentration \(C_i\) and volume \(V_i\), its amount is \(V_iC_i\). If the regions are disjoint and the measured pool combines them, total amount is the sum of regional amounts. Dividing by total volume gives the pooled concentration.

Result: For disjoint regions \(i=1,\ldots,m\) containing the same measured substance, with concentrations \(C_i\) and represented volumes \(V_i\), the pooled concentration is $$ C_{\text{pool}}=\frac{\sum_{i=1}^{m}V_iC_i}{\sum_{i=1}^{m}V_i}. $$ This mapping assumes that the regions do not overlap, the concentrations refer to compatible pools and units, and the pooled measurement corresponds to their combined amount and volume.

The expression follows from amount balance. Region \(i\) contains amount \(A_i=V_iC_i\). For disjoint regions, \(A_{\text{total}}=\sum_i A_i=\sum_i V_iC_i\), while \(V_{\text{total}}=\sum_i V_i\). By the definition \(C_{\text{pool}}=A_{\text{total}}/V_{\text{total}}\), the weighted-average formula follows. The weights are volumes because each concentration is amount per unit volume.

Worked Example: Pooling Concentrations From Two Tissues

Consider two hypothetical, non-overlapping tissue regions. Region A has volume 0.400 L and biomarker concentration 10.0 ng/mL; region B has volume 1.60 L and concentration 2.00 ng/mL. Calculate the concentration in the pooled 2.00 L sample.

Convert the volumes to milliliters so they match the concentration units: \(0.400\) L is \(400\) mL and \(1.60\) L is \(1{,}600\) mL. The amounts are

$$ A_A=400\ \text{mL}\times10.0\ \text{ng/mL}=4{,}000\ \text{ng}, \qquad A_B=1{,}600\ \text{mL}\times2.00\ \text{ng/mL}=3{,}200\ \text{ng}. $$

Thus the combined amount is \(7{,}200\) ng, and the pooled volume is \(2{,}000\) mL. The pooled concentration is

$$ C_{\text{pool}}=\frac{7{,}200\ \text{ng}}{2{,}000\ \text{mL}} =3.60\ \text{ng/mL}. $$

An unweighted mean would be \((10.0+2.00)/2=6.00\) ng/mL, which does not represent this pooled sample because the two regions have unequal volumes. The weighted result is lower because most of the pooled volume comes from region B, which has the lower concentration.

Diagnostics, Alternatives, and Limitations

A scale interface is ready for use when its units cancel or convert as intended, its quantities refer to compatible biological pools, and its mapping agrees with the model’s state definitions. Useful checks include tracing each state to its owner, checking whether mapped rates have the receiving state’s required units, and verifying that aggregation uses non-overlapping populations. For an amount balance, compare material entering and leaving the represented pool; for an aggregate measurement, recompute it from component amounts rather than averaging incompatible summaries.

Three common interface choices have different uses. Passing a derived quantity, such as occupancy, is appropriate when the receiving component needs that signal but should not independently track the underlying molecular pool. Passing a flux or rate is appropriate when material enters or leaves a separately represented compartment. Passing an empirical mapping may be useful when evidence supports a relationship between scales but does not support a detailed mechanistic conversion. The empirical option should not be presented as a mechanistic mechanism, and it must still specify input, output, units, and applicable context.

A mapping is not automatically valid outside the conditions used to define it. Uniform cell density, constant cell count, complete transfer, and compatible measurement pools are assumptions, not universal truths. If these assumptions fail, the mapped quantity may be biased or may not represent the intended biology. Also, matching units is necessary but not sufficient: two variables can share units while referring to different pools, times, or populations.

Common pitfall: A model tracks the amount of a biomarker in tissue and also adds a separate concentration state for that same tissue, assigning both independent rates of change. The reasoning is that the amount and concentration are “different variables.” The correction is to recognize that, for fixed volume, concentration is the derived quantity \(C=A/V\). Track one as the state and calculate the other from it, or define a clear conversion if volume changes. Otherwise, the model can imply two incompatible values for one biological pool.

A final audit is to ask whether every cross-scale quantity has a defined source, receiving role, transformation, and unit, and whether any quantity is represented as an independent state in more than one component. This does not determine whether the biological boundary is adequate for every perturbation; it makes the chosen boundary explicit enough to inspect and test.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. A cell component supplies an occupied fraction to a tissue component. What information should the interface specify, and why is the fraction not automatically a second state?
  2. A population has \(3.00\times10^7\) cells, each producing \(2.00\times10^{-17}\) mol/h. Calculate the total production rate.
  3. Two disjoint pools have equal volumes, with concentrations 4.00 and 8.00 mg/L. What is their pooled concentration? What would change if the second pool had three times the volume of the first?
  4. Explain why multiplying a concentration by a volume can help diagnose whether a pooled concentration has been calculated correctly.

Solutions. (1) Specify the source, receiving component, dimensionless fraction, mapping rule, and whether it is used as an input or derived quantity. It is not automatically a second state because it may be a value calculated from the molecular component’s state and passed to the tissue component. (2) \(3.00\times10^7\ \text{cells}\times2.00\times10^{-17}\ \text{mol}/(\text{cell}\cdot\text{h})=6.00\times10^{-10}\) mol/h. (3) Equal volumes give \((4.00+8.00)/2=6.00\) mg/L. If the second pool has three times the volume, the result is \((1\times4.00+3\times8.00)/(1+3)=7.00\) mg/L. (4) Concentration times volume gives amount. Summing amounts across disjoint pools and dividing by total volume checks that the aggregate is consistent with the represented material.

Clinical Biostats