Predicting Passive and Active Tissue:Plasma Partition Coefficients: Interindividual and Interspecies Variability

Unbound Fractions of PFAS in Human and Rodent Tissues: Rat Liver a Suitable Proxy for Evaluating Emerging PFAS?

Adverse health effects associated with exposures to perfluoroalkyl and polyfluoroalkyl substances (PFAS) are a concern for public health and are driven by their elimination half-lives and accumulation in specific tissues. However, data on PFAS binding in human tissues are limited. Accumulation of PFAS in human tissues has been linked to interactions with specific proteins and lipids in target organs. Additional data on PFAS binding and unbound fractions (funbound) in whole human tissues are urgently needed. Here, we address this gap by using rapid equilibrium dialysis to measure the binding and funbound of 16 PFAS with 3 to 13 perfluorinated carbon atoms (ηpfc = 3-13) and several functional headgroups in human liver, lung, kidney, heart, and brain tissue. We compare results to mouse (C57BL/6 and CD-1) and rat tissues. Results show that funbound decreases with increasing fluorinated carbon chain length and hydrophobicity. Among human tissues, PFAS binding was generally greatest in brain > liver ≈ kidneys ≈ heart > lungs. A correlation analysis among human and rodent tissues identified rat liver as a suitable surrogate for predicting funbound for PFAS in human tissues (R2 ≥ 0.98). The funbound data resulting from this work and the rat liver prediction method offer input parameters and tools for toxicokinetic models for legacy and emerging PFAS.

Application of partition coefficient methods to predict tissue:plasma affinities in common farm animals: Influence of ionisation state

Tissue affinities are conventionally determined from in vivo steady-state tissue and plasma or plasma-water chemical concentration data. In silico approaches were initially developed for preclinical species but standardly applied and tested in human physiologically-based kinetic (PBK) models. Recently, generic PBK models for farm animals have been made available and require partition coefficients as input parameters. In the current investigation, data for species-specific tissue compositions have been collected, and prediction of chemical distribution in various tissues of livestock species for cattle, chicken, sheep and swine have been performed. Overall, tissue composition was very similar across the four farm animal species. However, small differences were observed in moisture, fat and protein content in the various organs within each species. Such differences could be attributed to factors such as variations in age, breed, and weight of the animals and general conditions of the animal itself. With regards to the predictions of tissue:plasma partition coefficients, 80 %, 71 %, 77 % of the model predictions were within a factor 10 using the methods of Berezhkovskiy (2004), Rodgers and Rowland (2006) and Schmitt (2008). The method of Berezhkovskiy (2004) was often providing the most reliable predictions except for swine, where the method of Schmitt (2008) performed best. In addition, investigation of the impact of chemical classes on prediction performance, all methods had very similar reliability. Notwithstanding, no clear pattern regarding specific chemicals or tissues could be detected for the values predicted outside a 10-fold change in certain chemicals or specific tissues. This manuscript concludes with the need for future research, particularly focusing on lipophilicity and species differences in protein binding.

Interspecies scaling of toxicity reference values in human health versus ecological risk assessments: A critical review

Risk assessments that focus on anthropogenic chemicals in environmental media-whether considering human health or ecological effects-often rely on toxicity data from experimentally studied species to estimate safe exposures for species that lack similar data. Current default extrapolation approaches used in both human health risk assessments and ecological risk assessments (ERAs) account for differences in body weight between the test organisms and the species of interest, but the two default approaches differ in important ways. Human health risk assessments currently employ a default based on body weight raised to the three-quarters power. Ecological risk assessments for wildlife (i.e., mammals and birds) are typically based directly on body weight, as measured in the test organism and receptor species. This review describes differences in the experimental data underlying these default practices and discusses the many factors that affect interspecies variability in chemical exposures. The interplay of these different factors can lead to substantial departures from default expectations. Alternative methodologies for conducting more accurate interspecies extrapolations in ERAs for wildlife are discussed, including tissue-based toxicity reference values, physiologically based toxicokinetic and/or toxicodynamic modeling, chemical read-across, and a system of categorical defaults based on route of exposure and toxic mode of action. Integr Environ Assess Manag 2024;20:749-764. © 2023 SETAC.

A minimal PBPK model to accelerate preclinical development of drugs against tuberculosis

Introduction: Understanding drug exposure at disease target sites is pivotal to profiling new drug candidates in terms of tolerability and efficacy. Such quantification is particularly tedious for anti-tuberculosis (TB) compounds as the heterogeneous pulmonary microenvironment due to the infection may alter lung permeability and affect drug disposition. Murine models have been a longstanding support in TB research so far and are here used as human surrogates to unveil the distribution of several anti-TB compounds at the site-of-action via a novel and centralized PBPK design framework. Methods: As an intermediate approach between data-driven pharmacokinetic (PK) models and whole-body physiologically based (PB) PK models, we propose a parsimonious framework for PK investigation (minimal PBPK approach) that retains key physiological processes involved in TB disease, while reducing computational costs and prior knowledge requirements. By lumping together pulmonary TB-unessential organs, our minimal PBPK model counts 9 equations compared to the 36 of published full models, accelerating the simulation more than 3-folds in Matlab 2022b. Results: The model has been successfully tested and validated against 11 anti-TB compounds-rifampicin, rifapentine, pyrazinamide, ethambutol, isoniazid, moxifloxacin, delamanid, pretomanid, bedaquiline, OPC-167832, GSK2556286 - showing robust predictability power in recapitulating PK dynamics in mice. Structural inspections on the proposed design have ensured global identifiability and listed free fraction in plasma and blood-to-plasma ratio as top sensitive parameters for PK metrics. The platform-oriented implementation allows fast comparison of the compounds in terms of exposure and target attainment. Discrepancies in plasma and lung levels for the latest BPaMZ and HPMZ regimens have been analyzed in terms of their impact on preclinical experiment design and on PK/PD indices. Conclusion: The framework we developed requires limited drug- and species-specific information to reconstruct accurate PK dynamics, delivering a unified viewpoint on anti-TB drug distribution at the site-of-action and a flexible fit-for-purpose tool to accelerate model-informed drug design pipelines and facilitate translation into the clinic.

Paired Liver:Plasma PFAS Concentration Ratios from Adolescents in the Teen-LABS Study and Derivation of Empirical and Mass Balance Models to Predict and Explain Liver PFAS Accumulation

