Advancing Drug Discovery and Development through Clinical Pharmacological Approaches on Organs-on-Chips

2026 KBCS (Korean BioChip Society) · Medical Unmet Needs for BioChips — Invited Lecture

Sungpil Han, MD, PhD

Dept. of Pharmacology, College of Medicine, The Catholic University of Korea · Dept. of Clinical Pharmacology & Therapeutics, The Catholic University of Korea, Seoul St. Mary’s Hospital · PIPET

Roadmap

KBCS 2026 · Medical Unmet Needs for BioChips

  • Unmet need & regulatory shift — why preclinical systems fail, and why FDA is moving toward NAMs
  • Clinical pharmacology — PBPK & QSP as the bridge
  • Our engine — an LLM-augmented, open QSP library (258 diseases and growing)
  • Case studies — 12 vignettes, 6 deep dives, incl. a live app (Merigolix)
  • Our organ-chip work — Gut–Liver–Kidney MPS + PBPK, condensed
  • Where this goes — a browser-based platform, and a challenge to this community

The unmet need: why preclinical systems fail translation

Metric Value
Clinical-phase attrition ~90%
Time per approved drug 10–15 years
Cost per approved drug $1–2.6 billion
  • 2D cell cultures — no tissue architecture, no organ-level function
  • Animal models — persistent interspecies differences in drug disposition
  • Phase 1 trials characterize healthy-volunteer exposure well, but not the patients who need it most: IBD AKI Hepatic impairment Pregnancy Pediatric Rare disease

DiMasi et al., J Health Econ 2016; Kola & Landis, Nat Rev Drug Discov 2004.

The regulatory shift: from an animal mandate to NAMs

  • Since 1938, U.S. law required animal testing before human use
  • FDA Modernization Act 2.0 (2022) — replaced “animal tests” with “nonclinical tests,” explicitly opening the door to in vitro/in silico methods
  • FDA’s Roadmap to Reducing Animal Testing in Preclinical Safety Studies (2025) — phased rollout starting with monoclonal antibodies; goal: animal testing becomes the exception, not the default
  • Motivation: many drugs that pass animal testing still fail in humans

FDA Modernization Act of 2022, Pub. L. 117–328; FDA, Roadmap to Reducing Animal Testing in Preclinical Safety Studies (2025); Van Norman GA, JACC Basic Transl Sci 4(7):845–854 (2019).

Why this matters for organs-on-chips

  • FDA Modernization Act 2.0 opens the door for NAMs — including MPS — to replace some animal studies
  • MIDD already expects PBPK/QSP-grade quantitative packages; folding MPS data into that workflow is the credible path to regulatory acceptance
  • Organ-on-chip inside a clinical-pharmacology workflow = a translational bridge: less animal testing, faster candidate selection, more reliable human prediction

The NAM translation loop — where clinical pharmacology sits

PBPK and QSP are the mechanistic models this loop depends on — the qualification & translation step that AI/ML pattern-matching alone cannot supply.

Cao Y. New Approach Methodologies: What Clinical Pharmacologists Should Prepare For. Clin Pharmacol Ther 118(6):1269–1272 (2025).

From a chip curve to a dosing decision

  • What a chip gives you: a concentration–time curve · single- or multi-organ readouts · a mechanistically rich human-cell signal
  • What a dosing decision needs: clearance, tissue distribution, target engagement — for a whole human body, not a chip well
  • A chip curve is a descriptive observation until it passes through a mechanistic model

Clinical pharmacology — NCA, PBPK, QSP — is that bridge, and the subject of this talk.

What is QSP & PBPK?

PBPK and QSP describe the causal chain from drug to patient, mathematically:

𝐝𝐫𝐮𝐠𝐭𝐚𝐫𝐠𝐞𝐭𝐩𝐚𝐭𝐡𝐰𝐚𝐲𝐝𝐢𝐬𝐞𝐚𝐬𝐞𝐩𝐚𝐭𝐢𝐞𝐧𝐭\textbf{drug} \;\rightarrow\; \textbf{target} \;\rightarrow\; \textbf{pathway} \;\rightarrow\; \textbf{disease} \;\rightarrow\; \textbf{patient}

  • PBPK — where does the drug go, and at what concentration, in a whole human body?
  • QSP — given that concentration, what happens to the disease?
  • Together: dose selection, DDI risk, response in special populations, IVIVE of chip-derived parameters

