Prediction of Fraction Unbound in Microsomal and Hepatocyte Incubations: A Comparison of Methods across Industry Datasets

ADME characterization and PBK model development of 3 highly protein-bound UV filters through topical application

Estimating human exposure in the safety assessment of chemicals is crucial. Physiologically based kinetic (PBK) models which combine information on exposure, physiology, and chemical properties, describing the absorption, distribution, metabolism, and excretion (ADME) processes of a chemical, can be used to calculate internal exposure metrics such as maximum concentration and area under the concentration-time curve in plasma or tissues of a test chemical in next-generation risk assessment. This article demonstrates the development of PBK models for 3 UV filters, specifically octyl methoxycinnamate, octocrylene, and 4-methylbenzylidene camphor. The models were parameterized entirely based on data obtained from in vitro and/or in silico methods in a bottom-up modeling approach and then validated based on human dermal pharmacokinetic (PK) data. The 3 UV filters are "difficult to test" in in vitro test systems due to high lipophilicity, high binding affinity for proteins, and nonspecific binding, for example, toward plastic. This research work presents critical considerations in ADME data generation, interpretation, and parameterization to assure valid PBK model development to increase confidence in using PBK modeling to help make safety decisions in the absence of human PK data. The developed PBK models of the 3 chemicals successfully simulated the plasma concentration profiles of clinical PK data following dermal application, indicating the reliability of the ADME data generated and the parameters determined. The study also provides insights and lessons learned for characterizing ADME and developing PBK models for highly lipophilic and protein-bound chemicals in the future.

Plasma Protein-Mediated Uptake and Contradictions to the Free Drug Hypothesis: A Critical Review

According to the free drug hypothesis (FDH), only free, unbound drug is available to interact with biological targets. This hypothesis is the fundamental principle that continues to explain the vast majority of all pharmacokinetic and pharmacodynamic processes. Under the FDH, the free drug concentration at the target site is considered the driver of pharmacodynamic activity and pharmacokinetic processes. However, deviations from the FDH are observed in hepatic uptake and clearance predictions, where observed unbound intrinsic hepatic clearance (CL_int,u) is larger than expected. Such deviations are commonly observed when plasma proteins are present and form the basis of the so-called plasma protein-mediated uptake effect (PMUE). This review will discuss the basis of plasma protein binding as it pertains to hepatic clearance based on the FDH, as well as several hypotheses that may explain the underlying mechanisms of PMUE. Notably, some, but not all, potential mechanisms remained aligned with the FDH. Finally, we will outline possible experimental strategies to elucidate PMUE mechanisms. Understanding the mechanisms of PMUE and its potential contribution to clearance underprediction is vital to improving the drug development process.

Physiologically Based Pharmacokinetic Modeling To Guide Management of Drug Interactions between Elexacaftor-Tezacaftor-Ivacaftor and Antibiotics for the Treatment of Nontuberculous Mycobacteria

Nontuberculous mycobacteria (NTM) are the pathogens of concern in people with cystic fibrosis (pwCF) due to their association with deterioration of lung function. Treatment requires the use of a multidrug combination regimen, creating the potential for drug-drug interactions (DDIs) with cystic fibrosis transmembrane conductance regulator (CFTR)-modulating therapies, including elexacaftor, tezacaftor, and ivacaftor (ETI), which are eliminated mainly through cytochrome P450 (CYP) 3A-mediated metabolism. An assessment of the DDI risk for ETI coadministered with NTM treatments, including rifabutin, clofazimine, and clarithromycin, is needed to provide appropriate guidance on dosing.

Can We Predict Clinical Pharmacokinetics of Highly Lipophilic Compounds by Integration of Machine Learning or In Vitro Data into Physiologically Based Models? A Feasibility Study Based on 12 Development Compounds

