Model-informed Drug Discovery Research Articles

In-Depth Analysis of the Selection of PBPK Modeling Tools: Bibliometric and Social Network Analysis of the Open Systems Pharmacology Community.

Since the Open Source Initiative laid the foundation for the open source software environment in 1998, the popularity of free and open source software has been steadily increasing. Model-informed drug discovery and development (MID3), a key component of pharmaceutical research and development, heavily makes use of computational models which can be developed using various software including the Open Systems Pharmacology (OSP) software (PK-Sim/MoBi), a free and open source software tool for physiologically based pharmacokinetic (PBPK) modeling. In this study, we aimed to investigate the impact, application areas, and reach of the OSP software as well as the relationships and collaboration patterns between organizations having published OSP-related articles between 2017 and 2023. Therefore, we conducted a bibliometric analysis of OSP-related publications and a social network analysis of the organizations with which authors of OSP-related publications were affiliated. On several levels, we found evidence for a significant growth in the size of the OSP community as well as its visibility in the MID3 community since OSP's establishment in 2017. Specifically, the annual publication rate of PubMed-indexed PBPK-related articles using the OSP software outpaced that of PBPK-related articles using any software. Our bibliometric analysis and network analysis demonstrated that the expansion of the OSP community was predominantly driven by new authors and organizations without prior connections to the community involving the generation of research clusters de novo and an overall diversification of the network. These findings suggest an ongoing evolution of the OSP community toward a more segmented, diverse, and inclusive network.

Application of quantitative protein mass spectrometric data in the early predictive analysis of membrane-bound target engagement by monoclonal antibodies.

Model-informed drug discovery advocates the use of mathematical modeling and simulation for improved efficacy in drug discovery. In the case of monoclonal antibodies (mAbs) against cell membrane antigens, this requires quantitative insight into the target tissue concentration levels. Protein mass spectrometry data are often available but the values are expressed in relative, rather than in molar concentration units that are easier to incorporate into pharmacokinetic models. Here, we present an empirical correlation that converts the parts per million (ppm) concentrations in the PaxDb database to their molar equivalents that are more suitable for pharmacokinetic modeling. We evaluate the insight afforded to target tissue distribution by analyzing the likely tumor-targeting accuracy of mAbs recognizing either epidermal growth factor receptor or its homolog HER2. Surprisingly, the predicted tissue concentrations of both these targets exceed the Kd values of their respective therapeutic mAbs. Physiologically based pharmacokinetic (PBPK) modeling indicates that in these conditions only about 0.05% of the dosed mAb is likely to reach the solid tumor target cells. The rest of the dose is eliminated in healthy tissues via both nonspecific and target-mediated processes. The presented approach allows evaluation of the interplay between the target expression level in different tissues that determines the overall pharmacokinetic properties of the drug and the fraction that reaches the cells of interest. This methodology can help to evaluate the efficacy and safety properties of novel drugs, especially if the off-target cell degradation has cytotoxic outcomes, as in the case of antibody-drug conjugates.

Artificial Intelligence for Quantitative Modeling in Drug Discovery and Development: An Innovation and Quality Consortium Perspective on Use Cases and Best Practices.

Recent breakthroughs in artificial intelligence (AI) and machine learning (ML) have ushered in a new era of possibilities across various scientific domains. One area where these advancements hold significant promise is model-informed drug discovery and development (MID3). To foster a wider adoption and acceptance of these advanced algorithms, the Innovation and Quality (IQ) Consortium initiated the AI/ML working group in 2021 with the aim of promoting their acceptance among the broader scientific community as well as by regulatory agencies. By drawing insights from workshops organized by the working group and attended by key stakeholders across the biopharma industry, academia, and regulatory agencies, this white paper provides a perspective from the IQ Consortium. The range of applications covered in this white paper encompass the following thematic topics: (i) AI/ML-enabled Analytics for Pharmacometrics and Quantitative Systems Pharmacology (QSP) Workflows; (ii) Explainable Artificial Intelligence and its Applications in Disease Progression Modeling; (iii) Natural Language Processing (NLP) in Quantitative Pharmacology Modeling; and (iv) AI/ML Utilization in Drug Discovery. Additionally, the paper offers a set of best practices to ensure an effective and responsible use of AI, including considering the context of use, explainability and generalizability of models, and having human-in-the-loop. We believe that embracing the transformative power of AI in quantitative modeling while adopting a set of good practices can unlock new opportunities for innovation, increase efficiency, and ultimately bring benefits to patients.

