Ontogeny of Human Liver Aldehyde Oxidase: Developmental Changes and Implications for Drug Metabolism.

Despite the increasing importance of aldehyde oxidase (AO) in the drug metabolism of clinical candidates, ontogeny data for AO are limited. The objective of our study was to characterize the age-dependent AO content and activity in the human liver cytosolic fraction (HLC) and human hepatocytes (HH). HLC (n = 121 donors) and HH (n = 50 donors) were analyzed for (1) AO protein content by quantitative proteomics and (2) enzyme activity using carbazeran as a probe substrate. AO activity showed high technical variability and poor correlation with the content in HLC samples, whereas hepatocyte samples showed a strong correlation between the content and activity. Similarly, AO content and activity showed no significant age-dependent differences in HLC samples, whereas the average AO content and activity in hepatocytes increased significantly (∼20-40-fold) from the neonatal levels (0-28 days). Based on the hepatocyte data, the age at which 50% of the adult AO content is reached (age50) was 3.15 years (0.32-13.97 years, 95% CI). Metabolite profiling of carbazeran revealed age-dependent metabolic switching and the role of non-AO mechanisms (glucuronidation and desmethylation) in carbazeran elimination. The content-activity correlation in hepatocytes improved significantly (R2 = 0.95; p < 0.0001) in samples showing <10% contribution of glucuronidation toward the overall metabolism, confirming that AO-mediated oxidation and glucuronidation are the key routes of carbazeran metabolism. Considering the confounding effect of glucuronidation on AO activity, AO content-based ontogeny data are a more direct reflection of developmental changes in protein expression. The comprehensive ontogeny data of AO in HH samples are more reliable than HLC data, which are important for developing robust physiologically based pharmacokinetic models for predicting AO-mediated metabolism in children.

Lenacapavir: A Novel, Potent, and Selective First-in-Class Inhibitor of HIV-1 Capsid Function Exhibits Optimal Pharmacokinetic Properties for a Long-Acting Injectable Antiretroviral Agent.

Lenacapavir (LEN) is a picomolar first-in-class capsid inhibitor of human immunodeficiency virus type 1 (HIV-1) with a multistage mechanism of action and no known cross resistance to other existing antiretroviral (ARV) drug classes. LEN exhibits a low aqueous solubility and exceptionally low systemic clearance following intravenous (IV) administration in nonclinical species and humans. LEN formulated in an aqueous suspension or a PEG/water solution formulation showed sustained plasma exposure levels with no unintended rapid drug release following subcutaneous (SC) administration to rats and dogs. A high total fraction dose release was observed with both formulations. The long-acting pharmacokinetics (PK) were recapitulated in humans following SC administration of both formulations. The SC PK profiles displayed two-phase absorption kinetics in both animals and humans with an initial fast-release absorption phase, followed by a slow-release absorption phase. Noncompartmental and compartmental analyses informed the LEN systemic input rate from the SC depot and exit rate from the body. Modeling-enabled deconvolution of the input rates from two processes: absorption of the soluble fraction (minor) from a direct fast-release process leading to the early PK phase and absorption of the precipitated fraction (major) from an indirect slow-release process leading to the later PK phase. LEN SC PK showed flip-flop kinetics due to the input rate being substantially slower than the systemic exit rate. LEN input rates via the slow-release process in humans were slower than those in both rats and dogs. Overall, the combination of high potency, exceptional stability, and optimal release rate from the injection depot make LEN well suited for a parenteral long-acting formulation that can be administered once up to every 6 months in humans for the prevention and treatment of HIV-1.

Dissecting Parameters Contributing to the Underprediction of Aldehyde Oxidase-Mediated Metabolic Clearance of Drugs.

