Styrene lung cancer risk assessment: an alternative evaluation of human lung cancer risk assuming mouse lung tumors are potentially human relevant and operating by a threshold-based non-genotoxic mode of action.

Rodent inhalation studies indicate styrene is a mouse lung-specific carcinogen. Mode-of-action (MOA) analyses indicate that the lung tumors cannot be excluded as weakly quantitatively relevant to humans due to shared oxidative metabolites detected in rodents and humans. However, styrene also is not genotoxic following in vivo dosing. The objective of this review was to characterize occupational and general population cancer risks by conservatively assuming mouse lung tumors were relevant to humans but operating by a non-genotoxic MOA. Inhalation cancer values reference concentrations for respective occupational and general population exposures (RfCcar-occup and RfCcar-genpop) were derived from initial benchmark dose (BMD) modeling of mouse inhalation tumor dose-response data. An overall lowest BMDL10 of 4.7 ppm was modeled for lung tumors, which was further duration- and dose-adjusted by physiologically based pharmacokinetic (PBPK) modeling to derive RfCcar-occup/genpop values of 6.2 ppm and 0.8 ppm, respectively. With the exception of open-mold fiber reinforced composite workers not using personal protective equipment (PPE), the RfCcar-occup/genpop values are greater than typical occupational and general population human exposures, thus indicating styrene exposures represent a low potential for human lung cancer risk. Consistent with this conclusion, a review of styrene occupational epidemiology did not support a conclusion of an association between styrene exposure and lung cancer occurrence, and further supports a conclusion that the conservatively derived RfCcar-occup is lung cancer protective.

Translational two-pore PBPK model to characterize whole-body disposition of different-size endogenous and exogenous proteins.

Two-pore physiologically based pharmacokinetic (PBPK) modeling has demonstrated its potential in describing the pharmacokinetics (PK) of different-size proteins. However, all existing two-pore models lack either diverse proteins for validation or interspecies extrapolation. To fill the gap, here we have developed and optimized a translational two-pore PBPK model that can characterize plasma and tissue disposition of different-size proteins in mice, rats, monkeys, and humans. Datasets used for model development include more than 15 types of proteins: IgG (150 kDa), F(ab)2 (100 kDa), minibody (80 kDa), Fc-containing proteins (205, 200, 110, 105, 92, 84, 81, 65, or 60 kDa), albumin conjugate (85.7 kDa), albumin (67 kDa), Fab (50 kDa), diabody (50 kDa), scFv (27 kDa), dAb2 (23.5 kDa), proteins with an albumin-binding domain (26, 23.5, 22, 16, 14, or 13 kDa), nanobody (13 kDa), and other proteins (110, 65, or 60 kDa). The PBPK model incorporates: (i) molecular weight (MW)-dependent extravasation through large and small pores via diffusion and filtration, (ii) MW-dependent renal filtration, (iii) endosomal FcRn-mediated protection from catabolism for IgG and albumin-related modalities, and (iv) competition for FcRn binding from endogenous IgG and albumin. The finalized model can well characterize PK of most of these proteins, with area under the curve predicted within two-fold error. The model also provides insights into contribution of renal filtration and lysosomal degradation towards total elimination of proteins, and contribution of paracellular convection/diffusion and transcytosis towards extravasation. The PBPK model presented here represents a cross-modality, cross-species platform that can be used for development of novel biologics.

Translating dosimetry of Dibenzo[def,p]chrysene (DBC) and metabolites across dose and species using physiologically based pharmacokinetic (PBPK) modeling.