Pharmacokinetic analysis across studies to drive knowledge-integration: A tutorial on individual patient data meta-analysis (IPDMA).

Answering challenging questions in drug development sometimes requires pharmacokinetic (PK) data analysis across different studies, for example, to characterize PKs across diverse regions or populations, or to increase statistical power for subpopulations by combining smaller size trials. Given the growing interest in data sharing and advanced computational methods, knowledge integration based on multiple data sources is increasingly applied in the context of model-informed drug discovery and development. A powerful analysis method is the individual patient data meta-analysis (IPDMA), leveraging systematic review of databases and literature, with the most detailed data type of the individual patient, and quantitative modeling of the PK processes, including capturing heterogeneity of variance between studies. The methodology that should be used in IPDMA in the context of population PK analysis is summarized in this tutorial, highlighting areas of special attention compared to standard PK modeling, including hierarchical nested variability terms for interstudy variability, and handling between-assay differences in limits of quantification within a single analysis. This tutorial is intended for any pharmacological modeler who is interested in performing an integrated analysis of PK data across different studies in a systematic and thorough manner, to answer questions that transcend individual primary studies.

Deep Learning Methods Applied to Drug Concentration Prediction of Olanzapine.

Pharmacometrics and the utilization of population pharmacokinetics play an integral role in model-informed drug discovery and development (MIDD). Recently, there has been a growth in the application of deep learning approaches to aid in areas within MIDD. In this study, a deep learning model, LSTM-ANN, was developed to predict olanzapine drug concentrations from the CATIE study. A total of 1527 olanzapine drug concentrations from 523 individuals along with 11 patient-specific covariates were used in model development. The hyperparameters of the LSTM-ANN model were optimized through a Bayesian optimization algorithm. A population pharmacokinetic model using the NONMEM model was constructed as a reference to compare to the performance of the LSTM-ANN model. The RMSE of the LSTM-ANN model was 29.566 in the validation set, while the RMSE of the NONMEM model was 31.129. Permutation importance revealed that age, sex, and smoking were highly influential covariates in the LSTM-ANN model. The LSTM-ANN model showed potential in the application of drug concentration predictions as it was able to capture the relationships within a sparsely sampled pharmacokinetic dataset and perform comparably to the NONMEM model.

Physiologically based pharmacokinetic modeling of tacrolimus for food-drug and CYP3A drug-drug-gene interaction predictions.

The immunosuppressant and narrow therapeutic index drug tacrolimus is metabolized mainly via cytochrome P450 (CYP) 3A4 and CYP3A5. For its pharmacokinetics (PK), high inter- and intra-individual variability can be observed. Underlying causes include the effect of food intake on tacrolimus absorption as well as genetic polymorphism in the CYP3A5 gene. Furthermore, tacrolimus is highly susceptible to drug-drug interactions, acting as a victim drug when coadministered with CYP3A perpetrators. This work describes the development of a whole-body physiologically based pharmacokinetic model for tacrolimus as well as its application for investigation and prediction of (i) the impact of food intake on tacrolimus PK (food-drug interactions [FDIs]) and (ii) drug-drug(-gene) interactions (DD[G]Is) involving the CYP3A perpetrator drugs voriconazole, itraconazole, and rifampicin. The model was built in PK-Sim® Version 10 using a total of 37 whole blood concentration-time profiles of tacrolimus (training and test) compiled from 911 healthy individuals covering the administration of tacrolimus as intravenous infusions as well as immediate-release and extended-release capsules. Metabolism was incorporated via CYP3A4 and CYP3A5, with varying activities implemented for different CYP3A5 genotypes and study populations. The good predictive model performance is demonstrated for the examined food effect studies with 6/6 predicted FDI area under the curve determined between first and last concentration measurements (AUClast ) and 6/6 predicted FDI maximum whole blood concentration (Cmax ) ratios within twofold of the respective observed ratios. In addition, 7/7 predicted DD(G)I AUClast and 6/7 predicted DD(G)I Cmax ratios were within twofold of their observed values. Potential applications of the final model include model-informed drug discovery and development or the support of model-informed precision dosing.