We investigated the effect of variability and instability in aldehyde oxidase (AO) content and activity on the scaling of in vitro metabolism data. AO content and activity in human liver cytosol (HLC) and five recombinant human AO preparations (rAO) were determined using targeted proteomics and carbazeran oxidation assay, respectively. AO content was highly variable as indicated by the relative expression factor (REF; i.e., HLC to rAO content) ranging from 0.001 to 1.7 across different in vitro systems. The activity of AO in HLC degrades at a 10-fold higher rate in the presence of the substrate as compared with the activity performed after preincubation without substrate. To scale the metabolic activity from rAO to HLC, a protein-normalized activity factor (pnAF) was proposed wherein the activity was corrected by AO content, which revealed up to sixfold higher AO activity in HLC versus rAO systems. A similar value of pnAF was observed for another substrate, ripasudil. Physiologically based pharmacokinetic (PBPK) modeling revealed a significant additional clearance (CL; 66%), which allowed for the successful prediction of in vivo CL of four other substrates, i.e., O-benzyl guanine, BIBX1382, zaleplon, and zoniporide. For carbazeran, the metabolite identification study showed that the direct glucuronidation may be contributing to around 12% elimination. Taken together, this study identified differential protein content, instability of in vitro activity, role of additional AO clearance, and unaccounted metabolic pathways as plausible reasons for the underprediction of AO-mediated drug metabolism. Consideration of these factors and integration of REF and pnAF in PBPK models will allow better prediction of AO metabolism. SIGNIFICANCE STATEMENT: This study elucidated the plausible reasons for the underprediction of aldehyde oxidase (AO)-mediated drug metabolism and provided recommendations to address them. It demonstrated that integrating protein content and activity differences and accounting for the loss of AO activity, as well as consideration of extrahepatic clearance and additional pathways, would improve the in vitro to in vivo extrapolation of AO-mediated drug metabolism using physiologically based pharmacokinetic modeling.

Human and nonclinical disposition of [14C]bictegravir, a potent integrase strand-transfer inhibitor for the treatment of HIV-1 infection.

Bictegravir (BIC) is a potent small-molecule integrase strand-transfer inhibitor (INSTI) and a component of Biktarvy®, a single-tablet combination regimen that is currently approved for the treatment of human immunodeficiency virus type 1 (HIV-1) infection. The absorption, metabolism, distribution, and elimination (ADME) characteristics of BIC were determined through in vivo nonclinical and clinical studies (IND 121318).[14C]BIC was rapidly absorbed orally in mice, rats, monkeys and human. The cumulative dose recovery was high in nonclinical species (>80%) and humans (95.3%), with most of the excreted dose recovered in faeces. Quantifiable radioactivity with declining concentration was observed in rat tissues suggesting reversible binding. Unchanged BIC was the most abundant circulating component in all species along with two notable metabolites M20 (a sulphate conjugate of hydroxylated BIC) and M15 (a glucuronide conjugate of BIC). BIC was primarily eliminated by hepatic metabolism followed by excretion of the biotransformed products into faeces. In vitro drug-drug interaction (DDI) studies with M15 and M20 demonstrated that no clinically relevant interactions were expected.Overall, BIC is a novel and potent INSTI with a favourable resistance, PK, and ADME profile that provides important improvements over other currently available INSTIs for the treatment of HIV-1.

18O-Enabled High-Throughput Acyl Glucuronide Stability Assay.

Acyl glucuronides (AGs) are common metabolites of carboxylic acid-containing compounds. In some circumstances, AGs are suspected to be involved in drug toxicity due to formation of acyl migration products that bind covalently to cellular components. The risk of this adverse effect has been found to be correlated with the chemical stability of the AG, and assays have been described that monitor acyl migration by liquid chromatography coupled with mass spectrometry (LC-MS). This analysis can be challenging as it requires baseline chromatographic separation of the unmigrated 1-ß-acyl glucuronide from the migrated isomers and thus needs to be individually optimized for each aglycone. Therefore, a high-throughput assay that eliminates LC method development is desirable. Herein, we report an improved acyl glucuronide stability assay based on the rate of 18O-incorporation from [18O] water, which is compatible with high-throughput bioanalytical LC-MS workflows. Synthetic AGs with shorter migration half-lives showed faster incorporation of 18O. The level of differential incorporation of 18O following a 24 h incubation correlates well with the migration tendency of AGs. This assay was developed further, exploring in situ generation of AGs by human hepatic microsomal fraction. The results from 18 in situ-formed acyl glucuronides were similar to those obtained using authentic reference standards. In this format, this new 18O-labeling method offers a simplified workflow, requires no LC method development or AG reference standard, and thus facilitates AG liability assessment in early drug discovery.