Animal studies have pointed at the liver as a hotspot for per- and polyfluoroalkyl substances (PFAS) accumulation and toxicity; however, these findings have not been replicated in human populations. We measured concentrations of seven PFAS in matched liver and plasma samples collected at the time of bariatric surgery from 64 adolescents in the Teen-Longitudinal Assessment of Bariatric Surgery (Teen-LABS) study. Liver:plasma concentration ratios were perfectly explained (r2 > 0.99) in a multilinear regression (MLR) model based on toxicokinetic (TK) descriptors consisting of binding to tissue constituents and membrane permeabilities. Of the seven matched plasma and liver PFAS concentrations compared in this study, the liver:plasma concentration ratio of perfluoroheptanoic acid (PFHpA) was considerably higher than the liver:plasma concentration ratio of other PFAS congeners. Comparing the MLR model with an equilibrium mass balance model (MBM) suggested that complex kinetic transport processes are driving the unexpectedly high liver:plasma concentration ratio of PFHpA. Intratissue MBM modeling pointed to membrane lipids as the tissue constituents that drive the liver accumulation of long-chain, hydrophobic PFAS, whereas albumin binding of hydrophobic PFAS dominated PFAS distribution in plasma. The liver:plasma concentration data set, empirical MLR model, and mechanistic MBM modeling allow the prediction of liver from plasma concentrations measured in human cohort studies. Our study demonstrates that combining biomonitoring data with mechanistic modeling can identify underlying mechanisms of internal distribution and specific target organ toxicity of PFAS in humans.

Investigation of simplified physiologically-based pharmacokinetic models in rat and human

Whole-body physiologically-based pharmacokinetic (PBPK) models have many applications in drug research and development. It is often necessary to inform these models with animal or clinical data, updating model parameters, and making the model more predictive for future applications. This provides an opportunity and a challenge given the large number of parameters of such models. The aim of this work was to propose new mechanistic model structures with reduced complexity allowing for parameter optimization. These models were evaluated for their ability to estimate realistic values for unbound tissue to plasma partition coefficients (Kpu) and simulate observed pharmacokinetic (PK) data. Two approaches are presented: using either established kinetic lumping methods based on tissue time constants (drug-dependent) or a novel clustering analysis to identify tissues sharing common Kpu values or Kpu scalars based on similarities of tissue composition (drug-independent). PBPK models derived from these approaches were assessed using PK data of diazepam in rats and humans. Although the clustering analysis produced minor differences in tissue grouping depending on the method used, two larger groups of tissues emerged. One including the kidneys, liver, spleen, gut, heart, and lungs, and another including bone, brain, muscle, and pancreas whereas adipose and skin were generally considered distinct. Overall, a subdivision into four tissue groups appeared most physiologically relevant in terms of tissue composition. Several models were found to have similar abilities to describe diazepam i.v. data as empirical models. Comparability of estimated Kpus to experimental Kpu values for diazepam was one criterion for selecting the appropriate PK model structure.

Complex Mixtures: Array PBPK Modeling of Jet Fuel Components

An array physiologically-based pharmacokinetic (PBPK) model represents a streamlined method to simultaneously quantify dosimetry of multiple compounds. To predict internal dosimetry of jet fuel components simultaneously, an array PBPK model was coded to simulate inhalation exposures to one or more selected compounds: toluene, ethylbenzene, xylenes, n-nonane, n-decane, and naphthalene. The model structure accounts for metabolism of compounds in the lung and liver, as well as kinetics of each compound in multiple tissues, including the cochlea and brain regions associated with auditory signaling (brainstem and temporal lobe). The model can accommodate either diffusion-limited or flow-limited kinetics (or a combination), allowing the same structure to be utilized for compounds with different characteristics. The resulting model satisfactorily simulated blood concentration and tissue dosimetry data from multiple published single chemical rat studies. The model was then utilized to predict tissue kinetics for the jet fuel hearing loss study (JTEH A, 25:1-14). The model was also used to predict rat kinetic comparisons between hypothetical exposures to JP-8 or a Virent Synthesized Aromatic Kerosene (SAK):JP-8 50:50 blend at the occupational exposure limit (200 mg/m³). The array model has proven useful for comparing potential tissue burdens resulting from complex mixture exposures.

Prediction of drug pro-arrhythmic cardiotoxicity using a semi-physiologically based pharmacokinetic model linked to cardiac ionic currents inhibition

Drug-induced torsades de pointes (TdP) risks are responsible for the withdrawal of many drugs from the market. Nowadays, assessments of drug-induced TdP risks are mainly based on maximum effective free therapeutic plasma concentration (EFTPCmax) and cardiac ionic current inhibitions using the human ventricular myocytes model (Tor-ORd model). Myocytes are targets of drug-induced TdP. The TdP risks may be directly linked to myocyte drug concentrations. We aimed to develop a semi-physiologically based pharmacokinetic (Semi-PBPK) model linked to cardiac ionic current inhibition (pharmacodynamics, PD) (Semi-PBPK-PD) to simultaneously predict myocyte drug concentrations and their TdP risks in humans. Alterations in action potential duration (ΔAPD90) were simulated using the Tor-ORd model and ionic current inhibition parameters based on myocyte or plasma drug concentrations. The predicted ΔAPD90 values were translated into in vivo alterations in QT interval(ΔQTc) induced by moxifloxacin, dofetilide, or sotalol. Myocyte drug concentrations of moxifloxacin, dofetilide, and sotalol gave better predictions of ΔQTc than plasma. Following validating the developed semi-PBPK-PD model, TdP risks of 37 drugs were assessed using ΔAPD90 and early afterdepolarization occurrence, which were estimated based on 10 × EFTPCmax and 10 × EFTMCmax (maximum effective free therapeutic myocyte concentration). 10 × EFTMCmax gave more sensitive and accurate predictions of pro-arrhythmic cardiotoxicity and the predicted TdP risks were also closer to clinic practice than 10 × EFTPCmax. In conclusion, pharmacokinetics and TdP risks of 37 drugs were successfully predicted using the semi-PBPK-PD model. Myocyte drug concentrations gave better predictions of ΔQTc and TdP risks than plasma.

Case study on the impact of the source of metabolism parameters in next generation physiologically based pharmacokinetic models: Implications for occupational exposures to trimethylbenzenes

