Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Linear Concentration-Effect Models and Their Slope

Direct-Effect Models: Linear, Emax and Sigmoid Emax · Tutorial 15 of 598

Linear Concentration-Effect Models and Their Slope

Learn to derive and interpret a linear concentration-effect relationship, convert its slope across concentration scales, and calculate predicted responses within the model’s appropriate range.

Intermediate 9 min read

What You'll Learn

  • Derive the affine concentration-effect equation from a constant response change per concentration increment.
  • Interpret the intercept and slope, including the slope’s direction and units.
  • Convert a slope correctly when concentration units change.
  • Calculate expected responses for increasing and decreasing effects.
  • Identify diagnostics, assumptions, and limitations of a linear relationship.
  • Distinguish a direct concentration-effect model from cumulative exposure summaries and nonlinear alternatives.

Why Model a Linear Concentration-Effect Relationship?

A simple concentration-effect model can summarize how a pharmacodynamic response changes as drug concentration changes. Its central assumption is that, over the concentration range being considered, equal increments in concentration correspond to equal changes in the modeled response. This can be a useful approximation when the relationship is roughly straight over the observed range, even if it would not remain straight across all possible concentrations.

The model concerns a specified concentration and a specified effect measure. It does not say that dose, concentration, area under the concentration-time curve, and other exposure metrics are interchangeable. The concentration must be the metric relevant to the effect measurement, and the time or sampling window for each must be appropriate to the scientific question.

Definition: A linear concentration-effect model for the expected response at concentration \(C\) is
$$ E(C)=E_0+sC. $$
Here, \(E(C)\) is the model-predicted response at concentration \(C\), \(E_0\) is the predicted response at \(C=0\), and \(s\) is the slope: the change in predicted response per unit increase in concentration. The term “linear” is commonly used for this relationship, although mathematically it is an affine function when \(E_0\) is nonzero.

The units are part of the model’s interpretation. If response is measured in beats per minute and concentration in ng/mL, then \(s\) has units of beats per minute per ng/mL. The intercept \(E_0\) has the same units as the response. A positive slope means the predicted response increases with concentration; a negative slope means it decreases.

Why the Equation Has This Form

A constant slope means that a concentration increment \(\Delta C\) produces the same predicted response increment wherever that increment begins. In symbols, the response difference is proportional to the concentration difference:

$$ E(C+\Delta C)-E(C)=s\Delta C. $$

Set the starting concentration to zero and the increment to \(C\). The predicted response difference from zero concentration is then \(sC\), so the response at \(C\) is the response at zero concentration plus that difference:

$$ E(C)-E(0)=sC,\qquad E(C)=E(0)+sC. $$

Writing \(E(0)=E_0\) gives \(E(C)=E_0+sC\). Equivalently, for any two concentrations \(C_1\) and \(C_2\), the model implies

$$ E(C_2)-E(C_1)=s(C_2-C_1). $$

This identity is useful because the slope describes a predicted difference between concentrations, not merely a property of a plotted line. For example, if \(s\) is positive, raising concentration by \(5\) concentration units changes the predicted response by \(5s\) response units. The direction follows from the sign of \(s\); the magnitude depends on both \(s\) and the concentration increment.

The equation describes a model-predicted response, often understood as a conditional mean response at \(C\). Individual measurements may differ from that prediction because of measurement error, biological variation, or other sources of variation. The equation alone does not specify a residual error distribution or provide an uncertainty interval.

Assumptions and Interpretation

Use the model as a local or range-specific description, not as a universal law of pharmacology. Its main assumptions are:

  • Approximate straight-line behavior over the range used. The expected response changes by roughly the same amount for each equal concentration increment.
  • A meaningful direct concentration-effect pairing. The concentration metric and effect measurement refer to a scientifically appropriate time or window. A delay between concentration and effect, for example, can make a same-time pairing misleading.
  • A consistent response definition and concentration scale. The response units, concentration units, and measurement procedure do not change across the observations being described.
  • Other relevant conditions are comparable or accounted for. The model does not automatically separate a concentration relationship from changes in other factors that also affect response.
  • Predictions remain in a plausible domain. The unbounded line should not be assumed to respect biological limits unless its predictions are checked against the possible range of the response.

