Tutorials › Quantitative Systems Pharmacology (QSP) › State-Based Models of Multispecific Binding

Bispecific Antibodies and Antibody-Drug Conjugates · Tutorial 460 of 799

State-Based Models of Multispecific Binding

Learn to translate two-antigen binding into explicit antibody states and mass-action transitions, then use those transitions to calculate and check concentration changes.

Advanced 10 min read

What You'll Learn

  • Define free, singly bound, and doubly bound bispecific antibody states
  • Write both routes by which the doubly bound state can form
  • Translate association and dissociation reactions into differential equations
  • Calculate instantaneous species rates from hypothetical concentrations
  • Check a state model for consistent flows and identify its assumptions

Why Represent Binding as States?

A bispecific antibody can bind two different antigens, which we will call \(X\) and \(Y\). Tracking only “free antibody” and “bound antibody” loses important information: an antibody bound to \(X\) alone is not in the same state as one bound to both \(X\) and \(Y\). Those states have different available binding sites and different routes into and out of them.

A state-based model makes these distinctions explicit. It is a set of species and transitions: each association or dissociation moves antibody between states, while the concentrations of free antigens determine the association rates. The model below describes binding kinetics only; processes such as internalization or antibody distribution would require additional states or fluxes.

Definition: Let \(A\) denote free bispecific antibody, and let \(X\) and \(Y\) denote free antigens. The antibody states are \(A\), free antibody; \(AX\), antibody bound to \(X\) only; \(AY\), antibody bound to \(Y\) only; and \(AXY\), antibody bound to both antigens. Each symbol denotes a concentration when used in a rate equation.

The model assumes one available binding site for each antigen on each antibody, binding of \(X\) and \(Y\) to their respective sites, and no direct interaction between the antigens. It also assumes a well-mixed system in which mass-action rates are appropriate. For the binding-only equations, synthesis, degradation, internalization, and transport are omitted. These assumptions describe a minimal network, not every possible bispecific construct or biological environment.

Transitions Between the Four Antibody States

There are two ways to form the doubly bound state. Free antibody can bind \(X\) and then \(Y\), or it can bind \(Y\) and then \(X\). Dissociation reverses each of those steps. The reaction network is:

$$ \begin{aligned} A + X &\xrightleftharpoons[k_{X,\mathrm{off}}]{k_{X,\mathrm{on}}} AX, &\qquad A + Y &\xrightleftharpoons[k_{Y,\mathrm{off}}]{k_{Y,\mathrm{on}}} AY,\\ AX + Y &\xrightleftharpoons[k_{Y,\mathrm{off}}]{k_{Y,\mathrm{on}}} AXY, &\qquad AY + X &\xrightleftharpoons[k_{X,\mathrm{off}}]{k_{X,\mathrm{on}}} AXY . \end{aligned} $$

Here \(k_{X,\mathrm{on}}\) and \(k_{Y,\mathrm{on}}\) are association rate constants, with units of concentration\(^{-1}\) time\(^{-1}\); the dissociation rate constants \(k_{X,\mathrm{off}}\) and \(k_{Y,\mathrm{off}}\) have units of time\(^{-1}\). Thus an association flux is proportional to both reactant concentrations. For example, \(k_{X,\mathrm{on}}AX\) means \(k_{X,\mathrm{on}}\times[A]\times[X]\), not the product of a rate constant and the state label alone.

This simplest version uses the same \(X\)- and \(Y\)-specific rate constants whether the other site is unoccupied or occupied. That is a modeling assumption: if binding in the doubly bound context has different kinetics, separate state-specific rate constants are needed. The equations below use the simpler assumption so that each transition has a clearly defined flux.

Definition: Under mass action, the association flux for a reaction between two species is the association rate constant multiplied by the concentrations of both reactants. A dissociation flux is the dissociation rate constant multiplied by the concentration of the bound species.

From Transition Fluxes to Differential Equations

The concentration rate for a state is the sum of fluxes entering it minus the sum of fluxes leaving it. For instance, \(AX\) is formed when \(A\) binds \(X\), and is depleted when \(AX\) loses \(X\) or binds free \(Y\). It is also formed when \(AXY\) loses \(Y\). Applying the same accounting to all four antibody states gives:

$$ \begin{aligned} \frac{dA}{dt} &= -k_{X,\mathrm{on}}AX-k_{Y,\mathrm{on}}AY +k_{X,\mathrm{off}}AX+k_{Y,\mathrm{off}}AY,\\ \frac{d(AX)}{dt} &= k_{X,\mathrm{on}}AX-k_{X,\mathrm{off}}AX -k_{Y,\mathrm{on}}(AX)Y+k_{Y,\mathrm{off}}AXY,\\ \frac{d(AY)}{dt} &= k_{Y,\mathrm{on}}AY-k_{Y,\mathrm{off}}AY -k_{X,\mathrm{on}}(AY)X+k_{X,\mathrm{off}}AXY,\\ \frac{d(AXY)}{dt} &= k_{Y,\mathrm{on}}(AX)Y+k_{X,\mathrm{on}}(AY)X -(k_{X,\mathrm{off}}+k_{Y,\mathrm{off}})AXY. \end{aligned} $$