Physiologically based pharmacokinetic (PBPK) models are a means of making important linkages between exposure assessment and in vitro toxicity. A key constraint on rapid application of PBPK models in risk assessment is traditional reliance on substance-specific in vivo toxicokinetic data to evaluate model quality. Bounding conditions, in silico, in vitro, and chemical read-across approaches have been proposed as alternative sources for metabolic clearance estimates. A case study to test consistency of predictive ability across these approaches was conducted using trimethylbenzenes (TMB) as prototype chemicals. Substantial concordance was found among TMB isomers with respect to accuracy (or inaccuracy) of approaches to estimating metabolism; for example, the bounding conditions never reproduced the human in vivo toxicokinetic data within two-fold. Using only approaches that gave acceptable prediction of in vivo toxicokinetics for the source compound (1,2,4-TMB) substantially narrowed the range of plausible internal doses for a given external dose for occupational, emergency response, and environmental/community health risk assessment scenarios for TMB isomers. Thus, risk assessments developed using the target compound models with a constrained subset of metabolism estimates (determined for source chemical models) can be used with greater confidence that internal dosimetry will be estimated with accuracy sufficient for the purpose at hand.

Utility of preclinical species for uncertainty assessment and correction of prediction of human volume of distribution using the Rodgers-Lukacova model

Prediction of rat, dog, monkey, and human volume of distribution (VD_ss) by Rodgers-Lukacova model was evaluated using a data set of more than 100 compounds.The prediction accuracy was best for humans followed by monkeys and dogs with 59, 52, and 41% of compounds within 2-fold, respectively.The accuracy of predictions in preclinical species was indicative of the human situation. This was particularly true for monkeys, where 87% of the compounds that were predicted within 2-fold in monkeys were also predicted within 2-fold in humans.The model's tendency to underestimate VD_ss was higher in rats and dogs compared to humans and monkeys for all ion classes but zwitterions. Hence, correction of human predictions using prediction errors in rats and dogs resulted in overestimation of VD_ss.The model had a similar degree of underestimation in humans and monkeys. Correction using monkeys improved the accuracy of the human estimate, especially for basic and zwitterion compounds.A strategy is proposed based on the accuracy of prediction in monkey and monkey scalars for prediction and prospective assessment of the accuracy of human VD_ss.

Development and intercomparison of single and multicompartment physiologically-based toxicokinetic models: Implications for model selection and tiered modeling frameworks

This study describes the development and intercomparison of generic physiologically-based toxicokinetic (PBTK) models for humans comprised of internally consistent one-compartment (1Co-) and multi-compartment (MCo-) implementations (G-PBTK). The G-PBTK models were parameterized for an adult male (70 kg) using common physiological parameters and in vitro biotransformation rate estimates and subsequently evaluated using independent concentration versus time data (n = 6) and total elimination half-lives (n = 15) for diverse organic chemicals. The model performance is acceptable considering the inherent uncertainty in the biotransformation rate data and the absence of model calibration. The G-PBTK model was then applied using hypothetical neutral organics, acidic ionizable organics and basic ionizable organics (IOCs) to identify combinations of partitioning properties and biotransformation rates leading to substantial discrepancies between 1Co- and MCo-PBTK calculations for whole body concentrations and half-lives. The 1Co- and MCo-PBTK model calculations for key toxicokinetic parameters are broadly consistent unless biotransformation is rapid (e.g., half-life less than five days). When half-lives are relatively short, discrepancies are greatest for the neutral organics and least for the acidic IOCs which follows from the estimated volumes of distribution (e.g., VD_SS = 9.6-15.4 L/kg vs 0.3-1.6 L/kg for the neutral and acidic compounds respectively) and the related approach to internal chemical equilibrium. The model intercomparisons demonstrate that 1Co-PBTK models can be applied with confidence to many exposure scenarios, particularly those focused on chronic or repeat exposures and for prioritization and screening-level decision contexts. However, MCo-PBTK models may be necessary in certain contexts, particularly for intermittent, short-term and highly variable exposures. A key recommendation to guide model selection and the development of tiered PBTK modeling frameworks that emerges from this study is the need to harmonize models with respect to parameterization and process descriptions to the greatest extent possible when proceeding from the application of simpler to more complex modeling tools as part of chemical assessment activities.

Impact of Obesity on Brexpiprazole Pharmacokinetics: Proposal for Improved Initiation of Treatment

Brexpiprazole is an oral antipsychotic agent indicated for use in patients with schizophrenia or as adjunctive treatment for major depressive disorder. As obesity (body mass index ≥35 kg/m² ) has the potential to affect drug pharmacokinetics and is a common comorbidity of both schizophrenia and major depressive disorder, it is important to understand changes in brexpiprazole disposition in this population. This study uses a whole-body physiologically based pharmacokinetic model to compare the pharmacokinetics of brexpiprazole in obese and normal-weight (body mass index 18-25 kg/m² ) individuals known to be cytochrome P450 2D6 extensive metabolizers (EMs) and poor metabolizers (PMs). The physiologically based pharmacokinetic simulations demonstrated significant differences in the time to effective concentrations between obese and normal-weight individuals within metabolizer groups according to the label-recommended titration. Simulations using an alternative dosing strategy of 1 week of twice-daily dosing in obese EMs or 2 weeks of twice-daily dosing in obese poor metabolizers, followed by a return to once-daily dosing, yielded more consistent plasma concentrations between normal-weight and obese patients without exceeding the area under the plasma concentration-time curve observed in the normal-weight EMs. These alternative dosing strategies reduce the time to effective concentrations in obese patients and may improve clinical response to brexpiprazole.

Brexpiprazole Pharmacokinetics in CYP2D6 Poor Metabolizers: Using Physiologically Based Pharmacokinetic Modeling to Optimize Time to Effective Concentrations

