In Silico Prediction of Input Parameters for Simplified Physiologically Based Pharmacokinetic Models for Estimating Plasma, Liver, and Kidney Exposures in Rats after Oral Doses of 246 Disparate Chemicals

Formulation Strategy of BCS-II Drugs by Coupling Mechanistic In-Vitro and Nonclinical In-Vivo Data with PBPK: Fundamentals of Absorption-Dissolution to Parameterization of Modelling and Simulation

BCS class II candidates pose challenges in drug development due to their low solubility and permeability. Researchers have explored various techniques; co-amorphous and solid dispersion are major approaches to enhance in-vitro drug solubility and dissolution. However, in-vivo oral bioavailability remains challenging. Physiologically based pharmacokinetic (PBPK) modeling with a detailed understanding of drug absorption, distribution, metabolism, and excretion (ADME) using a mechanistic approach is emerging. This review summarizes the fundamentals of the PBPK, dissolution-absorption models, parameterization of oral absorption for BCS class II drugs, and provides information about newly emerging artificial intelligence/machine learning (AI/ML) linked PBPK approaches with their advantages, disadvantages, challenges and areas of further exploration. Additionally, the fully integrated workflow for formulation design for investigational new drugs (INDs) and virtual bioequivalence for generic molecules falling under BCS-II are discussed.

Predicting Pharmacokinetics in Rats Using Machine Learning: A Comparative Study Between Empirical, Compartmental, and PBPK-Based Approaches

A successful drug needs to combine several properties including high potency and good pharmacokinetic (PK) properties to sustain efficacious plasma concentration over time. To estimate required doses for preclinical animal efficacy models or for the clinics, in vivo PK studies need to be conducted. Although the prediction of ADME properties of compounds using machine learning (ML) models based on chemical structures is well established in drug discovery, the prediction of complete plasma concentration-time profiles has only recently gained attention. In this study, we systematically compare various approaches that integrate ML models with empiric or mechanistic PK models to predict PK profiles in rats after intravenous administration prior to synthesis. More specifically, we compare a standard noncompartmental analysis (NCA)-based approach (prediction of CL and Vss), a pure ML approach (non-mechanistic PK description), a compartmental modeling approach, and a physiologically based pharmacokinetic (PBPK) approach. Our study based on internal preclinical data shows that the latter three approaches yield PK profile predictions of comparable accuracy across a large data set (evaluated as geometric mean fold errors for each profile of over 1000 small molecules). In summary, we demonstrate the improved ability to prioritize drug candidates with desirable PK properties prior to synthesis with ML predictions.

Prediction of Multi-Pharmacokinetics Property in Multi-Species: Bayesian Neural Network Stacking Model with Uncertainty

Pharmacokinetic (PK) properties of a drug are vital attributes influencing its therapeutic effectiveness, playing an important role in the drug development process. Focusing on the difficult task of predicting PK parameters, we compiled an extensive data set comprising parameters across multiple species. Building upon this groundwork, we introduced the PKStack ensemble model to predict PK parameters across diverse species. PKStack integrates a variety of base models and includes uncertainty in its predictions. We also manually collected PK data from animals as an external test set. We predicted a total of 45 tasks for nine PK parameters in five species, and in general, the prediction accuracy was better for intravenous injections, including parameters such as human Vd (R2 = 0.72, RMSE = 0.31), human CL (R2 = 0.52, RMSE = 0.32), and others. In addition to predictive accuracy, we also considered the interpretability of the results and the definition of the model's application domain. Based on the findings, our model has great potential for practical applications in drug discovery.

A Practical In Silico Method for Predicting Compound Brain Concentration-Time Profiles: Combination of PK Modeling and Machine Learning

