Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Converting Ke0 to an Equilibration Half-Life

Effect-Compartment Models and Hysteresis · Tutorial 87 of 598

Converting Ke0 to an Equilibration Half-Life

Learn to convert \(k_{e0}\) into an interpretable equilibration timescale and avoid confusing it with plasma elimination half-life.

Advanced 9 min read

What You'll Learn

  • Derive the equilibration half-life from first-order effect-site dynamics.
  • Convert between \(k_{e0}\), equilibration half-life, and time constant.
  • Calculate the time to reach a specified fraction of a steady effect-site level.
  • Check rate-constant units and convert consistently between minutes and hours.
  • Distinguish effect-site equilibration half-life from plasma elimination half-life.
  • Recognize when a single equilibration half-life is not an adequate description.

Why Convert a Rate Constant to a Half-Life?

In an effect-compartment model, \(k_{e0}\) describes how quickly the effect-site concentration responds to plasma concentration. Its units are inverse time, such as \(\mathrm{h}^{-1}\) or \(\mathrm{min}^{-1}\). A half-life expresses the same rate in units of time, which can make the implied delay easier to interpret and communicate. The conversion is exact for the first-order equilibration model, not a separate biological process.

The crucial distinction is that the equilibration half-life describes how quickly the gap between plasma and effect-site concentrations closes when plasma concentration is held fixed. It is not, in general, the time for plasma drug concentration to halve. Plasma elimination and effect-site equilibration are different processes, each with its own rate parameter.

Definition: In the first-order effect-compartment model, \(C_p(t)\) is plasma concentration, \(C_e(t)\) is effect-site concentration, and \(k_{e0}>0\) is the effect-site equilibration rate constant. The model is $$ \frac{dC_e(t)}{dt}=k_{e0}\bigl(C_p(t)-C_e(t)\bigr). $$ For a fixed plasma concentration, the equilibration half-life, denoted \(t_{1/2,e}\), is the time required for the difference between the current effect-site concentration and its eventual equilibrium value to become half as large. Thus \(t_{1/2,e}=\ln(2)/k_{e0}\).

Deriving the Conversion

Suppose plasma concentration changes to a constant level \(C^\ast\) and remains there. Under the model, the effect-site concentration approaches \(C^\ast\). Define the remaining gap as \(D(t)=C^\ast-C_e(t)\). Differentiating this gap and substituting the model equation gives

$$ \frac{dD(t)}{dt} =-\frac{dC_e(t)}{dt} =-k_{e0}\bigl(C^\ast-C_e(t)\bigr) =-k_{e0}D(t). $$

This differential equation has the solution \(D(t)=D(0)e^{-k_{e0}t}\). The fraction of the original gap remaining after time \(t\) is therefore \(D(t)/D(0)=e^{-k_{e0}t}\). At the equilibration half-life, this fraction is \(1/2\). Solving that condition gives the conversion:

$$ e^{-k_{e0}t_{1/2,e}}=\frac{1}{2} \quad\Longrightarrow\quad -k_{e0}t_{1/2,e}=\ln\left(\frac{1}{2}\right) \quad\Longrightarrow\quad t_{1/2,e}=\frac{\ln(2)}{k_{e0}}. $$

Conversely, if the equilibration half-life is given, the rate constant is \(k_{e0}=\ln(2)/t_{1/2,e}\). Here \(\ln(2)\) is approximately 0.6931. The units follow directly: dividing a dimensionless number by a rate in \(\mathrm{h}^{-1}\) yields hours. Always express the time and rate in compatible units before calculating.

Equivalent parameterizations: The first-order equilibration process can be parameterized by \(k_{e0}\), by \(t_{1/2,e}=\ln(2)/k_{e0}\), or by the time constant \(\tau_e=1/k_{e0}\). The conversions are \(t_{1/2,e}=(\ln 2)\tau_e\) and \(\tau_e=t_{1/2,e}/\ln 2\). These are alternative descriptions of the same model rate, not independent parameters.

The time constant \(\tau_e\) is the time for the remaining gap to fall to \(e^{-1}\), or about 36.8%, of its initial value. Equivalently, about 63.2% of the change toward the fixed equilibrium level has occurred after one time constant. A half-life instead describes the time to close half the gap. Both are useful, but they are not numerically equal.