Brexpiprazole is an oral antipsychotic agent indicated for use in patients with schizophrenia, or as adjunctive treatment for major depressive disorder. As cytochrome P450 (CYP) 2D6 contributes significantly to brexpiprazole metabolism, there is a label-recommended 50% reduction in dose among patients with the CYP2D6 poor metabolizer phenotype. This study uses a whole-body physiologically based pharmacokinetic (PBPK) model to compare the pharmacokinetics of brexpiprazole in patients known to be extensive metabolizers (EMs) and poor metabolizers (PMs). A PBPK model was constructed, verified, and validated against brexpiprazole clinical data, and simulations of 500 subjects were performed to establish the median time to effective concentrations in EMs and PMs. The PBPK simulations captured brexpiprazole PK well and demonstrated significant differences in the time to effective concentrations between EMs and PMs according to the label-recommended titration. Additionally, these simulations suggest that CYP2D6 PMs consistently achieve lower minimum concentrations during the dosing interval than CYP2D6 EMs. Simulations using an alternative dosing strategy of twice-daily dosing (as opposed to once daily) in PMs during the first week of brexpiprazole dosing yielded more consistent plasma concentrations between EMs and PMs, without exceeding the area under the plasma concentration-time curve observed in the EMs. Taken together, the results of these PBPK simulations suggest that product labeling for brexpiprazole titration in CYP2D6 PMs likely overcompensates for the decreased clearance seen in this population. We propose an alternative dosing strategy that decreases the time to effective concentrations and recommend a reevaluation of steady-state PK in this population to potentially allow for higher daily doses in CYP2D6 PMs.

A hybrid modeling approach for assessing mechanistic models of small molecule partitioning in vivo using a machine learning-integrated modeling platform

Prediction of the first-in-human dosing regimens is a critical step in drug development and requires accurate quantitation of drug distribution. Traditional in vivo studies used to characterize clinical candidate’s volume of distribution are error-prone, time- and cost-intensive and lack reproducibility in clinical settings. The paper demonstrates how a computational platform integrating machine learning optimization with mechanistic modeling can be used to simulate compound plasma concentration profile and predict tissue-plasma partition coefficients with high accuracy by varying the lipophilicity descriptor logP. The approach applied to chemically diverse small molecules resulted in comparable geometric mean fold-errors of 1.50 and 1.63 in pharmacokinetic outputs for direct tissue:plasma partition and hybrid logP optimization, with the latter enabling prediction of tissue permeation that can be used to guide toxicity and efficacy dosing in human subjects. The optimization simulations required to achieve these results were parallelized on the AWS cloud and generated outputs in under 5 h. Accuracy, speed, and scalability of the framework indicate that it can be used to assess the relevance of other mechanistic relationships implicated in pharmacokinetic-pharmacodynamic phenomena with a lower risk of overfitting datasets and generate large database of physiologically-relevant drug disposition for further integration with machine learning models.

Comparison of Canine and Human Physiological Factors: Understanding Interspecies Differences that Impact Drug Pharmacokinetics

This review is a summary of factors affecting the drug pharmacokinetics (PK) of dogs versus humans. Identifying these interspecies differences can facilitate canine-human PK extrapolations while providing mechanistic insights into species-specific drug in vivo behavior. Such a cross-cutting perspective can be particularly useful when developing therapeutics targeting diseases shared between the two species such as cancer, diabetes, cognitive dysfunction, and inflammatory bowel disease. Furthermore, recognizing these differences also supports a reverse PK extrapolations from humans to dogs. To appreciate the canine-human differences that can affect drug absorption, distribution, metabolism, and elimination, this review provides a comparison of the physiology, drug transporter/enzyme location, abundance, activity, and specificity between dogs and humans. Supplemental material provides an in-depth discussion of certain topics, offering additional critical points to consider. Based upon an assessment of available state-of-the-art information, data gaps were identified. The hope is that this manuscript will encourage the research needed to support an understanding of similarities and differences in human versus canine drug PK.

Physiologically based metformin pharmacokinetics model of mice and scale-up to humans for the estimation of concentrations in various tissues

Metformin is the primary drug for type 2 diabetes treatment and a promising candidate for other disease treatment. It has significant deviations between individuals in therapy efficiency and pharmacokinetics, leading to the administration of an unnecessary overdose or an insufficient dose. There is a lack of data regarding the concentration-time profiles in various human tissues that limits the understanding of pharmacokinetics and hinders the development of precision therapies for individual patients. The physiologically based pharmacokinetic (PBPK) model developed in this study is based on humans’ known physiological parameters (blood flow, tissue volume, and others). The missing tissue-specific pharmacokinetics parameters are estimated by developing a PBPK model of metformin in mice where the concentration time series in various tissues have been measured. Some parameters are adapted from human intestine cell culture experiments. The resulting PBPK model for metformin in humans includes 21 tissues and body fluids compartments and can simulate metformin concentration in the stomach, small intestine, liver, kidney, heart, skeletal muscle adipose, and brain depending on the body weight, dose, and administration regimen. Simulations for humans with a bodyweight of 70kg have been analyzed for doses in the range of 500-1500mg. Most tissues have a half-life (T_1/2) similar to plasma (3.7h) except for the liver and intestine with shorter T_1/2 and muscle, kidney, and red blood cells that have longer T_1/2. The highest maximal concentrations (C_max) turned out to be in the intestine (absorption process) and kidney (excretion process), followed by the liver. The developed metformin PBPK model for mice does not have a compartment for red blood cells and consists of 20 compartments. The developed human model can be personalized by adapting measurable values (tissue volumes, blood flow) and measuring metformin concentration time-course in blood and urine after a single dose of metformin. The personalized model can be used as a decision support tool for precision therapy development for individuals.

Estimating the Equilibrium Distribution of Perfluoroalkyl Acids and 4 of Their Alternatives in Mammals

Perfluoroalkyl acids (PFAAs) mostly exist as ionic compounds that are of major concern because of their accumulative behavior. The discussion about their risk is ongoing considering the increasing production of structurally similar alternatives. We conducted model calculations based on equilibrium distribution coefficients that allow studying the distribution of PFAAs and their alternatives in various mammalian organs through comparison to empirical measurements in humans and rats. The calculations rely on experimentally determined distribution coefficients of a series of PFAAs and 4 of their alternatives to physiological matrices such as structural proteins, storage lipids, membrane lipids, albumin, and fatty acid binding protein (FABP). The relative sorption capacities in each organ were calculated from the combination of distribution coefficients and physiological data. The calculated distribution of PFAAs and alternatives within the organs showed that albumin and membrane lipids and, to a lesser extent, structural proteins have the highest relative sorption capacities for the compounds. Sorption to FABP is only relevant in the distribution of short-chain PFAAs. Storage lipids play a minor role in the distribution of all studied compounds. Our calculated distribution of PFAAs was evaluated by comparison to reported PFAA concentrations in various organs. Environ Toxicol Chem 2021;40:910-920. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Physiologically based pharmacokinetic (PBPK) modeling of RNAi therapeutics: Opportunities and challenges