Given the aging populations in advanced countries globally, many pharmaceutical companies have focused on developing central nervous system (CNS) drugs. However, due to the blood-brain barrier, drugs do not easily reach the target area in the brain. Although conventional screening methods for drug discovery involve the measurement of (unbound fraction of drug) brain-to-plasma partition coefficients, it is difficult to consider nonequilibrium between plasma and brain compound concentration-time profiles. To truly understand the pharmacokinetics/pharmacodynamics of CNS drugs, compound concentration-time profiles in the brain are necessary; however, such analyses are costly and time-consuming and require a significant number of animals. Therefore, in this study, we attempted to develop an in silico prediction method that does not require a large amount of experimental data by combining modeling and simulation (M&S) with machine learning (ML). First, we constructed a hybrid model linking plasma concentration-time profile to the brain compartment that takes into account the transit time and brain distribution of each compound. Using mouse plasma and brain time experimental values for 103 compounds, we determined the brain kinetic parameters of the hybrid model for each compound; this case was defined as scenario I (a positive control experiment) and included the full brain concentration-time profile data. Next, we built an ML model using chemical structure descriptors as explanatory variables and rate parameters as the target variable, and we then input the predicted values from 5-fold cross-validation (CV) into the hybrid model; this case was defined as scenario II, in which no brain compound concentration-time profile data exist. Finally, for scenario III, assuming that the brain concentration is obtained at only one time point, we used the brain kinetic parameters from the result of the 5-fold CV in scenario II as the initial values for the hybrid model and performed parameter refitting against the observed brain concentration at that time point. As a result, the RMSE/R2-values of the brain compound concentration-time profiles over time were 0.445/0.517 in scenario II and 0.246/0.805 in scenario III, indicating the method provides high accuracy and suggesting that it is a practical method for predicting brain compound concentration-time profiles.

Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans

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.

The Trends and Future Prospective of In Silico Models from the Viewpoint of ADME Evaluation in Drug Discovery

Drug discovery and development are aimed at identifying new chemical molecular entities (NCEs) with desirable pharmacokinetic profiles for high therapeutic efficacy. The plasma concentrations of NCEs are a biomarker of their efficacy and are governed by pharmacokinetic processes such as absorption, distribution, metabolism, and excretion (ADME). Poor ADME properties of NCEs are a major cause of attrition in drug development. ADME screening is used to identify and optimize lead compounds in the drug discovery process. Computational models predicting ADME properties have been developed with evolving model-building technologies from a simplified relationship between ADME endpoints and physicochemical properties to machine learning, including support vector machines, random forests, and convolution neural networks. Recently, in the field of in silico ADME research, there has been a shift toward evaluating the in vivo parameters or plasma concentrations of NCEs instead of using predictive results to guide chemical structure design. Another research hotspot is the establishment of a computational prediction platform to strengthen academic drug discovery. Bioinformatics projects have produced a series of in silico ADME models using free software and open-access databases. In this review, we introduce prediction models for various ADME parameters and discuss the currently available academic drug discovery platforms.

Ensemble Machine Learning Approaches Based on Molecular Descriptors and Graph Convolutional Networks for Predicting the Efflux Activities of MDR1 and BCRP Transporters

Multidrug resistance (MDR1) and breast cancer resistance protein (BCRP) play important roles in drug absorption and distribution. Computational prediction of substrates for both transporters can help reduce time in drug discovery. This study aimed to predict the efflux activity of MDR1 and BCRP using multiple machine learning approaches with molecular descriptors and graph convolutional networks (GCNs). In vitro efflux activity was determined using MDR1- and BCRP-expressing cells. Predictive performance was assessed using an in-house dataset with a chronological split and an external dataset. CatBoost and support vector regression showed the best predictive performance for MDR1 and BCRP efflux activities, respectively, of the 25 descriptor-based machine learning methods based on the coefficient of determination (R2). The single-task GCN showed a slightly lower performance than descriptor-based prediction in the in-house dataset. In both approaches, the percentage of compounds predicted within twofold of the observed values in the external dataset was lower than that in the in-house dataset. Multi-task GCN did not show any improvements, whereas multimodal GCN increased the predictive performance of BCRP efflux activity compared with single-task GCN. Furthermore, the ensemble approach of descriptor-based machine learning and GCN achieved the highest predictive performance with R2 values of 0.706 and 0.587 in MDR1 and BCRP, respectively, in time-split test sets. This result suggests that two different approaches to represent molecular structures complement each other in terms of molecular characteristics. Our study demonstrated that predictive models using advanced machine learning approaches are beneficial for identifying potential substrate liability of both MDR1 and BCRP.