While high lipophilicity tends to improve potency, its effects on pharmacokinetics (PK) are complex and often unfavorable. To predict clinical PK in early drug discovery, we built human physiologically based PK (PBPK) models integrating either (i) machine learning (ML)-predicted properties or (ii) discovery stage in vitro data. Our test set was composed of 12 challenging development compounds with high lipophilicity (mean calculated log P 4.2), low plasma-free fraction (50% of compounds with fu,p < 1%), and low aqueous solubility. Predictions focused on key human PK parameters, including plasma clearance (CL), volume of distribution at steady state (V_ss), and oral bioavailability (%F). For predictions of CL, the ML inputs showed acceptable accuracy and slight underprediction bias [an average absolute fold error (AAFE) of 3.55; an average fold error (AFE) of 0.95]. Surprisingly, use of measured data only slightly improved accuracy but introduced an overprediction bias (AAFE = 3.35; AFE = 2.63). Predictions of V_ss were more successful, with both ML (AAFE = 2.21; AFE = 0.90) and in vitro (AAFE = 2.24; AFE = 1.72) inputs showing good accuracy and moderate bias. The %F was poorly predicted using ML inputs [average absolute prediction error (AAPE) of 45%], and use of measured data for solubility and permeability improved this to 34%. Sensitivity analysis showed that predictions of CL limited the overall accuracy of human PK predictions, partly due to high nonspecific binding of lipophilic compounds, leading to uncertainty of unbound clearance. For accurate predictions of %F, solubility was the key factor. Despite current limitations, this work encourages further development of ML models and integration of their results within PBPK models to enable human PK prediction at the drug design stage, even before compounds are synthesized. Further evaluation of this approach with more diverse chemical types is warranted.

Non-specific binding of compounds in in vitro metabolism assays: a comparison of microsomal and hepatocyte binding in different species and an assessment of the accuracy of prediction models

Non-specific binding in in vitro metabolism systems leads to an underestimation of the true intrinsic metabolic clearance of compounds being studied. Therefore in vitro binding needs to be accounted for when extrapolating in vitro data to predict the in vivo metabolic clearance of a compound. While techniques exist for experimentally determining the fraction of a compound unbound in in vitro metabolism systems, early in drug discovery programmes computational approaches are often used to estimate the binding in the in vitro system.Experimental fraction unbound data (n = 60) were generated in liver microsomes (fu_mic) from five commonly used pre-clinical species (rat, mouse, dog, minipig, monkey) and humans. Unbound fraction in incubations with mouse, rat or human hepatocytes was determined for the same 60 compounds. These data were analysed to determine the relationship between experimentally determined binding in the different matrices and across different species. In hepatocytes there was a good correlation between fraction unbound in human and rat (r²=0.86) or mouse (r²=0.82) hepatocytes. Similar correlations were observed between binding in human liver microsomes and microsomes from rat, mouse, dog, Göttingen minipig or monkey liver microsomes (r² of >0.89, n = 51 - 52 measurements in different species). Physicochemical parameters (logP, pKa and logD) were predicted for all evaluated compounds. In addition, logP and/or logD were measured for a subset of compounds.Binding to human hepatocytes predicted using 5 different methods was compared to the measured data for a set of 59 compounds. The best methods evaluated used measured microsomal binding in human liver microsomes to predict hepatocyte binding. The collated physicochemical data were used to predict the human fu_mic using four different in silico models for a set of 53-60 compounds. The correlation (r²) and root mean square error between predicted and observed microsomal binding was 0.69 & 0.20, 0.47 & 0.23, 0.56 & 0.21 and 0.54 & 0.26 for the Turner-Simcyp, Austin, Hallifax-Houston and Poulin models, respectively. These analyses were extended to include measured literature values for binding in human liver microsomes for a larger set of compounds (n=697). For the larger dataset of compounds, microsomal binding was well predicted for neutral compounds (r²=0.67 - 0.70) using the Poulin, Austin, or Turner-Simcyp methods but not for acidic or basic compounds (r²<0.5) using any of the models. While the lipophilicity-based models can be used, the in vitro binding should be measured for compounds where more certainty is needed, using appropriately calibrated assays and possibly established weak, moderate, and strong binders as reference compounds to allow comparison across databases.

Drug clearance by aldehyde oxidase: can we avoid clinical failure?

Despite increased awareness of aldehyde oxidase (AO) as a major drug-metabolising enzyme, predicting the pharmacokinetics of its substrates remains challenging. Several drug candidates have been terminated due to high clearance, which were subsequently discovered to be AO substrates. Even retrospective extrapolation of human clearance, from models more sensitive to AO activity, often resulted in underprediction.The questions of the current work thus were: Is there an acceptable degree of in vitro AO metabolism that does not result in high in vivo human clearance? And, if so, how can this be predicted?We built an in vitro/in vivo correlation using known AO substrates, combining multiple in vitro parameters to calculate the blood metabolic clearance mediated by AO (CLb_AO). This value was compared with observed blood clearance (CLb_-obs), establishing cut-off CLb_AO values, to discriminate between low and high CLb_-obs. The model was validated using additional literature compounds, and CLb_-obs was predicted in the correct category.This simple, categorical, semi-quantitative yet multi-factorial model is readily applicable in drug discovery. Further, it is valuable for high-clearance compounds, as it predicts the CLb group, rather than an exact CLb value, for the substrates of this poorly-characterised enzyme.