Physiologically based pharmacokinetic (PBPK) modeling is a powerful tool with many demonstrated applications in various phases of drug development and regulatory review. RNA interference (RNAi)-based therapeutics are a class of drugs that have unique pharmacokinetic properties and mechanisms of action. With an increasing number of RNAi therapeutics in the pipeline and reaching the market, there is a considerable amount of active research in this area requiring a multidisciplinary approach. The application of PBPK models for RNAi therapeutics is in its infancy and its utility to facilitate the development of this new class of drugs is yet to be fully evaluated. From this perspective, we briefly discuss some of the current computational modeling approaches used in support of efficient development and approval of RNAi therapeutics. Considerations for PBPK model development are highlighted both in a relative context between small molecules and large molecules such as monoclonal antibodies and as it applies to RNAi therapeutics. In addition, the prospects for drawing upon other recognized avenues of PBPK modeling and some of the foreseeable challenges in PBPK model development for these chemical modalities are briefly discussed. Finally, an exploration of the potential application of PBPK model development for RNAi therapeutics is provided. We hope these preliminary thoughts will help initiate a dialogue between scientists in the relevant sectors to examine the value of PBPK modeling for RNAi therapeutics. Such evaluations could help standardize the practice in the future and support appropriate guidance development for strengthening the RNAi therapeutics development program.

Development of a physiologically based pharmacokinetic model of diisononyl phthalate (DiNP) in pregnant rat and human

A physiologically based pharmacokinetic (PBPK) model for di-isononyl phthalate (DiNP) was developed by adapting the existing models for di(2-ethylhexyl) phthalate (DEHP) and di-butylphthalate (DBP). Both pregnant rat and human time-course plasma and urine data were used to address the hydrolysis of DiNP in intestinal tract, plasma, and liver as well as hepatic oxidative metabolism and conjugation of the monoester and primary oxidative metabolites. Data in both rats and humans were available to inform the uptake and disposition of mono-isononyl phthalate (MiNP) as well as the three primary oxidative metabolites including hydroxy (7-OH)-, oxo (7-OXO)-, and carboxy (7-COX)-monoisononyl phthalate in plasma and urine. The DiNP model was reliable over a wide range of exposure levels in the pregnant rat as well as the two low exposure levels in humans including capturing the nonlinear behavior in the pregnant rat after repeated 750 mg/kg/day dosing. The presented DiNP PBPK model in pregnant rat and human, based upon an extensive kinetic dataset in both species, may provide a basis for assessing human equivalent exposures based upon either rodent or in vitro points of departure.

Interpretation of Drug Interaction Using Systemic and Local Tissue Exposure Changes

Systemic exposure of a drug is generally associated with its pharmacodynamic (PD) effect (e.g., efficacy and toxicity). In this regard, the change in area under the plasma concentration-time curve (AUC) of a drug, representing its systemic exposure, has been mainly considered in evaluation of drug-drug interactions (DDIs). Besides the systemic exposure, the drug concentration in the tissues has emerged as a factor to alter the PD effects. In this review, the status of systemic exposure, and/or tissue exposure changes in DDIs, were discussed based on the recent reports dealing with transporters and/or metabolic enzymes mediating DDIs. Particularly, the tissue concentration in the intestine, liver and kidney were referred to as important factors of PK-based DDIs.

Imbalance of Drug Transporter-CYP450s Interplay by Diabetes and Its Clinical Significance

The pharmacokinetics of a drug is dependent upon the coordinate work of influx transporters, enzymes and efflux transporters (i.e., transporter-enzyme interplay). The transporter-enzyme interplay may occur in liver, kidney and intestine. The influx transporters involving drug transport are organic anion transporting polypeptides (OATPs), peptide transporters (PepTs), organic anion transporters (OATs), monocarboxylate transporters (MCTs) and organic cation transporters (OCTs). The efflux transporters are P-glycoprotein (P-gp), multidrug/toxin extrusions (MATEs), multidrug resistance-associated proteins (MRPs) and breast cancer resistance protein (BCRP). The enzymes related to drug metabolism are mainly cytochrome P450 enzymes (CYP450s) and UDP-glucuronosyltransferases (UGTs). Accumulating evidence has demonstrated that diabetes alters the expression and functions of CYP450s and transporters in a different manner, disordering the transporter-enzyme interplay, in turn affecting the pharmacokinetics of some drugs. We aimed to focus on (1) the imbalance of transporter-CYP450 interplay in the liver, intestine and kidney due to altered expressions of influx transporters (OATPs, OCTs, OATs, PepTs and MCT6), efflux transporters (P-gp, BCRP and MRP2) and CYP450s (CYP3As, CYP1A2, CYP2E1 and CYP2Cs) under diabetic status; (2) the net contributions of these alterations in the expression and functions of transporters and CYP450s to drug disposition, therapeutic efficacy and drug toxicity; (3) application of a physiologically-based pharmacokinetic model in transporter-enzyme interplay.

Global Sensitivity Analysis of the Rodgers and Rowland Model for Prediction of Tissue: Plasma Partitioning Coefficients: Assessment of the Key Physiological and Physicochemical Factors That Determine Small-Molecule Tissue Distribution

In physiologically based pharmacokinetic (PBPK) modelling, the large number of input parameters, limited amount of available data and the structural model complexity generally hinder simultaneous estimation of uncertain and/or unknown parameters. These parameters are generally subject to estimation. However, the approaches taken for parameter estimation vary widely. Global sensitivity analyses are proposed as a method to systematically determine the most influential parameters that can be subject to estimation. Herein, a global sensitivity analysis was conducted to identify the key drug and physiological parameters influencing drug disposition in PBPK models and to potentially reduce the PBPK model dimensionality. The impact of these parameters was evaluated on the tissue-to-unbound plasma partition coefficients (Kpus) predicted by the Rodgers and Rowland model using Latin hypercube sampling combined to partial rank correlation coefficients (PRCC). For most drug classes, PRCC showed that LogP and fraction unbound in plasma (fu_p) were generally the most influential parameters for Kpu predictions. For strong bases, blood:plasma partitioning was one of the most influential parameter. Uncertainty in tissue composition parameters had a large impact on Kpu and Vss predictions for all classes. Among tissue composition parameters, changes in Kpu outputs were especially attributed to changes in tissue acidic phospholipid concentrations and extracellular protein tissue:plasma ratio values. In conclusion, this work demonstrates that for parameter estimation involving PBPK models and dimensionality reduction purposes, less influential parameters might be assigned fixed values depending on the parameter space, while influential parameters could be subject to parameters estimation.