Intravenous delivery of GS-441524 is efficacious in the African green monkey model of SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the COVID-19 pandemic, has infected over 260 million people over the past 2 years. Remdesivir (RDV, VEKLURY®) is currently the only antiviral therapy fully approved by the FDA for the treatment of COVID-19. The parent nucleoside of RDV, GS-441524, exhibits antiviral activity against numerous respiratory viruses including SARS-CoV-2, although at reduced in vitro potency compared to RDV in most assays. Here we find in both human alveolar and bronchial primary cells, GS-441524 is metabolized to the pharmacologically active GS-441524 triphosphate (TP) less efficiently than RDV, which correlates with a lower in vitro SARS-CoV-2 antiviral activity. In vivo, African green monkeys (AGM) orally dosed with GS-441524 yielded low plasma levels due to limited oral bioavailability of <10%. When GS-441524 was delivered via intravenous (IV) administration, although plasma concentrations of GS-441524 were significantly higher, lung TP levels were lower than observed from IV RDV. To determine the required systemic exposure of GS-441524 associated with in vivo antiviral efficacy, SARS-CoV-2 infected AGMs were treated with a once-daily IV dose of either 7.5 or 20 mg/kg GS-441524 or IV RDV for 5 days and compared to vehicle control. Despite the reduced lung TP formation compared to IV dosing of RDV, daily treatment with IV GS-441524 resulted in dose-dependent efficacy, with the 20 mg/kg GS-441524 treatment resulting in significant reductions of SARS-CoV-2 replication in the lower respiratory tract of infected animals. These findings demonstrate the in vivo SARS-CoV-2 antiviral efficacy of GS-441524 and support evaluation of its orally bioavailable prodrugs as potential therapies for COVID-19.

Comparison of Tissue Abundance of Non-Cytochrome P450 Drug-Metabolizing Enzymes by Quantitative Proteomics between Humans and Laboratory Animal Species.

The use of animal pharmacokinetic models as surrogates for humans relies on the assumption that the drug disposition mechanisms are similar between preclinical species and humans. However, significant cross-species differences exist in the tissue distribution and protein abundance of drug-metabolizing enzymes (DMEs) and transporters. We quantified non-cytochrome P450 (non-CYP) DMEs across commonly used preclinical species (cynomolgus and rhesus monkeys, beagle dog, Sprague Dawley and Wistar Han rats, and CD1 mouse) and compared these data with previously obtained human data. Aldehyde oxidase was abundant in humans and monkeys while poorly expressed in rodents, and not expressed in dogs. Carboxylesterase (CES) 1 abundance was highest in the liver while CES2 was primarily expressed in the intestine in all species with notable species differences. For example, hepatic CES1 was 3× higher in humans than in monkeys, but hepatic CES2 was 3-5× higher in monkeys than in humans. Hepatic UDP-glucuronosyltransferase (UGT) 1A2 abundance was ∼4× higher in dogs compared with rats, whereas UGT1A3 abundance was 3-5× higher in dog livers than its ortholog in human and monkey livers. UGT1A6 abundance was 5-6× higher in human livers compared with monkey and dog livers. Hepatic sulfotransferase 1B1 abundance was 5-7× higher in rats compared with the rest of the species. These quantitative non-CYP proteomics data can be used to explain unique toxicological profiles across species and can be integrated into physiologically based pharmacokinetic models for the mechanistic explanation of pharmacokinetics and tissue distribution of xenobiotics in animal species. SIGNIFICANCE STATEMENT: We characterized the quantitative differences in non-cytochrome P450 (non-CYP) drug-metabolizing enzymes across commonly used preclinical species (cynomolgus and rhesus monkeys, beagle dogs, Sprague Dawley and Wistar Han rats, and CD1 mice) and compared these data with previously obtained human data. Unique differences in non-CYP enzymes across species were observed, which can be used to explain significant pharmacokinetic and toxicokinetic differences between experimental animals and humans.

Human pharmacokinetics prediction with an in vitro-in vivo correction factor approach and in vitro drug-drug interaction profile of bictegravir, a potent integrase-strand transfer inhibitor component in approved biktarvy® for the treatment of HIV-1 infection.