Dibenzo[def,p]chrysene (DBC) is an environmental polycyclic aromatic hydrocarbon (PAH) that causes tumors in mice and has been classified as a probable human carcinogen by the International Agency for Research on Cancer. Animal toxicity studies often utilize higher doses than are found in relevant human exposures. Additionally, like many PAHs, DBC requires metabolic bioactivation to form the ultimate toxicant, and species differences in DBC and DBC metabolite metabolism have been observed. To understand the implications of dose and species differences, a physiologically based pharmacokinetic model (PBPK) for DBC and major metabolites was developed in mice and humans. Metabolism parameters used in the model were obtained from experimental in vitro metabolism assays using mice and human hepatic microsomes. PBPK model simulations were evaluated against mice dosed with 15 mg/kg DBC by oral gavage and human volunteers orally microdosed with 29 ng of DBC. DBC and its primary metabolite DBC-11,12-diol were measured in blood of mice and humans, while in urine, the majority of DBC metabolites were obeserved as conjugated DBC-11,12-diol, conjugated DBC tetrols, and unconjugated DBC tetrols. The PBPK model was able to predict the time course concentrations of DBC, DBC-11,12-diol, and other DBC metabolites in blood and urine of human volunteers and mice with reasonable accuracy. Agreement between model simulations and measured pharmacokinetic data in mice and human studies demonstrate the success and versatility of our model for interspecies extrapolation and applicability for different doses. Furthermore, our simulations show that internal dose metrics used for risk assessment do not necessarily scale allometrically, and that PBPK modeling provides a reliable approach to appropriately account for interspecies differences in metabolism and physiology.

Consideration of Fractional Distribution Parameter fd in the Chen and Gross Method for Tissue-to-Plasma Partition Coefficients: Comparison of Several Methods.

PURPOSE: The tissue-to-plasma partition coefficient (Kp) describes the extent of tissue distribution in physiologically-based pharmacokinetic (PBPK) models. Constant-rate infusion studies are common for experimental determination of the steady-state Kp,ss, while the tissue-plasma concentration ratio (CT/Cp) in the terminal phase after intravenous doses is often utilized. The Chen and Gross (C&G) method converts a terminal slope CT/Cp to Kp,ss based on assumptions of perfusion-limited distribution in tissue-plasma equilibration. However, considering blood flow (QT) and apparent tissue permeability (fupPSin) in the rate of tissue distribution, this report extends the C&G method by utilizing a fractional distribution parameter (fd). METHODS: Relevant PBPK equations for non-eliminating and eliminating organs along with lung and liver were derived for the conversion of CT/Cp values to Kp,ss. The relationships were demonstrated in rats with measured CT/Cp and Kp,ss values and the model-dependent fd for 8 compounds with a range of permeability coefficients. Several methods of assessing Kp were compared. RESULTS: Utilizing fd in an extended C&G method, our estimations of Kp,ss from CT/Cp were improved, particularly for lower permeability compounds. However, four in silico methods for estimating Kp performed poorly across tissues in comparison with measured Kp values. Mathematical relationships between Kp and Kp,ss that are generally applicable for eliminating organs with tissue permeability limitations necessitates inclusion of an extraction ratio (ER) and fd. CONCLUSION: Since many different types/sources of Kp are present in the literature and used in PBPK models, these perspectives and equations should provide better insights in measuring and interpreting Kp values in PBPK.

Integration of Toxicogenomics and Physiologically Based Pharmacokinetic Modeling in Human Health Risk Assessment of Perfluorooctane Sulfonate.