Lipid, water, and protein composition to facilitate kinetic modeling of the auditory pathway

Environments combining JP-8 jet fuel exposure with heightened ambient noise may accelerate hearing loss induced by noise. To reduce animal use and facilitate kinetic modeling of this military aviation fuel, tissue-specific parameters are required, including water, protein, and lipid content. However, tissues involved in hearing, including cochlea, brainstem, frontal, and temporal lobe, have not been characterized before. Therefore, water content was determined by lyophilization of rat auditory tissues and the protein of the freeze dried remainder was quantified using a bicinchoninic acid assay. Lipids were extracted from fresh-frozen rat auditory tissues and separated into neutral lipids, free fatty acids, neutral phospholipids, and acidic phospholipids using solid phase extraction. Phospholipid fractions were confirmed by ³¹P nuclear magnetic resonance analysis showing distinct phospholipid profiles. Lipid content in reference tissues, such as kidney and adipose, confirmed literature values. For the first time, lipid content in the rat auditory pathway was determined showing that total lipid content was lowest in cochlea and highest in brainstem compared with frontal and temporal lobes. Auditory tissues displayed distinct lipid fraction profiles. The information on water, protein, and lipid composition is necessary to validate algorithms used in mathematical models and predict partitioning of chemicals of future interest into these tissues. This research may reduce the use of animals to measure partition coefficients for prospective physiological models.

Prediction of Tissue-to-Plasma Ratios of Basic Compounds in Mice

BACKGROUND

Majority of reported studies so far developed correlation regression equations using the rat muscle-to-plasma drug concentration ratio (Kp-muscle) to predict tissue-to-plasma drug concentration ratios (Kp-tissues). Use of regression equations derived from rat Kp-muscle may not be ideal to predict the mice tissue-Kps as there are species differences.

OBJECTIVES

(i) To develop the linear regression equations using mouse tissue-Kps; (ii) to assess the correlation between organ blood flow and/or organ weight with tissue-Kps and (iii) compare the observed tissue-Kps from mice with corresponding predicted tissue-Kps using Richter's rat-Kp specific equations.

METHOD

Disposition of 12 small molecules were investigated extensively in mouse plasma and tissues after a single oral dose administration. Linear correlation was assessed for each of the tissue with rest of the other tissues, separately for weak and strong bases.

RESULT

Newly developed regression equations using mice tissue-Kps, predicted 79% data points within twofold. As observed correlation r ² range was 0.75-0.98 between Kp-muscle and Kp-brain, -spleen, -skin, -liver, -lung, suggesting superior correlation between the tissue-Kps. Order of tissue-Kps, showed that tissue concentrations were directly proportional to the organ blood flow and inversely to the organ weight. Further, the observed tissue-Kps from mice were compared with corresponding predicted tissue-Kps using Richter's rat-Kp specific equations. Overall, 46, 54 and 63% data points were under predicted (<0.5-fold) for liver, spleen and lung, respectively. Whereas 63 and 75% data points were over predicted (>twofold) for skin and brain, respectively. These findings suggest that cross species extrapolation predictability is poor.

CONCLUSION

All these findings together suggest that mouse specific regression equations developed under controlled experimental conditions could be most appropriate for predicting mouse tissue-Kps for compounds with wide range of volume of distribution.

Evaluation and calibration of high-throughput predictions of chemical distribution to tissues

Toxicokinetics (TK) provides critical information for integrating chemical toxicity and exposure assessments in order to determine potential chemical risk (i.e., the margin between toxic doses and plausible exposures). For thousands of chemicals that are present in our environment, in vivo TK data are lacking. The publicly available R package "httk" (version 1.8, named for "high throughput TK") draws from a database of in vitro data and physico-chemical properties in order to run physiologically-based TK (PBTK) models for 553 compounds. The PBTK model parameters include tissue:plasma partition coefficients (Kp) which the httk software predicts using the model of Schmitt (Toxicol In Vitro 22 (2):457-467, 2008). In this paper we evaluated and modified httk predictions, and quantified confidence using in vivo literature data. We used 964 rat Kp measured by in vivo experiments for 143 compounds. Initially, predicted Kp were significantly larger than measured Kp for many lipophilic compounds (log10 octanol:water partition coefficient > 3). Hence the approach for predicting Kp was revised to account for possible deficiencies in the in vitro protein binding assay, and the method for predicting membrane affinity was revised. These changes yielded improvements ranging from a factor of 10 to nearly a factor of 10,000 for 83 Kp across 23 compounds with only 3 Kp worsening by more than a factor of 10. The vast majority (92%) of Kp were predicted within a factor of 10 of the measured value (overall root mean squared error of 0.59 on log10-transformed scale). After applying the adjustments, regressions were performed to calibrate and evaluate the predictions for 12 tissues. Predictions for some tissues (e.g., spleen, bone, gut, lung) were observed to be better than predictions for other tissues (e.g., skin, brain, fat), indicating that confidence in the application of in silico tools to predict chemical partitioning varies depending upon the tissues involved. Our calibrated model was then evaluated using a second data set of human in vivo measurements of volume of distribution (Vss) for 498 compounds reviewed by Obach et al. (Drug Metab Dispos 36(7):1385-1405, 2008). We found that calibration of the model improved performance: a regression of the measured values as a function of the predictions has a slope of 1.03, intercept of - 0.04, and R2 of 0.43. Through careful evaluation of predictive methods for chemical partitioning into tissues, we have improved and calibrated these methods and quantified confidence for TK predictions in humans and rats.

httk: R Package for High-Throughput Toxicokinetics