Anatomy of a PBPK/QSP model

A coupled nonlinear dynamical system — an initial-value problem:

d𝐱dt=𝐟(𝐱(t),𝛉,𝐮(t)),𝐱(0)=𝐱0\frac{d\mathbf{x}}{dt} = \mathbf{f}\!\big(\mathbf{x}(t),\, \boldsymbol{\theta},\, \mathbf{u}(t)\big), \qquad \mathbf{x}(0)=\mathbf{x}_0

  • States 𝐱\mathbf{x} — concentrations, receptor occupancy, signaling species, cell populations, biomarkers, endpoints
  • Parameters 𝛉\boldsymbol{\theta} — rate constants, EC50EC_{50}/EmaxE_{max}, Hill nn, clearances, binding affinities
  • Inputs 𝐮(t)\mathbf{u}(t) — dosing regimen: the one controllable handle

Our library’s models: 15–35+ coupled ODEs, 70–100+ parameters per disease.

The recurring nonlinear structures

The same motifs recur across every disease — and every chip-to-human bridge:

Motif Governing idea Typical use
Hill / EmaxE_{max} E(C)=EmaxCn/(EC50n+Cn)E(C) = E_{max}C^n / (EC_{50}^n + C^n) any drug dose–response
Indirect response dR/dt=kinf()koutRdR/dt = k_{in}f(\cdot) - k_{out}R biomarker turnover
Mass-action / TMDD drug–target binding, target-mediated clearance mAbs, receptor-saturating drugs
Logistic growth dN/dt=rN(1N/K)kill(C)NdN/dt = rN(1-N/K) - kill(C)N tumor burden under therapy
Power-law feedback e.g. EPO(Hb0/Hb)γEPO \propto (Hb_0/Hb)^{\gamma} homeostatic feedback loops

Common grammar across a 258-model library and a chip–PBPK bridge.

Why this matters — Model-Informed Drug Development

Stage What PBPK/QSP contributes Effect on attrition
Target validation does target modulation reach the phenotype? kills implausible targets early
First-in-human dose exposure–response & safety margins from mechanism safer, leaner Phase 1
Trial design virtual patients → enrichment, endpoints, sample size higher trial success
Translation chip→animal→human, adult→pediatric extrapolation better external validity
Combinations synergy/antagonism in silico rational polypharmacy

FDA/EMA actively endorse MIDD — an established pathway, not a novelty.

The bottleneck: building one QSP model, traditionally

  • Deep literature synthesis across pathophysiology + pharmacology
  • Careful ODE formulation, parameterization, implementation, calibration
  • Cost: expert-months to years per disease

Consequences:

  • only a handful of diseases ever get a model
  • models are siloed, rarely open, rarely reproducible
  • the long tail — most chip programs’ target diseases — never gets one

One disease ≈ one expert-year+. Scarce · slow · closed.

Thesis — LLM-augmented model generation

An autonomous LLM coding agent builds complete, quality-gated, referenced QSP models — one disease per day — every artifact under version control.

  • 🗺️ Graphviz map · ⚙️ mrgsolve ODE model · 📊 Shiny dashboard · 📚 curated references
  • enforced quality gates, literature grounding for every parameter
  • standardized schema → comparable, auditable, reproducible

The Claude Code Routine — an autonomous daily loop

Each run completes one disease end-to-end and pushes it:

  1. Select an uncovered disease
  2. Research — ~50 PubMed sources
  3. Map — Graphviz network (≥100 nodes)
  4. Modelmrgsolve ODEs (≥15 compartments)
  5. Dashboard — Shiny (≥6 tabs)
  6. Ground — sectioned references (≥30)
  7. Commit & push — under a git stop-hook that won’t allow a half-done session

Guardrails — why trust code an LLM wrote

  • Quality gates: map ≥100 nodes/≥8 clusters · ODE ≥15 compartments/≥5 scenarios · dashboard ≥6 tabs · ≥30 references
  • Structural templates: reusable ODE motifs (Hill, turnover, TMDD) shrink the failure space
  • Literature grounding: every parameter set cites a calibration trial — the main anti-hallucination mechanism
  • Reproducibility: diffable, auditable, re-runnable code in git; human review before use