CYTOCON DB: A versatile database of human cell and molecule concentrations for accelerating model development

<scp>CYTOCON DB</scp>: A versatile database of human cell and molecule concentrations for accelerating model development

Application of in silico methods in clinical research and development of drugs and their formulation, a scoping review (Pub Med Central).

Abstract The research was conducted in March 2022, using the databases Latin American and Caribbean Literature on Health Sciences (Lilacs), National Library of Medicine (PubMed), MEDLINE, Web of Science and Scopus, in English. and Meta-Analyses extension for Scoping Reviews (PRSMA-ScR) Checklist (PRISMA, 2020). The criteria adopted for inclusion and exclusion and the specific keywords/descriptors: drug discovery; in silico clinical trial; computational methods; model-informed drug Discovery and development; virtual patients.

Application of in silico methods in clinical research and development of drugs and their formulation, a scoping review (BVS-Medline).

Abstract The research was conducted in March 2022, using the databases Latin American and Caribbean Literature on Health Sciences (Lilacs), National Library of Medicine (PubMed), MEDLINE, Web of Science and Scopus, in English. and Meta-Analyses extension for Scoping Reviews (PRSMA-ScR) Checklist (PRISMA, 2020). The criteria adopted for inclusion and exclusion and the specific keywords/descriptors: drug discovery; in silico clinical trial; computational methods; model-informed drug Discovery and development; virtual patients.

Early Feasibility Assessment: A Method for Accurately Predicting Biotherapeutic Dosing to Inform Early Drug Discovery Decisions.

The application of model-informed drug discovery and development (MID3) approaches in the early stages of drug discovery can help determine feasibility of drugging a target, prioritize between targets, or define optimal drug properties for a target product profile (TPP). However, applying MID3 in early discovery can be challenging due to the lack of pharmacokinetic (PK) and pharmacodynamic (PD) data at this stage. Early Feasibility Assessment (EFA) is the application of mechanistic PKPD models, built from first principles, and parameterized by data that is readily available early in drug discovery to make effective dose predictions. This manuscript demonstrates the ability of EFA to make accurate predictions of clinical effective doses for nine approved biotherapeutics and outlines the potential of extending this approach to novel therapeutics to impact early drug discovery decisions.

Quantitative system pharmacology as a legitimate approach to examine extrapolation strategies used to support pediatric drug development.

Extrapolation strategies from adult data for designing pediatric drug development programs are explored using the quantitative systems pharmacology (QSP) modeling approach, a mechanistic drug and disease modeling framework that can predict clinical response and guide pediatric drug development in general. This innovative model‐informed drug discovery and development approach can leverage adult‐pediatric pharmacology and disease similarity metrics to validate extrapolation assumptions. We describe the QSP model strategy and framework for extrapolation to design pediatric drug development programs by leveraging adult data across a wide range of therapeutic areas and illustrating stage‐gate decisions informed by such an approach.

Application of quantitative protein mass spectrometric data in the early predictive analysis of target engagement by monoclonal antibodies.

Model‐informed drug discovery is endorsed by the US Food and Drug Administration (FDA) to improve the flow of medicines from bench to bedside. In the case of monoclonal antibodies, this necessitates taking into account not only the pharmacokinetic (PK) properties of the drug, but also the tissue distribution, concentration, and turnover of the target to guide dose and affinity selection, as well as serve as a link to downstream pharmacology. Relevant information (e.g., tissue proteomic data from quantitative mass spectrometry), is increasingly available from public domain data repositories, although not necessarily in the form that is directly usable for the purpose of quantitative, predictive, and mechanistic PK/pharmacodynamic (PD) modeling based on molarity or similar frameworks instead. Using secreted plasma protein concentrations measured both by immunochemical methods and mass spectrometry, we addressed this gap and derived an optimized nonlinear empirical function that establishes the correlation between the two data sets and validated the approach taken using a wider data set of all proteins found in plasma. In addition, we present a semimechanistic framework for the plasma half‐life of soluble proteins where clearance is expressed as a nonlinear function of the molecular weight of the protein. Finally, we apply the approach to two established therapeutic antibody targets: complement factor C5 and PCSK9 to demonstrate how the described framework can be applied to predictive PK/PD modeling.