Bictegravir (BIC) is a potent small-molecule integrase strand-transfer inhibitor (INSTI) and a component of Biktarvy®, a single-tablet combination regimen that is currently approved for the treatment of human immunodeficiency virus type 1 (HIV-1) infection. The in vitro properties, pharmacokinetics (PK), and drug-drug interaction (DDI) profile of BIC were characterised in vitro and in vivo.BIC is a weakly acidic, ionisable, lipophilic, highly plasma protein-bound BCS class 2 molecule, which makes it difficult to predict human PK using standard methods. Its systemic plasma clearance is low, and the volume of distribution is approximately the volume of extracellular water in nonclinical species. BIC metabolism is predominantly mediated by cytochrome P450 enzyme (CYP) 3A and UDP-glucuronosyltransferase 1A1. BIC shows a low potential to perpetrate clinically meaningful DDIs via known drug metabolising enzymes or transporters.The human PK of BIC was predicted using a combination of bioavailability and volume of distribution scaled from nonclinical species and a modified in vitro-in vivo correlation (IVIVC) correction for clearance. Phase 1 studies in healthy subjects largely bore out the prediction and supported the methods used. The approach presented herein could be useful for other drug molecules where standard projections are not sufficiently accurate. .

Key Metabolic Enzymes Involved in Remdesivir Activation in Human Lung Cells.

Remdesivir (RDV; GS-5734, Veklury), the first FDA-approved antiviral to treat COVID-19, is a single-diastereomer monophosphoramidate prodrug of an adenosine analogue. RDV is taken up in the target cells and metabolized in multiple steps to form the active nucleoside triphosphate (TP) (GS-443902), which, in turn, acts as a potent and selective inhibitor of multiple viral RNA polymerases. In this report, we profiled the key enzymes involved in the RDV metabolic pathway with multiple parallel approaches: (i) bioinformatic analysis of nucleoside/nucleotide metabolic enzyme mRNA expression using public human tissue and lung single-cell bulk mRNA sequence (RNA-seq) data sets, (ii) protein and mRNA quantification of enzymes in human lung tissue and primary lung cells, (iii) biochemical studies on the catalytic rate of key enzymes, (iv) effects of specific enzyme inhibitors on the GS-443902 formation, and (v) the effects of these inhibitors on RDV antiviral activity against SARS-CoV-2 in cell culture. Our data collectively demonstrated that carboxylesterase 1 (CES1) and cathepsin A (CatA) are enzymes involved in hydrolyzing RDV to its alanine intermediate MetX, which is further hydrolyzed to the monophosphate form by histidine triad nucleotide-binding protein 1 (HINT1). The monophosphate is then consecutively phosphorylated to diphosphate and triphosphate by cellular phosphotransferases. Our data support the hypothesis that the unique properties of RDV prodrug not only allow lung-specific accumulation critical for the treatment of respiratory viral infection such as COVID-19 but also enable efficient intracellular metabolism of RDV and its MetX to monophosphate and successive phosphorylation to form the active TP in disease-relevant cells.

Characterization of Differential Tissue Abundance of Major Non-CYP Enzymes in Human.

The availability of assays that predict the contribution of cytochrome P450 (CYP) metabolism allows for the design of new chemical entities (NCEs) with minimal oxidative metabolism. These NCEs are often substrates of non-CYP drug-metabolizing enzymes (DMEs), such as UDP-glucuronosyltransferases (UGTs), sulfotransferases (SULTs), carboxylesterases (CESs), and aldehyde oxidase (AO). Nearly 30% of clinically approved drugs are metabolized by non-CYP enzymes. However, knowledge about the differential hepatic versus extrahepatic abundance of non-CYP DMEs is limited. In this study, we detected and quantified the protein abundance of eighteen non-CYP DMEs (AO, CES1 and 2, ten UGTs, and five SULTs) across five different human tissues. AO was most abundantly expressed in the liver and to a lesser extent in the kidney; however, it was not detected in the intestine, heart, or lung. CESs were ubiquitously expressed with CES1 being predominant in the liver, while CES2 was enriched in the small intestine. Consistent with the literature, UGT1A4, UGT2B4, and UGT2B15 demonstrated liver-specific expression, whereas UGT1A10 expression was specific to the intestine. UGT1A1 and UGT1A3 were expressed in both the liver and intestine; UGT1A9 was expressed in the liver and kidney; and UGT2B17 levels were significantly higher in the intestine than in the liver. All five SULTs were detected in the liver and intestine, and SULT1A1 and 1A3 were detected in the lung. Kidney abundance was the most variable among the studied tissues, and overall, high interindividual variability (>15-fold) was observed for UGT2B17, CES2 (intestine), SULT1A1 (liver), UGT1A9, UGT2B7, and CES1 (kidney). These differential tissue abundance data can be integrated into physiologically based pharmacokinetic (PBPK) models for the prediction of non-CYP drug metabolism and toxicity in hepatic and extrahepatic tissues.

