1
|
Geci R, Gadaleta D, de Lomana MG, Ortega-Vallbona R, Colombo E, Serrano-Candelas E, Paini A, Kuepfer L, Schaller S. Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans. Arch Toxicol 2024; 98:2659-2676. [PMID: 38722347 PMCID: PMC11272695 DOI: 10.1007/s00204-024-03764-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Accepted: 04/23/2024] [Indexed: 07/26/2024]
Abstract
Physiologically based kinetic (PBK) modelling offers a mechanistic basis for predicting the pharmaco-/toxicokinetics of compounds and thereby provides critical information for integrating toxicity and exposure data to replace animal testing with in vitro or in silico methods. However, traditional PBK modelling depends on animal and human data, which limits its usefulness for non-animal methods. To address this limitation, high-throughput PBK modelling aims to rely exclusively on in vitro and in silico data for model generation. Here, we evaluate a variety of in silico tools and different strategies to parameterise PBK models with input values from various sources in a high-throughput manner. We gather 2000 + publicly available human in vivo concentration-time profiles of 200 + compounds (IV and oral administration), as well as in silico, in vitro and in vivo determined compound-specific parameters required for the PBK modelling of these compounds. Then, we systematically evaluate all possible PBK model parametrisation strategies in PK-Sim and quantify their prediction accuracy against the collected in vivo concentration-time profiles. Our results show that even simple, generic high-throughput PBK modelling can provide accurate predictions of the pharmacokinetics of most compounds (87% of Cmax and 84% of AUC within tenfold). Nevertheless, we also observe major differences in prediction accuracies between the different parameterisation strategies, as well as between different compounds. Finally, we outline a strategy for high-throughput PBK modelling that relies exclusively on freely available tools. Our findings contribute to a more robust understanding of the reliability of high-throughput PBK modelling, which is essential to establish the confidence necessary for its utilisation in Next-Generation Risk Assessment.
Collapse
Affiliation(s)
- René Geci
- esqLABS GmbH, Saterland, Germany.
- Institute for Systems Medicine with Focus on Organ Interaction, University Hospital RWTH Aachen, Aachen, Germany.
| | | | - Marina García de Lomana
- Machine Learning Research, Research and Development, Pharmaceuticals, Bayer AG, Berlin, Germany
| | | | - Erika Colombo
- Istituto di Ricerche Farmacologiche Mario Negri IRCCS, Milan, Italy
| | | | | | - Lars Kuepfer
- Institute for Systems Medicine with Focus on Organ Interaction, University Hospital RWTH Aachen, Aachen, Germany
| | | |
Collapse
|
2
|
Li Y, Wang Z, Li Y, Du J, Gao X, Li Y, Lai L. A Combination of Machine Learning and PBPK Modeling Approach for Pharmacokinetics Prediction of Small Molecules in Humans. Pharm Res 2024; 41:1369-1379. [PMID: 38918309 DOI: 10.1007/s11095-024-03725-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2023] [Accepted: 06/02/2024] [Indexed: 06/27/2024]
Abstract
PURPOSE Recently, there has been rapid development in model-informed drug development, which has the potential to reduce animal experiments and accelerate drug discovery. Physiologically based pharmacokinetic (PBPK) and machine learning (ML) models are commonly used in early drug discovery to predict drug properties. However, basic PBPK models require a large number of molecule-specific inputs from in vitro experiments, which hinders the efficiency and accuracy of these models. To address this issue, this paper introduces a new computational platform that combines ML and PBPK models. The platform predicts molecule PK profiles with high accuracy and without the need for experimental data. METHODS This study developed a whole-body PBPK model and ML models of plasma protein fraction unbound ( f up ), Caco-2 cell permeability, and total plasma clearance to predict the PK of small molecules after intravenous administration. Pharmacokinetic profiles were simulated using a "bottom-up" PBPK modeling approach with ML inputs. Additionally, 40 compounds were used to evaluate the platform's accuracy. RESULTS Results showed that the ML-PBPK model predicted the area under the concentration-time curve (AUC) with 65.0 % accuracy within a 2-fold range, which was higher than using in vitro inputs with 47.5 % accuracy. CONCLUSION The ML-PBPK model platform provides high accuracy in prediction and reduces the number of experiments and time required compared to traditional PBPK approaches. The platform successfully predicts human PK parameters without in vitro and in vivo experiments and can potentially guide early drug discovery and development.
Collapse
Affiliation(s)
- Yuelin Li
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China
| | - Zonghu Wang
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China
| | - Yuru Li
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China
| | - Jiewen Du
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China
| | - Xiangrui Gao
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China
| | - Yuanpeng Li
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China
| | - Lipeng Lai
- XtalPi Innovation Center, XtalPi Inc., Beijing, 100080, China.
| |
Collapse
|
3
|
Najjar A, Grégoire S, Nicol B, Natsch A, Golbamaki N, Boisleve F, Irizar A, Wall B, Swinscoe A, Masini-Etévé V, Selechnik D, Api AM, Griem P, Hewitt N, Cardamone E. In vitro to in vivo extrapolation to derive a metabolism factor for estimating the aggregate exposure to salicylic acid after dermal exposure of its esters. Arch Toxicol 2024; 98:2199-2211. [PMID: 38658404 PMCID: PMC11169020 DOI: 10.1007/s00204-024-03749-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Accepted: 03/21/2024] [Indexed: 04/26/2024]
Abstract
As part of the safety assessment of salicylate esters in cosmetics, we developed a metabolism factor based on in vitro to in vivo extrapolation (IVIVE) to provide a better estimation of the aggregate internal exposure to the common metabolite, salicylic acid. Optimal incubation conditions using human liver S9 were identified before measuring salicylic acid formation from 31 substances. Four control substances, not defined as salicylic esters but which could be mistaken as such due to their nomenclature, did not form salicylic acid. For the remaining substances, higher in vitro intrinsic clearance (CLint, in vitro) values generally correlated with lower LogP values. A "High-Throughput Pharmacokinetic" (HTPK) model was used to extrapolate CLint, in vitro values to human in vivo clearance and half-lives. The latter were used to calculate the percentage of substance metabolised to salicylic acid in 24 h in vivo following human exposure to the ester, i.e. the "metabolism factor". The IVIVE model correctly reproduced the observed elimination rate of 3 substances using in silico or in vitro input parameters. For other substances, in silico only-based predictions generally resulted in lower metabolism factors than when in vitro values for plasma binding and liver S9 CLint, in vitro were used. Therefore, in vitro data input provides the more conservative metabolism factors compared to those derived using on in silico input. In conclusion, these results indicate that not all substances contribute equally (or at all) to the systemic exposure to salicylic acid. Therefore, we propose a realistic metabolism correction factor by which the potential contribution of salicylate esters to the aggregate consumer exposure to salicylic acid from cosmetic use can be estimated.
Collapse
Affiliation(s)
| | | | - Beate Nicol
- Safety and Environmental Assurance Centre, Unilever UK, Colworth Science Park, MK44 1LQ, Sharnbrook, United Kingdom
| | | | | | | | - Amaia Irizar
- The International Fragrance Association (IFRA), Geneva, Switzerland
| | - Brian Wall
- Colgate-Palmolive Company, Piscataway, NJ, 08854, USA
| | - Angus Swinscoe
- Whitman Laboratories, The Estée Lauder Companies, Petersfield, United Kingdom
| | | | - Dan Selechnik
- Research Institute for Fragrance Materials (RIFM), Inc., Woodcliff Lake, NJ, USA
| | - Anne Marie Api
- Research Institute for Fragrance Materials (RIFM), Inc., Woodcliff Lake, NJ, USA
| | | | | | | |
Collapse
|
4
|
Huang H, Zhao W, Qin N, Duan X. Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics. TOXICS 2024; 12:433. [PMID: 38922113 PMCID: PMC11209072 DOI: 10.3390/toxics12060433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Revised: 06/04/2024] [Accepted: 06/13/2024] [Indexed: 06/27/2024]
Abstract
Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.
Collapse
Affiliation(s)
| | | | - Ning Qin
- School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China; (H.H.); (W.Z.)
| | - Xiaoli Duan
- School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China; (H.H.); (W.Z.)
| |
Collapse
|
5
|
Bassani D, Parrott NJ, Manevski N, Zhang JD. Another string to your bow: machine learning prediction of the pharmacokinetic properties of small molecules. Expert Opin Drug Discov 2024; 19:683-698. [PMID: 38727016 DOI: 10.1080/17460441.2024.2348157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Accepted: 04/23/2024] [Indexed: 05/22/2024]
Abstract
INTRODUCTION Prediction of pharmacokinetic (PK) properties is crucial for drug discovery and development. Machine-learning (ML) models, which use statistical pattern recognition to learn correlations between input features (such as chemical structures) and target variables (such as PK parameters), are being increasingly used for this purpose. To embed ML models for PK prediction into workflows and to guide future development, a solid understanding of their applicability, advantages, limitations, and synergies with other approaches is necessary. AREAS COVERED This narrative review discusses the design and application of ML models to predict PK parameters of small molecules, especially in light of established approaches including in vitro-in vivo extrapolation (IVIVE) and physiologically based pharmacokinetic (PBPK) models. The authors illustrate scenarios in which the three approaches are used and emphasize how they enhance and complement each other. In particular, they highlight achievements, the state of the art and potentials of applying machine learning for PK prediction through a comphrehensive literature review. EXPERT OPINION ML models, when carefully crafted, regularly updated, and appropriately used, empower users to prioritize molecules with favorable PK properties. Informed practitioners can leverage these models to improve the efficiency of drug discovery and development process.
Collapse
Affiliation(s)
- Davide Bassani
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Basel, Switzerland
| | - Neil John Parrott
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Basel, Switzerland
| | - Nenad Manevski
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Basel, Switzerland
| | - Jitao David Zhang
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Basel, Switzerland
| |
Collapse
|
6
|
Haid RTU, Reichel A. Transforming the Discovery of Targeted Protein Degraders: The Translational Power of Predictive PK/PD Modeling. Clin Pharmacol Ther 2024. [PMID: 38708948 DOI: 10.1002/cpt.3273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Accepted: 03/21/2024] [Indexed: 05/07/2024]
Abstract
Targeted protein degraders (TPDs), an emerging therapeutic modality, are attracting considerable interest with the promise to address disease-related proteins that are not druggable with conventional small molecule inhibitors. Despite their novel mechanism of action, the PK/PD relationship of degraders is still approached with a mindset deeply rooted in inhibitor drugs. Here, we establish how predictive mechanistic modeling specifically tailored to TPDs can significantly enhance the value of the available information during lead generation and optimization. By integrating the results from in vitro assays with routinely collected PK data, modeling accurately predicts degradation in vivo. These predictions transform the prioritization of compounds for in vivo studies as well as the selection of optimal dose schedules and most informative measurement time points with the least number of animals. Moreover, the comprehensive modeling framework (1) identifies the PK/PD driver of targeted protein degradation and subsequent downstream pharmacodynamic effects, and (2) uncovers the fundamental difference between degrader and inhibitor PK/PD relationships. The practical utility of our predictive modeling is demonstrated with relevant use cases. This framework will allow researchers to transition from current, mostly serendipity-based approaches to more sound model-informed decision making. Going forward, the presented predictive PK/PD modeling framework lays out a rational path to incorporate inter-species differences in the pharmacology and thus promises to help with getting the dose right in clinical trials.
Collapse
Affiliation(s)
- Robin Thomas Ulrich Haid
- Preclinical Modeling & Simulation, Drug Metabolism & Pharmacokinetics, Preclinical Development, Bayer AG, Berlin, Germany
- Biopharmacy, Institute of Pharmaceutical Sciences, Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich, Switzerland
| | - Andreas Reichel
- Preclinical Modeling & Simulation, Drug Metabolism & Pharmacokinetics, Preclinical Development, Bayer AG, Berlin, Germany
| |
Collapse
|
7
|
Jones J, Clark RD, Lawless MS, Miller DW, Waldman M. The AI-driven Drug Design (AIDD) platform: an interactive multi-parameter optimization system integrating molecular evolution with physiologically based pharmacokinetic simulations. J Comput Aided Mol Des 2024; 38:14. [PMID: 38499823 DOI: 10.1007/s10822-024-00552-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 02/13/2024] [Indexed: 03/20/2024]
Abstract
Computer-aided drug design has advanced rapidly in recent years, and multiple instances of in silico designed molecules advancing to the clinic have demonstrated the contribution of this field to medicine. Properly designed and implemented platforms can drastically reduce drug development timelines and costs. While such efforts were initially focused primarily on target affinity/activity, it is now appreciated that other parameters are equally important in the successful development of a drug and its progression to the clinic, including pharmacokinetic properties as well as absorption, distribution, metabolic, excretion and toxicological (ADMET) properties. In the last decade, several programs have been developed that incorporate these properties into the drug design and optimization process and to varying degrees, allowing for multi-parameter optimization. Here, we introduce the Artificial Intelligence-driven Drug Design (AIDD) platform, which automates the drug design process by integrating high-throughput physiologically-based pharmacokinetic simulations (powered by GastroPlus) and ADMET predictions (powered by ADMET Predictor) with an advanced evolutionary algorithm that is quite different than current generative models. AIDD uses these and other estimates in iteratively performing multi-objective optimizations to produce novel molecules that are active and lead-like. Here we describe the AIDD workflow and details of the methodologies involved therein. We use a dataset of triazolopyrimidine inhibitors of the dihydroorotate dehydrogenase from Plasmodium falciparum to illustrate how AIDD generates novel sets of molecules.
Collapse
Affiliation(s)
- Jeremy Jones
- Simulations Plus, Inc., 42505 10th Street West, Lancaster, CA, 93534‑7059, USA.
| | - Robert D Clark
- The Indiana University Luddy School of Informatics, Computing and Engineering, 700 N. Woodlawn Avenue, Bloomington, IN, 47408, USA
| | - Michael S Lawless
- Simulations Plus, Inc., 42505 10th Street West, Lancaster, CA, 93534‑7059, USA
| | - David W Miller
- Simulations Plus, Inc., 42505 10th Street West, Lancaster, CA, 93534‑7059, USA
| | - Marvin Waldman
- Simulations Plus, Inc., 42505 10th Street West, Lancaster, CA, 93534‑7059, USA
| |
Collapse
|
8
|
Führer F, Gruber A, Diedam H, Göller AH, Menz S, Schneckener S. A deep neural network: mechanistic hybrid model to predict pharmacokinetics in rat. J Comput Aided Mol Des 2024; 38:7. [PMID: 38294570 DOI: 10.1007/s10822-023-00547-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Accepted: 12/21/2023] [Indexed: 02/01/2024]
Abstract
An important aspect in the development of small molecules as drugs or agrochemicals is their systemic availability after intravenous and oral administration. The prediction of the systemic availability from the chemical structure of a potential candidate is highly desirable, as it allows to focus the drug or agrochemical development on compounds with a favorable kinetic profile. However, such predictions are challenging as the availability is the result of the complex interplay between molecular properties, biology and physiology and training data is rare. In this work we improve the hybrid model developed earlier (Schneckener in J Chem Inf Model 59:4893-4905, 2019). We reduce the median fold change error for the total oral exposure from 2.85 to 2.35 and for intravenous administration from 1.95 to 1.62. This is achieved by training on a larger data set, improving the neural network architecture as well as the parametrization of mechanistic model. Further, we extend our approach to predict additional endpoints and to handle different covariates, like sex and dosage form. In contrast to a pure machine learning model, our model is able to predict new end points on which it has not been trained. We demonstrate this feature by predicting the exposure over the first 24 h, while the model has only been trained on the total exposure.
Collapse
Affiliation(s)
- Florian Führer
- Engineering & Technology, Applied Mathematics, Bayer AG, 51368, Leverkusen, Germany.
| | - Andrea Gruber
- Pharmaceuticals, R&D, Preclinical Modeling & Simulation, Bayer AG, 13353, Berlin, Germany
| | - Holger Diedam
- Crop Science, Product Supply, SC Simulation & Analysis, Bayer AG, 40789, Monheim, Germany
| | - Andreas H Göller
- Pharmaceuticals, R&D, Molecular Design, Bayer AG, 42096, Wuppertal, Germany
| | - Stephan Menz
- Pharmaceuticals, R&D, Preclinical Modeling & Simulation, Bayer AG, 13353, Berlin, Germany
| | | |
Collapse
|
9
|
Chou WC, Tell LA, Baynes RE, Davis JL, Cheng YH, Maunsell FP, Riviere JE, Lin Z. Development and application of an interactive generic physiologically based pharmacokinetic (igPBPK) model for adult beef cattle and lactating dairy cows to estimate tissue distribution and edible tissue and milk withdrawal intervals for per- and polyfluoroalkyl substances (PFAS). Food Chem Toxicol 2023; 181:114062. [PMID: 37769896 DOI: 10.1016/j.fct.2023.114062] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 09/20/2023] [Accepted: 09/23/2023] [Indexed: 10/03/2023]
Abstract
Humans can be exposed to per- and polyfluoroalkyl substances (PFAS) through dietary intake from milk and edible tissues from food animals. This study developed a physiologically based pharmacokinetic (PBPK) model to predict tissue and milk residues and estimate withdrawal intervals (WDIs) for multiple PFAS including PFOA, PFOS and PFHxS in beef cattle and lactating dairy cows. Results showed that model predictions were mostly within a two-fold factor of experimental data for plasma, tissues, and milk with an estimated coefficient of determination (R2) of >0.95. The predicted muscle WDIs for beef cattle were <1 day for PFOA, 449 days for PFOS, and 69 days for PFHxS, while the predicted milk WDIs in dairy cows were <1 day for PFOA, 1345 days for PFOS, and zero day for PFHxS following a high environmental exposure scenario (e.g., 49.3, 193, and 161 ng/kg/day for PFOA, PFOS, and PFHxS, respectively, for beef cattle for 2 years). The model was converted to a web-based interactive generic PBPK (igPBPK) platform to provide a user-friendly dashboard for predictions of tissue and milk WDIs for PFAS in cattle. This model serves as a foundation for extrapolation to other PFAS compounds to improve safety assessment of cattle-derived food products.
Collapse
Affiliation(s)
- Wei-Chun Chou
- Department of Environmental and Global Health, College of Public Health and Health Professions, University of Florida, Gainesville, FL, 32610, USA; Center for Environmental and Human Toxicology, University of Florida, FL, 32608, USA.
| | - Lisa A Tell
- Department of Medicine and Epidemiology, School of Veterinary Medicine, University of California-Davis, Davis, CA, 95616, USA.
| | - Ronald E Baynes
- Center for Chemical Toxicology Research and Pharmacokinetics, Department of Population Health and Pathobiology, College of Veterinary Medicine, North Carolina State University, Raleigh, NC, 27606, USA.
| | - Jennifer L Davis
- Department of Biomedical Sciences and Pathobiology, Virginia-Maryland College of Veterinary Medicine, Blacksburg, VA, 24060, USA.
| | - Yi-Hsien Cheng
- Institute of Computational Comparative Medicine, Kansas State University, Manhattan, KS, 66506, USA; Department of Anatomy and Physiology, College of Veterinary Medicine, Kansas State University, Manhattan, KS, 66506, USA.
| | - Fiona P Maunsell
- Department of Large Animal Clinical Sciences, College of Veterinary Medicine, University of Florida, Gainesville, FL, 32608, USA.
| | - Jim E Riviere
- Center for Chemical Toxicology Research and Pharmacokinetics, Department of Population Health and Pathobiology, College of Veterinary Medicine, North Carolina State University, Raleigh, NC, 27606, USA; 1Data Consortium, Kansas State University, Olathe, KS, 66061, USA.
| | - Zhoumeng Lin
- Department of Environmental and Global Health, College of Public Health and Health Professions, University of Florida, Gainesville, FL, 32610, USA; Center for Environmental and Human Toxicology, University of Florida, FL, 32608, USA.
| |
Collapse
|
10
|
Andrews-Morger A, Reutlinger M, Parrott N, Olivares-Morales A. A Machine Learning Framework to Improve Rat Clearance Predictions and Inform Physiologically Based Pharmacokinetic Modeling. Mol Pharm 2023; 20:5052-5065. [PMID: 37713584 DOI: 10.1021/acs.molpharmaceut.3c00374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/17/2023]
Abstract
During drug discovery and development, achieving appropriate pharmacokinetics is key to establishment of the efficacy and safety of new drugs. Physiologically based pharmacokinetic (PBPK) models integrating in vitro-to-in vivo extrapolation have become an essential in silico tool to achieve this goal. In this context, the most important and probably most challenging pharmacokinetic parameter to estimate is the clearance. Recent work on high-throughput PBPK modeling during drug discovery has shown that a good estimate of the unbound intrinsic clearance (CLint,u,) is the key factor for useful PBPK application. In this work, three different machine learning-based strategies were explored to predict the rat CLint,u as the input into PBPK. Therefore, in vivo and in vitro data was collected for a total of 2639 proprietary compounds. The strategies were compared to the standard in vitro bottom-up approach. Using the well-stirred liver model to back-calculate in vivo CLint,u from in vivo rat clearance and then training a machine learning model on this CLint,u led to more accurate clearance predictions (absolute average fold error (AAFE) 3.1 in temporal cross-validation) than the bottom-up approach (AAFE 3.6-16, depending on the scaling method) and has the advantage that no experimental in vitro data is needed. However, building a machine learning model on the bias between the back-calculated in vivo CLint,u and the bottom-up scaled in vitro CLint,u also performed well. For example, using unbound hepatocyte scaling, adding the bias prediction improved the AAFE in the temporal cross-validation from 16 for bottom-up to 2.9 together with the bias prediction. Similarly, the log Pearson r2 improved from 0.1 to 0.29. Although it would still require in vitro measurement of CLint,u., using unbound scaling for the bottom-up approach, the need for correction of the fu,inc by fu,p data is circumvented. While the above-described ML models were built on all data points available per approach, it is discussed that evaluation comparison across all approaches could only be performed on a subset because ca. 75% of the molecules had missing or unquantifiable measurements of the fraction unbound in plasma or in vitro unbound intrinsic clearance, or they dropped out due to the blood-flow limitation assumed by the well-stirred model. Advantageously, by predicting CLint,u as the input into PBPK, existing workflows can be reused and the prediction of the in vivo clearance and other PK parameters can be improved.
Collapse
Affiliation(s)
- Andrea Andrews-Morger
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, 4070 Basel, Switzerland
| | - Michael Reutlinger
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, 4070 Basel, Switzerland
| | - Neil Parrott
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, 4070 Basel, Switzerland
| | - Andrés Olivares-Morales
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, 4070 Basel, Switzerland
| |
Collapse
|
11
|
Wyszogrodzka-Gaweł G, Shuklinova O, Lisowski B, Wiśniowska B, Polak S. 3D printing combined with biopredictive dissolution and PBPK/PD modeling optimization and personalization of pharmacotherapy: Are we there yet? Drug Discov Today 2023; 28:103731. [PMID: 37541422 DOI: 10.1016/j.drudis.2023.103731] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 07/25/2023] [Accepted: 07/28/2023] [Indexed: 08/06/2023]
Abstract
Precision medicine requires selecting the appropriate dosage regimen for a patient using the right drug, at the right time. Model-Informed Precision Dosing (MIPD) is a concept suggesting utilization of model-based prediction methods for optimizing the treatment benefit-harm balance, based on individual characteristics of the patient, disease, treatment method, and other factors. Here, we discuss a theoretical workflow comprising several elements, beginning from the physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) models, through 3D printed tablets with the model proposed dose, information range and flow, and the patient themselves. We also describe each of these elements, and the connection between them, highlighting challenges and potential obstacles.
Collapse
Affiliation(s)
- Gabriela Wyszogrodzka-Gaweł
- Department of Social Pharmacy, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland.
| | - Olha Shuklinova
- Chair of Pharmaceutical Technology and Biopharmaceutics, Faculty of Pharmacy. Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
| | - Bartek Lisowski
- Chair of Pharmaceutical Technology and Biopharmaceutics, Faculty of Pharmacy. Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland.
| | - Barbara Wiśniowska
- Chair of Pharmaceutical Technology and Biopharmaceutics, Faculty of Pharmacy. Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland.
| | - Sebastian Polak
- Chair of Pharmaceutical Technology and Biopharmaceutics, Faculty of Pharmacy. Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland.
| |
Collapse
|
12
|
Karnati P, Murthy A, Gundeti M, Ahmed T. Modelling Based Approaches to Support Generic Drug Regulatory Submissions-Practical Considerations and Case Studies. AAPS J 2023; 25:63. [PMID: 37353655 DOI: 10.1208/s12248-023-00831-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 06/03/2023] [Indexed: 06/25/2023] Open
Abstract
Model informed drug development (MiDD) is useful to predict in vivo exposure of drugs during various stages of the drug development process. This approach employs a variety of quantitative tools to assess the risks during the drug development process. One important tool in the MiDD tool kit is the Physiologically Based Pharmacokinetic Modelling (PBPK). This tool is extensively used to reduce the development cost and to accelerate the access of medicines to the patients. In this work, we provide an overview of PBPK modelling approaches in the generic drug development process, with a special emphasis on the bio-waiver applications. We describe herein approaches and common pitfalls while submitting model based justifications as a response to the regulatory deficiencies during the generic drug development process. With some in-house case studies, we have attempted to provide a clear path for PBPK model based justifications for bio-waivers. With this review, the gap between theoretical knowledge and practical application of modelling and simulation tools for generic drug product development could be potentially reduced.
Collapse
Affiliation(s)
- Prajwala Karnati
- Biopharmaceutics Department, Global Clinical Management, Dr. Reddy's Laboratories Ltd., Hyderabad, India
| | - Aditya Murthy
- Biopharmaceutics Department, Global Clinical Management, Dr. Reddy's Laboratories Ltd., Hyderabad, India
| | - Manoj Gundeti
- Biopharmaceutics Department, Global Clinical Management, Dr. Reddy's Laboratories Ltd., Hyderabad, India
| | - Tausif Ahmed
- Biopharmaceutics Department, Global Clinical Management, Dr. Reddy's Laboratories Ltd., Hyderabad, India.
| |
Collapse
|
13
|
Manevski N, Umehara K, Parrott N. Drug Design and Success of Prospective Mouse In Vitro-In Vivo Extrapolation (IVIVE) for Predictions of Plasma Clearance (CL p) from Hepatocyte Intrinsic Clearance (CL int). Mol Pharm 2023. [PMID: 37235687 DOI: 10.1021/acs.molpharmaceut.2c01001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Hepatocyte intrinsic clearance (CLint) and methods of in vitro-in vivo extrapolation (IVIVE) are often used to predict plasma clearance (CLp) in drug discovery. While the prediction success of this approach is dependent on the chemotype, specific molecular properties and drug design features that govern these outcomes are poorly understood. To address this challenge, we investigated the success of prospective mouse CLp IVIVE across 2142 chemically diverse compounds. Dilution scaling, which assumes that the free fraction in hepatocyte incubations (fu,inc) is governed by binding to the 10% of serum in the incubation medium, was used as our default CLp IVIVE approach. Results show that predictions of CLp are better for smaller (molecular weight (MW) < 500 Da), less polar (total polar surface area (TPSA) < 100 Å2, hydrogen bond donor (HBD) ≤1, hydrogen bond acceptor (HBA) ≤ 6), lipophilic (log D > 3), and neutral compounds, with low HBD count playing the key role. If compounds are classified according to their chemical space, predictions were good for compounds resembling central nervous system (CNS) drugs [average absolute fold error (AAFE) of 2.05, average fold error (AFE) of 0.90], moderate for classical druglike compounds (according to Lipinski, Veber, and Ghose guidelines; AAFE of 2.55; AFE of 0.68), and poor for nonclassical "beyond the rule of 5" compounds (AAFE of 3.31; AFE of 0.41). From the perspective of measured druglike properties, predictions of CLp were better for compounds with moderate-to-high hepatocyte CLint (>10 μL/min/106 cells), high passive cellular permeability (Papp > 100 nm/s), and moderate observed CLp (5-50 mL/min/kg). Influences of plasma protein binding (fu,p) and P-glycoprotein (Pgp) apical efflux ratio (AP-ER) were less pronounced. If the extended clearance classification system (ECCS) is applied, predictions were good for class 2 (Papp > 50 nm/s; neutral or basic; AAFE of 2.35; AFE of 0.70) and acceptable for class 1A compounds (AAFE of 2.98; AFE of 0.70). Classes 1B, 3 A/B, and 4 showed poor outcomes (AAFE > 3.80; AFE < 0.60). Functional groups trending toward weaker CLp IVIVE were esters, carbamates, sulfonamides, carboxylic acids, ketones, primary and secondary amines, primary alcohols, oxetanes, and compounds liable to aldehyde oxidase metabolism, likely due to multifactorial reasons. Multivariate analysis showed that multiple properties are relevant, combining together to define the overall success of CLp IVIVE. Our results indicate that the current practice of prospective CLp IVIVE is suitable only for CNS-like compounds and well-behaved classical druglike space (e.g., high permeability or ECCS class 2) without challenging functional groups. Unfortunately, based on existing mouse data, prospective CLp IVIVE for complex and nonclassical chemotypes is poor and hardly better than random guessing. This is likely due to complexities such as extrahepatic metabolism and transporter-mediated disposition which are poorly captured by this methodology. With small-molecule drug discovery increasingly evolving toward nonclassical and complex chemotypes, existing CLp IVIVE methodology will require improvement. While empirical correction factors may bridge the gap in the near future, improved and new in vitro assays, data integration models, and machine learning (ML) methods are increasingly needed to address this challenge and reduce the number of nonclinical pharmacokinetic (PK) studies.
Collapse
Affiliation(s)
- Nenad Manevski
- Roche Pharmaceutical Research and Early Development (pRED), Roche Innovation Center Basel, 4070 Basel, Switzerland
| | - Kenichi Umehara
- Roche Pharmaceutical Research and Early Development (pRED), Roche Innovation Center Basel, 4070 Basel, Switzerland
| | - Neil Parrott
- Roche Pharmaceutical Research and Early Development (pRED), Roche Innovation Center Basel, 4070 Basel, Switzerland
| |
Collapse
|
14
|
Kollipara S, Ahmed T, Praveen S. Physiologically based pharmacokinetic modeling (PBPK) to predict drug-drug interactions for encorafenib. Part II. Prospective predictions in hepatic and renal impaired populations with clinical inhibitors and inducers. Xenobiotica 2023; 53:339-356. [PMID: 37584612 DOI: 10.1080/00498254.2023.2246153] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 08/06/2023] [Accepted: 08/06/2023] [Indexed: 08/17/2023]
Abstract
Encorafenib, a potent BRAF kinase inhibitor gets significantly metabolised by CYP3A4 (83%) and CYP2C19 (16%) and is a substrate for P-glycoprotein (P-gp). Due to significant metabolism by CYP3A4, encorafenib exposure can increase in hepatic and renal impairment and may lead to altered magnitude of drug-drug interactions (DDI). Hence, it is necessary to assess the exposures & DDI's in impaired population.Physiologically based pharmacokinetic modelling (PBPK) was utilised to determine the exposures of encorafenib in hepatic and renal impairment along with altered DDI's. Prospective DDI's were predicted with USFDA recommended clinical CYP3A4, CYP2C19, P-gp inhibitors and CYP3A4 inducers.PBPK models for encorafenib, perpetrators simulated PK parameters within 2-folds error. Encorafenib exposures significantly increased in hepatic as compared to renal impairment because of reduced CYP3A4 levels.Hepatic impairment caused changes in inhibition and induction DDI's, when compared to healthy population. Renal impairment did not cause significant changes in DDIs except for itraconazole. P-gp, CYP2C19 inhibitors did not result in altered DDI's. The DDI's were found to have insignificant correlation with relative exposure increase of perpetrators in case of impairment. Overall, this work signifies use of PBPK modelling for DDI's evaluations in hepatic and renal impairment populations.
Collapse
Affiliation(s)
- Sivacharan Kollipara
- KL College of Pharmacy, Koneru Lakshmaiah Education Foundation, Vaddeswaram, Guntur, Andhra Pradesh, India
| | - Tausif Ahmed
- Biopharmaceutics Group, Global Clinical Management, Dr. Reddy's Laboratories Ltd., Integrated Product Development Organization (IPDO), Hyderabad, Telangana, India
| | - Sivadasu Praveen
- KL College of Pharmacy, Koneru Lakshmaiah Education Foundation, Vaddeswaram, Guntur, Andhra Pradesh, India
| |
Collapse
|
15
|
Kollipara S, Ahmed T, Praveen S. Physiologically based pharmacokinetic modelling to predict drug-drug interactions for encorafenib. Part I. Model building, validation, and prospective predictions with enzyme inhibitors, inducers, and transporter inhibitors. Xenobiotica 2023; 53:366-381. [PMID: 37609899 DOI: 10.1080/00498254.2023.2250856] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 08/17/2023] [Accepted: 08/18/2023] [Indexed: 08/24/2023]
Abstract
Encorafenib, a potent BRAF kinase inhibitor undergoes significant metabolism by CYP3A4 (83%) and CYP2C19 (16%) and also a substrate of P-glycoprotein (P-gp). Because of this, encorafenib possesses potential for enzyme-transporter related interactions. Clinically, its drug-drug interactions (DDIs) with CYP3A4 inhibitors (posaconazole, diltiazem) were reported and hence there is a necessity to study DDIs with multiple enzyme inhibitors, inducers, and P-gp inhibitors.USFDA recommended clinical CYP3A4, CYP2C19, P-gp inhibitors, CYP3A4 inducers were selected and prospective DDIs were simulated using physiologically based pharmacokinetic modelling (PBPK). Impact of dose (50 mg vs. 300 mg) and staggering of administrations (0-10 h) on the DDIs were predicted.PBPK models for encorafenib, perpetrators simulated PK parameters within twofold prediction error. Clinically reported DDIs with posaconazole and diltiazem were successfully predicted.CYP2C19 inhibitors did not result in significant DDI whereas strong CYP3A4 inhibitors resulted in DDI ratio up to 4.5. Combining CYP3A4, CYP2C19 inhibitors yielded DDI equivalent CYP3A4 alone. Strong CYP3A4 inducers yielded DDI ratio up to 0.3 and no impact of P-gp inhibitors on DDIs was observed. The DDIs were not impacted by dose and staggering of administration. Overall, this work indicated significance of PBPK modelling for evaluating clinical DDIs with enzymes, transporters and interplay.
Collapse
Affiliation(s)
- Sivacharan Kollipara
- KL College of Pharmacy, Koneru Lakshmaiah Education Foundation, Vaddeswaram, Guntur, Andhra Pradesh, India
| | - Tausif Ahmed
- Biopharmaceutics Group, Global Clinical Management, Dr. Reddy's Laboratories Ltd., Integrated Product Development Organization (IPDO), Hyderabad, Telangana, India
| | - Sivadasu Praveen
- KL College of Pharmacy, Koneru Lakshmaiah Education Foundation, Vaddeswaram, Guntur, Andhra Pradesh, India
| |
Collapse
|
16
|
Stoyanova R, Katzberger PM, Komissarov L, Khadhraoui A, Sach-Peltason L, Groebke Zbinden K, Schindler T, Manevski N. Computational Predictions of Nonclinical Pharmacokinetics at the Drug Design Stage. J Chem Inf Model 2023; 63:442-458. [PMID: 36595708 DOI: 10.1021/acs.jcim.2c01134] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Although computational predictions of pharmacokinetics (PK) are desirable at the drug design stage, existing approaches are often limited by prediction accuracy and human interpretability. Using a discovery data set of mouse and rat PK studies at Roche (9,685 unique compounds), we performed a proof-of-concept study to predict key PK properties from chemical structure alone, including plasma clearance (CLp), volume of distribution at steady-state (Vss), and oral bioavailability (F). Ten machine learning (ML) models were evaluated, including Single-Task, Multitask, and transfer learning approaches (i.e., pretraining with in vitro data). In addition to prediction accuracy, we emphasized human interpretability of outcomes, especially the quantification of uncertainty, applicability domains, and explanations of predictions in terms of molecular features. Results show that intravenous (IV) PK properties (CLp and Vss) can be predicted with good precision (average absolute fold error, AAFE of 1.96-2.84 depending on data split) and low bias (average fold error, AFE of 0.98-1.36), with AutoGluon, Gaussian Process Regressor (GP), and ChemProp displaying the best performance. Driven by higher complexity of oral PK studies, predictions of F were more challenging, with the best AAFE values of 2.35-2.60 and higher overprediction bias (AFE of 1.45-1.62). Multi-Task approaches and pretraining of ChemProp neural networks with in vitro data showed similar precision to Single-Task models but helped reduce the bias and increase correlations between observations and predictions. A combination of GP-computed prediction variance, molecular clustering, and dimensionality-reduction provided valuable quantitative insights into prediction uncertainty and applicability domains. SHAPley Additive exPlanations (SHAPs) highlighted molecular features contributing to prediction outcomes of Vss, providing explanations that could aid drug design. Combined results show that computational predictions of PK are feasible at the drug design stage, with several ML technologies converging to successfully leverage historical PK data sets. Further studies are needed to unlock the full potential of this approach, especially with respect to data set sizes and quality, transfer learning between in vitro and in vivo data sets, model-independent quantification of uncertainty, and explainability of predictions.
Collapse
Affiliation(s)
- Raya Stoyanova
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Paul Maximilian Katzberger
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Leonid Komissarov
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Aous Khadhraoui
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Lisa Sach-Peltason
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Katrin Groebke Zbinden
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Torsten Schindler
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| | - Nenad Manevski
- Roche Pharmaceutical Research and Early Development, Roche Innovation Center Basel, 4070Basel, Switzerland
| |
Collapse
|
17
|
N-Derivatives of ( Z)-Methyl 3-(4-Oxo-2-thioxothiazolidin-5-ylidene)methyl)-1 H-indole-2-carboxylates as Antimicrobial Agents-In Silico and In Vitro Evaluation. Pharmaceuticals (Basel) 2023; 16:ph16010131. [PMID: 36678628 PMCID: PMC9865890 DOI: 10.3390/ph16010131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Revised: 01/09/2023] [Accepted: 01/10/2023] [Indexed: 01/18/2023] Open
Abstract
Herein, we report the experimental evaluation of the antimicrobial activity of seventeen new (Z)-methyl 3-(4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-1H-indole-2-carboxylate derivatives. All tested compounds exhibited antibacterial activity against eight Gram-positive and Gram-negative bacteria. Their activity exceeded those of ampicillin as well as streptomycin by 10-50 fold. The most sensitive bacterium was En. Cloacae, while E. coli was the most resistant one, followed by M. flavus. The most active compound appeared to be compound 8 with MIC at 0.004-0.03 mg/mL and MBC at 0.008-0.06 mg/mL. The antifungal activity of tested compounds was good to excellent with MIC in the range of 0.004-0.06 mg/mL, with compound 15 being the most potent. T. viride was the most sensitive fungal, while A. fumigatus was the most resistant one. Docking studies revealed that the inhibition of E. coli MurB is probably responsible for their antibacterial activity, while 14a-lanosterol demethylase of CYP51Ca is involved in the mechanism of antifungal activity. Furthermore, drug-likeness and ADMET profile prediction were performed. Finally, the cytotoxicity studies were performed for the most active compounds using MTT assay against normal MRC5 cells.
Collapse
|
18
|
In-Depth Analysis of Physiologically Based Pharmacokinetic (PBPK) Modeling Utilization in Different Application Fields Using Text Mining Tools. Pharmaceutics 2022; 15:pharmaceutics15010107. [PMID: 36678737 PMCID: PMC9860979 DOI: 10.3390/pharmaceutics15010107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2022] [Revised: 12/15/2022] [Accepted: 12/24/2022] [Indexed: 12/30/2022] Open
Abstract
In the past decade, only a small number of papers have elaborated on the application of physiologically based pharmacokinetic (PBPK) modeling across different areas. In this review, an in-depth analysis of the distribution of PBPK modeling in relation to its application in various research topics and model validation was conducted by text mining tools. Orange 3.32.0, an open-source data mining program was used for text mining. PubMed was used for data retrieval, and the collected articles were analyzed by several widgets. A total of 2699 articles related to PBPK modeling met the predefined criteria. The number of publications per year has been rising steadily. Regarding the application areas, the results revealed that 26% of the publications described the use of PBPK modeling in early drug development, risk assessment and toxicity assessment, followed by absorption/formulation modeling (25%), prediction of drug-disease interactions (20%), drug-drug interactions (DDIs) (17%) and pediatric drug development (12%). Furthermore, the analysis showed that only 12% of the publications mentioned model validation, of which 51% referred to literature-based validation and 26% to experimentally validated models. The obtained results present a valuable review of the state-of-the-art regarding PBPK modeling applications in drug discovery and development and related fields.
Collapse
|
19
|
Parrott N, Manevski N, Olivares-Morales A. Can We Predict Clinical Pharmacokinetics of Highly Lipophilic Compounds by Integration of Machine Learning or In Vitro Data into Physiologically Based Models? A Feasibility Study Based on 12 Development Compounds. Mol Pharm 2022; 19:3858-3868. [PMID: 36150125 DOI: 10.1021/acs.molpharmaceut.2c00350] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
While high lipophilicity tends to improve potency, its effects on pharmacokinetics (PK) are complex and often unfavorable. To predict clinical PK in early drug discovery, we built human physiologically based PK (PBPK) models integrating either (i) machine learning (ML)-predicted properties or (ii) discovery stage in vitro data. Our test set was composed of 12 challenging development compounds with high lipophilicity (mean calculated log P 4.2), low plasma-free fraction (50% of compounds with fu,p < 1%), and low aqueous solubility. Predictions focused on key human PK parameters, including plasma clearance (CL), volume of distribution at steady state (Vss), and oral bioavailability (%F). For predictions of CL, the ML inputs showed acceptable accuracy and slight underprediction bias [an average absolute fold error (AAFE) of 3.55; an average fold error (AFE) of 0.95]. Surprisingly, use of measured data only slightly improved accuracy but introduced an overprediction bias (AAFE = 3.35; AFE = 2.63). Predictions of Vss were more successful, with both ML (AAFE = 2.21; AFE = 0.90) and in vitro (AAFE = 2.24; AFE = 1.72) inputs showing good accuracy and moderate bias. The %F was poorly predicted using ML inputs [average absolute prediction error (AAPE) of 45%], and use of measured data for solubility and permeability improved this to 34%. Sensitivity analysis showed that predictions of CL limited the overall accuracy of human PK predictions, partly due to high nonspecific binding of lipophilic compounds, leading to uncertainty of unbound clearance. For accurate predictions of %F, solubility was the key factor. Despite current limitations, this work encourages further development of ML models and integration of their results within PBPK models to enable human PK prediction at the drug design stage, even before compounds are synthesized. Further evaluation of this approach with more diverse chemical types is warranted.
Collapse
Affiliation(s)
- Neil Parrott
- Pharmaceutical Sciences, Pharma Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070 Basel, Switzerland
| | - Nenad Manevski
- Pharmaceutical Sciences, Pharma Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070 Basel, Switzerland
| | - Andrés Olivares-Morales
- Pharmaceutical Sciences, Pharma Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070 Basel, Switzerland
| |
Collapse
|