Thousands of chemicals have been profiled by high-throughput screening programs such as ToxCast and Tox21; these chemicals are tested in part because most of them have limited or no data on hazard, exposure, or toxicokinetics. Toxicokinetic models aid in predicting tissue concentrations resulting from chemical exposure, and a "reverse dosimetry" approach can be used to predict exposure doses sufficient to cause tissue concentrations that have been identified as bioactive by high-throughput screening. We have created four toxicokinetic models within the R software package httk. These models are designed to be parameterized using high-throughput in vitro data (plasma protein binding and hepatic clearance), as well as structure-derived physicochemical properties and species-specific physiological data. The package contains tools for Monte Carlo sampling and reverse dosimetry along with functions for the analysis of concentration vs. time simulations. The package can currently use human in vitro data to make predictions for 553 chemicals in humans, rats, mice, dogs, and rabbits, including 94 pharmaceuticals and 415 ToxCast chemicals. For 67 of these chemicals, the package includes rat-specific in vitro data. This package is structured to be augmented with additional chemical data as they become available. Package httk enables the inclusion of toxicokinetics in the statistical analysis of chemicals undergoing high-throughput screening.

In Vitro and in Silico Tools To Assess Extent of Cellular Uptake and Lysosomal Sequestration of Respiratory Drugs in Human Alveolar Macrophages

Accumulation of respiratory drugs in human alveolar macrophages (AMs) has not been extensively studied in vitro and in silico despite its potential impact on therapeutic efficacy and/or occurrence of phospholipidosis. The current study aims to characterize the accumulation and subcellular distribution of drugs with respiratory indication in human AMs and to develop an in silico mechanistic AM model to predict lysosomal accumulation of investigated drugs. The data set included 9 drugs previously investigated in rat AM cell line NR8383. Cell-to-unbound medium concentration ratio (K_p,cell) of all drugs (5 μM) was determined to assess the magnitude of intracellular accumulation. The extent of lysosomal sequestration in freshly isolated human AMs from multiple donors (n = 5) was investigated for clarithromycin and imipramine (positive control) using an indirect in vitro method (±20 mM ammonium chloride, NH₄Cl). The AM cell parameters and drug physicochemical data were collated to develop an in silico mechanistic AM model. Three in silico models differing in their description of drug membrane partitioning were evaluated; model (1) relied on octanol-water partitioning of drugs, model (2) used in vitro data to account for this process, and model (3) predicted membrane partitioning by incorporating AM phospholipid fractions. In vitro K_p,cell ranged >200-fold for respiratory drugs, with the highest accumulation seen for clarithromycin. A good agreement in K_p,cell was observed between human AMs and NR8383 (2.45-fold bias), highlighting NR8383 as a potentially useful in vitro surrogate tool to characterize drug accumulation in AMs. The mean K_p,cell of clarithromycin (81, CV = 51%) and imipramine (963, CV = 54%) were reduced in the presence of NH₄Cl by up to 67% and 81%, respectively, suggesting substantial contribution of lysosomal sequestration and intracellular binding in the accumulation of these drugs in human AMs. The in vitro data showed variability in drug accumulation between individual human AM donors due to possible differences in lysosomal abundance, volume, and phospholipid content, which may have important clinical implications. Consideration of drug-acidic phospholipid interactions significantly improved the performance of the in silico models; use of in vitro K_p,cell obtained in the presence of NH₄Cl as a surrogate for membrane partitioning (model (2)) captured the variability in clarithromycin and imipramine K_p,cell observed in vitro and showed the best ability to predict correctly positive and negative lysosomotropic properties. The developed mechanistic AM model represents a useful in silico tool to predict lysosomal and cellular drug concentrations based on drug physicochemical data and system specific properties, with potential application to other cell types.

Prediction of Deoxypodophyllotoxin Disposition in Mouse, Rat, Monkey, and Dog by Physiologically Based Pharmacokinetic Model and the Extrapolation to Human

Deoxypodophyllotoxin (DPT) is a potential anti-tumor candidate prior to its clinical phase. The aim of the study was to develop a physiologically based pharmacokinetic (PBPK) model consisting of 13 tissue compartments to predict DPT disposition in mouse, rat, monkey, and dog based on in vitro and in silico inputs. Since large interspecies difference was found in unbound fraction of DPT in plasma, we assumed that K_t:pl,u (unbound tissue-to-plasma concentration ratio) was identical across species. The predictions of our model were then validated by in vivo data of corresponding preclinical species, along with visual predictive checks. Reasonable matches were found between observed and predicted plasma concentrations and pharmacokinetic parameters in all four animal species. The prediction in the related seven tissues of mouse was also desirable. We also attempted to predict human pharmacokinetic profile by both the developed PBPK model and interspecies allometric scaling across mouse, rat and monkey, while dog was excluded from the scaling. The two approaches reached similar results. We hope the study will help in the efficacy and safety assessment of DPT in future clinical studies and provide a reference to the preclinical screening of similar compounds by PBPK model.

Developing a Physiologically-Based Pharmacokinetic Model Knowledgebase in Support of Provisional Model Construction

Developing physiologically-based pharmacokinetic (PBPK) models for chemicals can be resource-intensive, as neither chemical-specific parameters nor in vivo pharmacokinetic data are easily available for model construction. Previously developed, well-parameterized, and thoroughly-vetted models can be a great resource for the construction of models pertaining to new chemicals. A PBPK knowledgebase was compiled and developed from existing PBPK-related articles and used to develop new models. From 2,039 PBPK-related articles published between 1977 and 2013, 307 unique chemicals were identified for use as the basis of our knowledgebase. Keywords related to species, gender, developmental stages, and organs were analyzed from the articles within the PBPK knowledgebase. A correlation matrix of the 307 chemicals in the PBPK knowledgebase was calculated based on pharmacokinetic-relevant molecular descriptors. Chemicals in the PBPK knowledgebase were ranked based on their correlation toward ethylbenzene and gefitinib. Next, multiple chemicals were selected to represent exact matches, close analogues, or non-analogues of the target case study chemicals. Parameters, equations, or experimental data relevant to existing models for these chemicals and their analogues were used to construct new models, and model predictions were compared to observed values. This compiled knowledgebase provides a chemical structure-based approach for identifying PBPK models relevant to other chemical entities. Using suitable correlation metrics, we demonstrated that models of chemical analogues in the PBPK knowledgebase can guide the construction of PBPK models for other chemicals.

The Role of Extracellular Binding Proteins in the Cellular Uptake of Drugs: Impact on Quantitative In Vitro-to-In Vivo Extrapolations of Toxicity and Efficacy in Physiologically Based Pharmacokinetic-Pharmacodynamic Research