Characterization of Seasonal Pharmacokinetic Variability in Woodchucks.

The eastern woodchuck (Marmota monax) is a hibernating species extensively used as an in vivo efficacy model for chronic human hepatitis B virus infection. Under laboratory conditions, woodchucks develop a pseudohibernation condition; thus, the pharmacokinetics (PK) of small-molecule therapeutics may be affected by the seasonal change. The seasonal PK of four probe compounds were characterized over 12 months in seven male and nine female laboratory-maintained woodchucks. These compounds were selected to study changes in oxidative metabolism [antipyrine (AP)], glucuronidation [raltegravir (RTG)], renal clearance [lamivudine (3TC)], and hepatic function [indocyanine green (ICG)]. Seasonal changes in physiologic parameters and PK were determined. Seasonal body weight increases were ≥30%. Seasonal changes in body temperature and heart rate were <10%. The mean AP exposure remained unchanged from April to August 2017, followed by a significant increase (≥1.0-fold) from August to December and subsequent decrease to baseline at the end of study. A similar trend was observed in RTG and 3TC exposures. The ICG exposure remained unchanged. No significant sex difference in PK was observed, although female woodchucks appeared to be less susceptible to seasonal PK and body weight changes. Significant seasonal PK changes for AP, RTG, and 3TC indicate decreases in oxidative metabolism, phase II glucuronidation, and renal clearance during pseudohibernation. The lack of seasonal change in ICG exposure suggests there are no significant changes in hepatic function. This information can be used to optimize the scheduling of woodchuck studies to avoid seasonally driven variation in drug PK. SIGNIFICANCE STATEMENT: Woodchuck is a hibernating species and is commonly used as a nonclinical model of hepatitis B infection. Investigation of seasonal PK changes is perhaps of greater interest to pharmaceutical industry scientists, who use the woodchuck model to optimize the scheduling of woodchuck studies to avoid seasonally driven variation in drug PK and/or toxicity. This information is also valuable to drug metabolism and veterinary scientists in understanding woodchuck's seasonal metabolism and behavior under the pseudohibernation condition.

Regional Proteomic Quantification of Clinically Relevant Non-Cytochrome P450 Enzymes along the Human Small Intestine.

Current challenges in accurately predicting intestinal metabolism arise from the complex nature of the intestine, leading to limited applicability of available in vitro tools as well as knowledge deficits in intestinal physiology, including enzyme abundance. In particular, information on regional enzyme abundance along the small intestine is lacking, especially for non-cytochrome P450 enzymes such as carboxylesterases (CESs), UDP-glucuronosyltransferases (UGTs), and sulfotransferases (SULTs). We used cryopreserved human intestinal mucosa samples from nine donors as an in vitro surrogate model for the small intestine and performed liquid chromatography tandem mass spectrometry-based quantitative proteomics for 17 non-cytochrome P450 enzymes using stable isotope-labeled peptides. Relative protein quantification was done by normalization with enterocyte marker proteins, i.e., villin-1, sucrase isomaltase, and fatty acid binding protein 2, and absolute protein quantification is reported as picomoles per milligram of protein. Activity assays in glucuronidations and sequential metabolisms were conducted to validate the proteomics findings. Relative or absolute quantifications are reported for CES1, CES2, five UGTs, and four SULTs along the small intestine: duodenum, jejunum, and ileum for six donors and in 10 segments along the entire small intestine (A-J) for three donors. Relative quantification using marker proteins may be beneficial in further controlling for technical variabilities. Absolute quantification data will allow for scaling factor generation and in vivo extrapolation of intestinal clearance using physiologically based pharmacokinetic modeling. SIGNIFICANCE STATEMENT: Current knowledge gaps exist in intestinal protein abundance of non-cytochrome P450 enzymes. Here, we employ quantitative proteomics to measure non-cytochrome P450 enzymes along the human small intestine in nine donors using cryopreserved human intestinal mucosa samples. Absolute and relative abundances reported here will allow better scaling of intestinal clearance.