Toxicogenomics and physiologically based pharmacokinetic (PBPK) models are useful approaches in chemical risk assessment, but the methodology to incorporate toxicogenomic data into a PBPK model to inform risk assessment remains to be developed. This study aimed to develop a probabilistic human health risk assessment approach by integrating toxicogenomic dose-response data and PBPK modeling using perfluorooctane sulfonate (PFOS) as a case study. Based on the available human in vitro and mouse in vivo toxicogenomic data, we identified the differentially expressed genes (DEGs) at each exposure paradigm/duration. Kyoto Encyclopedia of Genes and Genomes and disease ontology enrichment analyses were conducted on the DEGs to identify significantly enriched pathways and diseases. The dose-response data of DEGs were analyzed using the Bayesian benchmark dose (BMD) method. Using a previously published PBPK model, the gene BMDs were converted to human equivalent doses (HEDs), which were summarized to pathway and disease HEDs and then extrapolated to reference doses (RfDs) by considering an uncertainty factor of 30 for mouse in vivo data and 10 for human in vitro data. The results suggested that the median RfDs at different exposure paradigms were similar to the 2016 U.S. Environmental Protection Agency's recommended RfD, while the RfDs for the most sensitive pathways and diseases were closer to the recent European Food Safety Authority's guidance values. In conclusion, genomic dose-response data and PBPK modeling can be integrated to become a useful alternative approach in risk assessment of environmental chemicals. This approach considers multiple endpoints, provides toxicity mechanistic insights, and does not rely on apical toxicity endpoints.

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.

Development of a multi-route physiologically based pharmacokinetic (PBPK) model for nanomaterials: a comparison between a traditional versus a new route-specific approach using gold nanoparticles in rats.

BACKGROUND: Physiologically based pharmacokinetic (PBPK) modeling is an important tool in predicting target organ dosimetry and risk assessment of nanoparticles (NPs). The methodology of building a multi-route PBPK model for NPs has not been established, nor systematically evaluated. In this study, we hypothesized that the traditional route-to-route extrapolation approach of PBPK modeling that is typically used for small molecules may not be appropriate for NPs. To test this hypothesis, the objective of this study was to develop a multi-route PBPK model for different sizes (1.4-200 nm) of gold nanoparticles (AuNPs) in adult rats following different routes of administration (i.e., intravenous (IV), oral gavage, intratracheal instillation, and endotracheal inhalation) using two approaches: a traditional route-to-route extrapolation approach for small molecules and a new approach that is based on route-specific data that we propose to be applied generally to NPs. RESULTS: We found that the PBPK model using this new approach had superior performance than the traditional approach. The final PBPK model was optimized rigorously using a Bayesian hierarchical approach with Markov chain Monte Carlo simulations, and then converted to a web-based interface using R Shiny. In addition, quantitative structure-activity relationships (QSAR) based multivariate linear regressions were established to predict the route-specific key biodistribution parameters (e.g., maximum uptake rate) based on the physicochemical properties of AuNPs (e.g., size, surface area, dose, Zeta potential, and NP numbers). These results showed the size and surface area of AuNPs were the main determinants for endocytic/phagocytic uptake rates regardless of the route of administration, while Zeta potential was an important parameter for the estimation of the exocytic release rates following IV administration. CONCLUSIONS: This study suggests that traditional route-to-route extrapolation approaches for PBPK modeling of small molecules are not applicable to NPs. Therefore, multi-route PBPK models for NPs should be developed using route-specific data. This novel PBPK-based web interface serves as a foundation for extrapolating to other NPs and to humans to facilitate biodistribution estimation, safety, and risk assessment of NPs.

Prediction of drug-induced kidney injury in drug discovery.

Drug induced kidney injury is one of the leading causes of failure of drug development programs in the clinic. Early prediction of renal toxicity potential of drugs is crucial to the success of drug candidates in the clinic. The dynamic nature of the functioning of the kidney and the presence of drug uptake proteins introduce additional challenges in the prediction of renal injury caused by drugs. Renal injury due to drugs can be caused by a wide variety of mechanisms and can be broadly classified as toxic or obstructive. Several biomarkers are available for in vitro and in vivo detection of renal injury. In vitro static and dynamic (microfluidic) cellular models and preclinical models can provide valuable information regarding the toxicity potential of drugs. Differences in pharmacology and subsequent disconnect in biomarker response, differences in the expression of transporter and enzyme proteins between in vitro to in vivo systems and between preclinical species and humans are some of the limitations of current experimental models. The progress in microfluidic (kidney-on-chip) platforms in combination with the ability of 3-dimensional cell culture can help in addressing some of these issues in the future. Finally, newer in silico and computational techniques like physiologically based pharmacokinetic modeling and machine learning have demonstrated potential in assisting prediction of drug induced kidney injury.