In these equations, a juxtaposition such as \(AX\) in a rate term means the product of the concentrations of free \(A\) and free \(X\); the state \(AX\) in a derivative or dissociation term means the concentration of the \(X\)-bound antibody. To make that distinction explicit when calculating, one may write concentrations in square brackets, such as \([A]\), \([X]\), and \([AX]\). The equations use \(A\), \(X\), and the state labels as shorthand for those concentrations.

The terms involving \(AXY\) are especially informative. Its formation has two incoming fluxes, one from each singly bound state. Its loss has two outgoing fluxes: dissociation of \(X\) returns antibody to \(AY\), while dissociation of \(Y\) returns it to \(AX\). The destination depends on which antigen leaves.

Worked Example: Listing the States and Their Transitions

Worked Example: Listing the States and Their Transitions

Suppose a hypothetical bispecific antibody has one site for antigen \(X\) and one for antigen \(Y\). Write the antibody states and identify the transition that forms \(AXY\) from \(AX\).

The four antibody states are \(A\), \(AX\), \(AY\), and \(AXY\). The free antigens \(X\) and \(Y\) are also species in the model, but are not additional antibody states.

For formation from \(AX\), free \(Y\) associates with the available \(Y\)-binding site:

$$ AX + Y \xrightarrow{k_{Y,\mathrm{on}}} AXY, \qquad v_{\mathrm{form\ from\ }AX}=k_{Y,\mathrm{on}}[AX][Y]. $$

The reverse transition is \(AXY \xrightarrow{k_{Y,\mathrm{off}}} AX+Y\), with flux \(k_{Y,\mathrm{off}}[AXY]\). The other formation route is \(AY+X\to AXY\); it must also be included because the doubly bound state can be reached by either binding order.

Worked Example: Calculating Instantaneous Rates

Worked Example: Calculating Instantaneous Rates

Consider hypothetical concentrations, all in nM: \([A]=0.80\), \([X]=4.0\), \([Y]=3.0\), \([AX]=0.20\), \([AY]=0.10\), and \([AXY]=0.050\). Let \(k_{X,\mathrm{on}}=0.020\) nM\(^{-1}\) min\(^{-1}\), \(k_{Y,\mathrm{on}}=0.010\) nM\(^{-1}\) min\(^{-1}\), \(k_{X,\mathrm{off}}=0.040\) min\(^{-1}\), and \(k_{Y,\mathrm{off}}=0.030\) min\(^{-1}\). Calculate the instantaneous rate of change of \(AXY\).

The two formation fluxes are:

$$ \begin{aligned} v_{AX\to AXY} &= k_{Y,\mathrm{on}}[AX][Y] = (0.010)(0.20)(3.0) = 0.0060\ \mathrm{nM/min},\\ v_{AY\to AXY} &= k_{X,\mathrm{on}}[AY][X] = (0.020)(0.10)(4.0) = 0.0080\ \mathrm{nM/min}. \end{aligned} $$

The two dissociation fluxes sum to:

$$ (k_{X,\mathrm{off}}+k_{Y,\mathrm{off}})[AXY] =(0.040+0.030)(0.050) =0.0035\ \mathrm{nM/min}. $$

Therefore:

$$ \frac{d[AXY]}{dt} =0.0060+0.0080-0.0035 =0.0105\ \mathrm{nM/min}. $$

At these concentrations, \(AXY\) is increasing at an instantaneous rate of 0.0105 nM/min. This is a local rate at the stated moment, not a prediction that the concentration will rise by this amount every minute; the fluxes change as concentrations change.

A Short Derivation and a Bookkeeping Check

The equation for \(AXY\) follows directly by counting its incoming and outgoing transitions. Its incoming transitions are the association of \(Y\) with \(AX\) and the association of \(X\) with \(AY\). Its outgoing transitions are dissociation of \(X\) and dissociation of \(Y\). Under the mass-action assumptions, the four fluxes are respectively \(k_{Y,\mathrm{on}}[AX][Y]\), \(k_{X,\mathrm{on}}[AY][X]\), \(k_{X,\mathrm{off}}[AXY]\), and \(k_{Y,\mathrm{off}}[AXY]\). Subtracting outgoing from incoming fluxes gives:

$$ \frac{d[AXY]}{dt} = k_{Y,\mathrm{on}}[AX][Y] +k_{X,\mathrm{on}}[AY][X] -(k_{X,\mathrm{off}}+k_{Y,\mathrm{off}})[AXY]. $$