Why an LLM/AI is the right engine

  • Tireless & continuous — a full disease model nearly every day
  • Vast literature synthesis — ~12,800 citations (~50/model)
  • Breadth — 258 diseases, hundreds of drugs/targets, one framework
  • Consistency — identical schema and quality bar
  • Economics — expert-months → hours

The scarce resource is modeling expertise. The LLM scales it.

What’s inside each model

  • 🗺️ Mechanistic map — clustered directed graph, targets to disease
  • ⚙️ mrgsolve ODE model — PK coupled to PD; 15–35+ states, ≥5 scenarios
  • 📊 Shiny dashboard — patient profile, PK, endpoints, scenarios (6–8 tabs)
  • 📚 ~50 curated references per model

The exact structure a chip-to-human bridge needs to plug into.

Scale: a 258-model open QSP library

Metric Value
Disease QSP models 258
Therapeutic areas ~18
mrgsolve ODE systems 259
Pathway clusters ~3,100
PubMed references ~12,800
New models added +1 / day

Traditional QSP: a handful of diseases, expert-months each. This library: 258 diseases, hours each.

Breadth of coverage

Oncology Autoimmune/rheumatic Vasculitis Cardiovascular Respiratory Renal/urologic GI/hepatobiliary Endocrine/metabolic Neurologic Psychiatric Dermatologic Infectious Ophthalmic Rare/genetic

  • Cancers — breast, NSCLC, SCLC, glioblastoma, CML, multiple myeloma, melanoma…
  • Rare & genetic — Fabry, Gaucher, DMD, SMA, Huntington, ATTR amyloidosis…
  • Common chronic — diabetes, heart failure, COPD, CKD…
  • Immune / hematologic — lupus, RA, sickle cell, ITP, myelofibrosis…

Hundreds of drugs and molecular targets — very likely including yours.

Case studies: drugs · targets · modeled endpoints

Disease Drug (example) Target / mechanism Modeled endpoint
Rheumatoid arthritis tocilizumab IL-6 receptor DAS28 · CRP
Psoriasis secukinumab IL-17A PASI
Type 2 diabetes semaglutide GLP-1 receptor HbA1c · weight
Heart failure (rEF) sacubitril/valsartan neprilysin + AT₁ (ARNI) LVEF · NT-proBNP
COPD triple therapy (ICS/LABA/LAMA) β2/M3/glucocorticoid receptors FEV1 · exacerbation rate
Non-small-cell lung cancer osimertinib mutant EGFR tumor burden
Chronic myeloid leukemia imatinib BCR-ABL BCR-ABL transcript ratio
Pulmonary arterial HTN macitentan endothelin A/B receptor PVR · 6-min walk
IgA nephropathy sparsentan endothelin-A + AT₁ UPCR · eGFR
Sickle cell disease voxelotor HbS–O₂ affinity Hb · hemolysis
Multiple myeloma daratumumab CD38 (TMDD) M-protein
Endometriosis Merigolix (GnRH antagonist) GnRH receptor pain · lesion size

Every drug–target–endpoint link is grounded in the trial literature.

Case studies — one modeling insight each (1/4)

NSCLC — EGFR resistance

Osimertinib vs. T790M/C797S escape clones: timing → combination or sequencing.

HFrEF — GDMT sequencing

ARNI · β-blocker · MRA · SGLT2i order and up-titration speed vs. LVEF, NT-proBNP.

ATTR amyloidosis — gene silencing

TTR knockdown (patisiran/vutrisiran) → tetramer–monomer kinetics over years.

Case studies — one modeling insight each (2/4)

SLE — type-I interferon

Anifrolumab (anti-IFNAR): IFN → autoantibody → flare risk (SLEDAI).

T2DM — incretin axis

Tirzepatide (GLP-1/GIP): glucose–insulin–weight loop → HbA1c, weight.

ITP — platelet kinetics

TPO-RA vs. anti-platelet antibody: production vs. destruction balance.

Case studies — one modeling insight each (3/4)

Osteoporosis — transient uncoupling

Romosozumab: a ~12-month anabolic window before resorption rebounds.

Migraine — occupancy, not concentration

Erenumab: slow clearance sustains receptor occupancy between monthly doses.

