Tissue Profiles After an Intravenous Bolus

From a Bolus to a Tissue Concentration Profile

After an intravenous bolus, blood concentration commonly declines while drug enters tissues. A tissue concentration therefore need not decline at the same rate as blood: it can continue to rise after blood concentration has begun falling. The perfusion-limited model makes this lag quantitative. Given a blood concentration input, it predicts how tissue concentration changes and when the direction of net exchange reverses.

Here the blood concentration is specified as a declining exponential after the bolus. This is an input to the tissue calculation, not a complete model of the blood pool. In a fully coupled whole-body model, blood concentration would itself respond to tissue exchange. That distinction matters when interpreting the calculation.

Definition: In a perfusion-limited tissue, blood and tissue are assumed to equilibrate rapidly at the tissue interface. If $C_b(t)$ is blood concentration, $C_T(t)$ is tissue concentration, and $K_{T:B}$ is the tissue-to-blood partition coefficient, the equilibrium relation is $C_T=K_{T:B}C_b$. With tissue blood flow $Q_T$ and tissue volume $V_T$, the perfusion-limited tissue equation is $$ \frac{dC_T}{dt}=\frac{Q_T}{V_T}\left(C_b(t)-\frac{C_T(t)}{K_{T:B}}\right). $$

The concentration difference in parentheses determines the direction of net exchange. If $C_b>C_T/K_{T:B}$, drug enters the tissue; if $C_b<C_T/K_{T:B}$, the tissue returns drug to blood. The partition coefficient is dimensionless when both concentrations use the same units. The ratio $Q_T/V_T$ has units of inverse time, as does the tissue response rate $ \lambda_T=Q_T/(V_TK_{T:B}) $.

Deriving the Profile for an Exponential Blood Input

Assume that the post-bolus blood concentration is $C_b(t)=C_{b,0}e^{-k_bt}$, where $C_{b,0}$ is the blood concentration at time zero and $k_b>0$ is its decline rate. Take the tissue to be drug-free at time zero, so $C_T(0)=0$. For compactness, write $a=Q_T/V_T$ and $\lambda_T=a/K_{T:B}$. The tissue equation becomes

\frac{dC_T}{dt}+\lambda_T C_T=aC_{b,0}e^{-k_bt}.

Multiplying by the integrating factor $e^{\lambda_Tt}$ gives $ \frac{d}{dt}(e^{\lambda_Tt}C_T)=aC_{b,0}e^{(\lambda_T-k_b)t} $. Integrating from zero to $t$, and using $C_T(0)=0$, yields the profile below when $k_b\ne\lambda_T$.

C_T(t)=\frac{aC_{b,0}}{k_b-\lambda_T} \left(e^{-\lambda_Tt}-e^{-k_bt}\right).

The concentration is zero at time zero, as required. If the two rates are equal, the denominator in this expression is zero, but the model solution remains finite. Direct integration in that case gives $C_T(t)=aC_{b,0}t e^{-\lambda_Tt}$. This is the continuous limiting form of the unequal-rate solution.

Key result: For a drug-free tissue exposed to an exponential post-bolus blood concentration, the tissue profile is the difference of two exponential terms. Its shape depends on both the blood decline rate $k_b$ and the tissue response rate $\lambda_T$, not on either rate alone.

Worked Example: Calculating Uptake and Redistribution

Consider a hypothetical tissue with $Q_T/V_T=0.200\ \mathrm{h}^{-1}$ and $K_{T:B}=2.00$. After a hypothetical intravenous bolus, let $C_b(t)=12.0e^{-0.250t}\ \mathrm{mg/L}$, with $t$ in hours, and let $C_T(0)=0$. Calculate tissue concentration at 2 and 10 hours and determine the direction of net exchange at each time.

The tissue response rate is $ \lambda_T=(0.200\ \mathrm{h}^{-1})/2.00=0.100\ \mathrm{h}^{-1} $. Since $k_b=0.250\ \mathrm{h}^{-1}\ne\lambda_T$, the profile is

C_T(t)=\frac{(0.200)(12.0)}{0.250-0.100} \left(e^{-0.100t}-e^{-0.250t}\right) =16.0\left(e^{-0.100t}-e^{-0.250t}\right)\ \mathrm{mg/L}.

At 2 hours, $ C_T(2)=16.0(e^{-0.200}-e^{-0.500})=16.0(0.818731-0.606531)=3.395\ \mathrm{mg/L} $. The blood concentration is $ C_b(2)=12.0e^{-0.500}=7.278\ \mathrm{mg/L} $, and the blood-equivalent tissue concentration is $ C_T(2)/K_{T:B}=3.395/2.00=1.698\ \mathrm{mg/L} $. Because $7.278>1.698$, the concentration gradient is into the tissue and net uptake continues.