A useful check for a binding-only model is that each antibody transition transfers antibody between states rather than creating or destroying it. Adding the four antibody-state equations therefore cancels every transition flux:

$$ \frac{d}{dt}\bigl([A]+[AX]+[AY]+[AXY]\bigr)=0. $$

This check can reveal an omitted or incorrectly signed transition. It applies only to the stated binding-only system; adding processes such as internalization or degradation changes the total-antibody balance.

Worked Example: A Short-Euler Update

Worked Example: A Short-Euler Update

Use the concentrations and rate constants from the previous example. The corresponding instantaneous derivatives for the four antibody states are \(d[A]/dt=-0.0770\), \(d[AX]/dt=0.0515\), \(d[AY]/dt=0.0150\), and \(d[AXY]/dt=0.0105\) nM/min. Approximate the state concentrations after a short interval of 2.0 min using one forward-Euler step.

A forward-Euler update uses the initial derivative as a constant over the short interval: new concentration = initial concentration + (time step)(initial rate). Thus:

$$ \begin{aligned} [A]_{\mathrm{new}} &= 0.80+(2.0)(-0.0770)=0.646\ \mathrm{nM},\\ [AX]_{\mathrm{new}} &= 0.20+(2.0)(0.0515)=0.303\ \mathrm{nM},\\ [AY]_{\mathrm{new}} &= 0.10+(2.0)(0.0150)=0.130\ \mathrm{nM},\\ [AXY]_{\mathrm{new}} &= 0.050+(2.0)(0.0105)=0.0710\ \mathrm{nM}. \end{aligned} $$

The initial total antibody-state concentration is \(0.80+0.20+0.10+0.050=1.15\) nM. The updated sum is \(0.646+0.303+0.130+0.0710=1.15\) nM, consistent with the binding-only bookkeeping check. This is an approximation, not an exact solution over 2.0 min; a numerical solver recalculates rates as concentrations change.

Model Choice, Diagnostics, and Limitations

An explicit state model is useful when the distinction among free, singly bound, and doubly bound antibody matters to the question—for example, when a measurement distinguishes complexes or when a downstream process acts on a particular bound state. A lumped “bound antibody” variable can be simpler, but it combines \(AX\), \(AY\), and \(AXY\), obscuring which antigen is bound and which transitions are possible. An equilibrium occupancy description is another modeling choice, but it does not represent these time-dependent association and dissociation fluxes in the same explicit way.

Before interpreting model outputs, check that each transition has a source state, a destination state, a rate law with the right units, and a reverse transition if dissociation is part of the model. Useful warning signs include:

  • A concentration becomes negative during a simulation. This can indicate an integration step that is too large or an implementation problem.
  • The sum of antibody-state concentrations changes in a binding-only model. Recheck signs, missing transitions, and whether all state equations use the same flux definitions.
  • A second-antigen association term does not depend on the concentration of the free antigen. For a bimolecular mass-action association, it should.
  • The doubly bound state has only one formation route in a model that allows either antigen to bind first. The second route may have been omitted.

The model is deliberately limited. It does not account for target production or removal, antibody clearance, internalization, spatial separation, multiple copies of an antigen, or context-dependent binding rates. These processes should not be added implicitly by changing a binding constant; they require explicit model structure or justified parameter choices. A model should be no more detailed than the data and scientific question can support, but sufficiently detailed to preserve the states needed for the intended interpretation.

Key takeaway: For a two-antigen bispecific antibody, define the free, \(X\)-bound, \(Y\)-bound, and doubly bound states first. Then write each association and dissociation as a flux between specific states, and obtain each differential equation by adding incoming fluxes and subtracting outgoing fluxes.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. What is the association flux for \(AY+X\to AXY\), and what are its units?
  2. If \(AXY\) loses antigen \(X\), which antibody state is formed?
  3. Using \(k_{Y,\mathrm{on}}=0.015\) nM\(^{-1}\) min\(^{-1}\), \([AX]=0.12\) nM, and \([Y]=2.0\) nM, calculate the flux from \(AX\) to \(AXY\).
  4. Why does the \(AXY\) equation contain two formation terms?
  5. In the binding-only model, what should happen to the sum of the four antibody-state concentrations?

Solutions. (1) The flux is \(k_{X,\mathrm{on}}[AY][X]\), with units nM/min because nM\(^{-1}\) min\(^{-1}\) multiplied by two nM concentrations gives nM/min. (2) Dissociation of \(X\) from \(AXY\) leaves \(Y\) bound, forming \(AY\). (3) The flux is \((0.015)(0.12)(2.0)=0.0036\) nM/min. (4) \(AXY\) can form by binding \(Y\) to \(AX\) or by binding \(X\) to \(AY\). (5) It remains constant when binding and dissociation are the only processes, because every transition transfers antibody between the four states.

Clinical Biostats