Atopic dermatitis — dual blockade

Dupilumab: biomarkers (IgE, TARC) move early; EASI response lags.

Case studies — one modeling insight each (4/4)

Hemophilia A — a PD readout that lies

Emicizumab prolongs aPTT far past its real hemostatic benefit.

Myasthenia gravis — clearance vs. symptom lag

Efgartigimod: IgG drops ~60% in days; MG-ADL improvement lags.

Gout — when the drug beats itself

Pegloticase: anti-drug antibodies predict loss of response.

Deep dive — IgA nephropathy: the “four-hit” cascade

  • Galactose-deficient IgA1 → autoantibodies → immune complexes → mesangial deposition + complement → proteinuria, falling eGFR
  • Drugs modeled: RAAS blockade · budesonide-TRF · sparsentan (dual ETA+AT₁) · iptacopan (factor B) · sibeprenlimab (anti-APRIL)
  • Endpoints: UPCR, eGFR slope — NefIgArd, PROTECT, APPLAUSE-IgAN
  • 20 ODEs · 7 scenarios

Deep dive — sickle cell disease: polymerization to vaso-occlusion

  • Deoxy-HbS polymerizes → sickling & hemolysis → NO scavenging → endothelial activation → vaso-occlusion
  • Drugs modeled: hydroxyurea (HbF induction) · voxelotor (HbS–O₂ affinity) · crizanlizumab (anti-P-selectin) · L-glutamine
  • Endpoints: Hb, HbF%, vaso-occlusion rate, LDH — MSH, HOPE, SUSTAIN
  • 24 ODEs; power-law erythropoiesis feedback

Deep dive — multiple myeloma: an oncology exemplar

  • Plasma-cell clones secrete M-protein; IL-6 autocrine growth; RANKL/OPG/DKK1 imbalance → bone disease
  • Drugs modeled: bortezomib/carfilzomib (proteasome) · lenalidomide (cereblon) · daratumumab (anti-CD38, TMDD) · venetoclax · zoledronate
  • Endpoints: M-protein, free light chains — VRd/DRd/KRd/DVRd
  • Logistic growth + resistant clone + TMDD, coupled to bone-remodeling

Deep dive — rheumatoid arthritis: IL-6/TNF/JAK-driven synovitis

  • Synovial IL-6/TNF/JAK signaling → pannus → cartilage/bone erosion; CRP tracks IL-6 trans-signaling
  • Drugs modeled: methotrexate · tocilizumab (anti-IL-6R, TMDD) · adalimumab (anti-TNF) · baricitinib (JAK1/2)
  • Endpoints: DAS28, CRP, bone erosion
  • 19 ODEs · 7 scenarios

Deep dive — Duchenne muscular dystrophy: a genetic exemplar

  • Dystrophin loss → sarcolemmal fragility → chronic necrosis/regeneration → fibro-fatty replacement
  • Drugs modeled: deflazacort/prednisone · eteplirsen (exon-51-skipping ASO) · givinostat (HDAC inhibitor)
  • Endpoints: dystrophin (% of normal), functional decline
  • 22 ODEs · 6 scenarios

Deep dive — Alzheimer’s disease: amyloid to cognition

  • Amyloid-β production/clearance imbalance → plaques → synaptic/cholinergic loss → cognitive decline; APOE4 accelerates
  • Drugs modeled: donepezil (AChE inhibitor) · memantine (NMDA antagonist) · lecanemab (anti-amyloid mAb)
  • Endpoints: amyloid burden, MMSE-proxy cognition
  • 20 ODEs · 7 scenarios

Flagship case — Merigolix, a live drug–disease model

Merigolix — oral GnRH-receptor antagonist for endometriosis: ↓LH/FSH → ↓estradiol → lesion/pain regression.

  • 2-compartment PK + HPG-axis PD (E2/LH/FSH turnover, Hill inhibition)
  • Endpoints: pain (NRS), lesion size, hot flashes, BMD, endometrium

Try it live: pipetqsp.shinyapps.io/merigolix

Merigolix — the estrogen “threshold” trade-off

The therapeutic window

  • too little E2 suppression → no relief
  • too much → hot flashes, bone loss
  • target a partial E2 band (~20–40 pg/mL)