At 10 hours, $ C_T(10)=16.0(e^{-1.000}-e^{-2.500})=16.0(0.367879-0.082085)=4.573\ \mathrm{mg/L} $. Here $ C_b(10)=12.0e^{-2.500}=0.9850\ \mathrm{mg/L} $, while $ C_T(10)/K_{T:B}=4.573/2.00=2.286\ \mathrm{mg/L} $. Now the tissue-equivalent concentration exceeds blood concentration, so net exchange is out of the tissue. Tissue concentration is still above its value at 2 hours, but it is declining by 10 hours: redistribution has begun.

When Does Tissue Concentration Peak?

For the unequal-rate profile, a peak occurs when its time derivative is zero. Equivalently, the original mass-balance equation shows that the derivative is zero when $C_b=C_T/K_{T:B}$. This is the instant at which net exchange changes direction; it does not mean that blood and tissue concentrations themselves are equal.

Differentiating the profile and setting the derivative to zero gives $k_b e^{-k_bt}=\lambda_Te^{-\lambda_Tt}$. Solving for time gives the peak time

t_{\mathrm{peak}}= \frac{\ln(k_b/\lambda_T)}{k_b-\lambda_T}, \qquad k_b\ne\lambda_T.

For positive $k_b$ and $\lambda_T$, this time is positive: the logarithm and denominator have the same sign. The equal-rate case has $C_T(t)=aC_{b,0}t e^{-\lambda_Tt}$; differentiation gives a peak at $t_{\mathrm{peak}}=1/\lambda_T$.

Worked Example: Finding the Peak in the Bolus Profile

Continue with the hypothetical parameters in the previous example: $k_b=0.250\ \mathrm{h}^{-1}$, $\lambda_T=0.100\ \mathrm{h}^{-1}$, and $C_{b,0}=12.0\ \mathrm{mg/L}$. The peak time is

t_{\mathrm{peak}} =\frac{\ln(0.250/0.100)}{0.250-0.100} =\frac{\ln(2.50)}{0.150} =6.109\ \mathrm{h}.

Evaluating the profile at that time gives $ C_T(6.109)=16.0(e^{-0.6109}-e^{-1.5273})\approx5.212\ \mathrm{mg/L} $. At the exact peak, $C_T/K_{T:B}=C_b$, so both concentrations on the blood scale are approximately $5.212/2.00=2.606\ \mathrm{mg/L}$. This agrees with $C_b(6.109)=12.0e^{-0.250(6.109)}\approx2.606\ \mathrm{mg/L}$, allowing for rounding.

Before the peak, the tissue-equivalent concentration is below blood concentration and net uptake raises $C_T$. After the peak, it is above blood concentration and net return lowers $C_T$. A peak in tissue concentration therefore occurs even though the blood input has been declining since time zero.

How Flow and Rate Matching Shape the Profile

Flow affects the speed of tissue response through $a=Q_T/V_T$, while partitioning affects the equilibrium ratio and the response rate through $\lambda_T=a/K_{T:B}$. Holding $K_{T:B}$ fixed, a larger $Q_T/V_T$ makes the tissue respond more rapidly to the same blood input. It can produce greater early tissue uptake, but it does not change the equilibrium ratio itself.

Worked Example: Comparing Two Hypothetical Tissues

Use the same hypothetical blood input $C_b(t)=12.0e^{-0.250t}\ \mathrm{mg/L}$ and the same partition coefficient $K_{T:B}=2.00$, but compare a tissue with $Q_T/V_T=0.400\ \mathrm{h}^{-1}$ to one with $Q_T/V_T=0.100\ \mathrm{h}^{-1}$. Both start drug-free. At 2 hours, the higher-flow tissue has $\lambda_T=0.200\ \mathrm{h}^{-1}$, so $ C_T(2)=[(0.400)(12.0)/(0.250-0.200)](e^{-0.400}-e^{-0.500})=96.0(0.670320-0.606531)=6.124\ \mathrm{mg/L} $.

For the lower-flow tissue, $\lambda_T=0.0500\ \mathrm{h}^{-1}$, giving $ C_T(2)=[(0.100)(12.0)/(0.250-0.0500)](e^{-0.100}-e^{-0.500})=6.00(0.904837-0.606531)=1.790\ \mathrm{mg/L} $. Under these specified parameters, the higher-flow tissue has accumulated more drug by 2 hours. This comparison concerns early profile values; it does not imply that changing flow changes the partition coefficient or that flow alone determines every feature of a tissue profile.

Worked Example: Equal Blood and Tissue Response Rates