PBPK Absorption Modeling of Food Effect and Bioequivalence in Fed State for Two Formulations with Crystalline and Amorphous Forms of BCS 2 Class Drug in Generic Drug Development.

Prediction of the effect of food on drug's pharmacokinetics using modeling and simulation could cause difficulties due to complex in vivo processes. A generic formulation with amorphous form of BCS 2 class drug substance was developed and compared in vitro and in vivo to the reference drug product with drug substance in crystalline form. In order to approve generic formulation, some regulatory agencies are requesting to perform bioequivalence (BE) studies also in fed state. Food can have various effects on drug dissolution and absorption, depending also on drug's properties. A physiologically based pharmacokinetic (PBPK) absorption model was built in GastroPlus™ to predict the food effect on generic and reference formulation and to predict the fed BE study outcome. During model development, we were searching for model inputs that impact and describe in vivo behavior of amorphous and crystalline forms of active pharmaceutical ingredient (API) in fast and fed conditions. The developed model was able to predict the food effect with up to 10% prediction error (PE). Performed virtual BE trials confirmed the BE of drug products in fed state. Our model was able to capture the difference between the two drug products containing different forms of API (amorphous and crystalline) and predict the food effect on both formulations.

In Vitro, in Silico, and in Vivo Assessments of Intestinal Precipitation and Its Impact on Bioavailability of a BCS Class 2 Basic Compound.

In this study, a multipronged approach of in vitro experiments, in silico simulations, and in vivo studies was developed to evaluate the dissolution, supersaturation, precipitation, and absorption of three formulations of Compound-A, a BCS class 2 weak base with pH-dependent solubility. In in vitro 2-stage dissolution experiments, the solutions were highly supersaturated with no precipitation at the low dose but increasing precipitation at higher doses. No difference in precipitation was observed between the capsules and tablets. The in vitro precipitate was found to be noncrystalline with higher solubility than the crystalline API, and was readily soluble when the drug concentration was lowered by dilution. A gastric transit and biphasic dissolution (GTBD) model was developed to better mimic gastric transfer and intestinal absorption. Precipitation was also observed in GTBD, but the precipitate redissolved and partitioned into the organic phase. In vivo data from the phase 1 clinical trial showed linear and dose proportional PK for the formulations with no evidence of in vivo precipitation. While the in vitro precipitation observed in the 2-stage dissolution appeared to overestimate in vivo precipitation, the GTBD model provided absorption profiles consistent with in vivo data. In silico simulation of plasma concentrations by GastroPlus using biorelevant in vitro dissolution data from the tablets and capsules and assuming negligible precipitation was in line with the observed in vivo profiles of the two formulations. The totality of data generated with Compound-A indicated that the bioavailability differences among the three formulations were better explained by the differences in gastric dissolution than intestinal precipitation. The lack of intestinal precipitation was consistent with several other BCS class 2 basic compounds in the literature for which highly supersaturated concentrations and rapid absorption were also observed.

Physiologically based pharmacokinetic-quantitative systems toxicology and safety (PBPK-QSTS) modeling approach applied to predict the variability of amitriptyline pharmacokinetics and cardiac safety in populations and in individuals.

The physiologically based pharmacokinetic (PBPK) models allow for predictive assessment of variability in population of interest. One of the future application of PBPK modeling is in the field of precision dosing and personalized medicine. The aim of the study was to develop PBPK model for amitriptyline given orally, predict the variability of cardiac concentrations of amitriptyline and its main metabolite-nortriptyline in populations as well as individuals, and simulate the influence of those xenobiotics in therapeutic and supratherapeutic concentrations on human electrophysiology. The cardiac effect with regard to QT and RR interval lengths was assessed. The Emax model to describe the relationship between amitriptyline concentration and heart rate (RR) length was proposed. The developed PBPK model was used to mimic 29 clinical trials and 19 cases of amitriptyline intoxication. Three clinical trials and 18 cases were simulated with the use of PBPK-QSTS approach, confirming lack of cardiotoxic effect of amitriptyline in therapeutic doses and the increase in heart rate along with potential for arrhythmia development in case of amitriptyline overdose. The results of our study support the validity and feasibility of the PBPK-QSTS modeling development for personalized medicine.