Sagacious epitope selection for vaccines, and both antibody-based therapeutics and diagnostics: tips from virology and oncology.

The target of an antibody plays a significant role in the success of antibody-based therapeutics and diagnostics, and vaccine development. This importance is focused on the target binding site—epitope, where epitope selection as a part of design thinking beyond traditional antigen selection using whole cell or whole protein immunization can positively impact success. With purified recombinant protein production and peptide synthesis to display limited/selected epitopes, intrinsic factors that can affect the functioning of resulting antibodies can be more easily selected for. Many of these factors stem from the location of the epitope that can impact accessibility of the antibody to the epitope at a cellular or molecular level, direct inhibition of target antigen activity, conservation of function despite escape mutations, and even noncompetitive inhibition sites. By incorporating novel computational methods for predicting antigen changes to model-informed drug discovery and development, superior vaccines and antibody-based therapeutics or diagnostics can be easily designed to mitigate failures. With detailed examples, this review highlights the new opportunities, factors, and methods of predicting antigenic changes for consideration in sagacious epitope selection.

Physiologically-based pharmacokinetic modeling to evaluate in vitro-to-in vivo extrapolation for intestinal P-glycoprotein inhibition.

As one of the key components in model‐informed drug discovery and development, physiologically‐based pharmacokinetic (PBPK) modeling linked with in vitro‐to‐in vivo extrapolation (IVIVE) is widely applied to quantitatively predict drug–drug interactions (DDIs) on drug‐metabolizing enzymes and transporters. This study aimed to investigate an IVIVE for intestinal P‐glycoprotein (Pgp, ABCB1)‐mediated DDIs among three Pgp substrates, digoxin, dabigatran etexilate, and quinidine, and two Pgp inhibitors, itraconazole and verapamil, via PBPK modeling. For Pgp substrates, assuming unbound Michaelis‐Menten constant (K m) to be intrinsic, in vitro‐to‐in vivo scaling factors for maximal Pgp‐mediated efflux rate (J max) were optimized based on the clinically observed results without co‐administration of Pgp inhibitors. For Pgp inhibitors, PBPK models utilized the reported in vitro values of Pgp inhibition constants (K i), 1.0 μM for itraconazole and 2.0 μM for verapamil. Overall, the PBPK modeling sufficiently described Pgp‐mediated DDIs between these substrates and inhibitors with the prediction errors of less than or equal to ±25% in most cases, suggesting a reasonable IVIVE for Pgp kinetics in the clinical DDI results. The modeling results also suggest that Pgp kinetic parameters of both the substrates (K m and J max) and the inhibitors (K i) are sensitive to Pgp‐mediated DDIs, thus being key for successful DDI prediction. It would also be critical to incorporate appropriate unbound inhibitor concentrations at the site of action into PBPK models. The present results support a quantitative prediction of Pgp‐mediated DDIs using in vitro parameters, which will significantly increase the value of in vitro studies to design and run clinical DDI studies safely and effectively.

Physiologically Based Pharmacokinetic Modeling of Rosuvastatin to Predict Transporter-Mediated Drug-Drug Interactions