Worked Example: Converting an Hourly Rate to a Half-Life

Suppose a hypothetical model estimates \(k_{e0}=0.230\ \mathrm{h}^{-1}\). Calculate the equilibration half-life and time constant.

Half-life. Use \(t_{1/2,e}=\ln(2)/k_{e0}\), with the rate in inverse hours:

$$ t_{1/2,e} =\frac{0.6931}{0.230\ \mathrm{h}^{-1}} =3.014\ \mathrm{h}. $$

Thus, under the model's fixed-plasma condition, the remaining gap between effect-site concentration and its equilibrium value is halved in about 3.01 hours.

Time constant. Since \(\tau_e=1/k_{e0}\),

$$ \tau_e =\frac{1}{0.230\ \mathrm{h}^{-1}} =4.348\ \mathrm{h}. $$

The time constant is longer than the half-life, as expected: \(0.6931\times4.348\ \mathrm{h}=3.014\ \mathrm{h}\). Neither value describes how quickly plasma concentration is eliminated.

Other Useful Equilibration Times

A half-life does not mean that equilibrium is reached after one half-life. It means that half the original difference remains. For the fixed-plasma step described above, the fraction of the change completed by time \(t\) is \(1-e^{-k_{e0}t}\). More generally, the time to close a fraction \(q\) of the initial gap, where \(0<q<1\), is

$$ t_q=\frac{-\ln(1-q)}{k_{e0}}. $$

For example, closing 90% of the gap requires \(t_{90}=\ln(10)/k_{e0}\), which is about 3.322 equilibration half-lives. Such a time is a model-based measure of approach to equilibrium; it does not imply that the effect-site concentration becomes exactly equal to plasma concentration at a finite time.

Worked Example: Converting a Half-Life in Minutes to an Hourly Rate

In a hypothetical model, the reported equilibration half-life is 18.0 minutes. Find \(k_{e0}\) in both \(\mathrm{min}^{-1}\) and \(\mathrm{h}^{-1}\), and calculate the time constant.

First use minutes consistently:

$$ k_{e0} =\frac{0.6931}{18.0\ \mathrm{min}} =0.03851\ \mathrm{min}^{-1}. $$

There are 60 minutes in an hour, so convert the rate by multiplying by 60:

$$ 0.03851\ \mathrm{min}^{-1}\times60\ \frac{\mathrm{min}}{\mathrm{h}} =2.311\ \mathrm{h}^{-1}. $$

The time constant is \(1/k_{e0}\). Calculating in minutes avoids a unit mismatch:

$$ \tau_e =\frac{1}{0.03851\ \mathrm{min}^{-1}} =25.97\ \mathrm{min}. $$

As a check, \(0.6931\times25.97\ \mathrm{min}=18.00\ \mathrm{min}\). A rate expressed per hour is numerically larger than the rate per minute because one hour contains 60 minutes; the underlying equilibration speed has not changed.

Equilibration Half-Life Is Not Plasma Elimination Half-Life

For a one-compartment plasma concentration declining by first-order elimination after input stops, the elimination model is \(dC_p/dt=-k_{el}C_p\). Its plasma elimination half-life is \(t_{1/2,el}=\ln(2)/k_{el}\). Although this formula has the same mathematical form as the effect-site conversion, \(k_{el}\) and \(k_{e0}\) refer to different processes. The plasma half-life concerns the decline of plasma concentration; the equilibration half-life concerns the closing of the plasma-to-effect-site concentration gap under fixed plasma concentration.

There is no general rule that one of these half-lives must be longer than the other. They may be similar by coincidence, but that does not make them interchangeable. In more complex pharmacokinetic models, plasma decline can also have multiple phases, so a stated plasma half-life may be specific to a particular model or terminal phase. The effect-site equilibration parameter still describes the assumed effect-compartment process.

Worked Example: Comparing Two Distinct Half-Lives

Consider a hypothetical model with \(k_{e0}=0.120\ \mathrm{min}^{-1}\). Separately, suppose the plasma elimination half-life is 4.00 hours under a first-order plasma elimination model. Calculate both timescales and show why they must not be confused.

The effect-site equilibration half-life is

