Tutorials › Quantitative Systems Pharmacology (QSP) › Dimensional Audits of Composite Mechanistic Rate Expressions

Conservation Laws and Dimensional Consistency · Tutorial 65 of 799

Dimensional Audits of Composite Mechanistic Rate Expressions

You will learn to determine whether every term in a composite mechanistic rate has the required units, and to infer or correct parameter units when it does not.

Intermediate 9 min read

What You'll Learn

  • Distinguish concentration rates from amount rates when auditing model equations.
  • Infer parameter units for concentration products and saturable functions.
  • Check whether terms in a sum and factors in a product are dimensionally compatible.
  • Evaluate the effect of dimensionless and dimensional scaling factors.
  • Identify unit errors that numerical calculations alone may not reveal.
  • Apply a repeatable audit to composite rate expressions.

Why Audit a Composite Rate?

A mechanistic rate expression can look plausible and still be wrong if its terms have incompatible units. That error can change a predicted rate by orders of magnitude, or be hidden by a numerical implementation that silently treats every quantity as unitless. A dimensional audit checks whether the expression could represent the biological quantity assigned to it before parameter values are used.

This tutorial focuses on composite rates: expressions that combine concentration products, saturable terms, and scaling factors. The earlier material on deriving units for rates, clearances, and turnover parameters establishes the basic distinction between a concentration rate and an amount rate. Here, the key task is to carry those units through a whole expression, term by term.

Definition: A dimensional audit assigns units to every state variable, parameter, and factor in an expression, then uses multiplication, division, and addition rules to check whether the expression has the units required by its meaning. A concentration rate has units of concentration per time; an amount rate has units of amount per time. These are not interchangeable unless a volume conversion is included.

For examples, let concentration be measured in \(\mathrm{\mu M}\) and time in \(\mathrm{h}\). A rate of change of concentration must then have units \(\mathrm{\mu M\,h^{-1}}\). The numerical values and mechanisms below are hypothetical. Changing to another consistent unit system changes the numerical values of parameters, not the dimensional requirements.

Rules for Composite Expressions

Treat units algebraically. Products multiply units, quotients divide them, and a sum is valid only when every summand has the same units. A ratio such as \(C/(K_m+C)\) is dimensionless only if \(C\) and \(K_m\) are both concentrations. The numerator and denominator must also be commensurate: adding a concentration to a time or a rate constant is not defined.

A scaling factor must be audited like any other quantity. It may be dimensionless, or it may carry units needed to convert a rate per time into a concentration rate. Do not infer that a factor is dimensionless merely because it is called a “scale,” “multiplier,” or “gain.” State its units or derive them from the expression and the quantity it represents.

Audit principle: If an expression represents a concentration rate, its complete unit product must reduce to concentration per time. Each term being added must independently have those units. If it represents an amount rate, the target is amount per time instead.

A useful derivation is to start with a bimolecular concentration product \(A B\), where both \(A\) and \(B\) have units of concentration. If \(k\) is the parameter in \(kAB\) and the result is a concentration rate, then

$$ [k][A][B] = \frac{\mathrm{concentration}}{\mathrm{time}} \quad\Longrightarrow\quad [k] = \frac{\mathrm{concentration/time}}{\mathrm{concentration}^2} = \mathrm{concentration}^{-1}\mathrm{time}^{-1}. $$

Thus a second-order concentration product does not use a parameter with units of inverse time; it requires inverse concentration per time. More generally, if the product contains \(m\) concentration factors and no other dimensional factors, its rate coefficient has units \(\mathrm{concentration}^{1-m}\mathrm{time}^{-1}\). This follows by solving the unit equation for the coefficient, not by memorizing a rule detached from the expression.

Worked Example: Auditing a Reversible Binding Rate

Suppose a hypothetical binding reaction has forward and reverse contributions combined as \(v = k_{\mathrm{on}}RL-k_{\mathrm{off}}RL_{\mathrm{bound}}\). Here \(R\), \(L\), and \(RL_{\mathrm{bound}}\) are concentrations, and \(v\) is the net concentration rate of bound complex. Determine the required parameter units and calculate \(v\) for \(R=2.0\,\mathrm{\mu M}\), \(L=3.0\,\mathrm{\mu M}\), \(RL_{\mathrm{bound}}=1.0\,\mathrm{\mu M}\), \(k_{\mathrm{on}}=0.040\,\mathrm{\mu M^{-1}h^{-1}}\), and \(k_{\mathrm{off}}=0.10\,\mathrm{h^{-1}}\).