A physiologically based pharmacokinetic model to predict the pharmacokinetics of highly protein-bound drugs and the impact of errors in plasma protein binding.

Predicting the pharmacokinetics of highly protein-bound drugs is difficult. Also, since historical plasma protein binding data were often collected using unbuffered plasma, the resulting inaccurate binding data could contribute to incorrect predictions. This study uses a generic physiologically based pharmacokinetic (PBPK) model to predict human plasma concentration-time profiles for 22 highly protein-bound drugs. Tissue distribution was estimated from in vitro drug lipophilicity data, plasma protein binding and the blood: plasma ratio. Clearance was predicted with a well-stirred liver model. Underestimated hepatic clearance for acidic and neutral compounds was corrected by an empirical scaling factor. Predicted values (pharmacokinetic parameters, plasma concentration-time profile) were compared with observed data to evaluate the model accuracy. Of the 22 drugs, less than a 2-fold error was obtained for the terminal elimination half-life (t1/2 , 100% of drugs), peak plasma concentration (Cmax , 100%), area under the plasma concentration-time curve (AUC0-t , 95.4%), clearance (CLh , 95.4%), mean residence time (MRT, 95.4%) and steady state volume (Vss , 90.9%). The impact of fup errors on CLh and Vss prediction was evaluated. Errors in fup resulted in proportional errors in clearance prediction for low-clearance compounds, and in Vss prediction for high-volume neutral drugs. For high-volume basic drugs, errors in fup did not propagate to errors in Vss prediction. This is due to the cancellation of errors in the calculations for tissue partitioning of basic drugs. Overall, plasma profiles were well simulated with the present PBPK model. Copyright © 2016 John Wiley & Sons, Ltd.

Characterization of the human kinetic adjustment factor for the health risk assessment of environmental contaminants.

A default uncertainty factor of 3.16 (√10) is applied to account for interindividual variability in toxicokinetics when performing non-cancer risk assessments. Using relevant human data for specific chemicals, as WHO/IPCS suggests, it is possible to evaluate, and replace when appropriate, this default factor by quantifying chemical-specific adjustment factors for interindividual variability in toxicokinetics (also referred to as the human kinetic adjustment factor, HKAF). The HKAF has been determined based on the distributions of pharmacokinetic parameters (e.g., half-life, area under the curve, maximum blood concentration) in relevant populations. This article focuses on the current state of knowledge of the use of physiologically based algorithms and models in characterizing the HKAF for environmental contaminants. The recent modeling efforts on the computation of HKAF as a function of the characteristics of the population, chemical and its mode of action (dose metrics), as well as exposure scenario of relevance to the assessment are reviewed here. The results of these studies, taken together, suggest the HKAF varies as a function of the sensitive subpopulation and dose metrics of interest, exposure conditions considered (route, duration, and intensity), metabolic pathways involved and theoretical model underlying its computation. The HKAF seldom exceeded the default value of 3.16, except in very young children (i.e., <≈ 3 months) and when the parent compound is the toxic moiety. Overall, from a public health perspective, the current state of knowledge generally suggest that the default uncertainty factor is sufficient to account for human variability in non-cancer risk assessments of environmental contaminants.

Whole-Body Disposition and Physiologically Based Pharmacokinetic Modeling of Adeno-Associated Viruses and the Transgene Product.