PurposeTo build a physiologically based pharmacokinetic (PBPK) model of the clinical OATP1B1/OATP1B3/BCRP victim drug rosuvastatin for the investigation and prediction of its transporter-mediated drug-drug interactions (DDIs).MethodsThe Rosuvastatin model was developed using the open-source PBPK software PK-Sim®, following a middle-out approach. 42 clinical studies (dosing range 0.002–80.0 mg), providing rosuvastatin plasma, urine and feces data, positron emission tomography (PET) measurements of tissue concentrations and 7 different rosuvastatin DDI studies with rifampicin, gemfibrozil and probenecid as the perpetrator drugs, were included to build and qualify the model.ResultsThe carefully developed and thoroughly evaluated model adequately describes the analyzed clinical data, including blood, liver, feces and urine measurements. The processes implemented to describe the rosuvastatin pharmacokinetics and DDIs are active uptake by OATP2B1, OATP1B1/OATP1B3 and OAT3, active efflux by BCRP and Pgp, metabolism by CYP2C9 and passive glomerular filtration. The available clinical rifampicin, gemfibrozil and probenecid DDI studies were modeled using in vitro inhibition constants without adjustments. The good prediction of DDIs was demonstrated by simulated rosuvastatin plasma profiles, DDI AUClast ratios (AUClast during DDI/AUClast without co-administration) and DDI Cmax ratios (Cmax during DDI/Cmax without co-administration), with all simulated DDI ratios within 1.6-fold of the observed values.ConclusionsA whole-body PBPK model of rosuvastatin was built and qualified for the prediction of rosuvastatin pharmacokinetics and transporter-mediated DDIs. The model is freely available in the Open Systems Pharmacology model repository, to support future investigations of rosuvastatin pharmacokinetics, rosuvastatin therapy and DDI studies during model-informed drug discovery and development (MID3).

Bispecific antibodies: A guide to model informed drug discovery and development

Affinity (KD) optimization of monoclonal antibodies is one of the factors that impacts the stoichiometric binding and the corresponding efficacy of a drug. This impacts the dose and the dosing regimen, making the optimum KD a critical component of drug discovery and development. Its importance is further enhanced for bispecific antibodies, where affinity of the drug needs to be optimized with respect to two targets. Mathematical modeling can have critical impact on lead compound optimization. Here we build on previous work of using mathematical models to facilitate lead compound selection, expanding analysis from two membrane bound targets to soluble targets as well. Our analysis reveals the importance of three factors for lead compound optimization: drug affinity to both targets, target turnover rates, and target distribution throughout the body. We describe a method that leverages this information to help make early stage decisions on whether to optimize affinity, and if so, which arm of the bispecific should be optimized. We apply the proposed approach to a variety of scenarios and illustrate the ability to make improved decisions in each case. We integrate results to develop a bispecific antibody KD optimization guide that can be used to improve resource allocation for lead compound selection, accelerating advancement of better compounds. We conclude with a discussion of possible ways to assess the necessary levels of target engagement for affecting disease as part of an integrative approach for model-informed drug discovery and development.

In Silico Modeling and Simulation to Guide Bioequivalence Testing for Oral Drugs in a Virtual Population.

Model-informed drug discovery and development (MID3) shows great advantages in facilitating drug development. A physiologically based pharmacokinetic model is one of the powerful computational approaches of MID3, and the emerging field of virtual bioequivalence is well recognized to be the future of the physiologically based pharmacokinetic model. Based on the translational link between in vitro, in silico, and in vivo, virtual bioequivalence study can evaluate the similarity and potential difference of pharmacokinetic and clinical performance between test and reference formulations. With the aid of virtual bioequivalence study, the pivotal information of clinical trials can be provided to streamline the development for both new and generic drugs. However, a regulatory framework of virtual bioequivalence study has not reached its full maturity. Therefore, this article aims to present an overview of the current status of bioequivalence study, identify the framework of virtual bioequivalence studies for oral drugs, and also discuss the future opportunities of virtual bioequivalence in supporting the waiver and optimization of in vivo clinical trials.

Application of Machine Learning in Translational Medicine: Current Status and Future Opportunities

