Brivaracetam (BRV) is a new third-generation antiseizure medication for the treatment of focal epileptic seizures. Its use has been increasing among epileptic populations in recent years, but pharmacokinetic (PK) behavior may change in hepatic impairment and the elderly populations. Due to ethical constraints, clinical trials are difficult to conduct and data are limited. This study used PK-Sim® to develop a physiologically based pharmacokinetic (PBPK) model for adults and extrapolate it to hepatic impairment and the elderly populations. The model was evaluated with clinical PK data, and dosage explorations were conducted. For the adult population with mild hepatic impairment, the dose is recommended to be adjusted to 70% of the recommended dose, and to 60% for moderate and severe hepatic impairment. For the elderly population with mild hepatic impairment under 80 years old, it is recommended that the dose be adjusted to 60% of the recommended dose and to 50% for moderate and severe conditions. The elderly population with hepatic impairment over 80 years old is adjusted to 50% of the recommended dose for all stages. Healthy elderly do not need to adjust. The BRV PBPK model was successfully developed, studying exposure in hepatic impairment and elderly populations and optimizing dosing regimens.

Despite several screening levels for NDMA reported in water, soil, air, and drugs, the human risk assessment using biomonitoring concentrations has not been performed. In this study, gender-specific exposure guidance values were determined in humans, then biomonitoring measurements in healthy Korean subjects (32 men and 40 women) were compared to the exposure guidance values to evaluate the current exposure level to NDMA. For the human risk assessment of NDMA, the gender-specific physiologically based pharmacokinetic (PBPK) model was developed in humans using proper physiological parameters, partition coefficients, and biochemical parameters. Using the PBPK model, a Monte Carlo simulation was performed to describe the magnitudes of inter-individual variability and uncertainty on the single model predictions. The PBPK modeling and Monte Carlo simulation allowed the estimation of the relationship between external dose and blood concentration for the risk assessment. The procedure for the human risk assessment was summarized as follows: (1) estimating a steady-state blood concentration (Cavg) corresponding to the daily no observed adverse effect level (NOAEL) administration in rats; (2) applying uncertainty factors (UFs) for deriving the human Cavg; (3) determining the exposure guidance values as screening criteria; (4) interpreting the human biomonitoring measurements by forward and reverse dosimetry approaches. Using the biomonitoring concentrations, current daily exposures to NDMA were estimated to be 3.95 μg/day/kg for men and 10.60 μg/day/kg for women, respectively. The result of the study could be used as a basis for implementing further risk management and regulatory decision-making for NDMA.

Venlafaxine, a serotonin-norepinephrine reuptake inhibitor (SNRI), is indicated for the treatment of major depressive disorder, social anxiety disorder, generalized anxiety disorder, and panic disorder. Venlafaxine is metabolized to the active metabolite desvenlafaxine mainly by CYP2D6. Genetic polymorphism of CYP2D6 and coadministration with other medications can significantly affect the pharmacokinetics and/or pharmacodynamics of venlafaxine and its active metabolite. This study aimed to establish the PBPK models of venlafaxine and its active metabolite related to CYP2D6 genetic polymorphism and to predict drug-drug interactions (DDIs) with clarithromycin and paroxetine in different CYP2D6 genotypes. Clinical pharmacogenomic data for venlafaxine and desvenlafaxine were collected to build the PBPK model. Physicochemical and absorption, distribution, metabolism, and excretion (ADME) characteristics of respective compounds were obtained from previously reported data, predicted by the PK-Sim® software, or optimized to capture the plasma concentration-time profiles. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and plasma concentration-time profiles to the observed data. Predicted plasma concentration-time profiles of venlafaxine and its active metabolite were visually similar to the observed profiles and all predicted AUC and Cmax values for respective compounds were included in the twofold error range of observed values in non-genotyped populations and different CYP2D6 genotypes. When clarithromycin or clarithromycin plus paroxetine was concomitantly administered, predicted plasma concentration-time profiles of venlafaxine properly captured the observed profiles in two different CYP2D6 genotypes and all predicted DDI ratios for AUC and Cmax were included within the acceptance range. Consequently, the present model successfully captured the pharmacokinetic alterations of venlafaxine and its active metabolite according to CYP2D6 genetic polymorphism as well as the DDIs between venlafaxine and two CYP inhibitors. The present model can be used to predict the pharmacokinetics of venlafaxine and its active metabolite considering different races, ages, coadministered drugs, and CYP2D6 activity of individuals and it can contribute to individualized pharmacotherapy of venlafaxine.