To facilitate model-informed drug development (MIDD) of adeno-associated virus (AAV) therapy, here we have developed a physiologically based pharmacokinetic (PBPK) model for AAVs following preclinical investigation in mice. After 2E11 Vg/mouse dose of AAV8 and AAV9 encoding a monoclonal antibody (mAb) gene, whole-body disposition of both the vector and the transgene mAb was evaluated over 3 weeks. At steady-state, the following tissue-to-blood (T/B) concentration ratios were found for AAV8/9: ∼50 for liver; ∼10 for heart and muscle; ∼2 for brain, lung, kidney, adipose, and spleen; ≤1 for bone, skin, and pancreas. T/B values for mAb were compared with the antibody biodistribution coefficients, and five different clusters of organs were identified based on their transgene expression profile. All the biodistribution data were used to develop a novel AAV PBPK model that incorporates: (i) whole-body distribution of the vector; (ii) binding, internalization, and intracellular processing of the vector; (iii) transgene expression and secretion; and (iv) whole-body disposition of the secreted transgene product. The model was able to capture systemic and tissue PK of the vector and the transgene-produced mAb reasonably well. Pathway analysis of the PBPK model suggested that liver, muscle, and heart are the main contributors for the secreted transgene mAb. Unprecedented PK data and the novel PBPK model developed here provide the foundation for quantitative systems pharmacology (QSP) investigations of AAV-mediated gene therapies. The PBPK model can also serve as a quantitative tool for preclinical study design and preclinical-to-clinical translation of AAV-based gene therapies.

Preface to the special issue of Food and Chemical Toxicology on "New approach methodologies and machine learning in food safety and chemical risk assessment: Development of reproducible, open-source, and user-friendly tools for exposure, toxicokinetic, and toxicity assessments in the 21st century".

This Special Issue contains articles on applications of various new approach methodologies (NAMs) in the field of toxicology and risk assessment. These NAMs include in vitro high-throughput screening, quantitative structure-activity relationship (QSAR) modeling, physiologically based pharmacokinetic (PBPK) modeling, network toxicology analysis, molecular docking simulation, omics, machine learning, deep learning, and "template-and-anchor" multiscale computational modeling. These in vitro and in silico approaches complement each other and can be integrated together to support different applications of toxicology, including food safety assessment, dietary exposure assessment, chemical toxicity potency screening and ranking, chemical toxicity prediction, chemical toxicokinetic simulation, and to investigate the potential mechanisms of toxicities, as introduced further in selected articles in this Special Issue.

Applying Physiologically Based Pharmacokinetic Modeling to Interpret Carbamazepine's Nonlinear Pharmacokinetics and Its Induction Potential on Cytochrome P450 3A4 and Cytochrome P450 2C9 Enzymes.

Carbamazepine (CBZ) is commonly prescribed for epilepsy and frequently used in polypharmacy. However, concerns arise regarding its ability to induce the metabolism of other drugs, including itself, potentially leading to the undertreatment of co-administered drugs. Additionally, CBZ exhibits nonlinear pharmacokinetics (PK), but the root causes have not been fully studied. This study aims to investigate the mechanisms behind CBZ's nonlinear PK and its induction potential on CYP3A4 and CYP2C9 enzymes. To achieve this, we developed and validated a physiologically based pharmacokinetic (PBPK) parent-metabolite model of CBZ and its active metabolite Carbamazepine-10,11-epoxide in GastroPlus®. The model was utilized for Drug-Drug Interaction (DDI) prediction with CYP3A4 and CYP2C9 victim drugs and to further explore the underlying mechanisms behind CBZ's nonlinear PK. The model accurately recapitulated CBZ plasma PK. Good DDI performance was demonstrated by the prediction of CBZ DDIs with quinidine, dolutegravir, phenytoin, and tolbutamide; however, with midazolam, the predicted/observed DDI AUClast ratio was 0.49 (slightly outside of the two-fold range). CBZ's nonlinear PK can be attributed to its nonlinear metabolism caused by autoinduction, as well as nonlinear absorption due to poor solubility. In further applications, the model can help understand DDI potential when CBZ serves as a CYP3A4 and CYP2C9 inducer.