The Role of "Physiologically Based Pharmacokinetic Model (PBPK)" New Approach Methodology (NAM) in Pharmaceuticals and Environmental Chemical Risk Assessment

Physiologically Based Pharmacokinetic (PBPK) models are mechanistic tools generally employed in the pharmaceutical industry and environmental health risk assessment. These models are recognized by regulatory authorities for predicting organ concentration-time profiles, pharmacokinetics and daily intake dose of xenobiotics. The extension of PBPK models to capture sensitive populations such as pediatric, geriatric, pregnant females, fetus, etc., and diseased populations such as those with renal impairment, liver cirrhosis, etc., is a must. However, the current modelling practices and existing models are not mature enough to confidently predict the risk in these populations. A multidisciplinary collaboration between clinicians, experimental and modeler scientist is vital to improve the physiology and calculation of biochemical parameters for integrating knowledge and refining existing PBPK models. Specific PBPK covering compartments such as cerebrospinal fluid and the hippocampus are required to gain mechanistic understanding about xenobiotic disposition in these sub-parts. The PBPK model assists in building quantitative adverse outcome pathways (qAOPs) for several endpoints such as developmental neurotoxicity (DNT), hepatotoxicity and cardiotoxicity. Machine learning algorithms can predict physicochemical parameters required to develop in silico models where experimental data are unavailable. Integrating machine learning with PBPK carries the potential to revolutionize the field of drug discovery and development and environmental risk. Overall, this review tried to summarize the recent developments in the in-silico models, building of qAOPs and use of machine learning for improving existing models, along with a regulatory perspective. This review can act as a guide for toxicologists who wish to build their careers in kinetic modeling.

Machine learning and artificial intelligence in physiologically based pharmacokinetic modeling

Physiologically based pharmacokinetic (PBPK) models are useful tools in drug development and risk assessment of environmental chemicals. PBPK model development requires the collection of species-specific physiological, and chemical-specific absorption, distribution, metabolism, and excretion (ADME) parameters, which can be a time-consuming and expensive process. This raises a need to create computational models capable of predicting input parameter values for PBPK models, especially for new compounds. In this review, we summarize an emerging paradigm for integrating PBPK modeling with machine learning (ML) or artificial intelligence (AI)-based computational methods. This paradigm includes 3 steps (1) obtain time-concentration PK data and/or ADME parameters from publicly available databases, (2) develop ML/AI-based approaches to predict ADME parameters, and (3) incorporate the ML/AI models into PBPK models to predict PK summary statistics (eg, area under the curve and maximum plasma concentration). We also discuss a neural network architecture "neural ordinary differential equation (Neural-ODE)" that is capable of providing better predictive capabilities than other ML methods when used to directly predict time-series PK profiles. In order to support applications of ML/AI methods for PBPK model development, several challenges should be addressed (1) as more data become available, it is important to expand the training set by including the structural diversity of compounds to improve the prediction accuracy of ML/AI models; (2) due to the black box nature of many ML models, lack of sufficient interpretability is a limitation; (3) Neural-ODE has great potential to be used to generate time-series PK profiles for new compounds with limited ADME information, but its application remains to be explored. Despite existing challenges, ML/AI approaches will continue to facilitate the efficient development of robust PBPK models for a large number of chemicals.

Computational Predictions of Nonclinical Pharmacokinetics at the Drug Design Stage

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.

Updated in Silico Prediction Methods for Fractions Absorbed and Key Input Parameters of 355 Disparate Chemicals for Physiologically Based Pharmacokinetic Models for Time-Dependent Plasma Concentrations after Virtual Oral Doses in Humans

