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
  • 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
  • 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

Functional Characterization of 40 CYP3A4 Variants by Assessing Midazolam 1′-Hydroxylation and Testosterone 6β-Hydroxylation

Masaki Kumondai, Evelyn Marie Gutiérrez Rico, Eiji Hishinuma, Akiko Ueda, Sakae Saito, Daisuke Saigusa, Shu Tadaka, Kengo Kinoshita, Tomoki Nakayoshi, Akifumi Oda, Ai Abe, Masamitsu Maekawa, Nariyasu Mano, Noriyasu Hirasawa and Masahiro Hiratsuka
Drug Metabolism and Disposition March 2021, 49 (3) 212-220; DOI: https://doi.org/10.1124/dmd.120.000261
Masaki Kumondai
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Evelyn Marie Gutiérrez Rico
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eiji Hishinuma
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Akiko Ueda
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sakae Saito
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daisuke Saigusa
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shu Tadaka
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kengo Kinoshita
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tomoki Nakayoshi
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Akifumi Oda
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ai Abe
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Masamitsu Maekawa
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nariyasu Mano
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Noriyasu Hirasawa
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Masahiro Hiratsuka
Laboratory of Pharmacotherapy of Life-Style Related Diseases, Graduate School of Pharmaceutical Sciences (M.K., E.M.G.R., N.H., M.H.), Tohoku Medical Megabank Organization (E.H., S.S., D.S., S.T., K.K., M.H.), Advanced Research Center for Innovations in Next-Generation Medicine (E.H., A.U., N.H., M.H.), and Laboratory of Clinical Pharmacy, Faculty of Pharmaceutical Sciences (A.A., M.M., N.M.), Tohoku University, Sendai, Japan; Faculty of Pharmacy, Meijo University, Nagoya, Japan (T.N., A.O.); and Department of Pharmaceutical Sciences, Tohoku University Hospital, Sendai, Japan (M.M., N.M., N.H., M.H.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Masahiro Hiratsuka
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF
Loading

Visual Overview

Figure
  • Download figure
  • Open in new tab
  • Download powerpoint

Abstract

CYP3A4 is among the most abundant liver and intestinal drug-metabolizing cytochrome P450 enzymes, contributing to the metabolism of more than 30% of clinically used drugs. Therefore, interindividual variability in CYP3A4 activity is a frequent cause of reduced drug efficacy and adverse effects. In this study, we characterized wild-type CYP3A4 and 40 CYP3A4 variants, including 11 new variants, detected among 4773 Japanese individuals by assessing CYP3A4 enzymatic activities for two representative substrates (midazolam and testosterone). The reduced carbon monoxide–difference spectra of wild-type CYP3A4 and 31 CYP3A4 variants produced with our established mammalian cell expression system were determined by measuring the increase in maximum absorption at 450 nm after carbon monoxide treatment. The kinetic parameters of midazolam and testosterone hydroxylation by wild-type CYP3A4 and 29 CYP3A4 variants (Km, kcat, and catalytic efficiency) were determined, and the causes of their kinetic differences were evaluated by three-dimensional structural modeling. Our findings offer insight into the mechanism underlying interindividual differences in CYP3A4-dependent drug metabolism. Moreover, our results provide guidance for improving drug administration protocols by considering the information on CYP3A4 genetic polymorphisms.

SIGNIFICANCE STATEMENT CYP3A4 metabolizes more than 30% of clinically used drugs. Interindividual differences in drug efficacy and adverse-effect rates have been linked to ethnicity-specific differences in CYP3A4 gene variants in Asian populations, including Japanese individuals, indicating the presence of CYP3A4 polymorphisms resulting in the increased expression of loss-of-function variants. This study detected alterations in CYP3A4 activity due to amino acid substitutions by assessing the enzymatic activities of coding variants for two representative CYP3A4 substrates.

Introduction

CYP3A4 is the most abundant liver and intestinal drug-metabolizing cytochrome P450 enzyme, contributing to the metabolism of more than 30% of medications (Zanger and Schwab, 2013; Jackson et al., 2018; Niwa et al., 2020). Interindividual variability in CYP3A4 activity and expression influences drug responsiveness, efficacy, and adverse-effect rates (Lamba et al., 2002a). Besides the various factors that contribute to this variability, including CYP3A4 inducers and inhibitors and food-drug and drug-drug interactions, CYP3A4 genetic variations also affect drug responses (Racha et al., 2003; Huang et al., 2008; Elens et al., 2013; Hendriks et al., 2020). CYP3A5 genetic polymorphisms, including CYP3A5*3, have been widely investigated because of their high allele frequencies and clinical impact on the variability of CYP3A activity. CYP3A4 gene variations may be overlooked because of the higher impact of some CYP3A5 variants on activity. However, CYP3A4 gene variations should also be investigated to elucidate their influence on the combinational effects of CYP3A4 and CYP3A5 enzymatic activity and CYP3A in vivo activity (Fukushima-Uesaka et al., 2004; Lee et al., 2005; Zhou et al., 2017). For instance, several CYP3A4 variant alleles, including CYP3A4*1B (a promoter variant) and *22 (an intronic variant), reduce mRNA levels or cause splicing defects, although CYP3A4*6, *17, *20, and *26 are associated with diminished enzymatic activities by loss of function (Lamba et al., 2002a; Elens et al., 2013; Werk et al., 2014). The identification of these variants can improve patient safety and treatment efficacy; thus, increased dose-adjusted tacrolimus concentration levels are achieved in transplant patients carrying a loss-of-function genotype (CYP3A4*20) (Gómez-Bravo et al., 2018). Hence, CYP3A4 genotyping improves predictive outcomes.

The Human Cytochrome P450 Allele Nomenclature Database [https://www.pharmvar.org/gene/CYP3A4 (accessed September 21st, 2020)] contains records of CYP3A4 genetic polymorphisms (CYP3A4*2–*24, CYP3A4*26, and CYP3A4*28–*34). Moreover, several whole-genome sequencing (WGS) studies have identified CYP3A4 genetic variations, including c.1165C>T (Pro389Ser) and c.1423C>G (Leu475Val) (Apellaniz-Ruiz et al., 2015). A WGS project of the Tohoku Medical Megabank Organization, which targeted 4773 Japanese individuals, identified numerous genetic variations, including single nucleotide variations (SNVs), as shown in Table 1 [https://jmorp.megabank.tohoku.ac.jp (accessed August 30th, 2020)]. The identification of rare genetic variants, which may cause interindividual differences in drug efficacy and adverse reactions, facilitates the development of personalized medicine (Fujikura et al., 2015; Schärfe et al., 2017). CYP3A4 rare genetic variants have been implicated in high-dose simvastatin–related myopathy and affect serum voriconazole concentrations linked to adverse drug reactions (Ingelman-Sundberg et al., 2018). Although the CYP3A4 SNVs analyzed in this study are relatively rare (Table 1), their effects on CYP3A4 activity should be determined to improve treatment outcomes and manage adverse effects.

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

CYP3A4 allelic variants characterized in this study

Heterologous expression systems can be employed to evaluate the effects of CYP3A4 polymorphisms in vitro and complement in vivo analysis data. Performing in vivo assays is challenging because of low allele frequencies and increased patient burden. However, all CYP3A4 variant proteins should be evaluated using the same protocol and experimental conditions because methodological differences between studies may lead to inconsistent kinetic parameters for some variants (Hiratsuka, 2012). To date, characterization studies have used cDNA-based heterologous expression systems, including Escherichia coli, baculoviruses, and mammalian cell systems (Eiselt et al., 2001; Lee et al., 2005; Maekawa et al., 2010; Kumondai et al., 2018; Saito et al., 2018; Watanabe et al., 2018; Gutiérrez Rico et al., 2020). However, accurate assessments are limited by interspecies differences of post-translational modifications. Thus, mammalian cell expression systems should be preferred to evaluate the function of human CYP3A4 variants in drug metabolism (Clark and Pazdernik, 2016). Recently, an intronic CYP2D6 mutation was evaluated using a minigene expression system with transfected whole CYP2D6 sequences (<5 kb) (Zanger et al., 2020). Minigenes may aid the assessment of CYP3A4 levels. However, genomic DNA expression using minigenes, including the whole CYP3A4 gene (>20 kb), is challenging because transfection efficiency decreases with increased DNA length (Hornstein et al., 2016). Therefore, we used our heterologous expression system consisting of 293FT cells coexpressing CYP3A4, cytochrome P450 oxidoreductase (CPR), and cytochrome b5 (Kumondai et al., 2020).

Here, we functionally characterized wild-type CYP3A4 (CYP3A4.1) and 40 CYP3A4 allelic variants, including 11 novel variants identified in 4773 Japanese individuals, by using two representative CYP3A4 substrates (midazolam and testosterone). Notably, clinical and in vitro studies have assessed the functional differences and drug-drug interactions related to CYP3A4, including cytochrome P450 variants with substrate-specific enzymatic activities (Foti et al., 2010; Hiratsuka, 2016; Xiao et al., 2019). Moreover, carbon monoxide (CO)-difference spectroscopy and three-dimensional (3D) structure analysis were performed to gain insights into functional alterations caused by CYP3A4 variants.

Materials and Methods

Chemicals.

The following reagents were purchased from commercial sources: midazolam [PubChem Compound ID (CID) 4192] and flunitrazepam (Wako Pure Chemical Industries, Osaka, Japan), 1′-hydroxymidazolam (CID 107917) and 6β-hydroxytestosterone-d7 (Corning Incorporated, NY), testosterone (CID 6013) and 6β-testosterone (CID 65543) (Sigma-Aldrich, Steinheim, Germany), NADH and NADPH (Oriental Yeast, Tokyo, Japan), polyclonal anti-human CYP3A4 antibody (ab3572) (Abcam, Cambridge, UK), horseradish peroxidase–conjugated goat anti-rabbit IgG (ProteinSimple, Tokyo, Japan), dimanganese decacarbonyl (Sigma-Aldrich), and sodium cyanide and cytochrome c from horse heart (Nacalai Tesque, Kyoto, Japan). All other chemicals and reagents were of the highest quality commercially available.

Detection of CYP3A4 Sequence Alterations by Sanger Sequencing Analysis.

PCR amplification was conducted using peripheral blood leukocyte genome DNA isolated from whole blood using the Gentra Puregene Blood Kit (Qiagen, Hilden, Germany), as described previously (Nagasaki et al., 2015). Whole blood samples were obtained from Japanese subjects participating in the community-based cohort study conducted by the Tohoku Medical Megabank Organization, and written informed consent was obtained from all subjects before sample collection. Primers were used to amplify sequences containing each CYP3A4 exon (Supplemental Table 1). PCR amplification was performed with genomic DNA samples (10 ng), 2× AmpliTaq Gold 360 Master Mix (Applied Biosystems, Foster City, CA), and 0.5 µM of each primer in a total volume of 20 µl. Thermal cycling conditions included an initial denaturation step at 95°C for 10 minutes, followed by 30 cycles of denaturation at 95°C for 30 seconds, annealing at 60°C for 30 seconds, extension at 72°C for 30 seconds (for exons 2–12) or 1 minute (for exons 1 and 13), and a final extension at 72°C for 7 minutes. The PCR products were purified using columns and analyzed by Sanger sequencing utilizing the same primers for each exon as in the PCR except for the forward sequencing reaction of exon 1 (5′-CAG​GCG​TGG​AAA​CAC​AAT​GGT​GG-3′) and reverse sequencing reaction of exon 11 (5′-TGT​GGA​TGA​CTG​TAG​TTT​TC-3′). To determine the linkage between two SNVs identified in one subject, the respective PCR fragments were ligated into the pcDNA3.4 vector (ThermoFisher Scientific), and the ligation products were transformed into E. coli. Single colonies of the transformants were collected, and the plasmids were isolated and analyzed by Sanger sequencing. This analysis was performed in accordance with the Declaration of Helsinki, and prior written informed consent was obtained from all participants.

Construction of CYP3A4 cDNA Expression Vectors.

The wild-type CYP3A4 (CYP3A4*1) plasmid used in this study has been described in a previous report (Sakurai et al., 2005). A plasmid containing CYP3A4*1 cDNA was used as a template to generate CYP3A4 allelic variant constructs (CYP3A4*2–CYP3A4*5, CYP3A4*7–CYP3A4*21, CYP3A4*23, CYP3A4*24, Pro389Ser, Leu475Val, CYP3A4*28, CYP3A4*29, CYP3A4*31–CYP3A4*34, Thr138Ala, Met256Ile, Ile300Val, Thr323Ala, Met353Leu, Ile369Val, Ile369Asn, Ala370Val, Arg372Thr, Tyr399Ser, and Ala448Ser) using primer sets described in Supplemental Table 2 with the QuikChange Site-Directed Mutagenesis Kit (Stratagene, La Jolla, CA) following the manufacturer’s instructions. All prepared constructs were confirmed by Sanger sequencing. Wild-type and variant CYP3A4 cDNAs were subcloned into the mammalian expression vector pcDNA3.4. A mock plasmid, including the pcDNA3.4 vector, was prepared as previously reported (Kumondai et al., 2018; Saito et al., 2018).

Expression of CYP3A4 Variants in 293FT Cells.

The 293FT cells (ThermoFisher Scientific) were cultured in Dulbecco’s modified Eagle’s medium (Nacalai Tesque) containing 10% fetal bovine serum at 37°C under 5% CO2. Cells were plated at a density of 2.2 × 106 cells per 100-mm dish; 24 hours after plating, the cells were transfected with a plasmid carrying CYP3A4 (9.6 µg), CPR (0.2 µg), and cytochrome b5 (0.2 µg) cDNA using 30 μl of 1.0 mg/ml of Polyethylenimine Max (Polysciences, Inc., Warrington, PA), according to previously described methods (Kumondai et al., 2020). After incubation for 12 hours at 37°C, 0.25 mM of 5-aminolevulinic acid hydrochloride (Nacalai Tesque) and 0.25 mM iron (II) sulfate heptahydrate (Wako) were added to the medium. After incubation for 48 hours post-transfection at 37°C, the cells were scraped from the plates, and the microsomal fractions were prepared as previously described (Kumondai et al., 2020). Protein concentrations were determined using the BCA Protein Assay Kit (ThermoFisher Scientific).

Western Blotting.

CYP3A4 levels were quantified with an immunoassay using the Wes system (ProteinSimple). Each well was loaded with 25 μg/ml microsomes. CYP3A4 was detected using a polyclonal anti-human CYP3A4 antibody (diluted 1:100) and horseradish peroxidase–conjugated goat anti-rabbit IgG. A total protein assay was performed to normalize each signal using 25 μg/ml microsomes following the manufacturer’s instructions. Recombinant human CYP3A4 Supersomes (Corning Incorporated) were used as the standard to quantify the CYP3A4 protein level.

CO-Difference Spectroscopy.

The cytochrome P450 holoprotein content was spectrophotometrically measured according to a previously reported method based on the established protocol by Omura and Sato (Kumondai et al., 2020).

Measurement of NADPH-Cytochrome c Reduction Activity and Microsomal Cytochrome b5 Content Determination.

The activity of CPR and the content of cytochrome b5 were spectrophotometrically evaluated as previously described (Kumondai et al., 2020).

Midazolam 1′-Hydroxylation.

CYP3A4 midazolam 1′-hydroxylation activity was determined using a previously described method with several modifications (El Mahjoub and Staub, 2000). The reaction mixture, a total volume of 100 μl, contained the following components: the microsomal fraction (50 μg), midazolam (1, 2, 4, 6, 8, 10, and 15 μM), and 100 mM potassium phosphate buffer (pH 7.4). After preincubation at 37°C for 3 minutes, reactions were initiated by adding 10 mM NADPH, and incubation continued at 37°C for 10 minutes. Reactions were terminated by adding 100 μl of acetonitrile containing 500 nM flunitrazepam as an internal standard. After protein removal by centrifugation at 15,400g for 10 minutes, the supernatant was injected into a liquid chromatography-tandem mass spectrometry (LC-MS/MS) system, as previously described (Kumondai et al., 2020).

Testosterone 6β-Hydroxylation Assay.

CYP3A4 testosterone 6β-hydroxylation activity was determined using a previously described method with several modifications (Nguyen et al., 2020). The reaction mixture, in a total volume of 100 μl, contained the following components: the microsomal fraction (50 μg), testosterone (5, 10, 20, 50, 100, 200, or 500 μM), and 100 mM potassium phosphate buffer (pH 7.4). After preincubation at 37°C for 3 minutes, reactions were initiated by adding 10 mM NADPH, and the incubation was continued at 37°C for 10 minutes. The reactions were terminated by adding 100 μl of acetonitrile containing 10 μM 6β-hydroxytestosterone-d7 as an internal standard. After protein removal by centrifugation at 15,400g for 10 minutes, the supernatant was injected into an LC-MS/MS system.

The 6β-hydroxytestosterone content was measured using the LC-MS/MS system in the positive ion detection mode at the electrospray ionization interface (QTRAP 6500 LC-MS/MS system; AB Sciex, MA). Separation by ultra-high-performance liquid chromatography was conducted using the Nexera Ultra-High-Performance Liquid Chromatography system (Shimadzu, Kyoto, Japan). Chromatographic separation was performed with a Kinetex C18 column (2.1 × 50 mm, 5-µg particle size; Phenomenex, Shimadzu) maintained at 40°C. Mobile phases were prepared with deionized water containing 0.1% formic acid as eluent A and acetonitrile containing 0.1% formic acid as eluent B. The flow rate was 300 µl/min. The gradient program was as follows: elution was initiated with 10% B, followed by a linear gradient to 60% B from 1 to 7 minutes, held at 100% B for 2 minutes, and then immediately returned to the initial conditions, which were maintained for 3 minutes until the end of the run. Quantification analyses using Analyst Software (Sciex) were performed in the selected reaction monitoring mode in which ion transitions from the precursor into a product ion were monitored: mass-to-charge ratio (m/z) 289.1→97.0 for testosterone (collision energy, 14 V), m/z 305.0→269.1 for 6β-hydroxytestosterone (collision energy, 12 V), and m/z 312.1→276.1 for 6β-hydroxytestosterone-d6 (collision energy, 4 V). The optimized mass spectrometry settings were as follows: entrance potential, 2.0 V; curtain gas, 25.0 psi; ion transfer voltage, 5500.0 V; temperature, 750.0°C; gas 1, 70.0 psi; gas 2, 30.0 psi; and collision gas, 12.0. Standard curves for 6β-hydroxytestosterone were constructed within the range of 100–30,000 nM using metabolite standards, with a quantification limit of 100 nM.

Data Analysis.

Kinetic data, including the Michaelis constant (Km), turnover number (kcat), and catalytic efficiency (kcat/Km) values, were determined using the Enzyme Kinetics Module of SigmaPlot 12.5 (Systat Software, Inc., Chicago, IL), a curve-fitting program based on nonlinear regression analysis. The normality of our data sets was initially assessed using the Shapiro-Wilk test. Statistical analyses for multiple comparisons were performed through variance analysis by Dunnett’s T3 test or the Kruskal-Wallis method (IBM SPSS Statistics version 22; International Business Machines, Armonk, NY). All values are expressed as the means ± S.D. of experiments performed in triplicate. Differences with P < 0.05 were considered statistically significant. All assays and samples were prepared in triplicate to permit statistical analysis.

3D Structural Modeling of CYP3A4.

The 3D structural modeling of CYP3A4 was based on the CYP3A4 X-ray structure of Sevrioukova and Poulos (2017) (Protein Data Bank code 5te8). Swiss-model was used to complement a part of the lacking structure in CYP3A4. After removing the substrate, midazolam or testosterone was docked with the CYP3A4 X-ray structure according to the CDOCKER protocol of Discovery Studio 2.5 (BIOVIA, CA). Docking iterations were conducted, taking into consideration the binding orientation and binding energies under the condition that the volume of the space was defined as 9, and the heme iron charge status was Fe3+. After replacing each substitution, structural optimization was conducted as previously described (Oda et al., 2005).

Results

We performed resequencing for the subjects carrying novel SNVs, except for a subject carrying c.1342G>T because this DNA sample was not available for further sequencing. For SNVs located within exons, WGS identified 10 novel variants; these findings were identical with the results obtained using Sanger sequencing. Two SNVs (c.389G>A and c.412A>G) carried by one subject were located in different alleles (Supplemental Fig. 1). Thus, 11 novel CYP3A4 variants, including the c.1342C>T variant, were identified among 4773 Japanese individuals with low allele frequencies (0.01%–0.05%).

CYP3A4 protein levels normalized by total protein concentrations were assessed by Western blotting with a polyclonal CYP3A4 antibody that recognized all CYP3A4 variants (Supplemental Figs. 2, A and B). The protein levels of CYP3A4.1 and the CYP3A4 variants did not differ significantly (Supplemental Fig. 2C). Moreover, the aspect ratio of CPR activity and cytochrome b5 content did not significantly vary among the analyzed CYP3A4 variants (Supplemental Fig. 3). Thus, we used single microsomes expressed per CYP3A4 variant for functional characterization.

In our study, CYP3A4 exhibited high activity with no significant solvent effects. The holoprotein content of CYP3A4 variants expressed in 293FT cells was assessed by CO-difference spectroscopy (Fig. 1). The reduced CO-difference spectra of microsomal fractions of CYP3A4.1 and 31 CYP3A4 variants (CYP3A4.2–.5, .7, .9, .10, .14–.16, .18, .19, .23, .24, Pro389Ser, Leu475Val, .28, .29, .31–.34, Thr138Ala, Met256Ile, Ile300Val, Thr323Ala, Met353Leu, Ile369Val, Ile369Asn, Ala370Val, and Ala448Ser) were determined by measuring the increase in maximum absorption at 450 nm after CO treatment. The remaining nine variants (CYP3A4.8, .11–.13, .17, 488Frameshift, .21, Arg372Thr, and Tyr399Ser) had no detectable holoprotein content based on the lack of a significant increase in their maximum absorption at 450 nm.

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

Representative CO-difference spectra of CYP3A4 variant proteins expressed in 293FT cells. All assays and measurements were performed in triplicate using a single microsomal preparation. Mock 1 indicates transfection with 10 μg mock plasmids. Mock 2 indicates transfection with 9.6 μg mock plasmids, 0.2 μg CPR plasmids, and 0.2 μg cytochrome b5 plasmids.

We selected midazolam and testosterone as substrates because midazolam 1′-hydroxylation and testosterone 6β-hydroxylation have long been recognized as CYP3A member-mediated reactions. First, we confirmed that the formation of 1′-hydroxymidazolam or 6β-hydroxytestosterone was linear for up to 10 minutes of incubation by using 50 μg microsomal proteins and 15 µM midazolam or 500 µM testosterone, respectively. Both metabolites were linear in the presence of 0–50 μg microsomal proteins during 10 minutes of incubation. The kinetic parameters were calculated by Michaelis-Menten curves of midazolam 1′-hydroxylation and testosterone 6β-hydroxylation for wild-type CYP3A4 and CYP3A4 variants (Fig. 2; Table 2). However, because of the insufficient metabolite production, we could not determine the kinetic parameters of midazolam 1′-hydroxylation for 11 variants (CYP3A4.8, .11–.13, .17, 488Frameshift, .21, Ile369Asn, Ala370Val, Arg372Thr, and Tyr399Ser). We found that 15 variants had significantly lower catalytic efficiency values (P < 0.05) than CYP3A4.1, whereas there were no significant differences in Km and kcat values among these CYP3A4 variants. Conversely, the Km values of testosterone 6β-hydroxylation were significantly higher for CYP3A4.7, .16, .33, and Thr138Ala than for CYP3A4.1 (P < 0.05). Moreover, the kcat and catalytic efficiency values were significantly lower for CYP3A4.12 and Leu475Val than for CYP3A4.1 (P < 0.05). Insufficient metabolite production prevented us from determining the kinetic values for nine variants (CYP3A4.8, .11, .13, .17, 488Frameshift, .21, Ile369Asn, Ala370Val, and Arg372Thr). Lastly, we calculated the correlation between the midazolam 1′-hydroxylation and the testosterone 6β-hydroxylation catalytic efficiency values. As shown in Fig. 3, there was a significant correlation between these values among the different variants (R2 = 0.724, P < 0.001).

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

Michaelis-Menten curves for CYP3A4 variants. Determined kinetic parameters (Km, kcat, and catalytic efficiency) of midazolam 1′-hydroxylation and testosterone 6β-hydroxylation. All assays and measurements were performed in triplicate using a single microsomal preparation.

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

Kinetic parameters of midazolam 1′-hydroxylation and testosterone 6β-hydroxylation by microsomes from 293FT cells expressing wild-type CYP3A4 and variant CYP3A4 proteins

The data represent the means ± S.D. of the three independently performed catalytic assays. All assays and measurements were performed in triplicate using a single microsomal preparation. Mock 1 indicates transfection with 10 μg mock plasmids. Mock 2 indicates transfection with 9.6 μg mock plasmids, 0.2 μg CPR plasmids, and 0.2 μg cytochrome b5 plasmids.

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

Correlation between the catalytic efficiency ratios for midazolam 1′-hydroxylation and testosterone 6β-hydroxylation among CYP3A4 variants. Midazolam 1′-hydroxylation catalytic efficiency ratios are plotted on the horizontal axis, and testosterone 6β-hydroxylation catalytic efficiency ratios are plotted on the vertical axis. The numbers correspond to known CYP3A4 variants, and the amino acid substitutions indicate the newly identified variants.

We used 3D structural modeling to assess the molecular properties of CYP3A4 variants and their interactions with midazolam. As shown in Fig. 4, A and B, the residue substitution in CYP3A4.8 and CYP3A4.12 abolished the interaction with heme. The midazolam docking analysis indicated that the Thr138Ala substitution variant caused a hydrophobic interaction between Arg130 and Ser134 via the adjacently formed hydrophobic interaction between Ala138 and Ser134 (Fig. 4C). In CYP3A4.17, the substitution Ser189 located at the F-G loop lacked a hydrophobic interaction with Leu249 (Fig. 4D). CYP3A4.21 harbored its substitution in the I-helix, which formed a hydrophobic interaction between Cys319 and Val489 (Fig. 4E). Additionally, the substitution Ala448Ser abolished numerous hydrophobic interactions with heme (Fig. 4F). The substitution Tyr399Ser caused the formation of a hydrogen bond with Met371 (Fig. 4G). We further analyzed Ile369 and Val370, as shown in Fig. 4, H–M. The Ile369Val substitution had an additional hydrophobic interaction with midazolam and heme, affecting the distance between heme and the active site (Fig. 4H). Furthermore, differences in the midazolam and heme interactions due to the amino acid substitutions (Fig. 4, I, K, and L), as well as substrate-specific differences in the heme-surrounding interactions (Fig. 4, H, J, K, and M), were observed.

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

Diagram pairs showing the partial crystal structure of CYP3A4.1 (left image) and CYP3A4 variants (right image) for CYP3A4.8 (A), CYP3A4.12 (B), Thr138Ala (C), CYP3A4.17 (D), CYP3A4.21 (E), Ala448Ser (F), Tyr399Ser (G), Ile369Val coordinated with midazolam (H), Ile369Asn (I), Ile369Val coordinated with testosterone (J), CYP3A4.33 coordinated with midazolam (K), Ala370Val (L), and CYP3A4.33 coordinated with testosterone (M). Pink line, hydrophobic interactions. Green line, conventional hydrogen bonds. Gray line, carbon-hydrogen bonds. MDZ, midazolam; TST, testosterone.

Discussion

Several CYP3A4 genetic variations may contribute to interindividual differences in drug efficacy and adverse-effect rates (Lamba et al., 2002a). Several studies have identified CYP3A4 variants and characterized the associated in vitro enzyme kinetics, including previous studies of CYP3A4 variants among the Han Chinese population using numerous substrates (Supplemental Table 3) (Fang et al., 2017; Xu et al., 2018; Li et al., 2019a,b; Lin et al., 2019a,b; Yang et al., 2019; Zhou et al., 2019a; Chen et al., 2020). However, the underlying causes for functional differences produced by the amino acid alterations are not fully understood and require further research. This study characterized 40 CYP3A4 variants, including 11 novel variants identified among 4773 Japanese individuals, by assessing CO-difference spectra, enzymatic activities, and 3D structural modeling analysis.

A single individual had two allelic variants (CYP3A4.8 and Thr138Ala, Supplemental Fig. 1), each displaying significantly decreased CYP3A4 activity for midazolam and testosterone. These substitutions were located near the vital CPR-cytochrome P450 electron transfer site, where they could affect the enzymatic activities (Hasemann et al., 1995; Bridges et al., 1998). Accordingly, an increase in the maximum absorption at 450 nm after CO treatment was observed for the Thr138Ala variant (Fig. 1) as a result of interactions involving Ala138 that caused hydrogen bonds between Arg130 and Ser134, which could affect electron transfer because Arg130 located in the C-D loop plays an essential role in electrostatic interactions (Fig. 4C). Interestingly, CYP3A4.8 did not have an absorption peak at 450 nm after CO treatment (Fig. 1), reflecting a failed heme incorporation and resulting in zero CYP3A4.8 holoprotein content that could be based on disrupted interactions with heme (Fig. 4A). Although it may be an infrequent case, carriers of two allelic variants that are predicted to be poor metabolizers may be at an increased risk of adverse drug reactions when treated with drugs metabolized by CYP3A4.

The midazolam and testosterone hydroxylation activities were slightly decreased or absent in 11 variants (CYP3A4.8, .11–.13, .17, 488Frameshift, .21, Ile369Asn, Ala370Val, Arg372Thr, and Tyr399Ser). Excluding Ile369Asn and Ala370Val, the nine remaining variants did not exhibit absorption peaks at 450 nm after CO treatment (Fig. 1), indicating failed heme incorporation. These findings suggest that the remaining nine CYP3A4 variants lacked quantifiable enzymatic activity. Similarly, the zero holoprotein content for variants CYP3A4.8 and CYP3A4.12 may be explained by disrupted interactions with heme (Fig. 4, A and B). The lack of activity exhibited by 488Frameshift was due to extensive changes in its molecular conformation caused by the frameshift mutation. The substitutions in CYP3A4.11 and .13, which were located in the conserved glutamine-arginine-arginine motif and K region critical for maintaining the structure among cytochrome P450 isoforms, could affect holoprotein stability, as previously reported (Eiselt et al., 2001). The substitution in CYP3A4.17 located in the E-helix was associated with slightly decreased enzymatic activity. Ser189 lacked a hydrophobic interaction with Leu249 (Fig. 4D), resulting in conformational changes in the active site cavity around the F-G loop, affecting the main active site channel as well as holoprotein stability (Benkaidali et al., 2019). CYP3A4.21 harbored a substitution in the I-helix, which may affect the stability and flexibility of the β-sheet by forming a hydrophobic interaction between Cys319 and Val489 (Fig. 4E). However, other substitutions, such as Ala448Ser, which apparently affected the hydrophobic interactions with heme, did not result in a complete lack of an absorption peak at 450 nm (Fig. 4G); however, the Ala448Ser substitution decreased the enzymatic activities. Thus, it is important to evaluate CO-difference spectroscopy results in combination with the findings of 3D structural modeling to improve our understanding of the consequential impact of amino acid substitutions.

In addition, a novel substrate recognition site (SRS) map for mammalian cytochrome P450 enzymes was developed based on the 3D crystal structures and docking of 868 substrates (Zawaira et al., 2011). Among 40 variants characterized in this study, the following variants were classified as listed regions based on the SRS map: CYP3A4.4 and .29, SRS-1; CYP3A4.2 and .5, SRS-2; CYP3A4.18 and Ile300Val, SRS-4; CYP3A4.12, .33, Ile369Val, Ile369Asn, Ala370Val, and Arg372Thr, SRS-5; and Leu475Val, SRS-6. For instance, the substitution Ile369Val appeared to directly affect the enzyme’s interactions with midazolam, resulting in an increased distance between heme and the active site (Fig. 4H). Substitutions found in SRS-5, including Ile369 and Ala370, had the most significant effect on CYP3A4 activity (Sevrioukova, 2017; Sevrioukova and Poulos, 2017). Interestingly, interactions between the heme group and assayed substrates differed depending on amino acid substitutions, including Ile369Val, Ile369Asn, CYP3A4.33 (Ala370Ser), and Ala370Val, which ultimately could cause the observed differences in enzymatic activities (Fig. 4, H, I, K, and L). Interestingly, variants carrying amino acid substitutions around these sites showed a weaker correlation with the decreasing catalytic efficiencies for midazolam and testosterone than CYP3A4.1, as shown in Fig. 3. The differences in the interactions between amino acids (Ile369 and Ala370) and heme may play a critical role in the enzymatic activity (Fig. 4, H, J, K, and M). The Tyr399Ser substitution also affected the SRS-5 region by forming a hydrogen bond between Ser399 and Met371 (Fig. 4G). Thus, amino acid substitutions in SRS regions interfered with CYP3A4 activities, but the extent of the resulting differences varied among substrates. Although our study detected a correlation between midazolam and testosterone activities, in the context of previous data (Supplemental Table 3), the degree of CYP3A4 activity reduction is apparently highly substrate-specific.

The effects of CYP3A4*16 and CYP3A4*18 on several clinically used drugs have been evaluated because of their relatively high allele frequencies (1.4%–5.0% and 1.3%–10%, respectively) (Lamba et al., 2002b; Yamamoto et al., 2003; Fukushima-Uesaka et al., 2004). However, their enzymatic activities also appear to be substrate-dependent (Supplemental Table 3). In addition to SNVs in the coding region, intron variants, such as CYP3A4*16B and *18B, should also be evaluated in vivo because several promoter and intronic variants, including CYP3A4*1B and *22, exhibit decreased activities due to reduced mRNA expression (Murayama et al., 2002; Elens et al., 2013). Notably, CYP3A5*3 and other CYP3A5 allelic variants play important clinical roles because the effects of CYP3A4 and CYP3A5 on the metabolism of numerous drugs overlap (Fukushima-Uesaka et al., 2004). Moreover, microRNA-27b was identified as a potential regulator of CYP3A4 mRNA level, underscoring the fact that other regulatory factors, and not only CYP3A4 allelic variants, affect CYP3A4 mRNA and protein levels (Pan et al., 2009; Ekström et al., 2015). Thus, future studies should establish a dosing algorithm based on interindividual differences in the in vivo activity of CYP3A.

CYP3A4 can contribute to the catalysis of several sites in the same compound. For example, midazolam is metabolized by CYP3A4 to the main metabolite 1′-hydroxymidazolam and the secondary metabolite 4-hydroxymidazolam (Christensen et al., 2009). According to a previous report, the enzymatic activities of CYP3A4.16 and .18 for these two midazolam metabolites decreased to similar levels (Maekawa et al., 2010). Moreover, even with the amino acid sequence differences between CYP3A4 and CYP3A5 (83% homology), 4-hydroxymidazolam formation was several times lower than 1′-hydroxymidazolam, suggesting that single amino acid substitutions in CYP3A4 may not dramatically influence the catalysis of different primary products dependent on changing mutation sites (Niwa et al., 2014). However, a detailed study encompassing all midazolam metabolites is needed to fully understand the differences in the metabolic profiles associated with CYP3A4 genetic polymorphisms.

Cytochrome P450 variant characterization by assessing the enzymatic activities for all metabolites and using all possible substrates may prove essential to evaluate the impact caused by amino acid substitutions and to improve treatment outcomes. However, previously developed computational calculations have been successful in predicting the activities of over 90% of representative cytochrome P450 variants (Zhou et al., 2019b). Hence, using computational prediction to assess the effect of amino acid substitutions on enzymatic cytochrome P450 activities may be an acceptable surrogate approach, considering the challenging nature of in vitro and in vivo analyses.

In summary, we determined the in vitro metabolic activities of wild-type CYP3A4 and 40 CYP3A4 variant proteins using our recently established heterologous expression system. These data reveal the functional effects caused by CYP3A4 genetic polymorphisms, including 11 novel variants identified among 4773 Japanese individuals, and have been generated using a single microsomal protein preparation per CYP3A4 variant to reveal their qualitative features. We recognize the limitations associated with a single biologic replicate, and additional biologic replicates would provide greater confidence in estimates of intrinsic clearance for use in future applications. Reduced activity or inactive variants caused by CYP3A4 allelic variations located in the coding region may be the underlying cause of an altered therapeutic response leading to poor treatment outcomes. However, our assay that uses a cDNA-based expression system cannot evaluate the influence of splicing variants. Therefore, further studies including additional biologic replicates as well as the use of minigenes to assess the potential functional consequences of variants affecting splice sites on catalytic activity are required prior to further downstream applications such as in vitro-in vivo extrapolation. In the future, CYP3A4 genetic polymorphisms should be used in combination with other factors affecting the in vivo activity of CYP3A4, such as CYP3A4 inducers and inhibitors and CYP3A5 activity, to aid the implementation of personalized medicine, individualized dosing regimens, and predictive risk panels. Our current findings should provide guidance in advancing these efforts.

Acknowledgments

We thank the Biomedical Research Core at the Tohoku University Graduate School of Medicine for their technical support. We would like to thank Mizuguchi (Graduate School of Pharmaceutical Sciences, Osaka University) for gifting the plasmid vector. We would like to thank Editage (www.editage.com) for English language editing.

Authorship Contributions

Participated in research design: Kumondai, Hiratsuka.

Conducted experiments: Kumondai, Gutiérrez Rico, Hishinuma, Ueda, Saito, Abe, Maekawa.

Contributed new reagents or analytic tools: Saigusa, Kinoshita, Maekawa, Mano, Hirasawa, Hiratsuka.

Performed data analysis: Kumondai, Saito, Tadaka, Nakayoshi, Oda, Hiratsuka.

Wrote or contributed to the writing of the manuscript: Kumondai, Gutiérrez Rico, Hiratsuka.

Footnotes

    • Received September 28, 2020.
    • Accepted December 24, 2020.
  • M.H. was supported by grants from Japan Agency for Medical Research and Development (AMED) [Grant 20kk0305009], Takahashi Industrial and Economic Research Foundation, and Smoking Research Foundation. M.K. was supported by Japan Society for the Promotion of Science [Grant 19J10744] and the Pharmaceutical Society of Japan [Grant N-170603]. S.S., K.K., S.T., and D.S. were supported by grants from AMED [Grant JP20km0105001 and JP20km0105002]. This research was also supported in part by the Tohoku Medical Megabank Project: Promoting Public Utilization of Advanced Research Infrastructure, and the Sharing and Administrative Network for Research Equipment funded by the Ministry of Education, Culture, Sports, Science and Technology (MEXT). K.K. and S.T. were supported by grants for the Facilitation of R&D Platform for AMED Genome Medicine Support from AMED (Grant Number JP20km0405001).

  • The authors declare no conflict of interest.

  • https://doi.org/10.1124/dmd.120.000261.

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

Abbreviations

CO
carbon monoxide
CPR
cytochrome P450 oxidoreductase
3D
three-dimensional
LC-MS/MS
liquid chromatography-tandem mass spectrometry
m/z
mass-to-charge ratio
PCR
polymerase chain reaction
SNV
single nucleotide variation
SRS
substrate recognition site
WGS
whole-genome sequencing
  • Copyright © 2021 by The American Society for Pharmacology and Experimental Therapeutics

References

  1. ↵
    1. Apellániz-Ruiz M,
    2. Lee MY,
    3. Sánchez-Barroso L,
    4. Gutiérrez-Gutiérrez G,
    5. Calvo I,
    6. García-Estévez L,
    7. Sereno M,
    8. García-Donás J,
    9. Castelo B,
    10. Guerra E, et al.
    (2015) Whole-exome sequencing reveals defective CYP3A4 variants predictive of paclitaxel dose-limiting neuropathy. Clin Cancer Res 21:322–328.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Benkaidali L,
    2. André F,
    3. Moroy G,
    4. Tangour B,
    5. Maurel F, and
    6. Petitjean M
    (2019) Four major channels detected in the cytochrome P450 3A4: a Step toward understanding its multispecificity. Int J Mol Sci 20:987.
    OpenUrl
  3. ↵
    1. Bridges A,
    2. Gruenke L,
    3. Chang YT,
    4. Vakser IA,
    5. Loew G, and
    6. Waskell L
    (1998) Identification of the binding site on cytochrome P450 2B4 for cytochrome b5 and cytochrome P450 reductase. J Biol Chem 273:17036–17049.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Chen B,
    2. Zhang XD,
    3. Wen J,
    4. Zhang B,
    5. Chen D,
    6. Wang S,
    7. Cai JP, and
    8. Hu GX
    (2020) Effects of 26 recombinant CYP3A4 variants on brexpiprazole metabolism. Chem Res Toxicol 33:172–180.
    OpenUrl
  5. ↵
    1. Christensen H,
    2. Mathiesen L,
    3. Postvoll LW,
    4. Winther B, and
    5. Molden E
    (2009) Different enzyme kinetics of midazolam in recombinant CYP3A4 microsomes from human and insect sources. Drug Metab Pharmacokinet 24:261–268.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Clark DP and
    2. Pazdernik NJ
    (2016) Biotechnology, Elsevier/AP Cell, Amsterdam.
  7. ↵
    1. Eiselt R,
    2. Domanski TL,
    3. Zibat A,
    4. Mueller R,
    5. Presecan-Siedel E,
    6. Hustert E,
    7. Zanger UM,
    8. Brockmoller J,
    9. Klenk HP,
    10. Meyer UA, et al.
    (2001) Identification and functional characterization of eight CYP3A4 protein variants. Pharmacogenetics 11:447–458.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Ekström L,
    2. Skilving I,
    3. Ovesjö ML,
    4. Aklillu E,
    5. Nylén H,
    6. Rane A,
    7. Diczfalusy U, and
    8. Björkhem-Bergman L
    (2015) miRNA-27b levels are associated with CYP3A activity in vitro and in vivo. Pharmacol Res Perspect 3:e00192.
    OpenUrl
  9. ↵
    1. Elens L,
    2. van Gelder T,
    3. Hesselink DA,
    4. Haufroid V, and
    5. van Schaik RH
    (2013) CYP3A4*22: promising newly identified CYP3A4 variant allele for personalizing pharmacotherapy. Pharmacogenomics 14:47–62.
    OpenUrlCrossRefPubMed
  10. ↵
    1. El Mahjoub A and
    2. Staub C
    (2000) High-performance liquid chromatographic method for the determination of benzodiazepines in plasma or serum using the column-switching technique. J Chromatogr B Biomed Sci Appl 742:381–390.
    OpenUrlPubMed
  11. ↵
    1. Fang P,
    2. Tang PF,
    3. Xu RA,
    4. Zheng X,
    5. Wen J,
    6. Bao SS,
    7. Cai JP, and
    8. Hu GX
    (2017) Functional assessment of CYP3A4 allelic variants on lidocaine metabolism in vitro. Drug Des Devel Ther 11:3503–3510.
    OpenUrl
  12. ↵
    1. Foti RS,
    2. Rock DA,
    3. Wienkers LC, and
    4. Wahlstrom JL
    (2010) Selection of alternative CYP3A4 probe substrates for clinical drug interaction studies using in vitro data and in vivo simulation. Drug Metab Dispos 38:981–987.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Fujikura K,
    2. Ingelman-Sundberg M, and
    3. Lauschke VM
    (2015) Genetic variation in the human cytochrome P450 supergene family. Pharmacogenet Genomics 25:584–594.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Fukushima-Uesaka H,
    2. Saito Y,
    3. Watanabe H,
    4. Shiseki K,
    5. Saeki M,
    6. Nakamura T,
    7. Kurose K,
    8. Sai K,
    9. Komamura K,
    10. Ueno K, et al.
    (2004) Haplotypes of CYP3A4 and their close linkage with CYP3A5 haplotypes in a Japanese population. Hum Mutat 23:100.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Gómez-Bravo MA,
    2. Apellaniz-Ruiz M,
    3. Salcedo M,
    4. Fondevila C,
    5. Suarez F,
    6. Castellote J,
    7. Rufian S,
    8. Pons JA,
    9. Bilbao I,
    10. Alamo JM, et al.
    (2018) Influence of donor liver CYP3A4*20 loss-of-function genotype on tacrolimus pharmacokinetics in transplanted patients. Pharmacogenet Genomics 28:41–48.
    OpenUrl
  16. ↵
    1. Gutiérrez Rico EM,
    2. Kikuchi A,
    3. Saito T,
    4. Kumondai M,
    5. Hishinuma E,
    6. Kaneko A,
    7. Chan CW,
    8. Gitaka J,
    9. Nakayoshi T,
    10. Oda A, et al.
    (2020) CYP2D6 genotyping analysis and functional characterization of novel allelic variants in a Ni-Vanuatu and Kenyan population by assessing dextromethorphan O-demethylation activity. Drug Metab Pharmacokinet 35:89–101.
    OpenUrl
  17. ↵
    1. Hasemann CA,
    2. Kurumbail RG,
    3. Boddupalli SS,
    4. Peterson JA, and
    5. Deisenhofer J
    (1995) Structure and function of cytochromes P450: a comparative analysis of three crystal structures. Structure 3:41–62.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Hendriks DFG,
    2. Vorrink SU,
    3. Smutny T,
    4. Sim SC,
    5. Nordling Å,
    6. Ullah S,
    7. Kumondai M,
    8. Jones BC,
    9. Johansson I,
    10. Andersson TB, et al.
    (2020) Clinically relevant cytochrome P450 3A4 induction mechanisms and drug screening in three-dimensional spheroid cultures of primary human hepatocytes. Clin Pharmacol Ther 108:844–855.
    OpenUrl
  19. ↵
    1. Hiratsuka M
    (2012) In vitro assessment of the allelic variants of cytochrome P450. Drug Metab Pharmacokinet 27:68–84.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Hiratsuka M
    (2016) Genetic polymorphisms and in vitro functional characterization of CYP2C8, CYP2C9, and CYP2C19 allelic variants. Biol Pharm Bull 39:1748–1759.
    OpenUrl
  21. ↵
    1. Hornstein BD,
    2. Roman D,
    3. Arévalo-Soliz LM,
    4. Engevik MA, and
    5. Zechiedrich L
    (2016) Effects of circular DNA length on transfection efficiency by electroporation into HeLa cells. PLoS One 11:e0167537.
    OpenUrlCrossRef
  22. ↵
    1. Huang SM,
    2. Strong JM,
    3. Zhang L,
    4. Reynolds KS,
    5. Nallani S,
    6. Temple R,
    7. Abraham S,
    8. Habet SA,
    9. Baweja RK,
    10. Burckart GJ, et al.
    (2008) New era in drug interaction evaluation: US Food and Drug Administration update on CYP enzymes, transporters, and the guidance process. J Clin Pharmacol 48:662–670.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Ingelman-Sundberg M,
    2. Mkrtchian S,
    3. Zhou Y, and
    4. Lauschke VM
    (2018) Integrating rare genetic variants into pharmacogenetic drug response predictions. Hum Genomics 12:26.
    OpenUrlCrossRef
  24. ↵
    1. Jackson KD,
    2. Durandis R, and
    3. Vergne MJ
    (2018) Role of cytochrome P450 enzymes in the metabolic activation of tyrosine kinase inhibitors. Int J Mol Sci 19:2367.
    OpenUrl
  25. ↵
    1. Kumondai M,
    2. Hishinuma E,
    3. Gutiérrez Rico EM,
    4. Ito A,
    5. Nakanishi Y,
    6. Saigusa D,
    7. Hirasawa N, and
    8. Hiratsuka M
    (2020) Heterologous expression of high-activity cytochrome P450 in mammalian cells. Sci Rep 10:14193.
    OpenUrl
  26. ↵
    1. Kumondai M,
    2. Hosono H,
    3. Maekawa M,
    4. Yamaguchi H,
    5. Mano N,
    6. Oda A,
    7. Hirasawa N, and
    8. Hiratsuka M
    (2018) Functional characterization of 9 CYP2A13 allelic variants by assessment of nicotine C-oxidation and coumarin 7-hydroxylation. Drug Metab Pharmacokinet 33:82–89.
    OpenUrl
  27. ↵
    1. Lamba JK,
    2. Lin YS,
    3. Schuetz EG, and
    4. Thummel KE
    (2002a) Genetic contribution to variable human CYP3A-mediated metabolism. Adv Drug Deliv Rev 54:1271–1294.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Lamba JK,
    2. Lin YS,
    3. Thummel K,
    4. Daly A,
    5. Watkins PB,
    6. Strom S,
    7. Zhang J, and
    8. Schuetz EG
    (2002b) Common allelic variants of cytochrome P4503A4 and their prevalence in different populations. Pharmacogenetics 12:121–132.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Lee SJ,
    2. Bell DA,
    3. Coulter SJ,
    4. Ghanayem B, and
    5. Goldstein JA
    (2005) Recombinant CYP3A4*17 is defective in metabolizing the hypertensive drug nifedipine, and the CYP3A4*17 allele may occur on the same chromosome as CYP3A5*3, representing a new putative defective CYP3A haplotype. J Pharmacol Exp Ther 313:302–309.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Li YH,
    2. Lin QM,
    3. Pang NH,
    4. Zhang XD,
    5. Huang HL,
    6. Cai JP, and
    7. Hu GX
    (2019a) Functional characterization of 27 CYP3A4 protein variants to metabolize regorafenib in vitro. Basic Clin Pharmacol Toxicol 125:337–344.
    OpenUrl
  31. ↵
    1. Li YH,
    2. Lu XR,
    3. Lin QM,
    4. Huang HL,
    5. Liang XL,
    6. Cai JP,
    7. Cui J, and
    8. Hu GX
    (2019b) Functional characterization of 27 CYP3A4 variants on macitentan metabolism in vitro. J Pharm Pharmacol 71:1677–1683.
    OpenUrl
  32. ↵
    1. Lin QM,
    2. Li YH,
    3. Liu Q,
    4. Pang NH,
    5. Xu RA,
    6. Cai JP, and
    7. Hu GX
    (2019a) Functional characteristics of CYP3A4 allelic variants on the metabolism of loperamide in vitro. Infect Drug Resist 12:2809–2817.
    OpenUrl
  33. ↵
    1. Lin QM,
    2. Li YH,
    3. Lu XR,
    4. Wang R,
    5. Pang NH,
    6. Xu RA,
    7. Cai JP, and
    8. Hu GX
    (2019b) Characterization of genetic variation in CYP3A4 on the metabolism of cabozantinib in vitro. Chem Res Toxicol 32:1583–1590.
    OpenUrl
  34. ↵
    1. Maekawa K,
    2. Harakawa N,
    3. Yoshimura T,
    4. Kim SR,
    5. Fujimura Y,
    6. Aohara F,
    7. Sai K,
    8. Katori N,
    9. Tohkin M,
    10. Naito M, et al.
    (2010) CYP3A4*16 and CYP3A4*18 alleles found in East Asians exhibit differential catalytic activities for seven CYP3A4 substrate drugs. Drug Metab Dispos 38:2100–2104.
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Murayama N,
    2. Nakamura T,
    3. Saeki M,
    4. Soyama A,
    5. Saito Y,
    6. Sai K,
    7. Ishida S,
    8. Nakajima O,
    9. Itoda M,
    10. Ohno Y, et al.
    (2002) CYP3A4 gene polymorphisms influence testosterone 6beta-hydroxylation. Drug Metab Pharmacokinet 17:150–156.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Nagasaki M,
    2. Yasuda J,
    3. Katsuoka F,
    4. Nariai N,
    5. Kojima K,
    6. Kawai Y,
    7. Yamaguchi-Kabata Y,
    8. Yokozawa J,
    9. Danjoh I,
    10. Saito S, et al., and ToMMo Japanese Reference Panel Project
    (2015) Rare variant discovery by deep whole-genome sequencing of 1,070 Japanese individuals. Nat Commun 6:8018.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Nguyen V,
    2. Espiritu M, and
    3. Elbarbry F
    (2020) Development and validation of a sensitive and specific LC-MS/MS cocktail assay for CYP450 enzymes: application to study the effect of catechin on rat hepatic CYP activity. Biomed Chromatogr 34:e4789.
    OpenUrl
  38. ↵
    1. Niwa T,
    2. Okamoto A,
    3. Narita K,
    4. Toyota M,
    5. Kato K,
    6. Kobayashi K, and
    7. Sasaki S
    (2020) Comparison of steroid hormone hydroxylation mediated by cytochrome P450 3A subfamilies. Arch Biochem Biophys 682:108283.
    OpenUrl
  39. ↵
    1. Niwa T,
    2. Yasumura M,
    3. Murayama N, and
    4. Yamazaki H
    (2014) Comparison of catalytic properties of cytochromes P450 3A4 and 3A5 by molecular docking simulation. Drug Metab Lett 8:43–50.
    OpenUrl
  40. ↵
    1. Oda A,
    2. Yamaotsu N, and
    3. Hirono S
    (2005) New AMBER force field parameters of heme iron for cytochrome P450s determined by quantum chemical calculations of simplified models. J Comput Chem 26:818–826.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Pan YZ,
    2. Gao W, and
    3. Yu AM
    (2009) MicroRNAs regulate CYP3A4 expression via direct and indirect targeting. Drug Metab Dispos 37:2112–2117.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Racha JK,
    2. Zhao ZS,
    3. Olejnik N,
    4. Warner N,
    5. Chan R,
    6. Moore D, and
    7. Satoh H
    (2003) Substrate dependent inhibition profiles of fourteen drugs on CYP3A4 activity measured by a high throughput LCMS/MS method with four probe drugs, midazolam, testosterone, nifedipine and terfenadine. Drug Metab Pharmacokinet 18:128–138.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Saito T,
    2. Gutiérrez Rico EM,
    3. Kikuchi A,
    4. Kaneko A,
    5. Kumondai M,
    6. Akai F,
    7. Saigusa D,
    8. Oda A,
    9. Hirasawa N, and
    10. Hiratsuka M
    (2018) Functional characterization of 50 CYP2D6 allelic variants by assessing primaquine 5-hydroxylation. Drug Metab Pharmacokinet 33:250–257.
    OpenUrlCrossRef
  44. ↵
    1. Sakurai F,
    2. Kawabata K,
    3. Yamaguchi T,
    4. Hayakawa T, and
    5. Mizuguchi H
    (2005) Optimization of adenovirus serotype 35 vectors for efficient transduction in human hematopoietic progenitors: comparison of promoter activities. Gene Ther 12:1424–1433.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Schärfe CPI,
    2. Tremmel R,
    3. Schwab M,
    4. Kohlbacher O, and
    5. Marks DS
    (2017) Genetic variation in human drug-related genes. Genome Med 9:117.
    OpenUrl
  46. ↵
    1. Sevrioukova IF
    (2017) High-level production and properties of the cysteine-depleted cytochrome P450 3A4. Biochemistry 56:3058–3067.
    OpenUrl
  47. ↵
    1. Sevrioukova IF and
    2. Poulos TL
    (2017) Structural basis for regiospecific midazolam oxidation by human cytochrome P450 3A4. Proc Natl Acad Sci USA 114:486–491.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Watanabe T,
    2. Saito T,
    3. Rico EMG,
    4. Hishinuma E,
    5. Kumondai M,
    6. Maekawa M,
    7. Oda A,
    8. Saigusa D,
    9. Saito S,
    10. Yasuda J, et al.
    (2018) Functional characterization of 40 CYP2B6 allelic variants by assessing efavirenz 8-hydroxylation. Biochem Pharmacol 156:420–430.
    OpenUrl
  49. ↵
    1. Werk AN,
    2. Lefeldt S,
    3. Bruckmueller H,
    4. Hemmrich-Stanisak G,
    5. Franke A,
    6. Roos M,
    7. Küchle C,
    8. Steubl D,
    9. Schmaderer C,
    10. Bräsen JH, et al.
    (2014) Identification and characterization of a defective CYP3A4 genotype in a kidney transplant patient with severely diminished tacrolimus clearance. Clin Pharmacol Ther 95:416–422.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Xiao K,
    2. Gao J,
    3. Weng SJ,
    4. Fang Y,
    5. Gao N,
    6. Wen Q,
    7. Jin H, and
    8. Qiao HL
    (2019) CYP3A4/5 activity probed with testosterone and midazolam: correlation between two substrates at the microsomal and enzyme levels. Mol Pharm 16:382–392.
    OpenUrl
  51. ↵
    1. Xu RA,
    2. Wen J,
    3. Tang P,
    4. Wang C,
    5. Xie S,
    6. Zhang BW,
    7. Zhou Q,
    8. Cai JP, and
    9. Hu GX
    (2018) Functional characterization of 22 CYP3A4 protein variants to metabolize ibrutinib in vitro. Basic Clin Pharmacol Toxicol 122:383–387.
    OpenUrl
  52. ↵
    1. Yamamoto T,
    2. Nagafuchi N,
    3. Ozeki T,
    4. Kubota T,
    5. Ishikawa H,
    6. Ogawa S,
    7. Yamada Y,
    8. Hirai H, and
    9. Iga T
    (2003) CYP3A4*18: it is not rare allele in Japanese population. Drug Metab Pharmacokinet 18:267–268.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Yang CC,
    2. Zheng X,
    3. Liu TH,
    4. Wang CC,
    5. Tang PF,
    6. Chen Z,
    7. Zhang BW,
    8. Fang P,
    9. Hu GX, and
    10. Cai JP
    (2019) Functional characterization of 21 CYP3A4 variants on amiodarone metabolism in vitro. Xenobiotica 49:120–126.
    OpenUrl
  54. ↵
    1. Zanger UM,
    2. Momoi K,
    3. Hofmann U,
    4. Schwab M, and
    5. Klein K
    (2020) Tri-allelic haplotypes determine and differentiate functionally normal allele CYP2D6*2 and impaired allele CYP2D6*41. Clin Pharmacol Ther DOI: 10.1002/cpt.2078 [published ahead of print].
  55. ↵
    1. Zanger UM and
    2. Schwab M
    (2013) Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation. Pharmacol Ther 138:103–141.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Zawaira A,
    2. Ching LY,
    3. Coulson L,
    4. Blackburn J, and
    5. Wei YC
    (2011) An expanded, unified substrate recognition site map for mammalian cytochrome P450s: analysis of molecular interactions between 15 mammalian CYP450 isoforms and 868 substrates. Curr Drug Metab 12:684–700.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Zhou XY,
    2. Hu XX,
    3. Wang CC,
    4. Lu XR,
    5. Chen Z,
    6. Liu Q,
    7. Hu GX, and
    8. Cai JP
    (2019a) Enzymatic activities of CYP3A4 allelic variants on quinine 3-hydroxylation in vitro. Front Pharmacol 10:591.
    OpenUrl
  58. ↵
    1. Zhou Y,
    2. Ingelman-Sundberg M, and
    3. Lauschke VM
    (2017) Worldwide distribution of cytochrome P450 alleles: a meta-analysis of population-scale sequencing projects. Clin Pharmacol Ther 102:688–700.
    OpenUrlCrossRef
  59. ↵
    1. Zhou Y,
    2. Mkrtchian S,
    3. Kumondai M,
    4. Hiratsuka M, and
    5. Lauschke VM
    (2019b) An optimized prediction framework to assess the functional impact of pharmacogenetic variants. Pharmacogenomics J 19:115–126.
    OpenUrl
PreviousNext
Back to top

In this issue

Drug Metabolism and Disposition: 49 (3)
Drug Metabolism and Disposition
Vol. 49, Issue 3
1 Mar 2021
  • 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.
Functional Characterization of 40 CYP3A4 Variants by Assessing Midazolam 1′-Hydroxylation and Testosterone 6β-Hydroxylation
(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

Functional Characterization of 40 CYP3A4 Variants

Masaki Kumondai, Evelyn Marie Gutiérrez Rico, Eiji Hishinuma, Akiko Ueda, Sakae Saito, Daisuke Saigusa, Shu Tadaka, Kengo Kinoshita, Tomoki Nakayoshi, Akifumi Oda, Ai Abe, Masamitsu Maekawa, Nariyasu Mano, Noriyasu Hirasawa and Masahiro Hiratsuka
Drug Metabolism and Disposition March 1, 2021, 49 (3) 212-220; DOI: https://doi.org/10.1124/dmd.120.000261

Citation Manager Formats

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

Share
Research ArticleArticle

Functional Characterization of 40 CYP3A4 Variants

Masaki Kumondai, Evelyn Marie Gutiérrez Rico, Eiji Hishinuma, Akiko Ueda, Sakae Saito, Daisuke Saigusa, Shu Tadaka, Kengo Kinoshita, Tomoki Nakayoshi, Akifumi Oda, Ai Abe, Masamitsu Maekawa, Nariyasu Mano, Noriyasu Hirasawa and Masahiro Hiratsuka
Drug Metabolism and Disposition March 1, 2021, 49 (3) 212-220; DOI: https://doi.org/10.1124/dmd.120.000261
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Visual Overview
    • 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

  • Adipocyte PXR does not play an essential role in obesity.
  • CYP3A-mediated oxidation of DABE and BIBR0951
  • Biodistribution of Lipid in Rats
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