Development and Validation of Physiologically Based Pharmacokinetic Model of Levetiracetam to Predict Exposure and Dose Optimization in Pediatrics.

Levetiracetam (Lev) is an antiepileptic drug that has been increasingly used in the epilepsy pediatric population in recent years, but its pharmacokinetic behavior in pediatric population needs to be characterized clearly. Clinical trials for the pediatric drug remain difficult to conduct due to ethical and practical factors. The purpose of this study was to use the physiologically based pharmacokinetic (PBPK) model to predict changes in plasma exposure of Lev in pediatric patients and to provide recommendations for dose adjustment. A PBPK model of Lev in adults was developed using PK-Sim® software and extrapolated to the entire age range of the pediatric population. The model was evaluated using clinical pharmacokinetic data. The results showed the good fit between predictions and observations of the adult and pediatric models. The recommended doses for neonates, infants and children are 0.78, 1.67 and 1.22 times that of adults, respectively. Moreover, at the same dose, plasma exposure in adolescents was similar to that of adults. The PBPK models of Lev for adults and pediatrics were successfully developed and validated to provide a reference for the rational administration of drugs in the pediatric population.

A Literature Review of Changes in Phase II Drug-Metabolizing Enzyme and Drug Transporter Expression during Pregnancy.

The purpose of this literature review is to comprehensively summarize changes in the expression of phase II drug-metabolizing enzymes and drug transporters in both the pregnant woman and the placenta. Using PubMed®, a systematic search was conducted to identify literature relevant to drug metabolism and transport in pregnancy. PubMed was searched with pre-specified terms during the period of 26 May 2023 to 10 July 2023. The final dataset of 142 manuscripts was evaluated for evidence regarding the effect of gestational age and hormonal regulation on the expression of phase II enzymes (n = 16) and drug transporters (n = 38) in the pregnant woman and in the placenta. This comprehensive review exposes gaps in current knowledge of phase II enzyme and drug transporter localization, expression, and regulation during pregnancy, which emphasizes the need for further research. Moreover, the information collected in this review regarding phase II drug-metabolizing enzyme and drug transporter changes will aid in optimizing pregnancy physiologically based pharmacokinetic (PBPK) models to inform dose selection in the pregnant population.

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.

Physiologically Based Pharmacokinetic Modeling to Characterize the Effect of Molecular Charge on Whole-Body Disposition of Monoclonal Antibodies.

Motivated by a series of work demonstrating the effect of molecular charge on antibody pharmacokinetics (PK), physiological-based pharmacokinetic (PBPK) models are emerging that relate in silico calculated charge or in vitro measures of polyspecificity to antibody PK parameters. However, only plasma data has been used for model development in these studies, leading to unvalidated assumptions. Here, we present an extended platform PBPK model for antibodies that incorporate charge-dependent endothelial cell pinocytosis rate and nonspecific off-target binding in the interstitial space and on circulating blood cells, to simultaneously characterize whole-body disposition of three antibody charge variants. Predictive potential of various charge metrics was also explored, and the difference between positive charge patches and negative charge patches (i.e., PPC-PNC) was used as the charge parameter to establish quantitative relationships with nonspecific binding affinities and endothelial cell uptake rate. Whole-body disposition of these charge variants was captured well by the model, with less than 2-fold predictive error in area under the curve of most plasma and tissue PK data. The model also predicted that with greater positive charge, nonspecific binding was more substantial, and pinocytosis rate increased especially in brain, heart, kidney, liver, lung, and spleen, but remained unchanged in adipose, bone, muscle, and skin. The presented PBPK model contributes to our understanding of the mechanisms governing the disposition of charged antibodies and can be used as a platform to guide charge engineering based on desired plasma and tissue exposures.