Tutorials › Quantitative Systems Pharmacology (QSP) › Competing Reaction Pathways From a Shared Reactant

Reaction Kinetics and Mass Action · Tutorial 50 of 799

Competing Reaction Pathways From a Shared Reactant

Learn to calculate the rates and product fractions for two competing pathways that consume a shared reactant, including when the reactions are pseudo-first-order.

Beginner 12 min read

What You'll Learn

  • Write the mass-action rate laws for two pathways consuming the same reactant.
  • Derive the exponential depletion of the shared reactant and the corresponding product concentrations.
  • Distinguish instantaneous branching fractions from cumulative product fractions.
  • Convert bimolecular reactions to pseudo-first-order form when co-reactants remain effectively constant.
  • Check the assumptions behind constant branching fractions and identify signs they may fail.
  • Recover pathway-specific rate constants from depletion and product data.

Why a Shared Reactant Creates Branching

A reactant can be consumed by more than one reaction at the same time. For example, a compound may be converted into either product \(P_1\) or product \(P_2\). The amount of reactant available to each pathway is shared: as either reaction proceeds, it reduces the amount that remains for both. A useful first model asks how the competing rates determine the split of reactant between the products.

The key distinction is between the instantaneous rate split and the cumulative product split. Under constant first-order rate constants, these splits are equal and remain fixed over time. That convenient result depends on specific assumptions; it does not automatically apply when co-reactants change, products are removed, or pathways are reversible.

Definition: Consider two irreversible reactions, \(A\to P_1\) and \(A\to P_2\), each consuming one molecule of \(A\) and forming one molecule of its product. If both follow first-order mass-action kinetics, their rates are \(v_1=k_1[A]\) and \(v_2=k_2[A]\). The instantaneous branching fraction into pathway \(i\) is \(v_i/(v_1+v_2)\). The cumulative branching fraction is the amount of \(P_i\) formed divided by the total amount of \(A\) consumed over the stated time interval.

Here \([A]\), \([P_1]\), and \([P_2]\) denote concentrations, and \(k_1\) and \(k_2\) have units of inverse time. Assume a constant volume, constant rate constants, no input of \(A\), no reverse reactions, no other routes of consumption, and no loss or further reaction of either product. These conditions allow product formation to be counted directly from reactant consumption.

Deriving the First-Order Branching Equations

Both reactions remove \(A\), so their rates add in the balance for the shared reactant. Each reaction forms its own product at the rate at which it consumes \(A\). The resulting system is

$$ \frac{d[A]}{dt}=-(k_1+k_2)[A],\qquad \frac{d[P_1]}{dt}=k_1[A],\qquad \frac{d[P_2]}{dt}=k_2[A]. $$

Let \(K=k_1+k_2\), and let the initial concentrations be \([A](0)=[A]_0\), \([P_1](0)=[P_1]_0\), and \([P_2](0)=[P_2]_0\). Separating variables in the reactant equation gives

$$ \frac{d[A]}{[A]}=-K\,dt \quad\Longrightarrow\quad \ln\left(\frac{[A](t)}{[A]_0}\right)=-Kt \quad\Longrightarrow\quad [A](t)=[A]_0e^{-Kt}. $$

To derive product formation, substitute this reactant solution into each product equation and integrate from 0 to \(t\). For pathway 1,

$$ [P_1](t)-[P_1]_0 =\int_0^t k_1[A]_0e^{-Ks}\,ds =\frac{k_1}{K}[A]_0(1-e^{-Kt}). $$

The same calculation for pathway 2 replaces \(k_1\) with \(k_2\). Thus, the amount of \(A\) consumed by time \(t\) is \([A]_0(1-e^{-Kt})\), and the two product increments are fixed shares of that amount. Dividing the instantaneous rates gives the same shares:

$$ \frac{v_1}{v_1+v_2} =\frac{k_1[A]}{(k_1+k_2)[A]} =\frac{k_1}{k_1+k_2}, \qquad \frac{v_2}{v_1+v_2} =\frac{k_2}{k_1+k_2}. $$

The cancellation of \([A]\) is the central reason that the branching fractions do not change as \(A\) falls. The fractions are nonnegative and sum to one when \(k_1,k_2\geq 0\) and \(k_1+k_2>0\). If both rate constants are zero, no reaction occurs and branching fractions are undefined.

Worked Example: Product Formation Over Time

Suppose a hypothetical closed, constant-volume system starts with \([A]_0=100\ \mu\mathrm{M}\) and no products. The two first-order rate constants are \(k_1=0.12\ \mathrm{h}^{-1}\) and \(k_2=0.08\ \mathrm{h}^{-1}\). Find the reactant and product concentrations after 5 hours.

The total disappearance constant is \(K=0.12+0.08=0.20\ \mathrm{h}^{-1}\). The assumptions of constant first-order rates and no other routes mean the equations above apply. The remaining reactant is