This simulation run

  • E2 trough 2.6 pg/mL; pain 3.1 → 1.1; lesion 18.5 → 5.3 mm
  • BMD −1.22%; hot flashes 5.4 → 7.9/day
  • → motivates titration + hormonal add-back

QSP turns this trade-off into a tunable optimization problem — a chip program’s question too.

How this helps your chip program, directly

  • Working on a cytokine, transporter, receptor, or disease? A model likely already exists
  • Use it to generate a hypothesis before the chip experiment
  • Use it as a PD/disease layer on top of your own chip→PBPK bridge

Not “run all 258” — a platform to accelerate translation of what you already work on.

Fully open — every line of code

  • github.com/pipetcpt/qsp — every map, ODE model, dashboard, reference, openly versioned
  • Clone it, run a model, adapt it, or just browse for ideas
  • QR code, top-right corner of every slide, links to the repository

Where this goes: a browser-based PBPK/QSP platform

A web platform where any researcher can, from a browser alone —

  • search by disease, target, drug, or cytokine
  • browse the map and ODE structure with no install
  • run a simple simulation — a dose, a tweak, an instant plot
  • export parameters into your own workflow

Prerequisites for trusting a chip-derived number

Four checks before a chip number belongs inside a PBPK/QSP model:

  1. Chip-specific drug behavior — adsorption to plastics/membranes
  2. Scaling of cellular and fluidic dimensions — correct per-organ factors
  3. Verification against held-out clinical data
  4. Explicit, falsifiable correction rules — not curve-fitting

The next case study tests whether this discipline works in practice.

Our case study: Gut–Liver–Kidney MPS + PBPK

Metric Value
Drugs 6
Clinical observations 24
AUC GMFE (overall) 1.85×
Within 4-fold 100%
Clinical PK used to fit 0
Disease-state chip arms 2
  • Modular gut–liver–kidney MPS (hiEC–HepaRG–RPTEC), shared media recirculation
  • Chip-fitted parameters propagate to whole-body PBPK — no human PK enters the prediction
  • All 24 observations reserved for post-hoc validation

The model predicts human PK well

Predicted vs. observed AUC and Cmax across 24 observations, 6 drugs. Inner band = 2-fold; outer band = 4-fold.

How it works — chip kinetics to whole-body PBPK

  • Chip → kinetics: 33-state ODE fits chip data (parent + metabolite); Kpuu,hepK_{puu,hep} captures OATP-style hepatic uptake
  • Whole-body PBPK: 10-state perfusion-limited; Rodgers–Rowland KpK_p; Varma clearance; ACAT absorption
  • Inputs: physchem only (logP, pKapK_a, fu,pf_{u,p}, blood:plasma, MW, dose) — no clinical PK to fit

How it works — the falsifiable scaling rule

  • Apply a relative-activity factor (RAF) only if CYP3A4 ≥90%, biliary clearance <50%, hepatic uptake not already captured
  • Mechanism-based test, not curve-fitting
  • Uniform correction → GMFE 9.7×; selective rule → 1.85×
  • Disease-state extension: re-fit under DSS colitis / cisplatin injury; unaffected-organ parameters locked; same pipeline

Validation results — 24 observations, 6 drugs

Drug Mechanism AUC GMFE 2-fold 4-fold
Antipyrine CYP1A2, low E 1.59× 4/4 4/4
Benzylpenicillin OAT1/3 renal 1.23× 4/4 4/4
Testosterone CYP3A4 high-E + SHBG 1.98× 2/4 4/4
Crizotinib CYP3A4 high-E + biliary 2.61× 1/4 4/4
Gefitinib CYP3A4 + biliary-dominant 2.93× 0/4 4/4
Simvastatin OATP1B1 + CYP3A4 1.37× 4/4 4/4
Overall (n=24) 1.85× 15/24 24/24

Within the 2–5× GMFE benchmark of industry-standard PBPK platforms.

Disease-state extension: gut vs. hepatic localization

  • Gefitinib, DSS colitis chip: AUC 0.49× — gut-localized (Fg shifts, Fa/Fh unchanged)
  • Simvastatin, cisplatin-injury chip: AUC 2.19× — hepatic-localized (Kpuu,liverK_{puu,liver} halved, OATP1B1 suppression)
  • Same pipeline distinguishes gut vs. hepatic shifts — hypothesis-generating for IBD/AKI patients usually excluded from Phase 1

