A Fragment-Based Approach for the Computational Prediction of the Nonspecific Binding of Drugs to Hepatic Microsomes

Pharmacokinetics of Novel Furoxan/Coumarin Hybrids in Rats Using LC-MS/MS Method and Physiologically Based Pharmacokinetic Model

Novel furoxan/coumarin hybrids were synthesized, and pharmacologic studies showed that the compounds displayed potent antiproliferation activities via downregulating both the phosphatidylinositide 3-kinase (PI3K) pathway and the mitogen-activated protein kinase (MAPK) pathway. To investigate the preclinical pharmacokinetic (PK) properties of three candidate compounds (CY-14S-4A83, CY-16S-4A43, and CY-16S-4A93), liquid chromatography, in tandem with the mass spectrometry LC-MS/MS method, was developed and validated for the simultaneous determination of these compounds. The absorption, distribution, metabolism, and excretion (ADME) properties were investigated in in vitro studies and in rats. Meanwhile, physiologically based pharmacokinetic (PBPK) models were constructed using only in vitro data to obtain detailed PK information. Good linearity was observed over the concentration range of 0.01−1.0 μg/mL. The free drug fraction (fu) values of the compounds were less than 3%, and the clearance (CL) values were 414.5 ± 145.7 mL/h/kg, 2624.6 ± 648.4 mL/h/kg, and 500.6 ± 195.2 mL/h/kg, respectively. The predicted peak plasma concentration (Cmax) and the area under the concentration-time curve (AUC) were overestimated for the CY-16S-4A43 PBPK model compared with the experimental ones (fold error > 2), suggesting that tissue accumulation and additional elimination pathways may exist. In conclusion, the LC-MS/MS method was successively applied in the preclinical PK studies, and the detailed information from PBPK modeling may improve decision-making in subsequent new drug development.

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.

A Novel Method for Predicting the Human Inherent Clearance and Its Application in the Study of the Pharmacokinetics and Drug-Drug Interaction between Azidothymidine and Fluconazole Mediated by UGT Enzyme

In order to improve the benefit–risk ratio of pharmacokinetic (PK) research in the early development of new drugs, in silico and in vitro methods were constructed and improved. Models of intrinsic clearance rate (CL_int) were constructed based on the quantitative structure–activity relationship (QSAR) of 7882 collected compounds. Moreover, a novel in vitro metabolic method, the Bio-PK dynamic metabolic system, was constructed and combined with a physiology-based pharmacokinetic model (PBPK) model to predict the metabolism and the drug–drug interaction (DDI) of azidothymidine (AZT) and fluconazole (FCZ) mediated by the phase II metabolic enzyme UDP-glycosyltransferase (UGT) in humans. Compared with the QSAR models reported previously, the goodness of fit of our CL_int model was slightly improved (determination coefficient (R²) = 0.58 vs. 0.25–0.45). Meanwhile, compared with the predicted clearance of 61.96 L/h (fold error: 2.95–3.13) using CL_int (8 µL/min/mg) from traditional microsomal experiment, the predicted clearance using CL_int (25 μL/min/mg) from Bio-PK system was increased to 143.26 L/h (fold error: 1.27–1.36). The predicted C_max and AUC (the area under the concentration–time curve) ratio were 1.32 and 1.84 (fold error: 1.36 and 1.05) in a DDI study with an inhibition coefficient (Ki) of 13.97 μM from the Bio-PK system. The results indicate that the Bio-PK system more truly reflects the dynamic metabolism and DDI of AZT and FCZ in the body. In summary, the novel in silico and in vitro method may provide new ideas for the optimization of drug metabolism and DDI research methods in early drug development.

Successful and Unsuccessful Prediction of Human Hepatic Clearance for Lead Optimization

Development of new chemical entities is costly, time-consuming, and has a low success rate. Accurate prediction of pharmacokinetic properties is critical to progress compounds with favorable drug-like characteristics in lead optimization. Of particular importance is the prediction of hepatic clearance, which determines drug exposure and contributes to projection of dose, half-life, and bioavailability. The most commonly employed methodology to predict hepatic clearance is termed in vitro to in vivo extrapolation (IVIVE) that involves measuring drug metabolism in vitro, scaling-up this in vitro intrinsic clearance to a prediction of in vivo intrinsic clearance by reconciling the enzymatic content between the incubation and an average human liver, and applying a model of hepatic disposition to account for limitations of protein binding and blood flow to predict in vivo clearance. This manuscript reviews common in vitro techniques used to predict hepatic clearance as well as current challenges and recent theoretical advancements in IVIVE.

Evidence-based strategies for the characterisation of human drug and chemical glucuronidation in vitro and UDP-glucuronosyltransferase reaction phenotyping