For a separate hypothetical case, let $C_{b,0}=10.0\ \mathrm{mg/L}$, $k_b=0.100\ \mathrm{h}^{-1}$, $Q_T/V_T=0.200\ \mathrm{h}^{-1}$, and $K_{T:B}=2.00$. Then $\lambda_T=0.100\ \mathrm{h}^{-1}=k_b$, and the equal-rate solution is $C_T(t)=(0.200)(10.0)t e^{-0.100t}=2.00t e^{-0.100t}\ \mathrm{mg/L}$.

At 5 hours, $C_T(5)=10.0e^{-0.500}=6.065\ \mathrm{mg/L}$. Its peak occurs at $1/0.100=10.0\ \mathrm{h}$, where $C_T(10)=20.0e^{-1.000}=7.358\ \mathrm{mg/L}$. The blood concentration then is $10.0e^{-1.000}=3.679\ \mathrm{mg/L}$, and $C_T/K_{T:B}=7.358/2.00=3.679\ \mathrm{mg/L}$, as required at the peak.

Assumptions, Diagnostics, and Alternatives

These calculations assume constant tissue flow and volume, a constant partition coefficient, perfusion-limited exchange, a known monoexponential blood input, and a specified initial tissue concentration. The model treats concentrations within the tissue and its exchanging blood as adequately represented by the perfusion-limited description. It does not include a separate permeability barrier, time-varying binding, local metabolism, or a separately modeled blood input. The bolus profile above also assumes an initially drug-free tissue; with a nonzero initial concentration, the solution includes the additional term $C_T(0)e^{-\lambda_Tt}$.

Several checks help identify errors or inappropriate use. At $t=0$, the profile must equal the stated initial tissue concentration. For an initially drug-free tissue and positive rates, the tissue concentration must initially rise, since $dC_T/dt|_{t=0}=aC_{b,0}>0$. The computed derivative should be positive before the peak and negative afterward. At the peak, $C_T/K_{T:B}$ must equal $C_b$. Finally, concentrations must have units of concentration, while $k_b$, $Q_T/V_T$, and $\lambda_T$ must have units of inverse time.

A common mistake is to assume that falling blood concentration forces tissue concentration to fall immediately. The tissue equation shows otherwise: the sign of the tissue derivative depends on the difference $C_b-C_T/K_{T:B}$, not on whether $dC_b/dt$ is negative. Another error is to declare uptake based on comparing $C_T$ directly with $C_b$; the correct comparison is $C_T/K_{T:B}$ with $C_b$. These distinctions are especially important when $K_{T:B}\ne1$.

The analytic profile is convenient when the blood input is monoexponential and the assumptions hold. Numerical integration is an alternative when a specified blood input has a more complicated time course, although the same mass-balance equation and assumptions still need to be checked. An instantaneous-equilibrium approximation would instead set $C_T=K_{T:B}C_b$ at all times; it removes the finite response lag represented here and therefore cannot describe delayed tissue peaks. For a whole-body application, a coupled blood-and-tissue model may be needed rather than treating blood concentration as an imposed input. The present calculation cannot by itself establish that the assumed blood input or perfusion-limited representation is appropriate for a particular drug or tissue.

Key takeaway: A declining blood concentration can drive tissue concentration upward until $C_b=C_T/K_{T:B}$. After that crossing, net exchange reverses and tissue concentration falls. The lag and peak depend on the blood decline rate, tissue response rate, and partitioning.

Check Your Understanding

Answer each question using the model assumptions stated above. Solutions are included after the questions.

For an initially drug-free tissue with $k_b=0.300\ \mathrm{h}^{-1}$ and $\lambda_T=0.100\ \mathrm{h}^{-1}$, calculate the peak time.
At one time point, $C_b=3.0\ \mathrm{mg/L}$, $C_T=8.0\ \mathrm{mg/L}$, and $K_{T:B}=4.0$. Is net exchange into or out of the tissue?
When $k_b=\lambda_T$, what is the tissue profile for an initially drug-free tissue?
Why is a declining blood concentration alone insufficient to conclude that tissue concentration is declining?

Solutions. (1) $t_{\mathrm{peak}}=\ln(0.300/0.100)/(0.300-0.100)=\ln(3)/0.200=5.493\ \mathrm{h}$. (2) $C_T/K_{T:B}=8.0/4.0=2.0\ \mathrm{mg/L}$, which is less than $C_b=3.0\ \mathrm{mg/L}$, so net exchange is into the tissue. (3) The profile is $C_T(t)=aC_{b,0}t e^{-\lambda_Tt}$. (4) The tissue derivative is determined by $C_b-C_T/K_{T:B}$. Blood can decline while remaining above the blood-equivalent tissue concentration, in which case tissue concentration still rises.