Why Names and Units Are Part of the Model
A shared QSP model is interpretable only if its entries communicate what is represented and how quantities are measured. An identifier such as IL6 might refer to a molecule, a concentration in plasma, a concentration in tissue, or an assay result. A parameter called clearance might have units of inverse time or volume per time, depending on the equation in which it is used. Those are not cosmetic differences: confusing them can change the mathematical model while leaving its equations superficially plausible.
Naming and unit conventions should therefore make model meaning explicit, not attempt to replace a definition or equation. The conventions in this tutorial are practical choices for a project; no single identifier grammar is universal. The essential requirement is that a reader can determine each entry's biological referent, location, quantity type, unit and role without guessing.
These terms should not be collapsed into one label. For example, “IL-6” names an entity, “IL-6 concentration in plasma” names a state, and “serum IL-6 immunoassay result” names an observable that may be calculated from one or more states. A state is not automatically identical to what an assay measures.
A Practical Naming Convention
A useful convention combines a readable identifier with structured documentation. For state variables, a project might use the pattern entity_compartment_form, such as IL6_plasma_free. A complex can name both partners, as in IL6_IL6R_plasma_complex. The exact ordering is less important than applying it consistently and defining it in a project glossary.
Identifiers should be stable, unique, and specific enough to distinguish states that differ in location, molecular form, or biological meaning. Use a documented spelling and capitalization scheme; avoid spaces and punctuation that different software may parse differently. If an entity has species-of-origin or species-of-experiment distinctions, record them explicitly, for example in a species field or a controlled prefix such as human_IL6 and mouse_IL6. Do not rely on a reader to infer whether “species” means organism or molecular species.
For parameters, names should communicate the process and, where needed, the interacting entities or compartment. Examples include kon_IL6_IL6R_plasma and koff_IL6_IL6R_plasma. The identifier alone does not fully define a parameter. Its record should also state its description, units, scope, value or source, and any conditions or assumptions that delimit its use. Provenance belongs with that record; a plausible name is not evidence for a value.
Observable identifiers should distinguish the measured quantity from internal states. A name such as serum_IL6_total_pg_per_mL signals a matrix, quantity, and reporting unit. Documentation should additionally specify the assay or measurement definition, any transformation from model states, and whether the readout is total, free, or otherwise selective. If an observable uses a unit different from the underlying state, the conversion and its assumptions must be explicit.
Worked Example: Resolving Ambiguous State Names
A hypothetical shared model contains entries named R, IL6, and complex. The team says the model represents receptor binding in plasma and liver tissue. Determine what needs clarification and propose unambiguous state identifiers.
Clarify the entities and locations. “R” does not say which receptor is represented, and “IL6” does not say whether it is free or bound, or where it is located. “Complex” does not identify its components or location. Since the model has two compartments, a single unqualified name cannot distinguish plasma from liver tissue.
Make the distinctions explicit. Suppose the intended entities are IL-6 and IL-6 receptor, and all three states are concentrations. A clear set for the plasma binding system is IL6_plasma_free, IL6R_plasma_free, and IL6_IL6R_plasma_complex. The corresponding liver states use liver in place of plasma.
Record quantity type and units. Add metadata stating that each entry is a concentration in nM, rather than assuming that from the name. If the receptor state represents surface-bound receptor rather than a well-mixed concentration, revise its quantity definition and unit instead of retaining a misleading concentration label.
The names now distinguish location, free versus complexed form, and the components of the complex. They do not, by themselves, establish that plasma and liver binding have identical parameters or that the states are biologically interchangeable.
Units Must Follow the Equations
A unit convention is useful only when it is checked against the mathematical definition. For each state, record whether it is an amount, concentration, fraction, or another quantity, together with its unit. Record parameter units as well. State the model's base time unit, and keep conversions explicit when combining data or processes expressed on different scales.
Dimensional analysis provides a direct check: every term added or subtracted in an equation must have the same units. For a concentration state \(C\) measured in nM, its derivative with respect to time measured in h has units nM/h. In the equation below, \(k_{\mathrm{in}}\) must therefore be in nM/h, and \(k_{\mathrm{out}}\) must be in h\(^{-1}\).
Indeed, \(k_{\mathrm{out}}C\) has units \(\mathrm{h}^{-1}\times\mathrm{nM}=\mathrm{nM}/\mathrm{h}\), matching both the derivative and the input term. This check can expose an identifier that suggests one quantity while the equation uses another.
Worked Example: Diagnosing a Clearance Unit Mismatch
A hypothetical plasma drug concentration \(C\) is \(4.00\) nM, plasma volume \(V\) is \(1.50\) L, and systemic clearance \(CL\) is \(0.300\) L/h. The model uses \(dC/dt=-CL\,C/V\). Check the units and calculate the initial rate of change.
Check the units. The right-hand side has units \[ \frac{(\mathrm{L}/\mathrm{h})(\mathrm{nM})}{\mathrm{L}}=\mathrm{nM}/\mathrm{h}, \] which matches the derivative of a concentration measured in nM with time in hours.
Substitute the values. \[ \frac{dC}{dt}=-\frac{(0.300\ \mathrm{L}/\mathrm{h})(4.00\ \mathrm{nM})}{1.50\ \mathrm{L}} =-0.800\ \mathrm{nM}/\mathrm{h}. \]
The predicted initial concentration decrease is \(0.800\) nM per hour. The equivalent first-order rate constant is \(CL/V=0.300/1.50=0.200\) h\(^{-1}\), so the same equation can be written \(dC/dt=-k_{\mathrm{el}}C\) with \(k_{\mathrm{el}}=0.200\) h\(^{-1}\). Labelling \(CL\) as \(k_{\mathrm{el}}\) would be misleading: the numerical value and units of clearance differ from those of the first-order rate constant.
Justifying a Unit Conversion
A shared model may use concentrations in a mechanistic state while an assay reports mass per volume. A conversion requires a molecular mass and an explicit assumption about what molecular form the measurement represents. For a species with molecular mass \(M\) in g/mol and concentration \(C\) in mol/L, the mass concentration is \(CM\) in g/L. Unit conversion then gives the assay-scale value. This dimensional relation does not establish that the assay detects that species selectively; that is a separate measurement assumption.
For example, a molecular mass of 20 kDa is \(20{,}000\) g/mol. A concentration of 1 nM is \(10^{-9}\) mol/L, so its mass concentration is \(20{,}000\times10^{-9}=2.00\times10^{-5}\) g/L. Since \(1\) g/L equals \(10^6\) ng/mL, this is \(20.0\) ng/mL, or \(20{,}000\) pg/mL.
Worked Example: Mapping Model States to an Assay Observable
In a hypothetical model, free IL-6 in a sample is \(0.0180\) nM and complex-associated IL-6 is \(0.0120\) nM. A hypothetical assay reports total IL-6 in pg/mL. Assume a molecular mass of 20 kDa and, for this example, assume the assay responds equally to free and complex-associated IL-6.
Define the model-side quantity. Under the stated measurement assumption, total IL-6 is the sum of the two states: \[ C_{\mathrm{total}}=0.0180+0.0120=0.0300\ \mathrm{nM}. \]
Convert the unit. The molecular mass implies \(1\) nM \(=20{,}000\) pg/mL. Therefore, \[ 0.0300\ \mathrm{nM}\times20{,}000\ \frac{\mathrm{pg}/\mathrm{mL}}{\mathrm{nM}} =600\ \mathrm{pg}/\mathrm{mL}. \]
A suitable observable name is sample_IL6_total_pg_per_mL, with documentation stating that it is calculated from free plus complex-associated IL-6 and converted using the assumed molecular mass. If the assay does not detect complex-associated IL-6 equally, this mapping is not justified; the observable must instead represent the assay's response to the relevant forms.
Reviewing Shared Entries Systematically
A structured review helps identify problems before a collaborator reuses an entry. For each state, parameter and observable, ask what it represents, where it applies, what quantity it is, which units it uses, and how it connects to equations or measurements. When a model contains multiple organisms, tissues, molecular forms, or assay readouts, check that the distinctions are expressed consistently throughout equations, tables, plots, and documentation.
State the biological entity or process in words, including molecular form and organism when those distinctions matter.
Specify compartment and whether the entry is an amount, concentration, fraction, or another quantity.
Verify the dimensions of every equation term and ensure the time unit is consistent.
Document which states contribute, any transformations or conversions, and assumptions about assay response.
Common failure signals include identifiers reused for different entities, a compartment omitted where several are modelled, equations whose terms have incompatible dimensions, and plot labels that omit a unit or report a different unit from the underlying state. These are reasons to inspect the definition and model logic. They are not necessarily fixed by changing a display label.
A controlled glossary with stable identifiers and explicit units is generally easier to review than free-form names alone. Conversely, a rigid naming grammar can become unwieldy when it tries to encode every detail in a single string. A practical compromise is a concise, unique identifier supported by structured fields for entity, compartment, quantity type, unit, role, and description. Unit conversion can likewise be handled either by converting inputs to a common internal unit or by preserving source units and applying explicit conversions. In either case, conversions must be visible and testable.
Common Mistakes and Limitations
- Treating a familiar name as a definition. “Clearance” is not enough to determine whether a parameter is a volume-per-time clearance or a first-order rate constant. State its equation role and units.
- Using a molecule name for an assay output. An assay may measure total, free, or selectively detected material. Name and document the observable accordingly.
- Changing units without changing the numeric value. A concentration expressed in nM cannot be relabelled as pg/mL without a conversion and the assumptions needed for that conversion.
- Assuming matching names imply matching meaning. Two collaborators may use the same identifier for different compartments or molecular forms. Compare definitions and equations, not labels alone.
- Assuming unit consistency proves biological validity. Dimensional consistency is necessary for a coherent equation, but it does not establish that the chosen mechanism, parameter value, or measurement mapping is biologically appropriate.
Naming standards also have limits. They cannot resolve a genuinely ambiguous biological definition, establish the correct compartment structure, or guarantee that a parameter is supported by evidence. Those questions require explicit model documentation and scientific review. The purpose of consistent names and units is to make such questions visible and to let collaborators check whether they are using the same model concepts.
Check Your Understanding
Answer each question, then compare with the solutions below.
- A model has free drug concentrations in plasma and liver, both recorded in nM. Give two unambiguous state identifiers and name the additional metadata each needs.
- An equation is \(dC/dt=-kC\), where \(C\) is measured in mg/L and time in days. What units must \(k\) have?
- A parameter called CL is assigned units h\(^{-1}\) in an equation \(dC/dt=-CL\,C\). What is the likely terminology problem?
- A model predicts free biomarker concentration, while an assay reports total biomarker mass per volume. Name two requirements for a defensible observable mapping.
Solutions. (1) For example, drug_plasma_free and drug_liver_free. Each needs a definition of the drug form and quantity type, units of nM, and a clear compartment definition. (2) Since \(dC/dt\) has units mg/L/day and \(C\) has units mg/L, \(k\) must have units day\(^{-1}\). (3) In this equation the parameter is a first-order elimination rate constant, not clearance in the usual volume-per-time sense; name it accordingly, such as \(k_{\mathrm{el}}\), and document h\(^{-1}\). (4) Specify which model states contribute to the measured total and provide the unit conversion, including molecular-mass assumptions; also justify how the assay responds to free and any bound forms.