The intercept is a mathematical feature of the specified relationship: it is the prediction at \(C=0\). It is a baseline response only when that zero-concentration condition corresponds to the baseline of interest and the model is appropriate there. If all observed concentrations are far above zero, \(E_0\) may be an extrapolated intercept rather than a directly observed baseline.

Worked Example: Increasing Response With Concentration

Suppose a hypothetical direct-effect model describes heart rate in beats per minute using concentration in ng/mL. Let \(E_0=68\) beats per minute and \(s=0.42\) beats per minute per ng/mL. Calculate the predicted heart rate at \(C=15\) ng/mL and compare it with the prediction at \(C=35\) ng/mL.

The model is \(E(C)=E_0+sC\). At \(15\) ng/mL:

$$ E(15)=68\ \text{beats/min} +\left(0.42\ \frac{\text{beats/min}}{\text{ng/mL}}\right)(15\ \text{ng/mL}) =68+6.3=74.3\ \text{beats/min}. $$

At \(35\) ng/mL:

$$ E(35)=68+(0.42)(35)=68+14.7=82.7\ \text{beats/min}. $$

The concentration increase is \(35-15=20\) ng/mL, so the model predicts a response increase of \(0.42\times20=8.4\) beats per minute. This agrees with \(82.7-74.3=8.4\) beats per minute. The calculation is valid as a model prediction; whether the straight-line assumption is reasonable at both concentrations must be assessed using relevant data and scientific knowledge.

Changing Concentration Units Changes the Slope

A slope is a response change per concentration unit. Changing the concentration unit changes how many units represent the same physical concentration, so the numerical slope must change inversely. The predicted response itself does not change when the conversion is done correctly.

Suppose the original concentration is \(C\), and a new numerical concentration \(C^*\) is related by \(C^*=kC\). Substitution gives

$$ E(C)=E_0+sC =E_0+\frac{s}{k}C^*. $$

Thus the slope on the new scale is \(s^*=s/k\), with units of response per new concentration unit. This is a unit conversion, not a change in the underlying concentration-effect relationship.

Worked Example: Converting the Slope to a Different Concentration Scale

A hypothetical response model is \(E(C)=5+2.4C\), where \(E\) is measured in response units and \(C\) is measured in mg/L. Convert the model to use concentration in µg/L, then calculate the prediction at \(0.75\) mg/L.

Since \(1\) mg \(=1000\) µg, the numerical concentration in µg/L is \(C^*=1000C\). Therefore \(C=C^*/1000\), and

$$ E(C^*)=5+2.4\left(\frac{C^*}{1000}\right) =5+0.0024C^*. $$

The slope is now \(0.0024\) response units per µg/L, rather than \(2.4\) response units per mg/L. At \(0.75\) mg/L, the concentration is \(750\) µg/L. Either scale gives the same prediction:

$$ E=5+(2.4)(0.75)=5+1.8=6.8, \qquad E=5+(0.0024)(750)=5+1.8=6.8. $$

The equality checks the conversion. Keeping the slope at \(2.4\) after changing concentration to µg/L would incorrectly multiply the concentration contribution by \(1000\).

Negative Slopes and Decreasing Effects

A negative slope represents a decrease in the modeled response as concentration increases. The response could be a measurement such as enzyme activity, a biomarker, or a symptom score; the sign describes the direction of the chosen response measure, not whether the drug is beneficial or harmful. That judgment requires the scientific meaning of the response.

Worked Example: Predicting a Decreasing Response

Suppose a hypothetical model predicts percentage activity \(E\) from concentration \(C\) in mg/L: \(E(C)=100-1.8C\), where the slope has units of percentage points per mg/L. Calculate the prediction at \(C=12\) mg/L and the predicted change from \(C=4\) to \(C=12\) mg/L.

At \(12\) mg/L:

$$ E(12)=100\ \%-\left(1.8\ \frac{\text{percentage points}}{\text{mg/L}}\right)(12\ \text{mg/L}) =100\%-21.6\ \text{percentage points} =78.4\%. $$

At \(4\) mg/L, the prediction is \(E(4)=100-(1.8)(4)=92.8\%\). The predicted difference is

$$ E(12)-E(4)=78.4\%-92.8\%=-14.4\ \text{percentage points}. $$

The same result follows from the slope and concentration difference: \((-1.8)(12-4)=-14.4\) percentage points. The model predicts a decrease of \(14.4\) percentage points. This is not a relative decrease of \(14.4\%\); it is an absolute difference on the percentage scale. The predicted value at \(12\) mg/L is plausible for a percentage activity, but predictions at higher concentrations should also be checked for implausible values.

Diagnostics, Alternatives, and Limitations

A fitted straight line can look convincing over a narrow range even when it is not a good description more broadly. Diagnostics should therefore focus on whether the assumptions are reasonable in the concentration range used and whether predictions serve the intended purpose.

  • Look for systematic curvature. If residuals (observed responses minus model predictions) tend to be positive at low and high concentrations but negative in the middle, or show the reverse pattern, a straight line may be missing curvature.
  • Check for changes in spread. A much wider range of residuals at higher concentrations suggests that a constant-variance error model, if one is being used, may not describe the data adequately.
  • Inspect the concentration range and influential observations. A small number of observations at the extremes can strongly affect a fitted slope. Predictions beyond the observed range are extrapolations and need additional justification.
  • Check timing and repeated observations. If the response lags behind concentration, or if multiple measurements from one individual are treated as independent, a simple direct line may not represent the data structure.
  • Check whether predictions are possible. A line can predict negative values for a nonnegative response or values above a natural maximum. Such predictions indicate that the linear approximation is not suitable over that range.

A baseline-relative response can be written as a change from baseline rather than an absolute response. If baseline is \(E_0\), then the corresponding change is \(\Delta E(C)=E(C)-E_0=sC\). This re-expression clarifies the predicted change, but it does not alter the slope or remove the assumptions about linearity and concentration-effect pairing.

When the observed relationship bends or approaches a limit, a nonlinear concentration-effect model may be more appropriate. A saturating model, for instance, can represent responses that approach a maximum rather than increasing or decreasing indefinitely. A more complex model is not automatically better: the data must contain enough information to support its additional parameters, and its assumptions must match the question.

Common pitfall: “The concentration is now in µg/L, so I can keep the same numerical slope and just substitute the converted concentration.” This is wrong because the slope is expressed per concentration unit. If \(C^*=1000C\), then the new slope is \(s/1000\); changing the concentration scale without changing the slope changes the prediction, not merely its notation.

A second interpretive caution is to avoid treating the slope as a causal effect just because it is estimated from an observed association. The model specifies a predicted response change across concentrations. A causal interpretation requires a design and assumptions that support it; the equation alone does not establish those conditions.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. For \(E(C)=30-0.6C\), with \(C\) in mg/L and response in score units, calculate the prediction at \(C=10\) mg/L and interpret the slope.
  2. A model predicts an increase of \(12\) response units when concentration rises by \(4\) ng/mL. What is the slope, including units?
  3. A model has slope \(0.08\) response units per mg/L. What is the equivalent slope per µg/L?
  4. Why might \(E_0\) not be an observed baseline response when all measured concentrations are far above zero?

Solutions. (1) \(E(10)=30-(0.6)(10)=24\) score units. The slope is \(-0.6\) score units per mg/L, so each 1 mg/L increase corresponds to a predicted decrease of \(0.6\) score units within the model’s range. (2) \(s=12/4=3\) response units per ng/mL. (3) Since \(1\) mg/L \(=1000\) µg/L, divide the numerical slope by \(1000\): \(0.08/1000=0.000080\) response units per µg/L. (4) \(E_0\) is the model’s prediction at \(C=0\). If zero concentration is far outside the observed range, that value is an extrapolation and need not equal a measured baseline.

Clinical Biostats