$$ [A](5)=100e^{-(0.20)(5)} =100e^{-1} =36.79\ \mu\mathrm{M}. $$

Therefore, \(100-36.79=63.21\ \mu\mathrm{M}\) of \(A\) has been consumed. The branching fractions are \(0.12/0.20=0.600\) and \(0.08/0.20=0.400\). With no initial products,

$$ [P_1](5)=0.600(63.21)=37.93\ \mu\mathrm{M}, \qquad [P_2](5)=0.400(63.21)=25.28\ \mu\mathrm{M}. $$

The products sum to \(63.21\ \mu\mathrm{M}\), matching the reactant consumed. At the initial concentration, the two rates are \(v_1=0.12(100)=12\ \mu\mathrm{M}\,\mathrm{h}^{-1}\) and \(v_2=0.08(100)=8\ \mu\mathrm{M}\,\mathrm{h}^{-1}\); their ratio also gives the 60:40 split.

Pseudo-First-Order Competing Pathways

A pathway may be bimolecular rather than first-order in \(A\). Suppose the reactions are \(A+B\to P_1\) and \(A+C\to P_2\), with mass-action rates \(v_1=k_1[A][B]\) and \(v_2=k_2[A][C]\). If \([B]\) and \([C]\) are held constant, or are in such excess that their changes are negligible over the interval of interest, define \(k'_1=k_1[B]\) and \(k'_2=k_2[C]\). Then \(v_1=k'_1[A]\) and \(v_2=k'_2[A]\), which has the same mathematical form as the first-order model.

For this case, \(k_1\) and \(k_2\) have units of concentration\(^{-1}\) time\(^{-1}\), while the pseudo-first-order constants \(k'_1\) and \(k'_2\) have units of inverse time. The total disappearance constant is \(K'=k'_1+k'_2\). The reactant decays exponentially with \(K'\), and the product fractions are \(k'_1/K'\) and \(k'_2/K'\), provided the co-reactant concentrations stay effectively constant.

Worked Example: Branching With Constant Co-Reactants

In a hypothetical system, \(A+B\to P_1\) has \(k_1=0.020\ \mathrm{mM}^{-1}\mathrm{min}^{-1}\) while \([B]=10\ \mathrm{mM}\). The competing reaction \(A+C\to P_2\) has \(k_2=0.0050\ \mathrm{mM}^{-1}\mathrm{min}^{-1}\) while \([C]=20\ \mathrm{mM}\). Assume \(B\) and \(C\) remain effectively constant, \([A]_0=5.00\ \mathrm{mM}\), and both products initially equal zero. Calculate the initial rates and concentrations after 3.00 minutes.