An HPLC-UV validated bioanalytical method for measurement of in vitro phase 1 kinetics of α-synuclein binding bifunctional compounds

Aggregates of the protein α-synuclein are associated with pathophysiology of Parkinson's disease and are present in Lewy Bodies found in the brains of Parkinson's patients. We previously demonstrated that bifunctional compounds composed of caffeine linked via a six carbon chain to either 1-aminoindan (C₈-6-I) or nicotine (C₈-6-N) bind α-synuclein and protect yeast cells from α-synuclein mediated toxicity.A critical step in development of positron emission tomography (PET) probes for neurodegenerative diseases is evaluation of their metabolic stability. We determined that C₈-6-I, and C₈-6-N both undergo phase 1 P450 metabolism in mouse, rat, and human liver microsomes. We utilised this metabolic information to guide the design of fluorinated analogues for use as PET probes and determined that the fluorine in ¹⁹F-C₈-6-I and ¹⁹F-C₈-6-N is stable to P450 enzymes.We have developed and validated an analytical HPLC-UV method following FDA and EMA guidelines to measure in vitro phase 1 kinetics of these compounds and determine their V_max, K_M and CL_int,u in mouse liver microsomes. We found that C₈-6-I and ¹⁹F-C₈-6-I have a two- to fourfold lower CL_int,u than C₈-6-N, and ¹⁹F-C₈-6-N. Our approach shows a simple, specific, and effective system to design and develop compounds as PET probes.

Application of empirical scalars to enable early prediction of human hepatic clearance using IVIVE in drug discovery: an evaluation of 173 drugs

The utilization of in vitro data to predict drug pharmacokinetics (PK) in vivo has been a consistent practice in early drug discovery for decades. However, its success is hampered by mispredictions attributed to uncharacterized biological phenomena/experimental artifacts. Predicted drug clearance (CL) from experimental data (i.e. hepatocyte intrinsic clearance: CLint, fraction unbound in plasma: fu,p) is often systematically underpredicted using the well-stirred model (WSM). The objective of this study was to evaluate using empirical scalars in the WSM to correct for CL mispredictions. Drugs (N=28) were used to generate numerical scalars on CLint (α), and fu,p (β) to minimize the error (AAFE) for CL predictions. These scalars were validated using an additional dataset (N=28 drugs) and applied to a non-redundant AstraZeneca (AZ) dataset available in the literature (N=117 drugs) for a total of 173 compounds. CL predictions using the WSM were improved for most compounds using an α value of 3.66 (~64%<2-fold) compared to no scaling (~46%<2-fold). Similarly, using a β value of 0.55 or combination of α and β scalars (values of 1.74 and 0.66, respectively) resulted in a similar improvement in predictions (~64%<2-fold and ~65%<2-fold, respectively). For highly bound compounds (fu,p{less than or equal to}0.01), AAFE was substantially reduced across all scaling methods. Using the β scalar alone or a combination of α and β appeared optimal; and produce larger magnitude corrections for highly-bound compounds. Some drugs are still disproportionally mispredicted, however the improvements in prediction error and simplicity of applying these scalars suggests its utility for early-stage CL predictions. Significance Statement In early drug discovery, prediction of human clearance using in vitro experimental data plays an essential role in triaging compounds prior to in vivo studies. These predictions have been systematically underestimated. Here we introduce empirical scalars calibrated on the extent of plasma protein binding that appear to improve clearance prediction across multiple datasets. This approach can be used in early phases of drug discovery prior to the availability of pre-clinical data for early quantitative predictions of human clearance.

Opportunities and Considerations in the Application of Artificial Intelligence to Pharmacokinetic Prediction

The improvement in the ability of the pharmaceutical industry to predict human pharmacokinetic behavior are attributable to major technological shifts from 1990 to the present day. The opportunity for the application of AI/ML based approaches in the pharmaceutical industry is driven by the abundance of data sets that exist within individual pharmaceutical and biotech companies and the availability, within these environments, of abundant computing power. This chapter seeks to describe opportunities for artificial intelligence to contribute to the assessment and evaluation of the dug metabolism and pharmacokinetic (DMPK) properties of novel compounds across the drug discovery and development continuum. Many initiatives are already underway with respect to the application of AI/ML in predicting pharmacokinetic profiles so the question is not whether AI will influence pharmacokinetic prediction but rather how to best utilize and incorporate this and how to evaluate the value added from these applications. Since our understanding of the underlying biology of the in vitro and in vivo systems with respect to ADME, one of the key challenges to AI-based methods will be the ability to adapt to data sets that change in quality over time.

