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:
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$:
The slope of the concentration-effect curve at $EC_{50}$ is:
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:
For clinical safety, the certain safety factor (CSF) is more conservative:
Using the sigmoid Emax model, $ED_x$ (the dose producing x% of max effect) is:
Therefore:
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):
At a given concentration C, the differential activation is:
The concentration that maximizes the selectivity window (maximum $\Delta E$) is:
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):
Where $k_{e0}$ is the equilibration rate constant. The effect compartment has negligible volume (does not affect PK). After an IV bolus:
The time to peak effect site concentration:
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:
At baseline, $R_0 = k_{\text{in}}/k_{\text{out}}$. Four types exist based on whether the drug stimulates or inhibits production or elimination:
The maximum drug-induced change in the response at steady state (for Type I):
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
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server
Receptor Selectivity and PK/PD Effect Compartment Model
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server
Tolerance Development and Receptor Desensitization
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server
Direct vs Indirect Response PK/PD Models
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server