The forward term must have units \(\mathrm{\mu M\,h^{-1}}\). Since \(RL\) has units \(\mathrm{\mu M^2}\), \(k_{\mathrm{on}}\) must have units \(\mathrm{\mu M^{-1}h^{-1}}\). The reverse term contains one concentration, so \(k_{\mathrm{off}}\) must have units \(\mathrm{h^{-1}}\). Both terms are therefore concentration rates and may be subtracted.

Substitution gives \(v=(0.040\,\mathrm{\mu M^{-1}h^{-1}})(2.0\,\mathrm{\mu M})(3.0\,\mathrm{\mu M})-(0.10\,\mathrm{h^{-1}})(1.0\,\mathrm{\mu M})\). The forward contribution is \(0.24\,\mathrm{\mu M\,h^{-1}}\), and the reverse contribution is \(0.10\,\mathrm{\mu M\,h^{-1}}\). Therefore, \(v=0.24-0.10=0.14\,\mathrm{\mu M\,h^{-1}}\), a net positive rate of complex formation at the stated concentrations.

Saturable Terms and Sums

A standard saturable contribution has the form \(V_{\max}C/(K_m+C)\). If \(C\) and \(K_m\) are concentrations, their ratio is dimensionless, so \(V_{\max}\) must have the units of the complete rate. For a concentration rate, \(V_{\max}\) has units of concentration per time. For an amount rate, it has units of amount per time. The expression does not determine which kind of rate is intended; the model definition does.

When a saturable term is added to a product term, check each term separately. The fact that the total is intended to be a rate does not make unlike units compatible. A dimensionless multiplier preserves the units of the expression it multiplies. A dimensional multiplier changes them and must be included in the unit calculation.

Worked Example: Combining Saturation and a Concentration Product

Consider the hypothetical concentration rate \(r=s\left(V_{\max}C/(K_m+C)+kAB\right)\). Let \(A\), \(B\), \(C\), and \(K_m\) be concentrations, let \(V_{\max}\) be a maximum concentration rate, and let \(s\) be dimensionless. Find the units required for \(k\), then evaluate the rate for \(s=1.5\), \(V_{\max}=1.2\,\mathrm{\mu M\,h^{-1}}\), \(C=2.0\,\mathrm{\mu M}\), \(K_m=3.0\,\mathrm{\mu M}\), \(A=0.50\,\mathrm{\mu M}\), \(B=4.0\,\mathrm{\mu M}\), and \(k=0.10\,\mathrm{\mu M^{-1}h^{-1}}\).

The saturation ratio is dimensionless because \(C\) and \(K_m\) have the same units. Its contribution therefore has units \(\mathrm{\mu M\,h^{-1}}\). The second contribution must have the same units: \(kAB\) has units \([k]\mathrm{\mu M^2}\), so \(k\) must have units \(\mathrm{\mu M^{-1}h^{-1}}\). The sum is dimensionally valid, and the dimensionless \(s\) leaves its units unchanged.

The numerical calculation is \(r=1.5\left((1.2)(2.0)/(3.0+2.0)+(0.10)(0.50)(4.0)\right)\). The saturable contribution is \(1.2(2.0/5.0)=0.48\,\mathrm{\mu M\,h^{-1}}\). The product contribution is \(0.10(0.50)(4.0)=0.20\,\mathrm{\mu M\,h^{-1}}\). Thus \(r=1.5(0.48+0.20)=1.02\,\mathrm{\mu M\,h^{-1}}\).

Scaling Factors With Units

A factor need not be dimensionless. For example, multiplying a specific rate with units of inverse time by a concentration scale produces a concentration rate. Such a factor may represent a pool or abundance scale, but its biological interpretation must be specified separately from its dimensional role. If a model uses a conversion between amount and concentration, the relevant volume units must also appear; a bare numerical multiplier cannot perform that conversion consistently.

For a product \(qg\) to have target units \(\mathrm{\mu M\,h^{-1}}\), if \(g\) has units \(\mathrm{h^{-1}}\), then \(q\) must have units \(\mathrm{\mu M}\). Conversely, if \(q\) is dimensionless, \(g\) itself must already have units \(\mathrm{\mu M\,h^{-1}}\). These are alternative parameterizations; the model definition determines which is intended.

Worked Example: Inferring a Concentration Scale

A hypothetical rate is written \(r=q\left(k_{\mathrm{sat}}C/(K_m+C)+k_{\mathrm{loss}}\right)\), and \(r\) is a concentration rate. Suppose \(C\) and \(K_m\) are concentrations, \(k_{\mathrm{sat}}=0.30\,\mathrm{h^{-1}}\), and \(k_{\mathrm{loss}}=0.05\,\mathrm{h^{-1}}\). Infer the units of \(q\), then calculate \(r\) when \(q=4.0\,\mathrm{\mu M}\), \(C=2.0\,\mathrm{\mu M}\), and \(K_m=3.0\,\mathrm{\mu M}\).