HLM-beads: Rapid Assessment of Nonspecific Binding to Human Liver Microsomes Using Magnetizable Beads

In early drug development, drug-drug interaction risk is routinely assessed using human liver microsomes (HLMs). Nonspecific binding of drugs to HLMs can affect the determination of accurate enzyme parameters (K_m, K_i, K_I). Previously, we described a novel in vitro model consisting of HLMs bound to magnetizable beads [HLM-magnetizable-beads system (HLM-beads)]. The HLM-beads enable rapid separation of HLMs from incubation media by applying a magnetic field. Here, HLM-beads were further characterized and evaluated as a tool to assess HLM nonspecific binding of small molecules. The free fractions (f_u,mic) of 13 compounds (chosen based on their pKa values) were determined using HLM-beads under three HLM concentrations (0.025, 0.50, and 1.0 mg/ml) and compared with those determined by equilibrium dialysis. Most f_u,mic values obtained using HLM-beads were within 0.5- to 2-fold of the values determined using equilibrium dialysis. The highest fold difference were observed for high binders itraconazole and BIRT2584 (1.9- to 2.9-fold), as the pronounced adsorption of these compounds to the equilibrium dialysis apparatus interfered with their f_u,mic determination. Correlation and linear regression analysis of the f_u,mic values generated using HLM-beads and equilibrium dialysis was conducted. Overall, a good correlation of f_u,mic values obtained by the two methods were observed, as the r and R² values from correlational analysis and linear regression analysis were >0.9 and >0.89, respectively. These studies demonstrate that HLM-beads can produce comparable f_u,mic values as determined by equilibrium dialysis while reducing the time required for this type of study from hours to only 10 minutes and compound apparatus adsorption. SIGNIFICANCE STATEMENT: This work introduces a new method of rapidly assessing nonspecific microsomal binding using human liver microsomes bound to magnetizable beads.

Advances in Applying Computer-Aided Drug Design for Neurodegenerative Diseases

Neurodegenerative diseases (NDs) including Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, and Huntington's disease are incurable and affect millions of people worldwide. The development of treatments for this unmet clinical need is a major global research challenge. Computer-aided drug design (CADD) methods minimize the huge number of ligands that could be screened in biological assays, reducing the cost, time, and effort required to develop new drugs. In this review, we provide an introduction to CADD and examine the progress in applying CADD and other molecular docking studies to NDs. We provide an updated overview of potential therapeutic targets for various NDs and discuss some of the advantages and disadvantages of these tools.

Pharmacokinetics of 2-phenoxyethanol and its major metabolite, phenoxyacetic acid, after dermal and inhaled routes of exposure: application to development PBPK model in rats

2-Phenoxyethanol (PE), ethylene glycol monophenyl ether, is widely used as a preservative in cosmetic products as well as in non-cosmetics. Since PE has been used in many types of products, it can be absorbed via dermal or inhaled route for systemic exposures. In this study, the pharmacokinetic (PK) studies of PE and its major metabolite, phenoxyacetic acid (PAA), after dermal (30 mg and 100 mg) and inhaled administration (77 mg) of PE in rats were performed. PE was administered daily for 4 days and blood samples were collected at day 1 and day 4 for PK analysis. PE was rapidly absorbed and extensively metabolized to form PAA. After multiple dosing, the exposures of PE and PAA were decreased presumably due to the induction of metabolizing enzymes of PE and PAA. In dermal mass balance study using [¹⁴C]-phenoxyethanol ([¹⁴C]PE) as a microtracer, most of the PE and its derivatives were excreted in urine (73.03%) and rarely found in feces (0.66%). Based on these PK results, a whole-body physiologically-based pharmacokinetic (PBPK) model of PE and PAA after dermal application and inhalation in rats was successfully developed. Most of parameters were obtained from the literatures and experiments, and intrinsic clearance at steady-state (CL_int,ss) were optimized based on the observed multiple PK data. With the developed model, systemic exposures of PE and PAA after dermal application and inhalation were simulated following no-observed-adverse-effect level (NOAEL) of 500 mg/kg/day for dermal application and that of 12.7 mg/kg/day for inhalation provided by the Environmental Protection Agency. The area under the concentration-time curve at steady state (AUC_ss) in kidney and liver (and lung for inhalations), which are known target organs of exhibiting toxicity of PE, as well as AUC_ss in plasma of PE and PAA were obtained from the model.