The pseudo-first-order constants are \(k'_1=(0.020)(10)=0.200\ \mathrm{min}^{-1}\) and \(k'_2=(0.0050)(20)=0.100\ \mathrm{min}^{-1}\). Thus, \(K'=0.300\ \mathrm{min}^{-1}\), and the branch fractions are \(0.200/0.300=0.667\) and \(0.100/0.300=0.333\). The initial rates are \(v_1=(0.200)(5.00)=1.00\ \mathrm{mM}\,\mathrm{min}^{-1}\) and \(v_2=(0.100)(5.00)=0.500\ \mathrm{mM}\,\mathrm{min}^{-1}\).

At 3.00 minutes, the remaining reactant is

$$ [A](3.00)=5.00e^{-(0.300)(3.00)} =5.00e^{-0.900} =2.03\ \mathrm{mM}. $$

The consumed amount is \(5.00-2.03=2.97\ \mathrm{mM}\), so the product concentrations are approximately \([P_1]=0.667(2.97)=1.98\ \mathrm{mM}\) and \([P_2]=0.333(2.97)=0.989\ \mathrm{mM}\). Their sum, \(2.97\ \mathrm{mM}\), agrees with the amount of \(A\) consumed, subject to rounding.

When the Branching Fractions Can Change

If the co-reactants are not effectively constant, the pseudo-first-order constants vary with time. The instantaneous fraction into pathway 1 is then

$$ \frac{v_1(t)}{v_1(t)+v_2(t)} =\frac{k_1[B](t)}{k_1[B](t)+k_2[C](t)}, $$

assuming \([A](t)>0\). This expression can change as \(B\) or \(C\) is consumed or supplied. The cumulative fraction is not generally the instantaneous fraction at the endpoint; it depends on the integrated amounts through each pathway. In that situation, use the full coupled mass-action equations rather than assigning a fixed product fraction.

A useful check under the simple one-to-one model is the balance \([A]+[P_1]+[P_2]=[A]_0+[P_1]_0+[P_2]_0\). It follows by adding the three differential equations, whose sum is zero. If data substantially violate this balance, possible explanations include unmeasured pathways, product degradation, measurement error, or a mistaken stoichiometric model. Likewise, curvature in a plot of \(\ln([A](t)/[A]_0)\) against time is inconsistent with a single constant total first-order rate over that interval, though noisy measurements can also produce apparent curvature.

The mass-action model is mechanistic: it retains both pathways and predicts each product. A simpler lumped description may fit only the total disappearance of \(A\), using \(K=k_1+k_2\); that alone cannot identify the individual pathway constants or predict the product split. Conversely, if either pathway is reversible, saturable, or has changing co-reactants, the constant first-order branch formula is generally not adequate. Those cases require rate laws that represent the relevant mechanism rather than treating the fractions as fixed.

Worked Example: Inferring Pathway Rates From Depletion and Products

In a hypothetical experiment, \(100\ \mu\mathrm{M}\) of \(A\) is followed for \(6.931\) hours. At that time, \(50.0\ \mu\mathrm{M}\) of \(A\) remains, and the measured product increments are \(30.0\ \mu\mathrm{M}\) of \(P_1\) and \(20.0\ \mu\mathrm{M}\) of \(P_2\). Assume the products are stable and the two first-order pathways are the only routes. Estimate \(k_1\) and \(k_2\).

The total disappearance constant follows from \([A](t)=[A]_0e^{-Kt}\):

$$ K=-\frac{1}{t}\ln\left(\frac{[A](t)}{[A]_0}\right) =-\frac{1}{6.931}\ln\left(\frac{50.0}{100}\right) =0.1000\ \mathrm{h}^{-1}. $$

The consumed amount is \(100-50.0=50.0\ \mu\mathrm{M}\). The observed product fractions are \(30.0/50.0=0.600\) and \(20.0/50.0=0.400\); their sum is one, as required by the assumed balance. Since these equal \(k_1/K\) and \(k_2/K\), respectively,

$$ k_1=(0.600)(0.1000)=0.0600\ \mathrm{h}^{-1}, \qquad k_2=(0.400)(0.1000)=0.0400\ \mathrm{h}^{-1}. $$

The inferred constants add to \(0.1000\ \mathrm{h}^{-1}\), matching the depletion estimate. This agreement is a consistency check, not proof that the assumed mechanism is uniquely correct; other mechanisms can sometimes produce similar observations.

Common Pitfall and Practical Checks

A common error is to treat \(k_1\) and \(k_2\) themselves as product fractions. Rate constants have units of inverse time in the first-order model, whereas fractions are dimensionless. The correct comparison is \(k_1/(k_1+k_2)\), not \(k_1\) alone. For pseudo-first-order pathways, compare \(k_1[B]\) and \(k_2[C]\), not the underlying bimolecular constants without their co-reactant concentrations.

Before using a fixed branching fraction, check that the pathways consume the same shared reactant, follow the stated rate laws, and have constant rate coefficients over the analyzed interval. Confirm that the product stoichiometry supports the assumed one-to-one concentration balance. If measured product increments do not sum to reactant consumed, do not force them into fractions that sum to one without accounting for missing routes, losses, or stoichiometric differences.

1
Write each pathway rate.
Use mass action for each reaction, including any co-reactant concentrations that affect the rate.
2
Check whether the rates reduce to first order in the shared reactant.
For bimolecular reactions, this reduction requires co-reactants to remain effectively constant during the interval.
3
Add the consumption rates and calculate product formation.
The total disappearance rate determines the reactant trajectory; each pathway rate determines its product trajectory.
4
Check units and balances.
Confirm that rate constants and rates have compatible units and that predicted product formation is consistent with reactant consumption under the assumed stoichiometry.

Check Your Understanding

Answer each question using the stated assumptions. Brief solutions follow the questions.

  1. For first-order pathways with \(k_1=0.03\ \mathrm{h}^{-1}\) and \(k_2=0.07\ \mathrm{h}^{-1}\), what fraction of consumed \(A\) forms \(P_1\)?
  2. If \([A]_0=40\ \mu\mathrm{M}\) and \(K=0.20\ \mathrm{h}^{-1}\), what remains after \(5.00\) hours?
  3. For \(A+B\to P_1\) and \(A+C\to P_2\), what condition permits constant pseudo-first-order branching fractions?
  4. Why can total depletion data alone not determine the individual pathway constants in the two-pathway first-order model?

Solutions. (1) \(k_1/(k_1+k_2)=0.03/0.10=0.300\), or 30.0%. (2) \([A](5)=40e^{-(0.20)(5)}=40e^{-1}=14.7\ \mu\mathrm{M}\). (3) The concentrations of \(B\) and \(C\) must be constant or change negligibly, so \(k_1[B]\) and \(k_2[C]\) remain constant. (4) Depletion identifies only the sum \(K=k_1+k_2\); product-specific measurements or other information are needed to separate the two constants.

Clinical Biostats