Pharmacokinetics and Disposition of Momelotinib Revealed a Disproportionate Human Metabolite-Resolution for Clinical Development.

Momelotinib (MMB), a small-molecule inhibitor of Janus kinase (JAK)1/2 and of activin A receptor type 1 (ACVR1), is in clinical development for the treatment of myeloproliferative neoplasms. The pharmacokinetics and disposition of [14C]MMB were characterized in a single-dose, human mass-balance study. Metabolism and the pharmacologic activity of key metabolites were elucidated in multiple in vitro and in vivo experiments. MMB was rapidly absorbed following oral dosing with approximately 97% of the radioactivity recovered, primarily in feces with urine as a secondary route. Mean blood-to-plasma [14C] area under the plasma concentration-time curve ratio was 0.72, suggesting low association of MMB and metabolites with blood cells. [14C]MMB-derived radioactivity was detectable in blood for ≤48 hours, suggesting no irreversible binding of MMB or its metabolites. The major circulating human metabolite, M21 (a morpholino lactam), is a potent inhibitor of JAK1/2 and ACVR1 in vitro. Estimation of pharmacological activity index suggests M21 contributes significantly to the pharmacological activity of MMB for the inhibition of both JAK1/2 and ACVR1. M21 was observed in disproportionately higher amounts in human plasma than in rat or dog, the rodent and nonrodent species used for the general nonclinical safety assessment of this molecule. This discrepancy was resolved with additional nonclinical studies wherein the circulating metabolites and drug-drug interactions were further characterized. The human metabolism of MMB was mediated primarily by multiple cytochrome P450 enzymes, whereas M21 formation involved initial P450 oxidation of the morpholine ring followed by metabolism via aldehyde oxidase.

In Vitro Kinetic Characterization of Axitinib Metabolism.

N-Methyl-2-[3-((E)-2-pyridin-2-yl-vinyl)-1H-indazol-6-ylsulfanyl]-benzamide (axitinib) is an oral inhibitor of vascular endothelial growth factor receptors 1-3, which is approved for the treatment of advanced renal cell cancer. Human [(14)C]-labeled clinical studies indicate axitinib's primary route of clearance is metabolism. The aims of the in vitro experiments presented herein were to identify and characterize the enzymes involved in axitinib metabolic clearance. In vitro biotransformation studies of axitinib identified a number of metabolites including an axitinib sulfoxide, several less abundant oxidative metabolites, and glucuronide conjugates. The most abundant NADPH- and UDPGA-dependent metabolites, axitinib sulfoxide (M12) and axitinib N-glucuronide (M7) were selected for phenotyping and kinetic study. Phenotyping experiments with human liver microsomes (HLMs) using chemical inhibitors and recombinant human cytochrome P450s demonstrated axitinib was predominately metabolized by CYP3A4/5, with minor contributions from CYP2C19 and CYP1A2. The apparent substrate concentration at half-maximal velocity (Km) and Vmax values for the formation of axitinib sulfoxide by CYP3A4 or CYP3A5 were 4.0 or 1.9 µM and 9.6 or 1.4 pmol·min(-1)·pmol(-1), respectively. Using a CYP3A4-specific inhibitor (Cyp3cide) in liver microsomes expressing CYP3A5, 66% of the axitinib intrinsic clearance was attributable to CYP3A4 and 15% to CYP3A5. Axitinib N-glucuronidation was primarily catalyzed by UDP-glucuronosyltransferase (UGT) UGT1A1, which was verified by chemical inhibitors and UGT1A1 null expressers, with lesser contributions from UGTs 1A3, 1A9, and 1A4. The Km and Vmax values describing the formation of the N-glucuronide in HLM or rUGT1A1 were 2.7 µM or 0.75 µM and 8.9 or 8.3 pmol·min(-1)·mg(-1), respectively. In summary, CYP3A4 is the major enzyme involved in axitinib clearance with lesser contributions from CYP3A5, CYP2C19, CYP1A2, and UGT1A1.

Prediction of Drug-Drug Interactions with Crizotinib as the CYP3A Substrate Using a Physiologically Based Pharmacokinetic Model.