Primary Human Hepatocyte Spheroid Model as a 3D In Vitro Platform for Metabolism Studies

3D cultures of primary human hepatocytes (PHH) are emerging as a more in vivo-like culture system than previously available hepatic models. This work describes the characterisation of drug metabolism in 3D PHH spheroids. Spheroids were formed from three different donors of PHH and the expression and activities of important cytochrome P450 enzymes (CYP1A2, 2B6, 2C9, 2D6, and 3A4) were maintained for up to 21 days after seeding. The activity of CYP2B6 and 3A4 decreased, while the activity of CYP2C9 and 2D6 increased over time (P < 0.05). For six test compounds, that are metabolised by multiple enzymes, intrinsic clearance (CL_int) values were comparable to standard in vitro hepatic models and successfully predicted in vivo CL_int within 3-fold from observed values for low clearance compounds. Remarkably, the metabolic turnover of these low clearance compounds was reproducibly measured using only 1-3 spheroids, each composed of 2000 cells. Importantly, metabolites identified in the spheroid cultures reproduced the major metabolites observed in vivo, both primary and secondary metabolites were captured. In summary, the 3D PHH spheroid model shows promise to be used in drug discovery projects to study drug metabolism, including unknown mechanisms, over an extended period of time.

Direct Comparison of Total Clearance Prediction: Computational Machine Learning Model versus Bottom-Up Approach Using In Vitro Assay

The in vitro-in vivo extrapolation (IVIVE) approach for predicting total plasma clearance (CL_tot) has been widely used to rank order compounds early in discovery. More recently, a computational machine learning approach utilizing physicochemical descriptors and fingerprints calculated from chemical structure information has emerged, enabling virtual predictions even earlier in discovery. Previously, this approach focused more on in vitro intrinsic clearance (CL_int) prediction. Herein, we directly compare these two approaches for predicting CL_tot in rats. A structurally diverse set of 1114 compounds with known in vivo CL_tot, in vitro CL_int, and plasma protein binding was used as the basis for this evaluation. The machine learning models were assessed by validation approaches using the time- and cluster-split training and test sets, and five-fold cross validation. Assessed by five-fold validation, the random forest regression (RF) and radial basis function (RBF) models demonstrated better prediction performance in eight attempted machine learning models. The CL_tot values predicted by the RF and RBF models were within two-fold of the observed values for 67.7 and 71.9% of cluster-split test set compounds, respectively, while the predictivity was worse in the time-split dataset. The predictivity of both models tended to be improved by incorporating in vitro parameters, unbound fraction in plasma (f_u,p), and CL_int. CL_tot prediction utilizing in vitro CL_int and the well-stirred model, correcting for the fraction unbound in blood, was substantially worse compared to machine learning approaches for the same cluster-split test set. The reason that CL_tot is underestimated by IVIVE is not fully explained by considering the calculated microsomal unbound fraction (cf_u,mic), extended clearance classification system (ECCS), and omitting high clearance compounds in excess of hepatic blood flow. The analysis suggests that in silico machine learning models may have the power to reduce reliance on or replace in vitro and in vivo studies for chemical structure optimization in early drug discovery.

In vitro-in vivo extrapolation of metabolic clearance using human liver microsomes: factors showing variability and their normalization

In vitro-in vivo extrapolation (IVIVE) using human liver microsomes has been widely used to predict metabolic clearance, but some of the factors used in the process of prediction show variability for the same compound: notably, microsomal intrinsic clearance values corrected by the unbound fraction (CL_{int, u}), physiological parameters used for scale-up, and the source of in vivo clearance data.The purpose of this study was to assess the correlation between in vitro and in vivo CL_int with a focus on factors showing variability using four cytochrome P450 (CYP)3A substrates.We surveyed in vivo clearance values in literature and also determined the microsomal CL_{int, u} values. A scaling factor (SF_direct) was defined as in vivo CL_int divided by the microsomal CL_{int, u}, which ranged from 1190 to 2310 (mg protein per kg body weight). The application of a mean SF_direct of 1600 (mg protein per kg body weight) and further normalization by the microsomal CL_{int, u} values of midazolam, the most commonly used substrate, resulted in improved prediction accuracy for CL_{int, u} values from various microsomal batches.The results suggest the normalization of variability might be useful for predicting the in vivo CL_int.