A critical component in the development of physiologically based pharmacokinetic-pharmacodynamic (PBPK/PD) models for estimating target organ dosimetry in pharmacology and toxicology studies is the understanding of the uptake kinetics and accumulation of drugs and chemicals at the cellular level. Therefore, predicting free drug concentrations in intracellular fluid will contribute to our understanding of concentrations at the site of action in cells in PBPK/PD research. Some investigators believe that uptake of drugs in cells is solely driven by the unbound fraction; conversely, others argue that the protein-bound fraction contributes a significant portion of the total amount delivered to cells. Accordingly, the current literature suggests the existence of a so-called albumin-mediated uptake mechanism(s) for the protein-bound fraction (i.e., extracellular protein-facilitated uptake mechanisms) at least in hepatocytes and cardiac myocytes; however, such mechanism(s) and cells from other organs deserve further exploration. Therefore, the main objective of this present study was to discuss further the implication of potential protein-facilitated uptake mechanism(s) on drug distribution in cells under in vivo conditions. The interplay between the protein-facilitated uptake mechanism(s) and the effects of a pH gradient, metabolism, transport, and permeation limitation potentially occurring in cells was also discussed, as this should violate the basic assumption on similar free drug concentration in cells and plasma. This was made because the published equations used to calculate drug concentrations in cells in a PBPK/PD model did not consider potential protein-facilitated uptake mechanism(s). Consequently, we corrected some published equations for calculating the free drug concentrations in cells compared with plasma in PBPK/PD modeling studies, and we proposed a refined strategy for potentially performing more accurate quantitative in vitro-to-in vivo extrapolations (IVIVEs) of toxicity (efficacy) at the cellular level from data generated in cell assays. Overall, this present study may help to optimize the human dose prediction in preclinical and clinical studies, while prescribing drugs with narrow therapeutic windows that are highly bound to extracellular proteins and/or highly ionized at the physiological pH. This may facilitate building a more accurate safety (efficacy) profile for such drugs.

A paradigm shift in pharmacokinetic-pharmacodynamic (PKPD) modeling: rule of thumb for estimating free drug level in tissue compared with plasma to guide drug design

A basic assumption in pharmacokinetics-pharmacodynamics research is that the free drug concentration is similar in plasma and tissue, and, hence, in vitro plasma data can be used to estimate the in vivo condition in tissue. However, in a companion manuscript, it has been demonstrated that this assumption is violated for the ionized drugs. Nonetheless, these observations focus on in vitro static environments and do not challenge data with an in vivo dynamic system. Therefore, an extension from an in vitro to an in vivo system becomes the necessary next step. The objective of this study was to perform theoretical simulations of the free drug concentration in tissue and plasma by using a physiologically based pharmacokinetics (PBPK) model reproducing the in vivo conditions in human. Therefore, the effects of drug ionization, lipophilicity, and clearance have been taken into account in a dynamic system. This modeling exercise was performed as a proof of concept to demonstrate that free drug concentration in tissue and plasma may also differ in a dynamic system for passively permeable drugs that are ionized at the physiological pH. The PBPK model simulations indicated that free drug concentrations in tissue cells and plasma significantly differ for the ionized drugs because of the pH gradient effect between cells and interstitial space. Hence, a rule of thumb for potentially performing more accurate PBPK/PD modeling is suggested, which states that the free drug concentration in tissue and plasma will differ for the ionizable drugs in contrast to the neutral drugs. In addition to the pH gradient effect for the ionizable drugs, lipophilicity and clearance effects will increase or decrease the free drug concentration in tissue and plasma for each class of drugs; thus, higher will be the drug lipophilicity and clearance, lower would be the free drug concentration in plasma, and, hence, in tissue, in a dynamic in vivo system. Therefore, only considering the value of free fraction in plasma derived from a static in vitro environment might be biased to guide drug design (the old paradigm), and, hence, it is recommended to use a PBPK model to reproduce more accurately the in vivo condition in tissue (the new paradigm). This newly developed approach can be used to predict free drug concentration in diverse tissue compartments for small molecules in toxicology and pharmacology studies, which can be leveraged to optimize the pharmacokinetics drivers of tissue distribution based upon physicochemical and physiological input parameters in an attempt to optimize free drug level in tissue. Overall, this present study provides guidance on the application of plasma and tissue concentration information in PBPK/PD research in preclinical and clinical studies, which is in accordance with the recent literature.

Drug Distribution to Human Tissues: Prediction and Examination of the Basic Assumption in In Vivo Pharmacokinetics-Pharmacodynamics (PK/PD) Research

The tissue:plasma partition coefficients (Kp ) are good indicators of the extent of tissue distribution. Therefore, advanced tissue composition-based models were used to predict the Kp values of drugs under in vivo conditions on the basis of in vitro and physiological input data. These models, however, focus on animal tissues and do not challenge the predictions with human tissues for drugs. The first objective of this study was to predict the experimentally determined Kp values of seven human tissues for 26 drugs. In all, 95% of the predicted Kp values are within 2.5-fold error of the observed values in humans. Accordingly, these results suggest that the tissue composition-based model used in this study is able to provide accurate estimates of drug partitioning in the studied human tissues. Furthermore, as the Kp equals to the ratio of total concentration between tissue and plasma, or the ratio of unbound fraction between plasma (fup ) and tissue (fut ), this parameter Kp would deviate from the unity. Therefore, the second objective was to examine the corresponding relationships between fup and fut values experimentally determined in humans for several drugs. The results also indicate that fup may significantly deviate to fut ; the discrepancies are governed by the dissimilarities in the binding and ionization on both sides of the membrane, which were captured by the tissue composition-based model. Hence, this violated the basic assumption in in vivo pharmacokinetics-pharmacodynamics (PK/PD) research, since the free drug concentration in tissue and plasma was not equal particularly for the ionizable drugs due to the pH gradient effect on the fraction of unionized drug in plasma (fuip ) and tissue (fuit ) (i.e., fup × fuip × total plasma concentration = fut × fuit × total tissue concentration, and, hence, the free drug concentration in plasma and tissue differed by fuip/fuit). Therefore, this assumption should be adjusted for the ionized drugs, and, hence, a mathematical correction to the basic assumption of similar free drug concentration in plasma and tissues can be derived from the tissue composition-based model. Note that this assumption will be further challenged in a dynamic in vivo system in a companion manuscript. Overall, this study was a first attempt to predict the in vivo Kp values for specific human tissues by considering separately the effect of fup and fut , with the aim of facilitating the use of physiologically-based PK (PBPK) model in PK/PD studies.