Limitations, and what this means

  • No mucus/peristalsis layer — affects chip PappP_{app}
  • HepaRG 2D CYP3A4 ~10% of in vivo (RAF-corrected) · RPTEC OAT1/3 ~1% (PTC-anchor corrected)
  • Sparse clinical PK in disease-state patients — roadmap: prospective IBD/AKI cohorts
  • Open-source Python — full run reproduces in <5 minutes

Augments, does not replace, animal/clinical PK.

Clinical pharmacology’s unmet needs

Conclusion 1 — Unmet needs in our own field

  • A real shortage of clinical pharmacologists trained to interpret NAM signal through PBPK/QSP
  • No standardized workflow yet for FDA’s four NAM validation principles: context-of-use, biological relevance, technical characterization, fit-for-purpose
  • Few chip/organoid models validated head-to-head against human clinical data
  • The bottleneck isn’t technology — it’s people who can translate

Not a data problem, not a technology problem — a workforce problem.

A challenge to KBCS members

Conclusion 2 — A direct ask

  • Design chips with PBPK/QSP-ready outputs in mind — clearance, permeability, absorption
  • Let’s validate your chip data against held-out clinical data, together
  • Try the matching model in our 258-model library first — tell us where it’s wrong
  • Waiting to collaborate: github.com/pipetcpt/qsp · shan@catholic.ac.kr

If your chip already speaks a language PBPK/QSP can read, half the problem is solved before we meet.

Conclusion & an open invitation

  • Organs-on-chips generate human-relevant biology; clinical pharmacology turns it into decision-grade evidence
  • Our Gut–Liver–Kidney case study: achievable without clinical PK input, extends to disease-state chips
  • Our 258-model open QSP library exists so the NAMs/MPS community can adopt this now

This project is completely open. Building an organ-on-chip and want a PBPK/QSP bridge — or a disease model to build toward? I would love to collaborate.

Sungpil Han, MD, PhD shan@catholic.ac.kr

Selected references

Cao Y. New Approach Methodologies: What Clinical Pharmacologists Should Prepare For. Clin Pharmacol Ther 118(6), 1269–1272 (2025).

FDA. Roadmap to Reducing Animal Testing in Preclinical Safety Studies (2025).

FDA Modernization Act of 2022, Pub. L. 117–328.

Van Norman GA. Limitations of Animal Studies for Predicting Toxicity in Clinical Trials. JACC Basic Transl Sci 4(7), 845–854 (2019).

DiMasi JA et al. J Health Econ 47, 20–33 (2016).

Kola I & Landis J. Nat Rev Drug Discov 3, 711–715 (2004).

Huh D et al. Science 328, 1662–1668 (2010).

Bhatia SN & Ingber DE. Nat Biotechnol 32, 760–772 (2014).

Low LA et al. Nat Rev Drug Discov 20, 345–361 (2021).

Tsamandouras N et al. AAPS J 19, 1499–1512 (2017).

Edington CD et al. Sci Rep 8, 4530 (2018).

Kuepfer L et al. CPT Pharmacometrics Syst Pharmacol 5, 516–531 (2016).

Varma MV et al. Pharm Res 32, 3785–3802 (2015).

Guillouzo A et al. Chem Biol Interact 168, 66–73 (2007).

Mathialagan S et al. Drug Metab Dispos 45, 409–417 (2017).

Rodgers T & Rowland M. J Pharm Sci 95, 1238–1257 (2006).

Yang J et al. Curr Drug Metab 8, 676–684 (2007).

Contact

Sungpil Han, M.D./Ph.D. — Associate Professor

Dept. of Pharmacology, College of Medicine, The Catholic University of Korea
Dept. of Clinical Pharmacology & Therapeutics, The Catholic University of Korea, Seoul St. Mary’s Hospital
PIPET (Pharmacometrics Institute for Practical Education & Training)

222 Banpodaero, Seocho-gu, Seoul, Korea (06591)

Email: shan@catholic.ac.kr
Phone: +82-2-3147-8356 · Mobile: +82-10-6782-0522 · FAX: +82-2-2258-7876

github.com/pipetcpt/qsp

Thank you — questions and collaborators welcome.

open-source
github.com/pipetcpt/qsp