The saturation ratio is dimensionless, so \(k_{\mathrm{sat}}C/(K_m+C)\) has units \(\mathrm{h^{-1}}\). The loss term also has units \(\mathrm{h^{-1}}\); they can be added. Since the full bracket is a rate per time and the output is a concentration rate, \(q\) must have units of concentration.

Substituting the values gives \(r=(4.0\,\mathrm{\mu M})\left((0.30\,\mathrm{h^{-1}})(2.0/5.0)+0.05\,\mathrm{h^{-1}}\right)\). The saturation contribution inside the bracket is \(0.12\,\mathrm{h^{-1}}\), so the bracket is \(0.17\,\mathrm{h^{-1}}\). Hence \(r=4.0(0.17)=0.68\,\mathrm{\mu M\,h^{-1}}\).

If \(q\) had instead been defined as dimensionless, this same expression would have units \(\mathrm{h^{-1}}\), not \(\mathrm{\mu M\,h^{-1}}\). It would then fail to represent the stated concentration rate unless an additional concentration factor were included.

A Repeatable Audit and Its Limits

For a complicated expression, audit the algebraic structure before substituting numbers. The following sequence makes it less likely that an incompatible term is overlooked.

1
Name the target quantity.
Decide whether the expression is a concentration rate, amount rate, or another quantity, and write its required units.
2
List state and parameter units.
Include every scaling factor, conversion factor, and denominator quantity. Do not silently assume a multiplier is dimensionless.
3
Audit ratios and products.
Check that each sum in a numerator or denominator combines like units; reduce units through each product and quotient.
4
Check every additive term.
Each term in a sum must have the same units as the other terms and the required units of the complete rate.
5
Evaluate only after the audit.
Use values expressed in a consistent unit system, then report the numerical result with its units.

Dimensional analysis is a structural check, not a test of whether a mechanism is biologically correct. It cannot establish that the chosen pathway is present, that a parameter value is plausible, or that the model predicts observed data. Nor does a numerically reasonable output prove the units are correct: a program may calculate a number from incompatible quantities because it does not enforce unit algebra.

Two useful alternatives or complements are symbolic unit tracking and explicit conversion to a common unit system. Symbolic tracking is especially helpful for long products and nested expressions; conversion checks catch mistakes such as combining a concentration in \(\mathrm{mM}\) with one in \(\mathrm{\mu M}\) without converting. Neither replaces checking each additive term. If a unit mismatch appears, trace it back through the expression and correct a parameter definition, a missing factor, or the claimed target quantity rather than suppressing the mismatch.

A common mistake is to say, “The terms are both rates because they are added in the model.” That reasoning reverses the required check: addition is valid only after the units have been shown to match. For example, adding \(\mathrm{\mu M\,h^{-1}}\) to \(\mathrm{h^{-1}}\) is invalid even if both numbers happen to be around 0.1. A numerical result cannot repair a dimensional mismatch.

Key takeaway: Start with the units required of the rate, derive the units of each term and factor, and allow addition only between quantities with matching units. A valid audit is necessary for a coherent mechanistic expression, but it is not evidence by itself that the mechanism is biologically adequate.

Check Your Understanding

For each item, identify the unit requirement and compare your reasoning with the solutions below.

  1. If \(a\) and \(b\) are concentrations and \(kab\) is a concentration rate, what are the units of \(k\)?
  2. In \(V C/(K+C)\), what units must \(V\) have if \(C\) and \(K\) are concentrations and the expression is a concentration rate?
  3. A dimensionless factor multiplies a concentration rate. Does it change the rate units? What if the factor has units of concentration?
  4. Can a term in \(\mathrm{h^{-1}}\) be added to a term in \(\mathrm{\mu M\,h^{-1}}\) to form a concentration rate? Explain.

Solutions. (1) Since \([k]\mathrm{concentration}^2=\mathrm{concentration/time}\), \(k\) has units \(\mathrm{concentration^{-1}time^{-1}}\). (2) \(C/(K+C)\) is dimensionless, so \(V\) must have units of concentration per time. (3) A dimensionless factor preserves the rate units; a concentration factor multiplies them, producing concentration squared per time if applied directly to a concentration rate. (4) No. The terms have different units and cannot be added; the inverse-time term would need a concentration factor or another valid reformulation before it could contribute to a concentration rate.

Clinical Biostats