Abstract
Aldehyde oxidase (AO) efficiently metabolizes a range of compounds with N-containing heterocyclic aromatic rings and/or aldehydes. The limited knowledge of AO activity and abundance (in vitro and in vivo) has led to poor prediction of in vivo systemic clearance (CL) using in vitro–to–in vivo extrapolation approaches, which for drugs in development can lead to their discontinuation. We aimed to identify appropriate scaling factors to predict AO CL of future new chemical entities (NCEs). The metabolism of six AO substrates was measured in human liver cytosol (HLC) and S9 fractions. Measured blood-to-plasma ratios and free fractions (in the in vitro system and in plasma) were used to develop physiologically based pharmacokinetic models for each compound. The impact of extrahepatic metabolism was explored, and the intrinsic clearance required to recover in vivo profiles was estimated and compared with in vitro measurements. Using HLC data and assuming only hepatic metabolism, a systematic underprediction of clearance was observed (average fold underprediction was 3.8). Adding extrahepatic metabolism improved the accuracy of the results (average fold error of 1.9). A workflow for predicting metabolism of an NCE by AO is proposed, and an empirical (laboratory-specific) scaling factor of three on the predicted intravenous CL allows a reasonable prediction of the available clinical data. Alternatively, considering also extrahepatic metabolism, an scaling factor of 6.5 applied on the intrinsic clearance could be used. Future research should focus on the impact of the in vitro study designs and the contribution of extrahepatic metabolism to AO-mediated clearance to understand the mechanisms behind the systematic underprediction.
SIGNIFICANCE STATEMENT This works describes the development of scaling factors to allow in vitro–in vivo extrapolation of the clearance of compounds by aldehyde oxidase metabolism in humans. In addition, physiologically based pharmacokinetic models were developed for each of the aldehyde oxidase substrate compounds investigated.
Introduction
Aldehyde oxidase (AO) is a cytosolic molybdenum-containing enzyme that very efficiently oxidizes a range of N-containing heterocyclic aromatic rings and aldehydes (Montefiori et al., 2017). The limited knowledge about AO activity, abundance, and translation from in vitro to in vivo has led to poor prediction of in vivo clearance (CL) and consequently to clinical failure of some AO substrates (Fan et al., 2016; Jensen et al., 2017). The reasons for the poor prediction of CL of aldehyde oxidase substrates (CLAO) are multiple: Firstly, AO is only expressed in the cytosol, and therefore standard metabolism studies using human liver microsomes (HLMs) will overlook AO metabolism (Obach, 2011; Zientek and Youdim, 2015; Dalvie and Di, 2019). Secondly, scaling from animal data can be challenging because dogs have low AO expression, and the AO expression in different rat strains is very variable (Dalvie et al., 2013; Tanoue et al., 2013). Additionally, although monkey, rat, and rabbit AO have been investigated, the metabolism of compounds by AO in these species only moderately overlaps with human clearance by AO (Al salhen, 2014; Dick, 2018). Finally, the interindividual variability in AO expression in humans is high, and consequently, the risk of having data for a nonrepresentative individual is considerable (Hutzler et al., 2014). In 2010, Zientek et al. (2010) have predicted the in vivo intrinsic clearance (CLint) for AO from in vitro data and compared them to the CLint,AO estimated after intravenous administration for five compounds. The CLint,AO was underestimated by 13-fold (range: 5–32) when using human liver cytosol (HLC) data and by 15-fold (range: 3–52) when using human liver S9 (HLS9) data. Strategies have been formulated to handle AO-mediated clearance in drug discovery and development. However, understanding why human in vivo clearance is underpredicted using in vitro CLint,AO from HLC and HLS9 is necessary to understand the role of AO in the metabolism of new chemical entities in a wider chemical space.
In this work, we aimed 1) to assess the prediction of intravenous clearance (CLIV) of six AO substrates from in vitro data and then 2) to derive an empirical scaling factor that 3) could be used to predict the CL of future NCE using physiologically based pharmacokinetic (PBPK) modeling. Initially, the metabolism by AO was assumed to occur only in the liver, and results from HLC and HLS9 were compared. However, there is compelling evidence of extrahepatic metabolism (CLIV greater than the hepatic blood flow and AO expression data in extrahepatic tissues); hence, the impact of extrahepatic metabolism was also explored. Finally, scaling factors were estimated to optimally recover in vivo concentration-time profiles.
Materials and Methods
Pooled (150 donors; lot 38289) human liver cytosol, pooled (150 donors; lot 3829) human liver S9, and pooled (150 donors; equal sex mix) human liver microsomes were obtained from Corning Life Sciences (Woburn, MA). It is possible that the human donor liver tissues might have been perfused or preserved with University of Wisconsin solution or another allopurinol-containing buffer, which may exhibit aldehyde oxidase inhibition potential at high concentrations.
Frozen human plasma (pooled from 78 individuals, mixed sex) generated using K2-EDTA as an anticoagulant was purchased from BioreclamationIVT (Baltimore, MD). AO substrates O6-benzylguanine, zaleplon, zoniporide, and carbazeran were sourced from Sigma-Aldrich (Poole, UK). BIBX1382 was sourced from Santa Cruz Biotechnology (Dallas, TX). Ziprasidone was synthesized at AstraZeneca (Cambridge, UK). Formic acid, ammonium formate, and DMSO were purchased from Sigma-Aldrich. HPLC-grade methanol, water, and acetonitrile (ACN) were obtained from Thermo Fisher Scientific (Waltham, MA). All other solvents were HPLC-grade and, unless otherwise specified, all other reagents were purchased from Sigma-Aldrich.
Compound Selection.
The aim of this work was to assess the robustness of the prediction of CLAO. To remove the additional uncertainty associated with predicting processes influencing oral bioavailability, only drugs with reported intravenous clearance and a fraction metabolized (fm) by AO (fmAO) of greater than 5% were selected for inclusion in this exercise. Using these criteria, the involvement of AO in the metabolism of O6-benzylguanine, BIBX1382, carbazeran, zaleplon, ziprasidone, and zoniporide was investigated in this study (Figure 1). Physicochemical data, including molecular weight, logP and pKa, acid/base nature as well as blood binding properties, and information about the compound elimination, were compiled for all drugs (Table 1). When data from several reliable sources were available, a weighted mean value was used.
PBPK model input parameters and clinical data used for PBPK model verification. AGP, Alpha 1 Acid Glycoprotein; LOQ, limit of quantifiation; MW, molecular weight; Phys Chem, physicochemical properties; pka, negative log of the acid dissociation constant.
When several data were available, a weighted mean was calculated.
Determination of Aldehyde Oxidase Metabolic CLint.
The AO-mediated metabolism was measured in incubations containing either an HLC suspension at 1 mg protein/ml or HLS9 suspension at 2.5 mg protein/ml, both in phosphate buffer (100 mM), pH 7.4. The reactions were initiated by addition of prediluted compounds (2.5 µl from 100 µM in 100 mM Phosphate buffer/ACN/DMSO 90/9/1) to give a final nominal concentration of 1 µM. The solvent concentration did not exceed a total of 0.1%. The samples were then incubated at 37°C for either 120 minutes in HLC or 60 minutes in HLS9, with time points taken at 10, 30, 60, 90, and 120 minutes and 5, 10, 20, 40, and 60 minutes, respectively. The aliquots (25 µl) were precipitated with ACN (1 in 5 v/v) containing internal standard (historic AstraZeneca compound; AZ10024306) and centrifuged at 3500 rpm for 10 minutes, and the supernatant was diluted 1 in 7 (v/v) with ultra-pure HPLC water before analysis by liquid chromatography–tandem mass spectrometry (LC-MS/MS). All incubations were carried out in duplicate. The in vitro elimination rate constant corresponding to parent compound depletion was determined for each reaction using the first-order decay calculation in Microsoft Excel Sheet.
Determination of Unbound Fraction in Human Plasma.
The extent of binding of compounds to plasma proteins was determined by equilibrium dialysis at a compound concentration of 5 μM using the Rapid Equilibrium Device (Thermoscientific Pierce). Phosphate buffer (100 mM, pH 7.4) was added to the buffer chamber, and 300 μl of plasma was spiked with compound to the sample chamber. The unit was covered with a gas-permeable lid and incubated for 18 hours at 37°C at 300 rpm with 5% CO2. At the end of incubation, samples (50 μl) from both buffer and plasma chambers were removed for analysis. Samples and standards were matrix-matched and analyzed using LC-MS/MS. The unbound fraction in plasma (fu) was calculated as follows:(1)
Determination of Blood-to-Plasma Ratio.
A volume of plasma sufficient for the assay was obtained from whole human blood by centrifugation (3220 g for 10 minutes at 4°C). The test compound (10 µM) was added to 398 µl of the prewarmed human plasma and blood separately and incubated for 30 minutes. After incubation, the blood samples were centrifuged for 10 minutes at 3220 g (37°C), and the plasma samples were stored at 37°C. Aliquots (400 μl) of ice-cold acetonitrile containing internal standard were added to 100-μl samples of plasma separated from centrifuged whole blood and to reference plasma samples. These samples were then centrifuged, diluted with distilled water, and analyzed by LC-MS/MS to determine the compound concentration. Blood/plasma ratio (B/P) was calculated as follows:(2)
Determination of Unbound Fraction in HLM.
The extent of binding of compounds to HLM was determined by equilibrium dialysis using the HT Dialysis LLC device (Gales Ferry, CT) with HLM at a concentration of 1 mg protein/ml and a final compound concentration of 1 μM. PBS (150 μl) was added to the buffer well and 150 μl HLM containing the compound to the sample well and incubated at 37°C for 4 hours. After the incubation, 50-μl aliquots from both donor and receiver wells were removed for analysis. Samples and standards were matrix-matched and analyzed by LC-MS/MS. The unbound fraction in the incubation (fumic) was calculated as follows:(3)
LC-MS/MS Analysis.
The concentration of all compounds in the incubations was determined by LC-MS/MS. An Acquity ultra-performance liquid chromatography system (Waters, UK) coupled to a triple-quadrupole mass spectrometer (Xevo TQ-S; Waters, Milford, MA) was used to carry out the sample analysis. The details of quantification of analytes are described in Supplemental text. Detection of the ions was performed in the multiple reaction monitoring mode. Peak integration and calibrations were performed using TargetLynx software (Version 4.1; Waters).
Prediction of Intravenous Clearance Using PBPK Models.
The clinical trials providing the reference CLIV have been conducted in subjects with variable demographic characteristics (i.e., age range, proportion of females, healthy/patients with cancer). The specific demographics will influence some of the physiologic parameters (i.e., liver weight, plasma protein concentration) that in turn can impact the pharmokinetics (PK) parameters observed. Therefore, PBPK models were developed for each drug using the Simcyp Simulator V18R2, and the simulated trial designs and virtual population were selected accordingly to match the observed clinical trial (Table 1) (Kaye et al., 1984; Dolan et al., 1998; Rosen et al., 1999; Tserng et al., 2003; Miceli et al., 2005; Hutzler et al., 2012; Dalvie et al., 2013).
The CLint obtained in vitro from HLC and HLS9 fractions was corrected by the free fraction in the in vitro assay. The free fraction in HLM was measured using 1 mg/ml of microsomal protein. No clear trend concerning the difference in binding between HLC and HLM was observed, and therefore, the binding was assumed to stay the same in HLC and HLS9 fraction (Cubitt, 2009). When the protein concentrations used were different from the 1 mg/ml assessed in the binding experiments, the free fraction was extrapolated using the equation from Austin et al. (2002).
The B/P, fu and fumic, and in vitro metabolism data were used to develop PBPK models for each compound (Table 1). The main plasma-binding protein was assumed to be albumin for the acid, neutral, and ampholyte compounds and α-1 acid glycoprotein for the basic compounds. Physicochemical properties were gathered from literature sources and whole-body PBPK models with predicted volumes of distribution calculated using the Rodgers and Rowlands method were developed (Rodgers and Rowland, 2006) (). The contribution of microsomal metabolism and renal and biliary excretion to the clearance was added to the PBPK models when applicable (Fig. 1; Table 1).
Chemical structure of aldehyde oxidase substrate that have reported intravenous clearance.
Aldehyde oxidase is present in organs other than the liver (Moriwaki et al., 2001; Nishimura and Naito, 2006). To study the potential impact of extrahepatic metabolism, the activity per mg of cytosolic protein in the kidney was assumed to be the same as that of the liver, and the free intrinsic activity was scaled based on the human cytosolic protein per gram of kidney, kidney weight, and blood flow. A cytosolic protein per gram of kidney value of 40.6 mg/g was used (Scotcher, 2016). Similarly, lung metabolism was also explored, and the activity per mg of cytosolic protein was assumed to be the same as that of the liver. The IVIVE scaling approach (Fig. 2) using the well-stirred lung model (Yang, 2007) is integrated within the Simcyp Simulator and was simply entered as additional lung clearance; it was calculated with the following scaling parameters: cytosolic protein per gram of lung yield of 20 mg/g, a lung tissue weight (excluding blood) of 550 g, and a cardiac output of 386 l/h. Cytosolic protein per gram of lung was obtained using the S9 fraction in the lung of 28 mg/g tissue (Kozminski et al., 2019) and by assuming that the fraction of cytosolic protein to S9 protein is constant between the lung, liver, and kidney. Because of the limited expression of AO in the intestine (Moriwaki et al., 2001; Nishimura and Naito, 2006; Hutzler et al., 2012), intestinal metabolism was not considered in the current analysis.
In vitro–in vivo extrapolation for scaling AO metabolism in the lung. CLL, In-vivo plasma CLAO CL; CLuint, free intrinsic CL; CO, cardiac output; CPPGLu, cytosolic protein per gram of lung.
A linear regression between the predicted CLIV (dose/AUC0-infinity after a simulated single intravenous dose) and the observed CLIV was calculated using R (version 3.5.1, www.r-project.org) with a weighting option of 1/Y_pred to avoid bias toward the highest clearance value.
A sensitivity analysis was done on the AO intrinsic clearance of the six drugs to explore the impact of increasing the intrinsic clearance in the liver and in extrahepatic organs on the predicted CLIV. In the sensitivity analyses described here, the kidney was used as a surrogate organ to account for all the extrahepatic metabolism in the body. The kidney was chosen as the site of extrahepatic metabolism for practical reasons rather than splitting the clearance over several different organs. The ratio between the observed and predicted CLIV was calculated.
The intrinsic clearance of each compound was then optimized using available PK profiles except for BIBX1382, for which no concentration-time profile was available and only the clearance is reported in the literature. The following dosing regimens were used to simulate the PK profiles: O6-benzylguanine—bolus administration of 20 mg/m2 to seven patients with cancer aged 45–74 years (prop. of female = 0.42) (Tserng et al., 2003); carbazeran—10-minute infusion of 1.28 mg/kg of carbazeran to seven healthy male volunteers aged 20–50 years (Kaye et al., 1984); zaleplon—30-minute infusion of 5 mg to 10 healthy subjects aged 30–32 years (prop. of female = 0.5) (Rosen et al., 1999); ziprasidone—1-hour infusion of 5 mg to 13 male subjects aged 19–37 years (Miceli et al., 2005); and zoniporide—1-hour infusion of 80 mg to four male healthy subjects aged 18–55 years (Dalvie et al., 2010). Additionally, the Kp scalar was optimized for carbazeran (=0.13) and zoniporide (=0.45) to better fit the observed volume of distribution at steady state (Vss). The required (in silico) intrinsic clearance from the PBPK model was then compared with the measured (in vitro) intrinsic clearance to calculate a scaling factor for each compound.
Additionally, to verify the usefulness of this approach in predicting clearance of a new compound, an average scaling factor was applied to the in vitro intrinsic clearance of the studied drugs. The average scaling factor was calculated based on all the drugs except the one that was being predicted.
Results
The investigated compounds covered a wide range of fmAO, ranging from 0.064 (ziprasidone) to 0.98 (carbazeran). The log P values ranged from 1.04 (O6-benzylguanine) to 3.97 (BIBX1382), and there was one neutral compound, four basic compounds, and an ampholyte. Table 1 summarizes the physicochemical data and the measured CLint,u (free intrinsic clearance) values obtained in HLC, HLS9, and HLM. The results were overall in good agreement with literature data except for ziprasidone, in which the HLC CLint,u was more than 10-fold lower than previously reported values [33.4 vs. 410.3 µl/min per mg (Obach et al., 2012)] (Table 2).
Comparison of intrinsic clearance data obtained for AO from literature reports, in-house measurements, and retrograde scaling
Laboratory-specific scaling factors for six AO substrates.
Prediction of Intravenous Clearance Using PBPK Models.
Figure 3 compares the observed CLIV to the CLIV predicted from the PBPK model (Table 1).
Predicted CLIV (±S.D.) compared with mean observed CLIV (±S.D.). Results are compared with the unity line and 2- and 5-fold bias lines. HLS9, human liver S9; REF, reference.
HLC Liver Only.
Using HLC data, the average extent of underprediction of CLIV was 3.8 and ranged from 1.9 to 5.2 (Fig. 3). The best prediction was obtained for BIBX1382 (1.9-fold), and the biggest difference was for ziprasidone (5.2-fold). A coefficient of 2.77 was obtained using a weighted linear regression between the observed clearance and predicted clearance. For all the compounds, the average CV% for CL predicted by the PBPK models was 26% (range: 15.9%–37.5%). The predicted variability was in accordance with the mean observed CV of 21% (range: 5.7%–35.9%). Fig. 4 shows that no trend between the extent of underestimation and the predicted fm by AO could be observed.
Fold underprediction of clearance when considering only hepatic metabolism compared with the predicted fmAO.
HLS9 Liver Only.
Using HLS9 experiment, the CLIV was predicted to be lower than the CLIV obtained from HLC with an average extent of underprediction of 5.8, and it ranged from 2.9 (BIBX1382) to 10.6 (Zaleplon) (Fig. 3). A coefficient of 4.02 was obtained using a weighted linear regression between the observed clearance and predicted clearance. Simcyp workspaces for these six compounds using HLS9 data are provided as Supplemental Data Sets 1–6.
HLC Liver and Kidney.
When the metabolism in the liver and the kidney were included (assuming the same intrinsic clearance per mg of cytosolic protein), a small improvement of the prediction of the CLIV was observed. The extent of underprediction of CLIV was 3.5 and ranged from 1.4 (BIBX1382) to 5.2 (ziprasidone) (Fig. 3). A coefficient of 2.26 was obtained using a weighted linear regression between the observed clearance and predicted clearance.
HLC Liver, Kidney, and Lung.
When the metabolism in the liver, the kidney, and the lung were included (assuming the same activity per mg of cytosolic protein), a more pronounced increase in the accuracy of prediction of the CLIV was observed. The average extent of underprediction of CLIV was 3.2 and ranged from 0.9 (BIBX1382) to 5.2 (ziprasidone) (Fig. 3). A coefficient of 1.74 was obtained using a weighted linear regression between the observed clearance and predicted clearance.
Sensitivity Analyses.
Figure 5 shows the sensitivity analyses made on a scaling factor applied to the intrinsic clearance and considering the metabolism in liver and kidney. O6-benzylguanine and zaleplon have similar profiles and reach the observed CLIV with a scaling factor of around 15; BIXB1382 quickly reaches a clearance close to the observed clearance but also becomes quickly nonsensitive to any changes in CLint,u. Carbazeran reaches a plateau around 75% of the observed CLIV, and the blood flow is the limiting factor in this case. Ziprasidone CLint,u was a significantly lower value than reported in the literature (Table 2), and even with a 20-fold increase in the intrinsic clearance the observed CLIV is not attained. Zoniporide CLIV reaches the observed CLIV with a scaling factor of 3.5.
Sensitivity analyses on the intrinsic activity of AO in HLC assuming metabolism in the liver and the kidney.
PK Profiles.
Figure 6 shows the simulated profiles with and without optimization compared with the observed CLIV. The intrinsic clearance (CLint, AO) and scaling factor obtained are reported in Table 2. In this study, an average CLint, AO scaling factor of 6.5 was necessary to recover the PK profiles. Ziprasidone has a low fmAO, and the metabolic activity was lower in this study than previously reported, and for this reason ziprasidone was excluded from the calculation of the average scalar analysis.
Simulated mean profile with the default input in black (5th–95th percentile in gray) and after optimization in blue (5th–95th percentile in dashed light blue) assuming metabolism in the liver and kidney compared with observed concentrations (open circles). (A) Administration of a single 20 mg/m2–dose of O6-benzylguanine in intravenous bolus to seven patients with cancer aged 45–74 (prop. female = 0.42) (Tserng et al. 2003); (B) administration of a single 1.28-mg/kg i.v. infusion of carbazeran over 10 minutes to seven healthy male volunteers aged 20–50 years (Kaye et. al 1984); (C) administration of a single 5-mg i.v. infusion of zaleplon over 30 minutes to 10 healthy subjects aged 30–32 years (prop. female = 0.5) (Rosen et al. 1999); (D) administration of a single 5-mg i.v. infusion of ziprasidone over 1 hour to 13 male subjects aged 19–37 years (Miceli et al. 2005); (E) administration of a single 80-mg i.v. infusion of zoniporide over 1 hour to four male healthy subjects aged 18–55 years (Dalvie et al. 2010).
Table 3 shows that using the average scalar factor from the other drugs significantly improves the prediction of CLIV, with an average underprediction of 1.5 (range 0.98–1.96). All drugs were predicted within a 2-fold error, which is considered adequate for predictions early in drug discovery. Applying the scaling factor of 6.5 to ziprasidone, the CLIV was underpredicted by 3.9-fold, increasing the average fold underprediction to 1.9 for the six compounds.
Use of scaling factor from other compounds to predict the CL
The individual scaling factors are noted below the compounds; to predict the CL of a drug, the avg. scaling factor of the other drugs was used [i.e., the scalar used for 06-benzylguanine was 5.6 (avg. of 3.5, 4.35, 12, and 2.5)].
Discussion
Using IVIVE approaches to predict aldehyde oxidase–mediated clearance has typically resulted in a significant underprediction of the observed clearance, resulting in the clinical failure of multiple drugs that are metabolized by AO (Fan et al., 2016; Jensen et al., 2017). This study explored IVIVE of CLAO and aimed to develop a methodology to aid informed decision making on newly developed drug candidates that showed potential to be metabolized by AO, and a workflow for AO-mediated clearance prediction is proposed. Overall, using HLC, the intravenous CL continued to be underestimated, but HLC performed better than HLS9 in this study. Ideally the in vitro system of each laboratory should be characterized by measuring the CLint,u of the selected probe substrates present in this study, and then a scaling factor should be calculated using the information in Table 2 and applied on the intrinsic clearance assuming metabolism in the liver and the kidney as explained in the Materials and Methods section. The coefficient obtained with the linear regression could also be applied as an empirical scalar directly on the CLIV. If in-house probe substrate CLint data are not available, a scalar of 4.6 could be applied on the CLint,u based on the literature-published values of AO-mediated metabolism.
Measured B/P and fu were like values reported in the literature (Alousi et al., 2007; Zientek et al., 2010; Akabane et al., 2012). The average in vitro intrinsic clearances obtained from the literature were overall higher, especially for ziprasidone. One reason for this discrepancy in in vitro intrinsic clearances could be explained by the variability of activity in AO across human liver cytosolic fractions (Hutzler et al., 2014) or, alternatively, due to the nature of the experimental protocol (i.e., incubation and sampling time, buffer). It is possible that human donor liver tissues that were preserved in allopurinol-containing University of Wisconsin solution may exhibit weak AO inhibitory effect at higher concentrations and could affect the scaling factors calculated in this study. Nevertheless, Barr et al. (2014) have shown that small residual amounts of allopurinol or oxypurinol did not appear to impact AO activity. By analyzing the different protocols used in the literature, no clear association between experimental condition and measured CLint was observed, and a larger data set or new specific in vitro assays looking at the impact of the protocol would be required (Dick, 2018).
In this data set, the fraction metabolized by AO was estimated by including all known pathways; however, even if some inhibitors of aldehyde oxidase have been identified in vitro (Johns, 1967; Johnson et al., 1985; Obach, 2004), limited clinical drug interactions via inhibition of aldehyde oxidase have been recorded, thus restricting the additional validation of in vivo probe substrates or inhibitors for the verification of the fmAO of given substrate.
Some compounds have an intravenous clearance higher than the hepatic blood flow, suggesting extrahepatic elimination. Aldehyde oxidase is expressed in multiple tissues, including the kidney (Moriwaki et al., 2001; Nishimura and Naito, 2006). A scenario assuming an AO activity/expression per mg of cytosolic protein in the kidney identical to that in the liver was simulated in the PBPK models, and even though an improved prediction was observed, the underprediction of CL was still significant. Expression of AO in the intestine is limited, and incubations of AO substrate with human intestinal cytosol resulted in no measurable metabolism (Moriwaki et al., 2001; Nishimura and Naito, 2006; Hutzler et al., 2012); therefore, intestinal metabolism was not considered in the current analysis. The lung is a highly perfused organ with a significant tissue volume and an AO absolute abundance within 6-fold of that in the liver (Ezkurdia et al., 2015). Previous studies have attempted to incorporate lung metabolism into IVIVE approaches for AO-mediated clearance (Kozminski et al., 2019); however, in this study the reported intrinsic activity in the lung was almost 1000-fold lower than that in the liver, explaining the absence of significant impact of lung metabolism on the overall predicted clearance (Kozminski et al., 2019). In this study, even assuming a tissue activity per mg of lung cytosol identical to that in the liver did not explain the underestimation of clearance observed for all compounds (Fig. 3).
Aldehyde oxidase has been shown to have a limited stability, and freeze-thaw cycles might result in higher variability, and therefore, the metabolic activity might be underestimated in in vitro assays (Sherratt and Damani, 1989; Hutzler et al., 2012). The underestimation of AO clearance is likely to be due to an underestimation of the intrinsic clearance as well as extrahepatic metabolism.
An additional way to scale from in vitro to in vivo would be to use the absolute abundance rather than activity per mg of protein. This approach would allow the use of recombinant AO and, therefore, have an extremely specific system with less risk of contamination from other enzymes (i.e., xanthine oxidase). In addition, the absolute abundance in all of the different tissues in the body could be accounted for with a single in vitro metabolism measurement. So far, the absolute abundance in the liver has been measured (Barr et al., 2013; Fu et al., 2013; Ezkurdia et al., 2015; Wiśniewski et al., 2016), and recombinant AOs are available. Unfortunately, there is a lot of variability between the different laboratories (mean: 34.5 pmol/mg of cytosolic protein; range: 1.41–60.2 pmol/mg of cytosolic protein, four studies, a total of 30 livers), and the absolute abundance has not been measured in the recombinant systems.
An additional aim of this work was to gather the input information for PBPK models for the different AO substrate compounds so that the models could be available for use in future research efforts. The PBPK models could be used to investigate different aspects, such as interindividual differences in AO expression; to study the interaction between AO substrates and inhibitors [e.g., between zaleplon and cimetidine (Dalvie and Di, 2019)]; or to investigate the PK of these compounds in different populations of individuals.
Conclusion
A workflow for NCE metabolized by AO was suggested, and an empirical scaling factor of three on the predicted CLIV based on HLC data could be applied for NCEs that are significantly metabolized by AO when using PBPK models for predicting the exposure of NCEs in the human. Alternatively, a scaling factor of 6.5 could be applied to the AO intrinsic clearance in the liver and kidney. Ideally each laboratory should develop a correlation using a set of probe substrates under their own assay conditions; however, if the in vitro CLint,u values for probe substrates are not available in a given laboratory, an empirical scaling factor of 4.6 based on this work could be applied for CLint,u in HLC. Additional research on the impact of the in vitro study designs and extrahepatic metabolism is suggested to understand the mechanism behind the systematic underprediction observed for AO.
Authorship Contributions
Participated in research design: De Sousa Mendes, Gardner, Neuhoff, Pilla Reddy.
Conducted experiments: Orton, Jones.
Performed data analysis: De Sousa Mendes, Pilla Reddy.
Wrote or contributed to the writing of the manuscript: De Sousa Mendes, Orton, Humphries, Jones, Gardner, Neuhoff, Pilla Reddy.
Footnotes
- Received April 16, 2020.
- Accepted July 22, 2020.
↵1 M.D.S.M. and V.P.R. contributed equally as joint first authors.
No funding was received for this work.
V.P.R., A.O., and B.J. are full-time employees of AstraZeneca when this study was conducted and hold shares of AstraZeneca. M.D.S.M., H.E.H., I.G., and S.N. employees of Certara UK Limited when this study was conducted.
↵
This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- ACN
- acetonitrile
- AO
- aldehyde oxidase
- AUC
- area under the curve
- B/P
- blood/plasma ratio
- CL
- clearance
- CLAO
- CL of aldehyde oxidase substrates
- CLint
- intrinsic CL
- CLIV
- intravenous CL
- fm
- fraction metabolized
- fu
- unbound fraction in plasma
- HLC
- human liver cytosol
- HLM
- human liver microsome
- HLS9
- human liver S9
- HPLC
- high-pressure liquid chromatography
- IVIVE
- in vitro–to–in vivo extrapolation
- LC-MS/MS
- liquid chromatography–tandem mass spectrometry
- NCE
- new chemical entity
- PBPK
- physiologically based pharmacokinetic
- PK
- pharmacokinetics
- prop
- proportion
- Copyright © 2020 by The American Society for Pharmacology and Experimental Therapeutics