This study aims to explore and characterize the role of pediatric sedation via rectal route. A pediatric physiologically based pharmacokinetic-pharmacodynamic (PBPK/PD) model of midazolam gel was built and validated to support dose selection for pediatric clinical trials. Before developing the rectal PBPK model, an intravenous PBPK model was developed to determine drug disposition, specifically by describing the ontogeny model of the metabolic enzyme. Pediatric rectal absorption was developed based on the rectal PBPK model of adults. The improved Weibull function with permeability, surface area, and fluid volume parameters was used to extrapolate pediatric rectal absorption. A logistic regression model was used to characterize the relationship between the free concentrations of midazolam and the probability of sedation. All models successfully described the PK profiles with absolute average fold error (AAFE) < 2, especially our intravenous PBPK model that extended the predicted age to preterm. The simulation results of the PD model showed that when the free concentrations of midazolam ranged from 3.9 to 18.4 ng/mL, the probability of \"Sedation\" was greater than that of \"Not-sedation\" states. Combined with the rectal PBPK model, the recommended sedation doses were in the ranges of 0.44-2.08 mg/kg for children aged 2-3 years, 0.35-1.65 mg/kg for children aged 4-7 years, 0.24-1.27 mg/kg for children aged 8-12 years, and 0.20-1.10 mg/kg for adolescents aged 13-18 years. Overall, this model mechanistically quantified drug disposition and effect of midazolam gel in the pediatric population, accurately predicted the observed clinical data, and simulated the drug exposure for sedation that will inform dose selection for following pediatric clinical trials.

To ensure safety in chemical plants handling a wide variety of liquid and gaseous hazardous substances, it is necessary to carry out highly accurate risk assessments and take appropriate measures. In this study, a risk assessment method was developed for the problem of the leakage of liquid hazardous substances. The risk assessment of toxic liquid leaks must consider the exposure of workers to the liquid and toxic gases produced by vaporization. The absorption and subsequent metabolism of hazardous substances in the body via multiple pathways after exposure to liquids and gases was calculated using a pharmacokinetic model. Estimation of exposure concentrations of toxic gases volatilized from leaked liquids was reproduced by computational fluid dynamics simulation. In this study, ethanol was selected as the hazardous substance and the risk of hazardous liquid leakage was assessed. The results of the analysis, which considered liquid and gas exposure under the conditions of the assumed scenario, showed that the maximum blood concentration of ethanol was 1640 µmol/L, which is sufficiently low compared to the concentration of 10,900 µmol/L at which acute toxic effects become apparent. These results suggest that work can be carried out safely under the conditions of the assumed scenario. The risk assessment methodology for liquid spills in this study confirms that risk assessment is possible under multiple scenarios, including individual differences, activity conditions, and the use of protective equipment.

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

Pantoprazole is used to treat gastroesophageal reflux disease (GERD), maintain healing of erosive esophagitis (EE), and control symptoms related to Zollinger-Ellison syndrome (ZES). Pantoprazole is mainly metabolized by cytochrome P450 (CYP) 2C19, converting to 4\'-demethyl pantoprazole. CYP2C19 is a genetically polymorphic enzyme, and the genetic polymorphism affects the pharmacokinetics and/or pharmacodynamics of pantoprazole. In this study, we aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of pantoprazole in populations with various CYP2C19 metabolic activities. A comprehensive investigation of previous reports and drug databases was conducted to collect the clinical pharmacogenomic data, physicochemical data, and disposition properties of pantoprazole, and the collected data were used for model establishment. The model was evaluated by comparing the predicted plasma concentration-time profiles and/or pharmacokinetic parameters (AUC and Cmax) with the clinical observation results. The predicted plasma concentration-time profiles in different CYP2C19 phenotypes properly captured the observed profiles. All fold error values for AUC and Cmax were included in the two-fold range. Consequently, the minimal PBPK model for pantoprazole related to CYP2C19 genetic polymorphism was properly established and it can predict the pharmacokinetics of pantoprazole in different CYP2C19 phenotypes. The present model can broaden the insight into the individualized pharmacotherapy for pantoprazole.

