Hybrid Mechanistic-Neural Modeling of Concentration-Time Dynamics Generalizes Human Pharmacokinetics Prediction Across Unseen Chemical Space

Abstract

The pharmacokinetics concentration-time profile encodes ADME dynamics. Its optimization in drug discovery curbs late-stage attrition, yet forecasting human pharmacokinetics from chemical structure remains difficult. Machine learning models ignore continuous ADME dynamics, whereas physiologically based pharmacokinetics (PBPK) models demand unavailable parameters and impose rigid, misspecified compartmental structures. We present PK-MUSE (PharmacoKinetic hybrid Mechanistic Universal System Equation), augmenting a two-compartment model with bounded, state- and time-dependent neural corrections to elimination and intercompartmental distribution. Molecular structure parameterizes compound-specific behavior, while learned corrections refine dynamic equations. Across scaffold and temporal out-of-distribution evaluations, PK-MUSE achieved the lowest mean exposure-endpoint MdFE, reducing MdFE by 7.2%, 45.2%, and 26.2% versus state-of-the-art machine learning (HML- RF), deep learning compartmental (DeepCt), and PBPK (PK-Sim) models. Interrogating the corrections localized misspecification to distribution rather than elimination. By extending mechanistic simulators with data-driven constrained neural dynamics that preserve interpretability and generalize to unseen chemotypes, PK-MUSE advances digital twins for drug discovery.