Human metabolic profiles for substances such as toxic food-derived compounds are usually allometrically extrapolated from traditionally determined in vivo rat concentration profiles. To evaluate internal exposures in humans without any reference to experimental data, physiologically based pharmacokinetic (PBPK) modeling could be used if the model input parameters could be estimated in silico. This approach would simplify the use of PBPK models for forward dosimetry after oral doses. In this study, the in silico estimation of input parameters for PBPK models (i.e., fraction absorbed × intestinal availability, absorption rate constants, and volumes of the systemic circulation) was updated for an panel of 355 chemicals (212 previously analyzed and 143 additional substances) using a light gradient boosting machine learning algorithms (LightGBM) based on between 11 and 29 in silico-calculated chemical descriptors. Simplified human PBPK models were then used to calculate virtual maximum plasma concentrations (C_max) and areas under the concentration-time curve (AUC) based on two sets of input parameters, i.e., traditionally derived values from in vivo data and those calculated in silico using the current updated systems. Both sets of C_max and AUC data were well correlated (r = 0.87 and r = 0.73, respectively; p < 0.01, n = 355). Therefore, input parameters for human PBPK models for a diverse range of compounds could be successfully estimated using chemical descriptors and in silico tools. This approach to pharmacokinetic modeling has potential for application in computational toxicology and in the clinical setting for assessing the potential risk of general chemicals.

Advances in computational methods along the exposure to toxicological response paradigm

Human health risk assessment is a function of chemical toxicity, bioavailability to reach target biological tissues, and potential environmental exposure. These factors are complicated by many physiological, biochemical, physical and lifestyle factors. Furthermore, chemical health risk assessment is challenging in view of the large, and continually increasing, number of chemicals found in the environment. These challenges highlight the need to prioritize resources for the efficient and timely assessment of those environmental chemicals that pose greatest health risks. Computational methods, either predictive or investigative, are designed to assist in this prioritization in view of the lack of cost prohibitive in vivo experimental data. Computational methods provide specific and focused toxicity information using in vitro high throughput screening (HTS) assays. Information from the HTS assays can be converted to in vivo estimates of chemical levels in blood or target tissue, which in turn are converted to in vivo dose estimates that can be compared to exposure levels of the screened chemicals. This manuscript provides a review for the landscape of computational methods developed and used at the U.S. Environmental Protection Agency (EPA) highlighting their potentials and challenges.

Machine learning and artificial intelligence in toxicological sciences

Machine learning and artificial intelligence approaches have revolutionized multiple disciplines, including toxicology. This review summarizes representative recent applications of machine learning and artificial intelligence approaches in different areas of toxicology, including physiologically based pharmacokinetic (PBPK) modeling, quantitative structure-activity relationship modeling for toxicity prediction, adverse outcome pathway analysis, high-throughput screening, toxicogenomics, big data and toxicological databases. By leveraging machine learning and artificial intelligence approaches, now it is possible to develop PBPK models for hundreds of chemicals efficiently, to create in silico models to predict toxicity for a large number of chemicals with similar accuracies compared to in vivo animal experiments, and to analyze a large amount of different types of data (toxicogenomics, high-content image data, etc.) to generate new insights into toxicity mechanisms rapidly, which was impossible by manual approaches in the past. To continue advancing the field of toxicological sciences, several challenges should be considered: (1) not all machine learning models are equally useful for a particular type of toxicology data, and thus it is important to test different methods to determine the optimal approach; (2) current toxicity prediction is mainly on bioactivity classification (yes/no), so additional studies are needed to predict the intensity of effect or dose-response relationship; (3) as more data become available, it is crucial to perform rigorous data quality check and develop infrastructure to store, share, analyze, evaluate, and manage big data; and (4) it is important to convert machine learning models to user-friendly interfaces to facilitate their applications by both computational and bench scientists.

Low cerebrospinal fluid-to-plasma ratios of orally administered lenalidomide mediated by its low cell membrane permeability in patients with hematologic malignancies