Enzymes of the UDP-glucuronosyltransferase (UGT) superfamily contribute to the elimination of drugs from almost all therapeutic classes. Awareness of the importance of glucuronidation as a drug clearance mechanism along with increased knowledge of the enzymology of drug and chemical metabolism has stimulated interest in the development and application of approaches for the characterisation of human drug glucuronidation in vitro, in particular reaction phenotyping (the fractional contribution of the individual UGT enzymes responsible for the glucuronidation of a given drug), assessment of metabolic stability, and UGT enzyme inhibition by drugs and other xenobiotics. In turn, this has permitted the implementation of in vitro - in vivo extrapolation approaches for the prediction of drug metabolic clearance, intestinal availability, and drug-drug interaction liability, all of which are of considerable importance in pre-clinical drug development. Indeed, regulatory agencies (FDA and EMA) require UGT reaction phenotyping for new chemical entities if glucuronidation accounts for ≥25% of total metabolism. In vitro studies are most commonly performed with recombinant UGT enzymes and human liver microsomes (HLM) as the enzyme sources. Despite the widespread use of in vitro approaches for the characterisation of drug and chemical glucuronidation by HLM and recombinant enzymes, evidence-based guidelines relating to experimental approaches are lacking. Here we present evidence-based strategies for the characterisation of drug and chemical glucuronidation in vitro, and for UGT reaction phenotyping. We anticipate that the strategies will inform practice, encourage development of standardised experimental procedures where feasible, and guide ongoing research in the field.

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

The fraction unbound in the incubation, f_u,inc, is an important parameter to consider in the evaluation of intrinsic clearance measurements performed in vitro in hepatocytes or microsomes. Reliable estimates of f_u,inc based on a compound's structure have the potential to positively impact the screening timelines in drug discovery. Previous works suggested that f_u,inc is primarily driven by passive processes and can be described using physicochemical properties such as lipophilicity and the protonation state of the molecule. While models based on these principles proved predictive in relatively small datasets that included marketed drugs, their applicability domain has not been extensively explored. The work presented here from the in silico ADME discussion group (part of the International Consortium for Innovation through Quality in Pharmaceutical Development, the IQ consortium) describes the accuracy of these models in large proprietary datasets that include several thousand of compounds across chemical space. Overall, the models do well for compounds with low lipophilicity. In other words, the equations correctly predict that f_u,inc is, in general, above 0.5 for compounds with a calculated logP of less than 3. When applied to lipophilic compounds, the models failed to produce quantitatively accurate predictions of f_u,inc, with a high risk of underestimating binding properties. These models can, therefore, be used quantitatively for less lipophilic compounds. On the other hand, internal machine-learning models using a company's own proprietary dataset also predict compounds with higher lipophilicity reasonably well. Additionally, the data shown indicate that microsomal binding is, in general, a good proxy for hepatocyte binding.

Fraction Unbound for Liver Microsome and Hepatocyte Incubations for All Major Species Can Be Approximated Using a Single-Species Surrogate

It is well recognized that nonspecific binding of a drug within an in vitro assay (f u) can have a large impact on in vitro to in vivo correlations of intrinsic clearance. Typically, this value is determined experimentally across multiple species in the drug-discovery stage. Herein we examine the feasibility of using a single species (rat) as a surrogate for other species using a panel of small molecules representing highly diverse structures and physiochemical classes. The study demonstrated that 86% and 92% of the tested compounds measured in the mouse, dog, monkey, and human were within 2-fold of rat values for f u in microsomes and hepatocytes, respectively. One compound, amiodarone, exhibited unique species-dependent binding where the f u was approximately 10-fold higher in human microsomes and 20-fold higher in human hepatocytes compared with the average of the other species tested. Overall, these data indicate that using a single species (rat) f u as a surrogate for other major species, including humans, is a means to increase the throughput of measuring nonspecific binding in vitro.

Intrinsic Clearance Assay Incubational Binding: A Method Comparison

The fraction of unbound drug (fuinc) in in vitro intrinsic clearance (CLint) incubation is an important parameter in the pursuit of accurate clearance predictions and is often predicted using algorithms based on drug lipophilicity measures. However, analysis of an AstraZeneca database suggests that simple lipophilicity alone is a relatively poor predictor of fuinc measured using equilibrium dialysis. He fuinc value can also be measured directly in CLint assays using multiple concentrations of hepatocytes or microsomal protein. Since this approach informs of the unbound drug concentration in the assay used to predict in vivo clearance, it should be considered the gold standard method. As a starting point for building better predictive algorithms we aimed to determine if equilibrium dialysis really is an appropriate assay for assessing fuinc Employing a large number of compounds with a wide range of lipophilicities, experiments were performed to measure fuinc using rat hepatocytes (RH) and human liver microsomes (HLM) in both assay formats. A high percentage (94% and 93% for HLM and RH, respectively) of the fuinc values were within 2-fold when the compound distribution coefficient describing the ratio of compound concentration in octanol and pH 7.4 buffer when the test system is at equilibrium (lipophilicity measure) (logD7.4) values were less than 3.5. However, with logD7.4 values greater than these, the agreement was considerably worse. Additional experimental data generated indicated that this discrepancy was likely due to failings in the direct method when drug binding is high. Thus, we conclude that unbound CLint can be indeed calculated indirectly by incorporating equilibrium dialysis data with measured CLint but that simple lipophilicity descriptors alone may be inadequate for predicting fuinc.