Drugs acting on the heart, the blood vessels and circulating lipids. The large classes are developed against continuous biomarkers, blood pressure and LDL cholesterol [1,2], which is what makes dense exposure-response work practical here.
In July 1973, after screening thousands of fungal broths at Sankyo in Tokyo, Akira Endo isolated the first statin, ML-236B, from a Penicillium citrinum strain found on rice from a Kyoto grain shop [3].
Modeling notes
Endpoint. Continuous surrogates measured repeatedly: blood pressure, LDL cholesterol [4,5]. Both lag plasma concentration, because the drug acts on the turnover of the measured quantity [6,7].
Typical structure. Indirect response models [8,9], near-universal for LDL [10,11] and common for blood pressure. Effect-compartment models only where the delay is genuinely short, as for the acute heart rate effect of a beta blocker [12,13].
Counter-regulation (C02, C03, C09). Baroreflex activation, renin release and renal sodium handling push back against blood pressure effects [14], which is why feedback and tolerance terms appear in antihypertensive models [15].
Baseline is the trap. Blood pressure has a circadian rhythm, large measurement error and regression to the mean on enrolment [16,17], and placebo arms show genuine falls [18,19]. How baseline is handled changes bias and precision in the estimated drug effect, and normalising each observation by its own baseline is the worst common option [20,21].
Covariates worth pre-specifying. Renal function for diuretics and renally cleared renin-angiotensin agents [22,23], CYP2D6 for metoprolol and flecainide [24,25], SLCO1B1 for statins [26]. A data-driven search is underpowered at realistic dataset sizes [27,28].
Safety. A concentration-QTc analysis is expected under ICH E14 for essentially every new drug, as a separate pre-specified analysis [29,30] (see C-QTc modeling). Class-matched effects such as hypotension and bradycardia are exaggerated on-target pharmacology, and keep climbing after the efficacy curve flattens [31,32].
Approval rests on the surrogate. Approval often rests on the surrogate, with the outcome trial reporting years later [33,34]. The pharmacometric job is linking biomarker change to event hazard [35,36], and handling the adherence and dropout endemic to chronic outpatient dosing [37,38].
References
[1]
Kim CO, Jeon S, Han S, Hong T, Park MS, Yoon Y-R, et al. Decreased potency of fimasartan in liver cirrhosis was quantified using mixed-effects analysis. Translational and Clinical Pharmacology 2017;25:43.
https://doi.org/10.12793/tcp.2017.25.1.43.
[2]
Jadhav SB, Amore BM, Bockbrader H, Crass RL, Chapel S, Sasiela WJ, et al. Population pharmacokinetic and pharmacokinetic-pharmacodynamic modeling of bempedoic acid and low-density lipoprotein cholesterol in healthy subjects and patients with dyslipidemia. Journal of Pharmacokinetics and Pharmacodynamics 2023;50:351–64.
https://doi.org/10.1007/s10928-023-09864-w.
[4]
Biomarkers Definitions Working Group. Biomarkers and surrogate endpoints: Preferred definitions and conceptual framework. Clinical Pharmacology & Therapeutics 2001;69:89–95.
https://doi.org/10.1067/mcp.2001.113989.
[5]
Liaigre L, Guigui A, Manceau M, Cracowski J-L, Khouri C, Roustit M. Trial-level surrogacy of non-high-density and low-density lipoprotein cholesterol reduction on the clinical efficacy of statins. European Heart Journal - Cardiovascular Pharmacotherapy 2025;11:387–92.
https://doi.org/10.1093/ehjcvp/pvaf016.
[7]
Nicolas X, Djebli N, Rauch C, Brunet A, Hurbin F, Martinez J-M, et al. Population pharmacokinetic/pharmacodynamic analysis of alirocumab in healthy volunteers or hypercholesterolemic subjects using an indirect response model to predict low-density lipoprotein cholesterol lowering: Support for a biologics license application submission: Part II. Clinical Pharmacokinetics 2019;58:115–30.
https://doi.org/10.1007/s40262-018-0670-5.
[8]
Dayneka NL, Garg V, Jusko WJ. Comparison of four basic models of indirect pharmacodynamic responses. Journal of Pharmacokinetics and Biopharmaceutics 1993;21:457–78.
https://doi.org/10.1007/BF01061691.
[9]
Jusko WJ, Ko HC. Physiologic indirect response models characterize diverse types of pharmacodynamic effects. Clinical Pharmacology & Therapeutics 1994;56:406–19.
https://doi.org/10.1038/clpt.1994.155.
[10]
Kakara M, Nomura H, Fukae M, Gotanda K, Hirota T, Matsubayashi S, et al. Population pharmacodynamic analysis of
LDL-cholesterol lowering effects by statins and co-medications based on electronic medical records. British Journal of Clinical Pharmacology 2014.
https://doi.org/10.1111/bcp.12405.
[12]
Francheteau P, Steimer J-L, Merdjan H, Guerret M, Dubray C. A mathematical model for dynamics of cardiovascular drug action: Application to intravenous dihydropyridines in healthy volunteers. Journal of Pharmacokinetics and Biopharmaceutics 1993.
https://doi.org/10.1007/BF01059111.
[13]
Belz GG, Essig J, Kleinbloesem CH, Hoogkamer JF, Wiegand UW, Wellstein A. Interactions between cilazapril and propranolol in man; plasma drug concentrations, hormone and enzyme responses, haemodynamics, agonist dose-effect curves and baroreceptor reflex. British Journal of Clinical Pharmacology 1988.
https://doi.org/10.1111/j.1365-2125.1988.tb05294.x.
[15]
Brynne L, McNay JL, Schaefer HG, Swedberg K, Wiltse CG, Karlsson MO. Pharmacodynamic models for the cardiovascular effects of moxonidine in patients with congestive heart failure. British Journal of Clinical Pharmacology 2001.
https://doi.org/10.1046/j.1365-2125.2001.01320.x.
[16]
Hempel G, Karlsson MO, Alwis DP de, Toublanc N, McNay J, Schaefer H-G. Population pharmacokinetic-pharmacodynamic modeling of moxonidine using 24-hour ambulatory blood pressure measurements. Clinical Pharmacology and Therapeutics 1998.
https://doi.org/10.1016/S0009-9236(98)90053-4.
[18]
Preston RA, Materson BJ, Reda DJ, Williams DW. Placebo-associated blood pressure response and adverse effects in the treatment of hypertension: Observations from a
Department of Veterans Affairs cooperative study. Archives of Internal Medicine 2000;160:1449–54.
https://doi.org/10.1001/archinte.160.10.1449.
[19]
Patel HC, Hayward C, Ozdemir BA, Rosen SD, Krum H, Lyon AR, et al. Magnitude of blood pressure reduction in the placebo arms of modern hypertension trials: Implications for trials of renal denervation. Hypertension 2015;65:401–6.
https://doi.org/10.1161/HYPERTENSIONAHA.114.04640.
[21]
Sun YN, Jusko WJ. Role of baseline parameters in determining indirect pharmacodynamic responses. Journal of Pharmaceutical Sciences 1999.
https://doi.org/10.1021/js9901155.
[22]
Hammarlund-Udenaes M, Benet LZ. Furosemide pharmacokinetics and pharmacodynamics in health and disease: An update. Journal of Pharmacokinetics and Biopharmaceutics 1989;17:1–46.
https://doi.org/10.1007/BF01059086.
[23]
Elung-Jensen T, Heisterberg J, Kamper A-L, Sonne J, Strandgaard S. Blood pressure response to conventional and low-dose enalapril in chronic renal failure. British Journal of Clinical Pharmacology 2003;55:139–46.
https://doi.org/10.1046/j.1365-2125.2003.01764.x.
[24]
Meloche M, Khazaka M, Kassem I, Barhdadi A, Dubé M-P, Denus S de.
CYP2D6 polymorphism and its impact on the clinical response to metoprolol: A systematic review and meta-analysis. British Journal of Clinical Pharmacology 2020;86:1015–33.
https://doi.org/10.1111/bcp.14247.
[25]
Thomas CD, Mosley SA, Kim S, Lingineni K, El Rouby N, Langaee TY, et al. Examination of metoprolol pharmacokinetics and pharmacodynamics across
CYP2D6 genotype-derived activity scores. CPT: Pharmacometrics & Systems Pharmacology 2020;9:678–85.
https://doi.org/10.1002/psp4.12563.
[26]
Türkmen D, Masoli JAH, Kuo C-L, Bowden J, Melzer D, Pilling LC. Statin treatment effectiveness and the
SLCO1B1*5 reduced function genotype: Long-term outcomes in women and men. British Journal of Clinical Pharmacology 2022;88:3230–40.
https://doi.org/10.1111/bcp.15245.
[28]
Sanghavi K, Ribbing J, Rogers JA, Ahmed MA, Karlsson MO, Holford N, et al. Covariate modeling in pharmacometrics: General points for consideration. CPT: Pharmacometrics & Systems Pharmacology 2024;13:710–28.
https://doi.org/10.1002/psp4.13115.
[29]
Garnett C, Bonate PL, Dang Q, Ferber G, Huang D, Liu J, et al. Scientific white paper on concentration-
QTc modeling. Journal of Pharmacokinetics and Pharmacodynamics 2018;45:383–97.
https://doi.org/10.1007/s10928-017-9558-5.
[30]
Darpo B, Benson C, Dota C, Ferber G, Garnett C, Green CL, et al. Results from the
IQ-CSRC prospective study support replacement of the thorough
QT study by
QT assessment in the early clinical phase. Clinical Pharmacology & Therapeutics 2015;97:326–35.
https://doi.org/10.1002/cpt.60.
[32]
Overgaard RV, Ingwersen SH, Tornøe CW. Establishing good practices for exposure-response analysis of clinical endpoints in drug development. CPT: Pharmacometrics & Systems Pharmacology 2015;4:565–75.
https://doi.org/10.1002/psp4.12015.
[33]
Sabatine MS, Giugliano RP, Keech AC, Honarpour N, Wiviott SD, Murphy SA, et al. Evolocumab and clinical outcomes in patients with cardiovascular disease. New England Journal of Medicine 2017;376:1713–22.
https://doi.org/10.1056/NEJMoa1615664.
[34]
Nissen SE, Lincoff AM, Brennan D, Ray KK, Mason D, Kastelein JJP, et al. Bempedoic acid and cardiovascular outcomes in statin-intolerant patients. New England Journal of Medicine 2023.
https://doi.org/10.1056/NEJMoa2215024.
[36]
Zhudenkov K, Gavrilov S, Sofronova A, Stepanov O, Kudryashova N, Helmlinger G, et al. A workflow for the joint modeling of longitudinal and event data in the development of therapeutics: Tools, statistical methods, and diagnostics. CPT: Pharmacometrics & Systems Pharmacology 2022;11:425–37.
https://doi.org/10.1002/psp4.12763.
[37]
Vrijens B, Urquhart J. Methods for measuring, enhancing, and accounting for medication adherence in clinical trials. Clinical Pharmacology & Therapeutics 2014;95:617–26.
https://doi.org/10.1038/clpt.2014.59.