Lenalidomide is a synthetic analog of thalidomide formed by the removal of one keto group (plus the addition of an amino group); it has anti-tumor activities beneficial for the treatment of hematologic malignancies. However, lenalidomide distribution to brain in animal models is reportedly low compared with that of thalidomide. The aim of this study was to evaluate plasma and cerebrospinal fluid concentrations of lenalidomide in three patients with malignant hematologic malignancies. Lenalidomide was detected in plasma from the three Japanese patients 1.5 h following oral administration of 20 mg lenalidomide using liquid chromatography/mass spectrometry, despite the in vitro gastrointestinal permeability of lenalidomide being low. Clinically observed cerebrospinal fluid-to-plasma ratios of lenalidomide were low (1.3-2.4%). Observed influx permeability values for lenalidomide in monkey blood-brain barrier model and human placental cell systems were one order of magnitude lower than those of thalidomide and another second-generation drug, pomalidomide along with a positive permeability control, caffeine. Because of the low cell-barrier permeability of lenalidomide demonstrated in in vitro assays, clinically relevant pharmacokinetic profiles of lenalidomide resulted in low penetrability from plasma into cerebrospinal fluid in patients with hematologic malignancies. Lenalidomide is conclusively suggested to expert its favorable immunomodulatory effects via systemic exposures in the patients.

Pharmacokinetics of primary metabolites 5-hydroxythalidomide and 5'-hydroxythalidomide formed after oral administration of thalidomide in the rabbit, a thalidomide-sensitive species

The teratogenicity of the chemotherapeutic drug thalidomide is species-specific and affects humans, non-human primates, and rabbits. The primary oxidation of thalidomide in previously investigated rodents predominantly resulted in the formation of deactivated 5'-hydroxythalidomide. In the current study, similar in vivo biotransformations to 5-hydroxythalidomide and 5'-hydroxythalidomide were confirmed by the analysis of blood plasma from male rabbits, a thalidomide-sensitive species, after oral administration of thalidomide (2.0 mg/kg). Similar levels of thalidomide in seminal plasma and in blood plasma were detected using liquid chromatography-tandem mass spectrometry at 4 hr and 7 hr after oral doses in male rabbits. Seminal plasma concentrations of 5-hydroxythalidomide and 5'-hydroxythalidomide were also seen in male rabbits in a roughly similar time-dependent manner to those in the blood plasma after oral doses of thalidomide (2.0 mg/kg). Furthermore, the values generated by a simplified physiologically based pharmacokinetic rabbit model were in agreement with the measured in vivo blood plasma data under metabolic ratios of 0.01 for the hepatic intrinsic clearance of thalidomide to both unconjugated 5-hydroxythalidomide and 5'-hydroxythalidomide. These results suggest that metabolic activation of thalidomide may be dependent on rabbit liver enzymes just it was for cytochrome P450 enzymes in humanized-liver mice; in contrast, rodent livers predominantly mediate biotransformation of thalidomide to 5'-hydroxythalidomide. A developmental toxicity test system with experimental animals that involves intravaginal exposures to the chemotherapeutic drug thalidomide via semen should be considered in the future.

Machine Learning Prediction of the Three Main Input Parameters of a Simplified Physiologically Based Pharmacokinetic Model Subsequently Used to Generate Time-Dependent Plasma Concentration Data in Humans after Oral Doses of 212 Disparate Chemicals

Physiologically based pharmacokinetic (PBPK) modeling has the potential to play significant roles in estimating internal chemical exposures. The three major PBPK model input parameters (i.e., absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearances) were generated in silico for 212 chemicals using machine learning algorithms. These input parameters were calculated based on sets of between 17 and 65 chemical properties that were generated by in silico prediction tools before being processed by machine learning algorithms. The resulting simplified PBPK models were used to estimate plasma concentrations after virtual oral administrations in humans. The estimated absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearance values for the 212 test compounds determined traditionally (i.e., based on fitting to measured concentration profiles) and newly estimated had correlation coefficients of 0.65, 0.68, and 0.77 (p < 0.01, n = 212), respectively. When human plasma concentrations were modeled using traditionally determined input parameters and again using in silico estimated input parameters, the two sets of maximum plasma concentrations (r = 0.85, p < 0.01, n = 212) and areas under the curve (r = 0.80, p < 0.01, n = 212) were correlated. Virtual chemical exposure levels in liver and kidney were also estimated using these simplified PBPK models along with human plasma levels. These results indicate that the PBPK model input parameters for humans of a diverse set of compounds can be reliability estimated using chemical descriptors calculated using in silico tools.

