3. Pharmacodynamics

Pharmacodynamics (PD) describes what the drug does to the body -- the relationship between drug concentration at the site of action and the resulting pharmacological effect. It encompasses dose-response relationships, therapeutic indices, selectivity, and receptor-mediated signaling.

Historical Context

The concept of dose-response was first quantified by A.J. Clark (1926) using acetylcholine on frog heart preparations. J.H. Gaddum (1937) introduced the concept of competitive antagonism, while Arunlakshana and Schild (1959) developed the graphical method for quantifying antagonist potency that bears Schild's name.

The integration of PK and PD into unified models began with Sheiner and colleagues in the 1970s. The effect compartment concept (Sheiner 1981) resolved the apparent disconnect between plasma concentration and effect timing. Modern PK/PD modeling uses population approaches implemented in software like NONMEM, Monolix, and nlmixr.

Derivation 1: The Emax Model

The simplest pharmacodynamic model assumes a hyperbolic relationship between drug concentration and effect, directly analogous to Michaelis-Menten enzyme kinetics:

$$E = \frac{E_{\max} \cdot C}{EC_{50} + C}$$

Where $E_{\max}$ is the maximum achievable effect and $EC_{50}$is the concentration producing 50% of maximum effect. The sigmoid Emax model introduces the Hill coefficient $\gamma$:

$$E = \frac{E_{\max} \cdot C^{\gamma}}{EC_{50}^{\gamma} + C^{\gamma}}$$

The slope of the concentration-effect curve at $EC_{50}$ is:

$$\left.\frac{dE}{dC}\right|_{C=EC_{50}} = \frac{\gamma \cdot E_{\max}}{4 \cdot EC_{50}}$$

When $\gamma = 1$, the effect goes from 20% to 80% of maximum over an 81-fold concentration range. When $\gamma = 3$, this range narrows to only 4.3-fold, creating a steep "switch-like" response that can be clinically dangerous if the dose varies near $EC_{50}$.

Derivation 2: Therapeutic Index and Safety Margins

The therapeutic index (TI) quantifies the margin between therapeutic and toxic doses:

$$TI = \frac{TD_{50}}{ED_{50}}$$

For clinical safety, the certain safety factor (CSF) is more conservative:

$$CSF = \frac{TD_1}{ED_{99}}$$

Using the sigmoid Emax model, $ED_x$ (the dose producing x% of max effect) is:

$$ED_x = EC_{50} \cdot \left(\frac{x}{100 - x}\right)^{1/\gamma}$$

Therefore:

$$CSF = \frac{TD_{50}}{ED_{50}} \cdot \left(\frac{1/99}{99/1}\right)^{1/\gamma_{\text{tox}} - 1/\gamma_{\text{eff}}}$$

Drugs with TI less than 2 (warfarin, lithium, digoxin, theophylline, phenytoin) require therapeutic drug monitoring. The narrower the TI, the more critical are factors affecting drug levels: renal/hepatic impairment, drug interactions, genetic polymorphisms, and adherence.

Derivation 3: Receptor Selectivity

Selectivity is the ratio of potencies at different receptor subtypes. For a drug acting at receptor A (desired) and receptor B (undesired):

$$\text{Selectivity ratio} = \frac{EC_{50,B}}{EC_{50,A}}$$

At a given concentration C, the differential activation is:

$$\Delta E = \frac{C}{C + EC_{50,A}} - \frac{C}{C + EC_{50,B}}$$

The concentration that maximizes the selectivity window (maximum $\Delta E$) is:

$$C_{\text{opt}} = \sqrt{EC_{50,A} \cdot EC_{50,B}}$$

At this optimal concentration, the maximum selectivity window is:$\Delta E_{\max} = \left(\frac{\sqrt{S} - 1}{\sqrt{S} + 1}\right)^2$ where$S = EC_{50,B}/EC_{50,A}$. A 100-fold selectivity ratio gives only 67% differential activation, explaining why even "selective" drugs have off-target effects at high doses.

