How effective are ionization state-based QSPKR models at predicting pharmacokinetic parameters in humans?

Optimizing the pharmacokinetics (PK) of a drug candidate to support oral dosing is a key challenge in drug development. PK parameters are usually estimated from the concentration-time profile following intravenous administration; however, traditional methods are time-consuming and expensive. In recent years, quantitative structure-pharmacokinetic relationship (QSPKR), an in silico tool that aims to develop a mathematical relationship between the structure of a molecule and its PK properties, has emerged as a useful alternative to experimental testing. Due to the complex nature of the various processes involved in dictating the fate of a drug, the development of adequate QSPKR models that can be used in real-world pre-screening situations has proved challenging. Given the crucial role played by a molecule's ionization state in determining its PK properties, this work aims to build predictive QSPKR models for PK parameters in humans using an ionization state-based strategy. We divide a high-quality dataset into clusters based on ionization state at physiological pH and build global and ion subset-based 'local' models for three major PK parameters: plasma clearance (CL), steady-state volume of distribution (VD_ss), and half-life (t_1/2). We use a robust methodology developed in our lab entitled 'EigenValue ANalySis' that accounts for the stereospecificity in drug disposition and use the support vector machine algorithm for model building. Our findings suggest that categorizing compounds in accordance with ionization state does not result in improved QSPKR models. The narrow ranges in the endpoints along with redundancies in the data adversely affect the ion subset-based QSPKR models. We suggest alternative approaches such as elimination route-based models that account for drug-transporter interactions for CL and chemotype-specific QSPKR for VD_ss.