An Updated In Silico Prediction Method for Volumes of Systemic Circulation of 323 Disparate Chemicals for Use in Physiologically Based Pharmacokinetic Models to Estimate Plasma and Tissue Concentrations after Oral Doses in Rats

Updated algorithms for predicting the volumes of systemic circulation (V₁), along with absorption rate constants and hepatic intrinsic clearances, as input parameters for physiologically based pharmacokinetic (PBPK) models were established to improve the accuracy of estimated plasma and tissue concentrations of 323 chemicals after virtual oral administrations in rats. Using ridge regression with an enlarged set of chemical descriptors (up to 99), the estimated input V₁ values resulted in an improved correlation coefficient (from 246 compounds) with the traditionally determined values. The PBPK model input parameters for rats of diverse compounds can be precisely estimated by increasing the number of descriptors.

Metabolic profiles for the pyrrolizidine alkaloid neopetasitenine and its metabolite petasitenine in humans extrapolated from rat in vivo and in vitro data sets using a simplified physiologically based pharmacokinetic model

Naturally occurring food substances may constitute safety hazards. The risks associated with plant-derived pyrrolizidine alkaloids have been extensively evaluated. Petasites japonicus (common Japanese name, fuki) is a widely consumed water-soluble pyrrolizidine alkaloid-producing plant. In this study, neopetasitenine (acetylfukinotoxin) was selected as a model food substrate (for which human pharmacokinetics were estimated) because of its high concentration in fuki, along with petasitenine (fukinotoxin), its carcinogenic deacetylated metabolite. Although neopetasitenine was rapidly absorbed and converted to petasitenine after oral administration of 1.0 mg/kg in rats, petasitenine was slowly cleared from plasma. Forward dosimetry was conducted using in silico simplified physiologically based pharmacokinetic (PBPK) modeling formulated on experimental pharmacokinetic rat data. From ~2 hr after the oral administration of neopetasitenine in rats, the plasma concentrations of petasitenine were higher than those of neopetasitenine under the present conditions. A human PBPK model was established following an allometric scaling approach applied to rat parameters (without considering interspecies factors) to estimate human intrinsic hepatic clearances from empirical rat values. Human in silico neopetasitenine and petasitenine plasma concentration curves were simulated after daily oral administrations of 3.0 and 1.3 mg/kg neopetasitenine. These doses were taken from reported acute/short-term cases of pyrrolizidine alkaloid toxicity. In vitro hepatotoxicity of neopetasitenine and petasitenine was caused by their high concentrations in the medium for human hepatocyte-like cell line HepaRG cells as an index of lactate dehydrogenase leakage. Neopetasitenine was estimated to be rapidly absorbed and converted to deacetylated carcinogenic petasitenine, even after hepatotoxic doses of 1.0 mg/kg in humans. If the water-soluble pyrrolizidine alkaloid-producing plant P. japonicus were daily consumed as food, current simulation results suggest that dangerous amounts of deacetylated petasitenine could be continuously present in human plasma.

Prediction of permeability across intestinal cell monolayers for 219 disparate chemicals using in vitro experimental coefficients in a pH gradient system and in silico analyses by trivariate linear regressions and machine learning