The exponential increase in our ability to harness multi-dimensional biological and clinical data from experimental to real-world settings has transformed pharmaceutical research and development in recent years, with increasing applications of artificial intelligence (AI) and machine learning (ML). Patient-centered iterative forward and reverse translation is at the heart of precision medicine discovery and development across the continuum from target validation to optimization of pharmacotherapy. Integration of advanced analytics into the practice of Translational Medicine is now a fundamental enabler to fully exploit information contained in diverse sources of big data sets such as “omics” data, as illustrated by deep characterizations of the genome, transcriptome, proteome, metabolome, microbiome, and exposome. In this commentary, we provide an overview of ML applications in drug discovery and development, aligned with the three strategic pillars of Translational Medicine (target, patient, dose) and offer perspectives on their potential to transform the science and practice of the discipline. Opportunities for integrating ML approaches into the discipline of Pharmacometrics are discussed and will revolutionize the practice of model-informed drug discovery and development. Finally, we posit that joint efforts of Clinical Pharmacology, Bioinformatics, and Biomarker Technology experts are vital in cross-functional team settings to realize the promise of AI/ML-enabled Translational and Precision Medicine.

Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling to Support Waivers of In Vivo Clinical Studies: Current Status, Challenges, and Opportunities.

Physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling has been extensively applied to quantitatively translate invitro data, predict the invivo performance, and ultimately support waivers of invivo clinical studies. In the area of biopharmaceutics and within the context of model-informed drug discovery and development (MID3), there is a rapidly growing interest in applying verified and validated mechanistic PBPK models to waive invivo clinical studies. However, the regulatory acceptance of PBPK analyses for biopharmaceutics and oral drug absorption applications, which is also referred to variously as "PBPK absorption modeling" [Zhang et al. CPT: Pharmacometrics Syst. Pharmacol. 2017, 6, 492], "physiologically based absorption modeling", or "physiologically based biopharmaceutics modeling" (PBBM), remains rather low [Kesisoglou et al. J. Pharm. Sci. 2016, 105, 2723] [Heimbach et al. AAPS J. 2019, 21, 29]. Despite considerable progress in the understanding of gastrointestinal (GI) physiology, invitro biopharmaceutic and in silico tools, PBPK models for oral absorption often suffer from an incomplete understanding of the physiology, overparameterization, and insufficient model validation and/or platform verification, all of which can represent limitations to their translatability and predictive performance. The complex interactions of drug substances and (bioenabling) formulations with the highly dynamic and heterogeneous environment of the GI tract in different age, ethnic, and genetic groups as well as disease states have not been yet fully elucidated, and they deserve further research. Along with advancements in the understanding of GI physiology and refinement of current or development of fully mechanistic in silico tools, we strongly believe that harmonization, interdisciplinary interaction, and enhancement of the translational link between invitro, in silico, and invivo will determine the future of PBBM. This Perspective provides an overview of the current status of PBBM, reflects on challenges and knowledge gaps, and discusses future opportunities around PBPK/PD models for oral absorption of small and large molecules to waive invivo clinical studies.

Model-Informed Drug Discovery and Development in Pulmonary Delivery: Biopharmaceutical Pharmacometric Modeling for Formulation Evaluation of Pulmonary Suspensions.

For respiratory conditions, targeted drug delivery to the lungs could produce higher local concentrations with reduced risk of adverse events compared to systemic administration. Despite the increasing interest in pulmonary delivery, the pharmacokinetics (PK) of drugs following pulmonary administration remains to be elucidated. In this context, the application of modeling and simulation methodologies to characterize PK properties of compounds following pulmonary administration remains a scarcity. Pseudomonas aeruginosa (PA) lung infections are resistant to many of the current antibiotic therapies. Targeted treatments for pulmonary delivery could be particularly beneficial for these local conditions. In this study, we report the application of biopharmaceutical pharmacometrics (BPMX) for the analysis of PK data from three investigational antimicrobial agents following pulmonary administration of a suspension formulation. The observed drug concentration–time profiles in lungs and plasma of the compound series were combined for simultaneous analysis and modeling. The developed model describes the PK data, taking into account formulation properties, and provides a mechanism to predict dissolved drug concentrations in the lungs available for activity. The model was then used to evaluate formulation effects and the impact of variability on total and dissolved drug concentrations in lungs and plasma. The predictions suggest that these therapies for lung delivery should ideally be delivered in a sustained release formulation with high solubility for maximum local exposure in lungs for efficacy, with rapid systemic clearance in plasma for reduced risk of unwanted systemic adverse effects. This work shows the potential benefits of BPMX and the role it can play to support drug discovery and development in pulmonary delivery.