Irbesartan, a potent and selective angiotensin II type-1 (AT1) receptor blocker (ARB), is one of the representative medications for the treatment of hypertension. Cytochrome P450 (CYP) 2C9 is primarily involved in the oxidation of irbesartan. CYP2C9 is highly polymorphic, and genetic polymorphism of this enzyme is the leading cause of significant alterations in the pharmacokinetics of irbesartan. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes. The irbesartan PBPK model was established using the PK-Sim® software. Our previously reported pharmacogenomic data for irbesartan was leveraged in the development of the PBPK model and collected clinical pharmacokinetic data for irbesartan was used for the validation of the model. Physicochemical and ADME properties of irbesartan were obtained from previously reported data, predicted by the modeling software, or optimized to fit the observed plasma concentration-time profiles. Model evaluation was performed by comparing the predicted plasma concentration-time profiles and pharmacokinetic parameters to the observed results. Predicted plasma concentration-time profiles were visually similar to observed profiles. Predicted AUCinf in CYP2C9*1/*3 and CYP2C9*1/*13 genotypes were increased by 1.54- and 1.62-fold compared to CYP2C9*1/*1 genotype, respectively. All fold error values for AUC and Cmax in non-genotyped and CYP2C9 genotyped models were within the two-fold error criterion. We properly established the PBPK model of irbesartan in different CYP2C9 genotypes. It can be used to predict the pharmacokinetics of irbesartan for personalized pharmacotherapy in individuals of various races, ages, and CYP2C9 genotypes.

A permeability-limited physiologically based pharmacokinetic (PBPK) model featuring four subcompartments (corresponding to the intracellular and extracellular water of the tissue, the residual plasma, and blood cells) for each tissue has been developed in MATLAB/SimBiology and applied to various what-if scenario simulations. This model allowed us to explore the complex interplay of passive permeability, metabolism in tissue or residual blood, active uptake or efflux transporters, and different dosing routes (intravenous (IV) or oral (PO)) in determining the dynamics of the tissue/plasma partition coefficient (Kp) and volume of distribution (Vd) within a realistic pseudo-steady state. Based on the modeling exercise, the permeability, metabolism, and transporters demonstrated significant effects on the dynamics of the Kp and Vd for IV bolus administration and PO fast absorption, but these effects were not as pronounced for IV infusion or PO slow absorption. Especially for low-permeability compounds, uptake transporters were found to increase both the Kp and Vd at the pseudo-steady state (Vdss), while efflux transporters had the opposite effect of decreasing the Kp and Vdss. For IV bolus administration and PO fast absorption, increasing tissue metabolism was predicted to elevate the Kp and Vdss, which contrasted with the traditional derivation from the steady-state perfusion-limited PBPK model. Moreover, metabolism in the residual blood had more impact on the Kp and Vdss compared to metabolism in tissue. Due to its ability to offer a more realistic description of tissue dynamics, the permeability-limited PBPK model is expected to gain broader acceptance in describing clinical PK and observed Kp and Vdss, even for certain small molecules like cyclosporine, which are currently treated as perfusion-limited in commercial PBPK platforms.

Withdrawal periods for diclazuril in broilers have traditionally been determined through regression analysis. However, over the last two decades, the physiologically based pharmacokinetic (PBPK) model has gained prominence as a predictive tool for veterinary drug residues, which offers an alternative method for establishing appropriate withdrawal periods for veterinary drugs. In this current study, a flow-limited PBPK model was developed to predict diclazuril concentrations in broilers following long-duration administration via medicated feed and water. This model consists of nine compartments, including arterial and venous plasma, lung, muscle, skin + fat, kidney, liver, intestine contents, and the rest of the body compartment. Physiological parameters such as tissue weights (Vcxx) and blood flow (Qcxx) were gathered from published studies, and tissue/plasma partition coefficients (Pxx) were calculated through the area method or parameter optimization. Published diclazuril concentrations were compared to the predicted values, indicating the accuracy and validity of the model. The sensitivity analysis showed that parameters associated with cardiac output, drug absorption, and elimination significantly affected diclazuril concentrations in the muscle. Finally, a Monte Carlo analysis, consisting of 1000 iterations, was conducted to calculate the withdrawal period. Based on the Chinese MRL values, we calculated a withdrawal period of 0 days for both recommended dosing regimens (through mediated water and feed at concentrations of 0.5-1 mg/L and 1 mg/kg, respectively). However, based on the European MRLs, longer periods were determined for the mediated feed dosing route. Our model provides a foundation for scaling other coccidiostats and poultry species.