An orally available multiple tyrosine kinase inhibitor, crizotinib (Xalkori), is a CYP3A substrate, moderate time-dependent inhibitor, and weak inducer. The main objectives of the present study were to: 1) develop and refine a physiologically based pharmacokinetic (PBPK) model of crizotinib on the basis of clinical single- and multiple-dose results, 2) verify the crizotinib PBPK model from crizotinib single-dose drug-drug interaction (DDI) results with multiple-dose coadministration of ketoconazole or rifampin, and 3) apply the crizotinib PBPK model to predict crizotinib multiple-dose DDI outcomes. We also focused on gaining insights into the underlying mechanisms mediating crizotinib DDIs using a dynamic PBPK model, the Simcyp population-based simulator. First, PBPK model-predicted crizotinib exposures adequately matched clinically observed results in the single- and multiple-dose studies. Second, the model-predicted crizotinib exposures sufficiently matched clinically observed results in the crizotinib single-dose DDI studies with ketoconazole or rifampin, resulting in the reasonably predicted fold-increases in crizotinib exposures. Finally, the predicted fold-increases in crizotinib exposures in the multiple-dose DDI studies were roughly comparable to those in the single-dose DDI studies, suggesting that the effects of crizotinib CYP3A time-dependent inhibition (net inhibition) on the multiple-dose DDI outcomes would be negligible. Therefore, crizotinib dose-adjustment in the multiple-dose DDI studies could be made on the basis of currently available single-dose results. Overall, we believe that the crizotinib PBPK model developed, refined, and verified in the present study would adequately predict crizotinib oral exposures in other clinical studies, such as DDIs with weak/moderate CYP3A inhibitors/inducers and drug-disease interactions in patients with hepatic or renal impairment.

Metabolism, excretion and pharmacokinetics of [14C]crizotinib following oral administration to healthy subjects.

1. Crizotinib (XALKORI®), an oral inhibitor of anaplastic lymphoma kinase (ALK) and mesenchymal-epithelial transition factor kinase (c-Met), is currently approved for the treatment of patients with non-small cell lung cancer that is ALK-positive. 2. The metabolism, excretion and pharmacokinetics of crizotinib were investigated following administration of a single oral dose of 250 mg/100 µCi [(14)C]crizotinib to six healthy male subjects. 3. Mean recovery of [(14)C]crizotinib-related radioactivity in excreta samples was 85% of the dose (63% in feces and 22% in urine). 4. Crizotinib and its metabolite, crizotinib lactam, were the major components circulating in plasma, accounting for 33% and 10%, respectively, of the 0-96 h plasma radioactivity. Unchanged crizotinib was the major excreted component in feces (∼ 53% of the dose). In urine, crizotinib and O-desalkyl crizotinib lactam accounted for ∼ 2% and 5% of the dose, respectively. Collectively, these data indicate that the primary clearance pathway for crizotinib in humans is oxidative metabolism/hepatic elimination. 5. Based on plasma exposure in healthy subjects following a single dose of crizotinib and in vitro potency against ALK and c-Met, the crizotinib lactam diastereomers are not anticipated to contribute significantly to in vivo activity; however, additional assessment in cancer patients is warranted.

Pharmacokinetics, metabolism, and excretion of [14C]axitinib, a vascular endothelial growth factor receptor tyrosine kinase inhibitor, in humans.

The disposition of a single oral dose of 5 mg (100 µCi) of [(14)C]axitinib was investigated in fasted healthy human subjects (N = 8). Axitinib was rapidly absorbed, with a median plasma Tmax of 2.2 hours and a geometric mean Cmax and half-life of 29.2 ng/ml and 10.6 hours, respectively. The plasma total radioactivity-time profile was similar to that of axitinib but the AUC was greater, suggesting the presence of metabolites. The major metabolites in human plasma (0-12 hours), identified as axitinib N-glucuronide (M7) and axitinib sulfoxide (M12), were pharmacologically inactive, and with axitinib comprised 50.4%, 16.2%, and 22.5% of the radioactivity, respectively. In excreta, the majority of radioactivity was recovered in most subjects by 48 hours postdose. The median radioactivity excreted in urine, feces, and total recovery was 22.7%, 37.0%, and 59.7%, respectively. The recovery from feces was variable across subjects (range, 2.5%-60.2%). The metabolites identified in urine were M5 (carboxylic acid), M12 (sulfoxide), M7 (N-glucuronide), M9 (sulfoxide/N-oxide), and M8a (methylhydroxy glucuronide), accounting for 5.7%, 3.5%, 2.6%, 1.7%, and 1.3% of the dose, respectively. The drug-related products identified in feces were unchanged axitinib, M14/15 (mono-oxidation/sulfone), M12a (epoxide), and an unidentified metabolite, comprising 12%, 5.7%, 5.1%, and 5.0% of the dose, respectively. The proposed mechanism to form M5 involved a carbon-carbon bond cleavage via M12a, followed by rearrangement to a ketone intermediate and subsequent Baeyer-Villiger rearrangement, possibly through a peroxide intermediate. In summary, the study characterized axitinib metabolites in circulation and primary elimination pathways of the drug, which were mainly oxidative in nature.

In vitro characterization of axitinib interactions with human efflux and hepatic uptake transporters: implications for disposition and drug interactions.

Axitinib is an inhibitor of tyrosine kinase vascular endothelin growth factor receptors 1, 2, and 3. The ATP-binding cassette (ABC) and solute carrier (SLC) transport properties of axitinib were determined in selected cellular systems. Axitinib exhibited high passive permeability in all cell lines evaluated (Papp ≥ 6 × 10(-6) cm/s). Active efflux was observed in Caco-2 cells, and further evaluation in multidrug resistance gene 1 (MDR1) or breast cancer resistance protein (BCRP) transfected Madin-Darby canine kidney cells type 2 (MDCK) cells indicated that axitinib is at most only a weak substrate for P-glycoprotein (P-gp) but not BCRP. Axitinib showed incomplete inhibition of P-gp-mediated transport of digoxin in Caco-2 cells and BCRP transport of topotecan in BCRP-transfected MDCK cells with IC50 values of 3 µM and 4.4 µM, respectively. Axitinib (10 mg) did not pose a risk for systemic drug interactions with P-gp or BCRP per regulatory guidance. A potential risk for drug interactions through inhibition of P-gp and BCRP in the gastrointestinal tract was identified because an axitinib dose of 10 mg divided by 250 mL was greater than 10-fold the IC50 for each transporter. However, a GastroPlus simulation that considered the low solubility of axitinib resulted in lower intestinal concentrations and suggested a low potential for gastrointestinal interactions with P-gp and BCRP substrates. Organic anion transporting polypeptide 1B1 (OATP1B1) and OATP1B3 transfected human embryonic kidney 293 (HEK293) cells transported axitinib to a minor extent but uptake into suspended hepatocytes was not inhibited by rifamycin SV suggesting that high passive permeability predominates. Mouse whole-body autoradiography revealed that [(14)C]axitinib-equivalents showed rapid absorption and distribution to all tissues except the brain. This suggests that efflux transport of axitinib may occur at the mouse blood-brain barrier.

Application of target-mediated drug disposition model to small molecule heat shock protein 90 inhibitors.

Replacement of hydrogen with fluorine within three pairs of structurally similar small molecule inhibitors of heat shock protein 90 (HSP90) resulted in differences in inhibition constants (K(i)) in vitro as well as marked differences in rat intravenous pharmacokinetic profiles. The difference in pharmacokinetic profiles between lower and higher affinity inhibitors (LAIs and HAIs, respectively) was characterized by remarkably different estimates for steady-state volumes of distribution (V(ss): 1.8-2.0 versus 10-13 l/kg) with comparable clearance estimates (3.2-3.5 l/h per kilogram). When the observed V(ss) estimates were compared with the values predicted with the tissue-composition-based model, the observed V(ss) estimates for HAIs were 4- to 8-fold larger than the predicted values, whereas the V(ss) values for LAIs were comparable. Accordingly, a negative relationship between in vitro HSP90 K(i) versus in vivo V(ss) estimates was observed among these inhibitors. We therefore hypothesized that pharmacokinetic profiles of these inhibitors could be characterized by a target-mediated drug disposition (TMDD) model. In vivo equilibrium dissociation constant (K(D)) estimates for HAIs due to target binding by TMDD model with rapid binding approximation were 1-6 nM (equivalent to 0.3-2 nM free drug), which appeared comparable to the in vitro K(i) estimates (2-3 nM). In vivo KD values of LAIs were not accurately determined by the TMDD model, likely due to nonspecific binding-dependent tissue distribution obscuring TMDD profiles. Overall, these results suggest that the observed large Vss estimates for potent HSP90 inhibitors are likely due to pharmacological target binding.