$$ t_{1/2,e} =\frac{0.6931}{0.120\ \mathrm{min}^{-1}} =5.776\ \mathrm{min}. $$

The stated plasma elimination half-life is 4.00 hours, or \(4.00\times60=240\) minutes. Its corresponding elimination rate is

$$ k_{el} =\frac{0.6931}{240\ \mathrm{min}} =0.002888\ \mathrm{min}^{-1} =0.1733\ \mathrm{h}^{-1}. $$

The effect-site rate, \(0.120\ \mathrm{min}^{-1}\), and plasma elimination rate, \(0.002888\ \mathrm{min}^{-1}\), differ substantially. The model therefore implies an effect-site equilibration half-life of about 5.78 minutes and a plasma elimination half-life of 240 minutes. Substituting the plasma half-life into the \(k_{e0}\) conversion would answer a different question.

Assumptions, Checks, and Limitations

The conversion is valid when the effect-site dynamics are represented by a constant, positive first-order rate \(k_{e0}\). For its half-life interpretation, plasma concentration is held at a fixed level so that there is a well-defined equilibrium value toward which \(C_e\) approaches. Concentration is assumed to be represented consistently in the effect-compartment equation, and time units must match the units used for \(k_{e0}\). If plasma concentration varies over time, the same parameter remains a characteristic rate of the specified model, but the effect-site trajectory need not have a single observed “halfway time” that equals \(t_{1/2,e}\).

A practical check is to convert any reported rate to an intuitive time before interpreting it. A value reported in \(\mathrm{min}^{-1}\) may imply a timescale of only minutes, whereas the same numerical value in \(\mathrm{h}^{-1}\) implies a timescale of hours. Inspect the original model specification and units rather than relying on a parameter label alone. Also check that \(k_{e0}>0\): zero would imply no movement toward equilibrium, and a negative rate would not describe stable first-order equilibration.

The three parameterizations are equivalent only within the same first-order model. A pure fixed delay, for instance, means no response-site change until a specified time has passed; first-order equilibration instead begins changing immediately and approaches its target gradually. More elaborate delay or turnover structures may have different timescales and cannot generally be reduced to one \(k_{e0}\) without additional assumptions.

Common pitfall: “The plasma half-life is 4 hours, so \(k_{e0}=0.693/4\ \mathrm{h}^{-1}\).” This incorrectly treats elimination and equilibration as the same process. Use \(k_{e0}=\ln(2)/t_{1/2,e}\) only when the given time is the effect-site equilibration half-life; use the plasma elimination half-life only to calculate the corresponding plasma elimination rate under its stated model.

When an inferred or reported half-life seems implausible, first check units and conversion direction: \(k_{e0}=\ln(2)/t_{1/2,e}\), not \(t_{1/2,e}/\ln(2)\). Then check whether the quoted half-life refers to the effect site, the plasma, or another process. Even with correct arithmetic, the simple conversion cannot validate whether a first-order effect-compartment model adequately represents the data; that requires assessing the model itself.

Check Your Understanding

Use the first-order equilibration model and keep units explicit. Solutions follow the questions.

  1. A hypothetical model has \(k_{e0}=0.350\ \mathrm{h}^{-1}\). What is its equilibration half-life in hours?
  2. A reported equilibration half-life is 30.0 minutes. What is \(k_{e0}\) in \(\mathrm{h}^{-1}\)?
  3. After one equilibration half-life under fixed plasma concentration, what fraction of the original gap remains, and what fraction of the change toward equilibrium has occurred?
  4. Why can a 6-hour plasma elimination half-life not, by itself, determine \(k_{e0}\)?

Solutions. (1) \(t_{1/2,e}=0.6931/(0.350\ \mathrm{h}^{-1})=1.980\ \mathrm{h}\). (2) Convert 30.0 minutes to 0.500 hours, then \(k_{e0}=0.6931/(0.500\ \mathrm{h})=1.386\ \mathrm{h}^{-1}\). (3) Half the original gap remains, so half the change toward equilibrium has occurred. (4) The plasma elimination half-life describes plasma decline and determines \(k_{el}\) under its model; \(k_{e0}\) describes effect-site equilibration. They are distinct parameters, so the elimination half-life does not specify \(k_{e0}\).

Clinical Biostats