For medicines, the apparent membrane permeability coefficients (P_app) across human colorectal carcinoma cell line (Caco-2) monolayers under a pH gradient generally correlate with the fraction absorbed after oral intake. Furthermore, the in vitro P_app values of 29 industrial chemicals were found to have an inverse association with their reported no-observed effect levels for hepatotoxicity in rats. In the current study, we expanded our influx permeability predictions for the 90 previously investigated chemicals to both influx and efflux permeability predictions for 207 diverse primary compounds, along with those for 23 secondary compounds. Trivariate linear regression analysis found that the observed influx and efflux logP_app values determined by in vitro experiments significantly correlated with molecular weights and the octanol-water distribution coefficients at apical and basal pH levels (pH 6.0 and 7.4, respectively) (apical to basal, r = 0.76, n = 198; and basal to apical, r = 0.77, n = 202); the distribution coefficients were estimated in silico. Further, prediction accuracy was enhanced by applying a light gradient boosting machine learning system (LightGBM) to estimate influx and efflux logP_app values that incorporated 17 and 19 in silico chemical descriptors (r = 0.83-0.84, p < 0.001). The determination in vitro and/or prediction in silico of permeability coefficients across intestinal cell monolayers of a diverse range of industrial chemicals/food components/medicines could contribute to the safety evaluations of oral intakes of general chemicals in humans. Such new alternative methods could also reduce the need for animal testing during toxicity assessment.

Hepatotoxicological potential of P-toluic acid in humanised-liver mice investigated using simplified physiologically based pharmacokinetic models

p-Toluic acid, a metabolite of organic solvent xylene, has a high reported no-observed-effect level (NOEL, 1000 mg/kg) in rats, possibly because of direct glycine conjugation to methylhippuric acid. In this study, plasma levels of p-toluic acid and its glycine conjugate in mice and humanised-liver mice were evaluated after oral administrations.Although rapid conversion of p-toluic acid to its glycine conjugate was evident from mouse plasma concentrations, the biotransformation of p-toluic acid was slower in humanised-liver mice. The input parameters for physiologically based pharmacokinetic (PBPK) models were determined using fitting procedures to create PBPK-generated plasma concentration curves.The PBPK-modelled hepatic concentrations of p-toluic acid in humanised-liver mice were higher than those observed in plasma. PBPK-modelled hepatic and plasma concentrations of p-toluic acid also indicated slow elimination in humans.These results suggest that rapid conjugations of p-toluic acid reportedly observed in rats could result in overestimation of NOELs for conjugatable chemicals when extrapolated to humanised-liver mice or humans.

Pharmacokinetics of primary oxidative metabolites of thalidomide in rats and in chimeric mice humanized with different human hepatocytes

The approved drug thalidomide is teratogenic in humans, nonhuman primates, and rabbits but not in rodents. The extensive biotransformation of 5'-hydroxythalidomide after oral administration of thalidomide (250 mg/kg) in rats was investigated in detail using liquid chromatography-tandem mass spectrometry. Probable metabolites 5'-hydroxythalidomide sulfate and glucuronide were extensively formed, with approximately tenfold and onefold peak areas, respectively, to the primary 5'-hydroxythalidomide measured using authentic standards. As a minor metabolite, 5-hydroxythalidomide was also detected. The output of simplified physiologically based pharmacokinetic rat models was consistent with the observed in vivo data under a metabolic ratio of 0.05 for the hepatic intrinsic clearance of thalidomide to unconjugated 5'-hydroxythalidomide. The aggregate of unconjugated and sulfate/glucuronide conjugated 5'-hydroxythalidomide forms appear to be the predominant metabolites in rats. Two hours after oral administration of thalidomide (100 mg/kg) to chimeric mice humanized with four different batches of genotyped human hepatocytes, the plasma concentration ratios of 5-hydroxythalidomide to 5'-hydroxythalidomide were correlated with replacement indexes of human liver cells previously transplanted in immunodeficient mice. These results indicate that rodent livers mediate thalidomide primary oxidation, leading to extensive deactivation in vivo to unconjugated/conjugated 5'-hydroxythalidomide and suggest that thalidomide activation might be dependent on the humanized livers in mice transplanted with human hepatocytes.