Derivation 4: Effect Compartment Model

The effect compartment resolves PK/PD hysteresis by introducing a hypothetical compartment representing the biophase (site of action):

$$\frac{dC_e}{dt} = k_{e0}(C_p - C_e)$$

Where $k_{e0}$ is the equilibration rate constant. The effect compartment has negligible volume (does not affect PK). After an IV bolus:

$$C_e(t) = \frac{D \cdot k_{e0}}{V_d(k_{e0} - k_e)}\left(e^{-k_e t} - e^{-k_{e0} t}\right)$$

The time to peak effect site concentration:

$$t_{\text{peak,e}} = \frac{\ln(k_{e0}/k_e)}{k_{e0} - k_e}$$

The $t_{1/2,k_{e0}}$ characterizes the delay between plasma concentration changes and effect. For propofol, $t_{1/2,k_{e0}} \approx 2.6$ min (rapid equilibration with brain). For warfarin, $t_{1/2,k_{e0}} \approx 24$ hours (slow clotting factor turnover).

Derivation 5: Indirect Response Models

Many drugs act by modifying the production or elimination of an endogenous mediator. The indirect response model describes the response variable R:

$$\frac{dR}{dt} = k_{\text{in}} - k_{\text{out}} \cdot R$$

At baseline, $R_0 = k_{\text{in}}/k_{\text{out}}$. Four types exist based on whether the drug stimulates or inhibits production or elimination:

$$\text{Type I:} \quad \frac{dR}{dt} = k_{\text{in}}\left(1 - \frac{I_{\max} \cdot C}{IC_{50} + C}\right) - k_{\text{out}} \cdot R$$
$$\text{Type II:} \quad \frac{dR}{dt} = k_{\text{in}} - k_{\text{out}}\left(1 + \frac{S_{\max} \cdot C}{SC_{50} + C}\right) R$$

The maximum drug-induced change in the response at steady state (for Type I):

$$R_{\min} = R_0 \cdot (1 - I_{\max})$$

Type I: inhibition of production (warfarin inhibiting clotting factor synthesis, corticosteroids suppressing cortisol). Type II: stimulation of elimination (diuretics enhancing sodium excretion). Type III: stimulation of production (erythropoietin). Type IV: inhibition of elimination (probenecid blocking uric acid excretion).

Clinical Applications

Warfarin PK/PD

Warfarin exemplifies indirect response pharmacodynamics. Despite a plasma half-life of ~40 hours, the anticoagulant effect takes 3-5 days to develop because the drug inhibits synthesis of clotting factors whose own half-lives are 6-60 hours (factor VII ~6h, factor II ~60h). The INR response follows a Type I indirect model.

Anesthetic Depth

The effect compartment concept is critical for titrating IV anesthetics. Propofol's rapid brain equilibration ($k_{e0} = 0.26$ /min) allows precise control of anesthetic depth. Target-controlled infusion (TCI) pumps use PK/PD models to maintain desired effect-site concentrations automatically.

Antimicrobial PK/PD

Three PK/PD indices predict antimicrobial efficacy: $C_{\max}/MIC$ for aminoglycosides (concentration-dependent killing), $T > MIC$ for beta-lactams (time-dependent killing), and $AUC/MIC$ for fluoroquinolones (mixed killing pattern).

Disease Progression Models

Modern PD models incorporate disease progression as a baseline drift:$R(t) = R_0 + \alpha t + E(t)$. This is critical for Alzheimer's trials (ADAS-cog worsening) and oncology (tumor growth kinetics) where placebo patients worsen over time.

Python Simulations

Dose-Response Curves: Potency, Efficacy, and Therapeutic Index

Python
script.py84 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Receptor Selectivity and PK/PD Effect Compartment Model

Python
script.py98 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Tolerance Development and Receptor Desensitization

Python
script.py104 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Direct vs Indirect Response PK/PD Models

Python
script.py90 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Rate this chapter: