Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • Log out
  • My Cart

Search

  • Advanced search
Drug Metabolism & Disposition
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • Log out
  • My Cart
Drug Metabolism & Disposition

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Visit dmd on Facebook
  • Follow dmd on Twitter
  • Follow ASPET on LinkedIn
Research ArticleArticle

An Investigation into the Prediction of the Plasma Concentration-Time Profile and Its Interindividual Variability for a Range of Flavin-Containing Monooxygenase Substrates Using a Physiologically Based Pharmacokinetic Modeling Approach

Venkatesh Pilla Reddy, Barry C. Jones, Nicola Colclough, Abhishek Srivastava, Joanne Wilson and Danxi Li
Drug Metabolism and Disposition September 2018, 46 (9) 1259-1267; DOI: https://doi.org/10.1124/dmd.118.080648
Venkatesh Pilla Reddy
Departments of Modelling and Simulation, Oncology Drug Metabolism and Pharmacokinetics (V.P.R.), Departments of Drug Metabolism and Pharmacokinetics and Oncology (B.C.J., N.C., J.W.), and Department of Drug Safety and Metabolism (A.S.), IMED Biotech Unit, AstraZeneca, Cambridge, United Kingdom; and Pharmaron, Beijing, China (D.L.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Barry C. Jones
Departments of Modelling and Simulation, Oncology Drug Metabolism and Pharmacokinetics (V.P.R.), Departments of Drug Metabolism and Pharmacokinetics and Oncology (B.C.J., N.C., J.W.), and Department of Drug Safety and Metabolism (A.S.), IMED Biotech Unit, AstraZeneca, Cambridge, United Kingdom; and Pharmaron, Beijing, China (D.L.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicola Colclough
Departments of Modelling and Simulation, Oncology Drug Metabolism and Pharmacokinetics (V.P.R.), Departments of Drug Metabolism and Pharmacokinetics and Oncology (B.C.J., N.C., J.W.), and Department of Drug Safety and Metabolism (A.S.), IMED Biotech Unit, AstraZeneca, Cambridge, United Kingdom; and Pharmaron, Beijing, China (D.L.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Abhishek Srivastava
Departments of Modelling and Simulation, Oncology Drug Metabolism and Pharmacokinetics (V.P.R.), Departments of Drug Metabolism and Pharmacokinetics and Oncology (B.C.J., N.C., J.W.), and Department of Drug Safety and Metabolism (A.S.), IMED Biotech Unit, AstraZeneca, Cambridge, United Kingdom; and Pharmaron, Beijing, China (D.L.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joanne Wilson
Departments of Modelling and Simulation, Oncology Drug Metabolism and Pharmacokinetics (V.P.R.), Departments of Drug Metabolism and Pharmacokinetics and Oncology (B.C.J., N.C., J.W.), and Department of Drug Safety and Metabolism (A.S.), IMED Biotech Unit, AstraZeneca, Cambridge, United Kingdom; and Pharmaron, Beijing, China (D.L.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Danxi Li
Departments of Modelling and Simulation, Oncology Drug Metabolism and Pharmacokinetics (V.P.R.), Departments of Drug Metabolism and Pharmacokinetics and Oncology (B.C.J., N.C., J.W.), and Department of Drug Safety and Metabolism (A.S.), IMED Biotech Unit, AstraZeneca, Cambridge, United Kingdom; and Pharmaron, Beijing, China (D.L.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF
Loading

Abstract

Our recent paper demonstrated the ability to predict in vivo clearance of flavin-containing monooxygenase (FMO) drug substrates using in vitro human hepatocyte and human liver microsomal intrinsic clearance with standard scaling approaches. In this paper, we apply a physiologically based pharmacokinetic (PBPK) modeling and simulation approach (M&S) to predict the clearance, area under the curve (AUC), and Cmax values together with the plasma profile of a range of drugs from the original study. The human physiologic parameters for FMO, such as enzyme abundance in liver, kidney, and gut, were derived from in vitro data and clinical pharmacogenetics studies. The drugs investigated include itopride, benzydamine, tozasertib, tamoxifen, moclobemide, imipramine, clozapine, ranitidine, and olanzapine. The fraction metabolized by FMO for these drugs ranged from 21% to 96%. The developed PBPK models were verified with data from multiple clinical studies. An attempt was made to estimate the scaling factor for recombinant FMO (rFMO) using a parameter estimation approach and automated sensitivity analysis within the PBPK platform. Simulated oral clearance using in vitro hepatocyte data and associated extrahepatic FMO data predicts the observed in vivo plasma concentration profile reasonably well and predicts the AUC for all of the FMO substrates within 2-fold of the observed clinical data; seven of the nine compounds fell within 2-fold when human liver microsomal data were used. rFMO overpredicted the AUC by approximately 2.5-fold for three of the nine compounds. Applying a calculated intersystem extrapolation scalar or tissue-specific scalar for the rFMO data resulted in better prediction of clinical data. The PBPK M&S results from this study demonstrate that human hepatocytes and human liver microsomes can be used along with our standard scaling approaches to predict human in vivo pharmacokinetic parameters for FMO substrates.

Introduction

Flavin-containing monooxygenase (FMO) belongs to a family of enzymes that oxygenate a wide variety of nucleophilic heteroatom-containing chemicals and drugs. FMO3 is the predominant functional isoform in adult human liver, while FMO1 and 5 are expressed in adult kidney and liver, respectively (Yeung et al., 2000). FMO enzymes can be differentiated from cytochrome P450 (P450) enzymes via heat inactivation and selective substrate inhibition by methimazole (Cashman and Zhang, 2006; Taniguchi-Takizawa et al., 2015). There is a lack of validated methods for extrapolating in vitro hepatic intrinsic clearance (CLint) for FMO to in vitro/in vivo extrapolation (IVIVE) clearance (CL). Our earlier paper (Jones et al., 2017) investigated a set of 10 literature compounds with varying degrees of FMO involvement. The compounds were profiled in in vitro assays, which defined the extent of FMO involvement based on heat inactivation and inhibition by the selective substrate methimazole. The FMO contribution varied from 21% for imipramine to 96% for itropride. In addition, these data were used to predict the unbound intrinsic clearance and compared the results with human clearance data obtained from the literature. Using standard scaling scaled methods, 70% of the compounds studied had predicted unbound intrinsic clearance within 2-fold of the observed unbound intrinsic clearance.

In this paper, we apply a mechanistic, physiologically based pharmacokinetic (PBPK) modeling and simulation approach (M&S) to predict the observed clinical data of drugs that are metabolized by FMO. The utility of PBPK modeling via a top-down approach [typically involving the estimation of the model parameters using observed clinical pharmacokinetic (PK) profiles] for FMO was shown for itopride (Zhou et al., 2017), which evaluated the impact of FMO3 polymorphism on exposure in Asian subjects. This study showed that genotyping for FMO3 to exclude subjects with homozygous Lys158/Gly308 alleles is not clinically relevant in terms of exposure, and thus is not required. However, the elimination kinetics of itopride from in vitro studies using human liver microsome (HLM) data resulted in about 8-fold underprediction of itopride clearance (Obach et al., 1997; Humphries et al., 2015; Zhou et al., 2017). This was similar to the FMO substrate benzydamine, where underprediction was reported to be about 1.5- to 11-fold (Obach et al., 1997; Humphries et al., 2015; Zhou et al., 2017). Thus, an additional investigation to achieve a better in vitro/in vivo extrapolation was warranted. As a result of the aforementioned observations, Jones et al. (2017) generated in vitro data under well-defined experimental conditions. To date, there has not been a comprehensive assessment of FMO scaling using in vitro data and a bottom-up PBPK M&S approach (which typically uses in vitro data along with mechanistic understanding of absorption, distribution, metabolism, and excretion processes of a drug to prospectively simulate the PK profiles) for substrates with differing FMO contributions.

A human PBPK model, once established, can be a valuable tool to inform the clinical study design of an investigational drug under development, and help to link the exposure from healthy subjects to patients by taking into account pharmacogenetics and extrinsic factors or simulating the PK profiles in special populations such as pediatric subjects. Hence, the aim of this study was to investigate the performance of PBPK M&S to predict the concentration-time profiles, PK parameters [area under the curve (AUC), Cmax and CL], and interindividual variability for a set of compounds for which all or part of the in vivo clearance has been determined to be mediated by FMO.

Materials and Methods

Compound Selection and Clinical Data Collection

Nine orally administered drugs, namely, itopride, benzydamine, tozasertib, tamoxifen, moclobemide, imipramine, clozapine, olanzapine, and ranitidine were selected for PBPK modeling based on the availability of the in vitro and in vivo data within AstraZeneca or in the public domain. Clinical plasma profiles were digitized using an AstraZeneca internal tool namely, GI-SIM (Sjögren et al., 2013), if literature data were used. The detailed clinical trial information for each of the nine drugs is described in Table 1. The PK data and parameters for these drugs were either extracted from the literature or obtained from internal clinical studies (e.g., tamoxifen).

View this table:
  • View inline
  • View popup
TABLE 1

Clinical data used for PBPK model verification

PBPK Modeling Approach

Whole-body PBPK M&S of clinical PK data was performed using the population-based absorption, distribution, metabolism, and excretion simulator, Simcyp version 16, release 1 (Certara, Sheffield, United Kingdom) (Jamei et al., 2013). The modeling assumptions and parameters were verified using clinical data for the nine FMO substrates.

PBPK Modeling Input Parameters

System-Related Parameters.

Among the five FMOs, only FMO1, 3, and 5 are functional in human. FMO1 is a major form in human kidney (Yeung et al., 2000), while FMO3 and FMO5 are expressed mainly in liver (Dolphin et al., 1996; Koukouritaki and Hines, 2005). The reported FMO1 (kidney), FMO3 (liver), and FMO5 (gut) protein abundance values of 47, 71, and 22 pmol/mg protein were used in this analysis along with the reported coefficient of variance (Haining et al., 1997; Overby et al., 1997; Yeung et al., 2000; Chen et al., 2016). Enzyme abundance, phenotype frequencies, and relative activities of FMO3 were accounted for within the virtual population of the simulator. The intrinsic catalytic activity per unit amount of FMO1 and FMO3 enzyme was assumed to be the same in healthy Caucasian, Japanese, and Chinese subjects since the available expression data pool broadly represents the Caucasian and Asian populations (Cashman and Zhang, 2006).

The Simcyp platform provides an enzyme kinetics pane for assessing the contribution to metabolic clearance in the liver, kidney, or intestine for P450 and UDP‐glucuronosyltransferase (UGT) enzymes, but not yet for FMOs. However, a user-defined enzyme option, or unused UGT enzyme, within the UGT enzyme kinetic pane allows entry of recombinant FMO (rFMO) CLint data, along with expression levels, frequency of polymorphism, and turnover rate (Reddy et al., 2010), to define the FMO enzyme to enable IVIVE.

Drug-Related Parameters.

The PBPK model input parameters used are listed in Table 2. The in vitro parameters such as CLint, plasma protein binding (unbound fraction in plasma), and fraction unbound in hepatocyte incubation were experimentally measured (Jones et al., 2017); in some instances, the PBPK platform’s prediction toolbox was employed to predict the drug-related properties such as the effective human permeability (Peff) value using the mechanistic permeability model or the polar surface area/hydrogen bond donor model if permeability data were not available. The first-order absorption rate constant and fraction absorbed were estimated from clinical data or predicted based on the Peff value within the Simcyp simulator. Unless otherwise specified, the volume of distribution at steady state for all compounds was predicted using the Rodgers and Rowland method (Rodgers et al., 2005; Rodgers and Rowland, 2006) within the PBPK simulator. The fraction unbound in microsomes and hepatocytes (and fraction unbound in hepatocyte incubation) was also included in the model.

View this table:
  • View inline
  • View popup
TABLE 2

PBPK model input parameters

Intrinsic Clearance Measurement Using Hepatocytes and Human Liver Microsomes.

The measured in vitro data are detailed in Supplemental Table 1, and for methodology refer to Jones et al. (2017). The assignment of FMO metabolic intrinsic clearance and methodology employed to assess drug elimination are given Supplemental Table 2.

Intrinsic Clearance Measurement in rFMO1, rFMO3, and rFMO5.

The compounds (final concentration 1 μM) were incubated with rFMO1, rFMO3, and rFMO5 (final concentration 0.5 mg/ml protein). Incubations were carried out in duplicate in 100 mM potassium phosphate buffer (pH 7.4) containing NADPH (1 mM) in a total volume of 0.2 ml. Incubations were commenced with the addition of NADPH and carried out in a 96-well plate at a temperature of 37°C. At times ranging between 0 and 50 minutes, aliquots (0.02 ml) were removed and added to 0.08 ml of termination solution (acetonitrile + formic acid) containing buspirone (4 nM) as an internal standard. Samples were vortexed and centrifuged at 3500 rpm for 5 minutes, and the resulting supernatant was collected for liquid chromatography–tandem mass spectrometry analysis as shown in the Supplemental Material and Supplemental Table 5.

Estimating the Tissue-Specific Scalar or Intersystem Extrapolation Factor for Extrapolation Using rFMO Data.

Human PK predictions using PBPK models and in vitro assays such as hepatocytes or HLM data combined with the appropriate scalars appears to be robust for P450 enzymes (Wagner et al., 2015; Luzon et al., 2017). However, very limited knowledge of scaling factors exists for drugs that are metabolized by non-P450s such as FMOs. A drug-specific intersystem extrapolation factor (ISEF) value can be derived using the following formula (see Supplemental Table 1):Embedded Image(1)When developing a model based on the rFMO data, the elimination via FMO1, FMO3, and FMO5 enzymes was considered separately along with respective liver, kidney, and intestine expression levels. The intrinsic catalytic activity per unit amounts of FMO1, FMO3, and FMO5 enzymes were 71 pmol/mg protein in liver, 47 pmol/mg protein in kidney, and 22 pmol/mg protein in intestine to reflect the abundances of FMO3, FMO1, and FMO5, respectively. These abundance values were included in the population file using unused UGT enzyme with the user-defined option within Simcyp (Supplemental Table 2). The impact of FMO3 polymorphism was included in the final PBPK model. The remaining metabolic CL was accounted for via major metabolizing enzymes of the drug via HLM in enzyme kinetics or via recombinant P450 if P450 contribution was known, as shown in Table 2.

The automated sensitivity analysis (ASA) or parameter estimation (PE) option within Simcyp was used to determine the tissue-specific rFMO scalar in liver, intestine, and kidney tissue or to determine the ISEF value. To this end, we used clinical data of the drugs that are mainly metabolized by FMO such as itopride. When the PBPK model was optimized to the clinical data of itopride using the ASA approach, the lower and upper bounds provided for the rFMO scalar parameter were 0–100 and the used log-distributed step size with 10 steps (increasing with increased parameter value) within the Simcyp PBPK platform. The effect of the rFMO scalar on the Cmax, AUC, CL, and PK profiles were reviewed to understand the elimination trends of a drug. In addition to ASA, PE was performed using the weighted least-squares objective function to allow fitting to the mean profile. Weighting by the reciprocal of the prediction was used since this provides a good fit for both Cmax and the terminal phase when oral data are fitted to the model. The Nelder-Mead fitting algorithm was used with the default settings (reflection coefficient = 1, expansion coefficient = 2, contraction coefficient = 0.5, and shrink coefficient = 0.5).

PBPK Model Verification.

PK simulations were conducted during the PBPK model development as a verification step to ensure appropriate model input parameters were used. The predicted PK parameters and simulated concentration-time profiles were compared with those observed in clinical trials (Table 1). The same trial designs as the clinical studies (e.g., dose, number of subjects, age, gender, population, etc.) were used in the simulations. When a specific population was not available in the Simcyp simulator, for example, the Korean population, an alternative validated population such as the Japanese population was used instead. In the case of olanzapine and clozapine, we assumed no difference between healthy subjects and patients since no validated schizophrenic population was available.

In addition to verification of the PK parameters and PK profiles, we also investigated the nature of the correlation between human in vitro and in vivo clearance predicted by the PBPK approach using hepatocyte, HLM, and rFMO data with ISEF and tissue-specific scalars.

Results

PBPK Model Development and Optimization of PBPK Models.

The plasma concentration profiles for nine orally or intravenously administrated FMO substrates were simulated by PBPK modeling using the hepatocyte CLint data and were compared with those observed clinically, as shown in Fig. 1. By visual inspection, the developed PBPK models adequately simulated the observed clinical data, and thus it was assumed that the models reflect the metabolic routes of elimination based on in vitro hepatocyte or HLM data (Fig. 2), along with any known additional pathways such as renal or biliary clearance. The observed and model-simulated PK parameters such as AUC, Cmax, and CL are shown in Table 3. Table 3 shows that using hepatocyte CLint data the PBPK model successfully predict the AUC and Cmax for all nine drugs within a factor of 2 compared with the clinical data. When HLM CLint data were used seven of the nine drugs had these parameters correctly predicted within a factor of 2 (Fig. 2; Supplemental Table 3). Use of rFMO CLint data over- or underpredicted the observed clinical data for five of the nine compounds as depicted in the forest plot using the AUC as an end point (Fig. 2) when no tissue-specific factor (zero value) or ISEF option was used. The drug-specific ISEF values for each of the rFMO isoforms determined using eq. 1 are shown in Supplemental Table 1.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Mean simulated using hepatocyte CLint data (solid line) and observed (data points) concentrations of FMO substrates after administration of a single dose to humans for nine FMO substrates. The gray lines around the solid black line represent simulated individual trials matched to a clinical study. Lower and upper gray lines represent the 90% confidence interval of the respective simulations.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Forest plot showing the PBPK modeling performance of FMO substrates. Predictions are expressed as ratios of predicted over observed AUC. The dashed line represents the identity (predicted/observed ratio), the gray shaded represents the 0.5–2.0 ratio window, and the blue shaded area represents the 0.67- to 1.50-fold ratio. Mean AUC using hepatocytes (blue squares); mean AUC using HLM (red circles); mean AUC using rFMO with no scaling (green triangles), with error bars represented as percentile range (fifth and 95th percentiles).

View this table:
  • View inline
  • View popup
TABLE 3

Observed mean PK parameters for FMO substrates using hepatocyte data

Based on the ASA and PE analyses using itopride clinical data, an rFMO tissue scalar value of 0.2 or a compound-specific ISEF appears to be a reasonable value to obtain a good IVIVE, as shown in Fig. 3. Performance of the tozasertib PBPK model with an rFMO scalar value is shown as an example in Fig. 3. Figure 3, A–C shows the model fit to observed data without and with the ISEF or tissue scalar for itopride, while Fig. 3, D–F shows the model fit without and with ISEF or rFMO scalar for tozasertib, respectively. This generic scalar value of 0.2 was applied for each of the FMO1, FMO3, and FMO5 isoforms and for all other compounds. This optimized rFMO scalar also resulted in better predictions of the observed PK profiles (Fig. 3) with an increased correlation value (r2) for the correlation of predicted versus observed clearance, as shown in Fig. 4. The forest plot showing the model performance using rFMO data of FMO substrates with and without scalar for all compounds is summarized in Supplemental Fig. 1. The results of the PK parameters for these two options are given in Supplemental Table 3.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Determination of rFMO scalar for itopride and its application to tozasertib as an example. (A–C) Model fit without rFMO scalar, with rFMO tissue scalar, and with rFMO ISEF scalar for itopride, respectively. (D–F) Tozasertib model without rFMO scalar, with rFMO tissue scalar, and with rFMO ISEF scalar, respectively.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Correlations for the PBPK model–predicted clearance using human hepatocytes, HLM, and rFMO with ISEF scalar and tissue-specific scalar.

Discussion

The prediction of human PK properties can be critical in decision making when ranking short listed compounds, the assessment of lead series, and for clinical drug candidates. A bottom-up PBPK modeling approach involves modeling of the mechanisms of the drug disposition processes and simulating the human PK profiles by accounting for variability arising from enzyme expression at the patient level. The application of bottom-up models depends on the availability of quality data and established IVIVE factors and scalars. Human PK predictions using PBPK models and data from in vitro assays such as hepatocytes or HLM CLint combined with the appropriate scalars appears to be robust for P450 enzymes (Wagner et al., 2015; Luzon et al., 2017). However, very limited knowledge of scaling factors exists for drugs that are metabolized by non-P450s such as FMOs. As a result, early bottom-up models will require verification with in vivo clinical data, and in some cases calibration of parameters including scaling factors through a middle-out approach (i.e., combining clinical and in vitro data) to be considered as more impactful in the drug development process. A method that many consider would provide greater certainty is to compare the predictions from two or more methods. An example would be comparing two predictions based on the same input data set; in this case, comparing the results shown in our earlier study (Jones et al., 2017) versus the PBPK modeling approach.

In this paper, we use a PBPK M&S approach to predict the PK parameters, to verify and simulate the observed clinical data of FMO substrates, and to compare the results with our earlier IVIVE results (Jones et al., 2017). Human hepatocytes produced good concordance (R2 = 0.72) between the predicted and observed in vivo clearance values (Fig. 4), while all drugs fell within 2-fold. The correlation between the predicted and observed in vivo clearance was weaker for HLM and rFMO, which is consistent with our earlier study (Jones et al., 2017).

The PBPK results using rFMO in vitro data deviated from the observed exposure by approximately 2.5-fold, when a scalar of zero was used. This biased prediction may not be surprising since individual rFMO and HLM scaling factors are required to compensate for differences in the intrinsic activity per unit enzyme of rFMO relative to HLM. The observed clinical data can be used to refine and/or improve the model performance via a middle-out approach (combination of top-down and bottom-up approaches). We estimated a scaling factor to convert enzyme kinetic values obtained in different expression systems to the equivalent in HLM.

This rFMO scalar concept is well established for recombinant P450s within PBPK platforms such as Simcyp and has been shown to be applicable for recombinant UGTs as well (Lin, 2015). In addition to a compound-specific ISEF (Supplemental Table 1), an attempt was made to generate a tissue-specific scalar similar to that which exists for UGT enzymes within the simulator. To determine the scalar value associated with a particular rFMO isoform we performed sensitivity and PE analyses using itopride clinical data (Fig. 3). The ASA approach suggested that an rFMO scalar value of around 0.2 is required to recover the itopride clinical data. This rFMO scalar of 0.2 also resulted in better predictions for other compounds, as shown in Fig. 3 and Supplemental Fig. 1, than when the rFMO scalar of zero was used. This analysis shows that the observed data can be used to refine the model performance via fitting the clinical data to estimate a few of the model parameters (in this case, the rFMO scalar) to improve the model’s predictive performance of other FMO substrates using data of a corresponding rFMO in vitro assay.

When developing the model based on the rFMO data the elimination via FMO1, FMO3, and FMO5 was considered separately along with respective liver, kidney, and intestine expression levels. The intrinsic catalytic activity per unit amounts of FMO1, FMO3, and FMO5 enzymes were 71 pmol/mg protein in liver, 47 pmol/mg protein in kidney, and 22 pmol/mg protein in intestine to reflect the abundances of FMO3, FMO1 and FMO5, respectively, these amounts were included in the population file. The impact of FMO3 polymorphisms was included in the final PBPK model.

PBPK model approaches allow the simulation of PK profiles by taking into account parameter uncertainty. This uncertainty reflects both uncertainty in the experimental values of input parameters in a model along with variability arising from FMO expression and demographics. This is unlike simple IVIVE methods that do not take variability arising from expression, demographics, and experimental errors into account (Jones et al., 2017). Thus, the correlation values from PBPK modeling give more confidence in the scaling approach on top of simple IVIVE methods.

Once the PBPK model is verified it can be applied to predict the drug exposure differences in subjects with FMO3 wild-type and homozygous FMO Lys158/Gly308 mutant genotypes, as shown in Zhou et al. (2017) for itopride. In addition, we can apply the developed PBPK model to special populations, such as pediatrics, since FMO abundance data recently quantified will allow accounting for an ontogeny.

The ontogeny of hepatic FMO exhibits a similar trend to the CYP3A family, where FMO1 resembles the pattern of CYP3A7 and FMO1 and FMO3 expression data have been quantified using a HLM system (Koukouritaki et al., 2002). Moreover, it was shown that FMO3 matures to the adult level by 10 years of age (Hines and McCarver, 2002; Koukouritaki et al., 2002), requiring a dose adjustment because of this ontogeny process. Currently, the pediatric clinical data for FMO substrates is lacking, but a prospective simulation for itopride across the age range from 6 months to 12 years of age is shown in Supplemental Table 4.

Our results confirm that clearance by FMO in human can be assessed and derisked using a set of in vitro drug metabolism and PK assays. When combined with PBPK modeling, it allows for more robust prediction of AUC and other PK parameters further supporting derisking during the drug development process.

Acknowledgments

We thank Andy Sykes, Helen Humphries, Iain Gardner, and Sibylle Neuhofff for scientific input in carrying out the PBPK analyses described in this paper.

Authorship Contributions

Participated in research design: Pilla Reddy, Jones, Colclough.

Conducted experiments: Srivastava, Wilson, Li.

Performed data analysis: Pilla Reddy, Jones.

Wrote or contributed to the writing of the manuscript: Pilla Reddy, Jones, Srivastava, Colclough, Wilson, Li.

Footnotes

    • Received January 25, 2018.
    • Accepted June 7, 2018.
  • https://doi.org/10.1124/dmd.118.080648.

  • ↵Embedded ImageThis article has supplemental material available at dmd.aspetjournals.org.

Abbreviations

ASA
automated sensitivity analysis
AUC
area under the curve
CL
clearance
CLint
intrinsic clearance
FMO
flavin-containing monooxygenase
HLM
human liver microsomes
ISEF
intersystem extrapolation factor
IVIVE
in vitro/in vivo extrapolation
M&S
modeling and simulation approach
P450
cytochrome P450
PBPK
physiologically based pharmacokinetic
PE
parameter estimation
PK
pharmacokinetic
rFMO
recombinant FMO
UGT
UDP‐glucuronosyltransferase
  • Copyright © 2018 by The American Society for Pharmacology and Experimental Therapeutics

References

    1. Albers LJ,
    2. Reist C,
    3. Vu RL,
    4. Fujimoto K,
    5. Ozdemir V,
    6. Helmeste D,
    7. Poland R, and
    8. Tang SW
    (2000) Effect of venlafaxine on imipramine metabolism. Psychiatry Res 96:235–243.
    OpenUrlCrossRefPubMed
    1. Baldock GA,
    2. Brodie RR,
    3. Chasseaud LF,
    4. Taylor T,
    5. Walmsley LM, and
    6. Catanese B
    (1991) Pharmacokinetics of benzydamine after intravenous, oral, and topical doses to human subjects. Biopharm Drug Dispos 12:481–492.
    OpenUrlPubMed
    1. Brøsen K and
    2. Gram LF
    (1988) First-pass metabolism of imipramine and desipramine: impact of the sparteine oxidation phenotype. Clin Pharmacol Ther 43:400–406.
    OpenUrlCrossRefPubMed
  1. ↵
    1. Cashman JR and
    2. Zhang J
    (2006) Human flavin-containing monooxygenases. Annu Rev Pharmacol Toxicol 46:65–100.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Chen Y,
    2. Zane NR,
    3. Thakker DR, and
    4. Wang MZ
    (2016) Quantification of flavin-containing monooxygenases 1, 3, and 5 in human liver microsomes by UPLC-MRM-based targeted quantitative proteomics and its application to the study of ontogeny. Drug Metab Dispos 44:975–983.
    OpenUrlAbstract/FREE Full Text
    1. Chiu CC,
    2. Lane HY,
    3. Huang MC,
    4. Liu HC,
    5. Jann MW,
    6. Hon YY,
    7. Chang WH, and
    8. Lu ML
    (2004) Dose-dependent alternations in the pharmacokinetics of olanzapine during coadministration of fluvoxamine in patients with schizophrenia. J Clin Pharmacol 44:1385–1390.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Dolphin CT,
    2. Cullingford TE,
    3. Shephard EA,
    4. Smith RL, and
    5. Phillips IR
    (1996) Differential developmental and tissue-specific regulation of expression of the genes encoding three members of the flavin-containing monooxygenase family of man, FMO1, FMO3 and FM04. Eur J Biochem 235:683–689.
    OpenUrlPubMed
    1. Elshafeey AH,
    2. Elsherbiny MA, and
    3. Fathallah MM
    (2009) A single-dose, randomized, two-way crossover study comparing two olanzapine tablet products in healthy adult male volunteers under fasting conditions. Clin Ther 31:600–608.
    OpenUrlCrossRefPubMed
    1. Fuchs WS,
    2. Leary WP,
    3. van der Meer MJ,
    4. Gay S,
    5. Witschital K, and
    6. von Nieciecki A
    (1996) Pharmacokinetics and bioavailability of tamoxifen in postmenopausal healthy women. Arzneimittelforschung 46:418–422.
    OpenUrlPubMed
    1. Garg DC,
    2. Eshelman FN, and
    3. Weidler DJ
    (1985) Pharmacokinetics of ranitidine following oral administration with ascending doses and with multiple-fixed doses. J Clin Pharmacol 25:437–443.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Haining RL,
    2. Hunter AP,
    3. Sadeque AJ,
    4. Philpot RM, and
    5. Rettie AE
    (1997) Baculovirus-mediated expression and purification of human FMO3: catalytic, immunochemical, and structural characterization. Drug Metab Dispos 25:790–797.
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Hines RN and
    2. McCarver DG
    (2002) The ontogeny of human drug-metabolizing enzymes: phase I oxidative enzymes. J Pharmacol Exp Ther 300:355–360.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Humphries H,
    2. Neuhoff S, and
    3. Gardner I
    (2015) Application of PBPK modelling for prediction of FMO metabolism using benzydamine as a probe for FMO3, in Proceedings of the 44th DMDG Open Meeting, Cambridge UK September 16th to 18th 2015.
  7. ↵
    1. Jamei M,
    2. Marciniak S,
    3. Edwards D,
    4. Wragg K,
    5. Feng K,
    6. Barnett A, and
    7. Rostami-Hodjegan A
    (2013) The Simcyp population based simulator: architecture, implementation, and quality assurance. In Silico Pharmacol 1:9.
  8. ↵
    1. Jones BC,
    2. Srivastava A,
    3. Colclough N,
    4. Wilson J,
    5. Reddy VP,
    6. Amberntsson S, and
    7. Li D
    (2017) An investigation into the prediction of in vivo clearance for a range of flavin-containing monooxygenase substrates. Drug Metab Dispos 45:1060–1067.
    OpenUrlAbstract/FREE Full Text
    1. Katagiri F,
    2. Shiga T,
    3. Inoue S,
    4. Sato Y,
    5. Itoh H, and
    6. Takeyama M
    (2006) Effects of itopride hydrochloride on plasma gut-regulatory peptide and stress-related hormone levels in healthy human subjects. Pharmacology 77:115–121.
    OpenUrlPubMed
  9. ↵
    1. Koukouritaki SB and
    2. Hines RN
    (2005) Flavin-containing monooxygenase genetic polymorphism: impact on chemical metabolism and drug development. Pharmacogenomics 6:807–822.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Koukouritaki SB,
    2. Simpson P,
    3. Yeung CK,
    4. Rettie AE, and
    5. Hines RN
    (2002) Human hepatic flavin-containing monooxygenases 1 (FMO1) and 3 (FMO3) developmental expression. Pediatr Res 51:236–243.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Lin J
    (2015) In-vitro-in-vivo extraploation of UGT1A9-mediated clerance and application in a PBPK model for dapagliflozin, in Proceedings of the 44th DMDG Open Meeting, Cambridge UK September 16th to 18th 2015.
  12. ↵
    1. Luzon E,
    2. Blake K,
    3. Cole S,
    4. Nordmark A,
    5. Versantvoort C, and
    6. Berglund EG
    (2017) Physiologically based pharmacokinetic modeling in regulatory decision-making at the European Medicines Agency. Clin Pharmacol Ther 102:98–105.
    OpenUrl
    1. McNeil JJ,
    2. Mihaly GW,
    3. Anderson A,
    4. Marshall AW,
    5. Smallwood RA, and
    6. Louis WJ
    (1981) Pharmacokinetics of the H2-receptor antagonist ranitidine in man. Br J Clin Pharmacol 12:411–415.
    OpenUrlCrossRefPubMed
    1. Nguyen HQ,
    2. Callegari E, and
    3. Obach RS
    (2016) The use of in vitro data and physiologically-based pharmacokinetic modeling to predict drug metabolite exposure: desipramine exposure in cytochrome P4502D6 extensive and poor metabolizers following administration of imipramine. Drug Metab Dispos 44:1569–1578.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Obach RS,
    2. Baxter JG,
    3. Liston TE,
    4. Silber BM,
    5. Jones BC,
    6. MacIntyre F,
    7. Rance DJ, and
    8. Wastall P
    (1997) The prediction of human pharmacokinetic parameters from preclinical and in vitro metabolism data. J Pharmacol Exp Ther 283:46–58.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Overby LH,
    2. Carver GC, and
    3. Philpot RM
    (1997) Quantitation and kinetic properties of hepatic microsomal and recombinant flavin-containing monooxygenases 3 and 5 from humans. Chem Biol Interact 106:29–45.
    OpenUrlCrossRefPubMed
    1. Raaflaub J,
    2. Haefelfinger P, and
    3. Trautmann KH
    (1984) Single-dose pharmacokinetics of the MAO-inhibitor moclobemide in man. Arzneimittelforschung 34:80–82.
    OpenUrlPubMed
  15. ↵
    1. Reddy RR,
    2. Ralph EC,
    3. Motika MS,
    4. Zhang J, and
    5. Cashman JR
    (2010) Characterization of human flavin-containing monooxygenase (FMO) 3 and FMO5 expressed as maltose-binding protein fusions. Drug Metab Dispos 38:2239–2245.
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Rodgers T,
    2. Leahy D, and
    3. Rowland M
    (2005) Physiologically based pharmacokinetic modeling 1: predicting the tissue distribution of moderate-to-strong bases. J Pharm Sci 94:1259–1276.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Rodgers T and
    2. Rowland M
    (2006) Physiologically based pharmacokinetic modelling 2: predicting the tissue distribution of acids, very weak bases, neutrals and zwitterions. J Pharm Sci 95:1238–1257.
    OpenUrlCrossRefPubMed
    1. Schoerlin MP,
    2. Mayersohn M,
    3. Korn A, and
    4. Eggers H
    (1987) Disposition kinetics of moclobemide, a monoamine oxidase-A enzyme inhibitor: single and multiple dosing in normal subjects. Clin Pharmacol Ther 42:395–404.
    OpenUrlPubMed
  18. ↵
    1. Sjögren E,
    2. Westergren J,
    3. Grant I,
    4. Hanisch G,
    5. Lindfors L,
    6. Lennernäs H,
    7. Abrahamsson B, and
    8. Tannergren C
    (2013) In silico predictions of gastrointestinal drug absorption in pharmaceutical product development: application of the mechanistic absorption model GI-Sim. Eur J Pharm Sci 49:679–698.
    OpenUrlCrossRefPubMed
    1. Takano A,
    2. Suhara T,
    3. Kusumi I,
    4. Takahashi Y,
    5. Asai Y,
    6. Yasuno F,
    7. Ichimiya T,
    8. Inoue M,
    9. Sudo Y, and
    10. Koyama T
    (2006) Time course of dopamine D2 receptor occupancy by clozapine with medium and high plasma concentrations. Prog Neuropsychopharmacol Biol Psychiatry 30:75–81.
    OpenUrlPubMed
  19. ↵
    1. Taniguchi-Takizawa T,
    2. Shimizu M,
    3. Kume T, and
    4. Yamazaki H
    (2015) Benzydamine N-oxygenation as an index for flavin-containing monooxygenase activity and benzydamine N-demethylation by cytochrome P450 enzymes in liver microsomes from rats, dogs, monkeys, and humans. Drug Metab Pharmacokinet 30:64–69.
    OpenUrlCrossRefPubMed
    1. Tassaneeyakul W,
    2. Kittiwattanagul K,
    3. Vannaprasaht S,
    4. Kampan J,
    5. Tawalee A,
    6. Puapairoj P,
    7. Tiamkao S,
    8. Juthagridsada S,
    9. Kukongviriyapan V, and
    10. Tassaneeyakul W
    (2005) Steady-state bioequivalence study of clozapine tablet in schizophrenic patients. J Pharm Pharm Sci 8:47–53.
    OpenUrlPubMed
    1. Traynor AM,
    2. Hewitt M,
    3. Liu G,
    4. Flaherty KT,
    5. Clark J,
    6. Freedman SJ,
    7. Scott BB,
    8. Leighton AM,
    9. Watson PA,
    10. Zhao B, et al.
    (2011) Phase I dose escalation study of MK-0457, a novel Aurora kinase inhibitor, in adult patients with advanced solid tumors. Cancer Chemother Pharmacol 67:305–314.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Wagner C,
    2. Zhao P,
    3. Pan Y,
    4. Hsu V,
    5. Grillo J,
    6. Huang SM, and
    7. Sinha V
    (2015) Application of physiologically based pharmacokinetic (PBPK) modeling to support dose selection: report of an FDA public workshop on PBPK. CPT Pharmacometrics Syst Pharmacol 4:226–230.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Yeung CK,
    2. Lang DH,
    3. Thummel KE, and
    4. Rettie AE
    (2000) Immunoquantitation of FMO1 in human liver, kidney, and intestine. Drug Metab Dispos 28:1107–1111.
    OpenUrlAbstract/FREE Full Text
    1. Yoon S,
    2. Lee H,
    3. Kim TE,
    4. Lee S,
    5. Chee DH,
    6. Cho JY,
    7. Yu KS, and
    8. Jang IJ
    (2014) Comparative steady-state pharmacokinetic study of an extended-release formulation of itopride and its immediate-release reference formulation in healthy volunteers. Drug Des Devel Ther 8:123–128.
    OpenUrl
  22. ↵
    1. Zhou W,
    2. Humphries H,
    3. Neuhoff S,
    4. Gardner I,
    5. Masson E,
    6. Al-Huniti N, and
    7. Zhou D
    (2017) Development of a physiologically based pharmacokinetic model to predict the effects of flavin-containing monooxygenase 3 (FMO3) polymorphisms on itopride exposure. Biopharm Drug Dispos 38:389–393.
    OpenUrl
PreviousNext
Back to top

In this issue

Drug Metabolism and Disposition: 46 (9)
Drug Metabolism and Disposition
Vol. 46, Issue 9
1 Sep 2018
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Drug Metabolism & Disposition article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
An Investigation into the Prediction of the Plasma Concentration-Time Profile and Its Interindividual Variability for a Range of Flavin-Containing Monooxygenase Substrates Using a Physiologically Based Pharmacokinetic Modeling Approach
(Your Name) has forwarded a page to you from Drug Metabolism & Disposition
(Your Name) thought you would be interested in this article in Drug Metabolism & Disposition.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Research ArticleArticle

PBPK Modeling of FMO Substrates

Venkatesh Pilla Reddy, Barry C. Jones, Nicola Colclough, Abhishek Srivastava, Joanne Wilson and Danxi Li
Drug Metabolism and Disposition September 1, 2018, 46 (9) 1259-1267; DOI: https://doi.org/10.1124/dmd.118.080648

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Share
Research ArticleArticle

PBPK Modeling of FMO Substrates

Venkatesh Pilla Reddy, Barry C. Jones, Nicola Colclough, Abhishek Srivastava, Joanne Wilson and Danxi Li
Drug Metabolism and Disposition September 1, 2018, 46 (9) 1259-1267; DOI: https://doi.org/10.1124/dmd.118.080648
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Sex- and lifestyle-related factors affect hepatic CYP levels
  • Adipocyte PXR does not play an essential role in obesity.
  • CYP3A-mediated oxidation of DABE and BIBR0951
Show more Articles

Similar Articles

Advertisement
  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About DMD
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Journal of Pharmacology and Experimental Therapeutics
  • Molecular Pharmacology
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-009X (Online)

Copyright © 2023 by the American Society for Pharmacology and Experimental Therapeutics