Abstract
The major mode of metabolism of nicotine is by hydroxylation via cytochrome P450 (CYP) 2A6, but it can also undergo glucuronidation by UDP-glucuronosyltransferases and oxidation by flavin monooxygenases (FMO). The goal of this study was to examine the potential importance of FMOs in nicotine metabolism and assess the potential impact of missense polymorphisms in active FMOs on nicotine-N′-oxide (NOX) formation.
Urine samples from 106 current Chinese smokers were analyzed for nicotine metabolites by mass spectrometry. Wild-type FMOs 1–5 and their most prevalent nonsynonymous variants were cloned and overexpressed in HEK293 cells, and were tested in oxidation reactions against nicotine.
A strong inverse correlation was observed between the ratio of urinary 3′-hydroxycotinine/cotinine, a measure of CYP2A6 activity, and the urinary levels of NOX alone (r = −0.383; P < 0.001) or NOX measured as a ratio of total nicotine metabolites (r = −0.414; P < 0.001) in smokers. In addition to FMO1 and FMO3, the functional FMO2427Q isoform was active against nicotine, whereas FMO4 and FMO5 exhibited low activity against nicotine (Km > 5.0 mmol/L). Significant (P < 0.05) decreases in N′-oxidation activity (Vmax/Km) were observed for the FMO1I303V, FMO3N61S, FMO3D132H, FMO3V257M, and FMO3E308G variants in vitro when compared with their respective wild-type isoforms; the truncated FMO2Q472stop isoform exhibited no enzyme activity.
These data indicate that increases in nicotine-N′-oxidation occur in subjects with deficient CYP2A6 activity, and that several FMO enzymes are active in nicotine-N′-oxidation.
Several common missense FMO variants are associated with altered enzyme activity against nicotine and may play an important role in nicotine metabolism in low-CYP2A6 activity subjects.
Introduction
Currently, nearly 480,000 deaths in the United States are attributable to cigarette smoking and second-hand smoke exposure annually, and tobacco use is the leading cause of preventable deaths worldwide (1, 2). Approximately 90% of lung cancer in humans is attributed to cigarette smoking and lung cancer results in more deaths each year than any other cancer type (3, 4).
Nicotine is abundant in tobacco and tobacco smoke, activating the nicotinic acetylcholine receptors in the brain, resulting in tobacco addiction (5). Nicotine is rapidly metabolized in vivo, with a half-life of approximately 2 hours (6). The major metabolic pathway for nicotine is hydroxylation by cytochrome P450 (CYP) 2A6. Nicotine also exhibits metabolism through either glucuronidation, catalyzed by the UDP-glucuronosyltransferases (UGT), or oxygenation by flavin monooxygenases (FMO; see Fig. 1). Besides nicotine, the FMO family catalyzes the oxidation of heteroatoms (such as nitrogen, sulfur, phosphorus, and selenium), xenobiotics (such as drugs and pesticides; ref. 7), and endogenous compounds (such as cysteamine, ref. 8; trimethylamine, ref. 9; and phenethylamine and tyramine, ref. 10). The human family of FMOs includes 5 known functional genes—FMO1, FMO2, FMO3, FMO4, and FMO5. The FMO1, FMO2, FMO3, and FMO4 genes are clustered in chromosome 1, region q24.3, while the FMO5 gene is located at q21 of the same chromosome (11). These enzymes exhibit a tissue-specific expression pattern in adults. FMO1 is mainly expressed in the kidney and FMO3-FMO5 in the liver, with FMO4 and FMO5 expressed at lower levels than FMO3 (12). FMO2 is mainly expressed in the lung (13), but its expression in other organ tissues has not been thoroughly investigated. The most common FMO2*1 allele (rs6661174) acts as a pseudogene, coding for an early stop codon (Q472Stop) resulting in a truncated protein without the last 64 C-terminal amino acids (13). This common variant accounts for >95% of all FMO2 alleles in all populations combined. The remaining 5% of FMO2 alleles encode a full-length, active FMO2 protein, with African-Americans exhibiting a minor allele frequency (MAF) of 14% for these functional alleles (14). The high overall allelic prevalence corresponding to the truncated form of this protein has likely resulted in the null expression of this enzyme observed in most human studies.
It has been reported that FMOs are the enzymes responsible for the stereoselective N′-oxidation of nicotine (15). Nicotine-N′-oxide (NOX) formation accounts for 4%–7% of total urinary nicotine metabolites in humans under normal circumstances (16, 17). In a small study of 5 subjects, when the CYP2A6-mediated hydroxylation pathway is diminished or impaired, urinary levels of NOX are increased to up to 31% of total nicotine metabolites (18). While some studies have shown that FMO1 and FMO3 may also play a role in NOX formation (15, 19), a comprehensive analysis of all five FMO enzymes on the impact of nicotine oxidation has not yet been performed.
Studies have shown a marked difference in tobacco addiction profiles among different populations (20). This could be due to several factors including genetic differences in the enzymes that metabolize nicotine (21–23). Functional polymorphisms in CYP2A6 have been correlated with differences in cotinine and 3-hydroxycotinine formation (22). Similarly, genetic variants in several UGTs are correlated with differences in the glucuronidation of urinary nicotine, cotinine, and 3-hydroxycotinine (24, 25). Genome-wide association studies (GWAS) have identified associations between FMO1 and FMO3 SNPs and nicotine dependence (7, 19, 26–28). Previous studies have shown that genetic variations in FMO3 in different smoker cohorts were associated with nicotine dependence (28), while a recent in vivo study demonstrated an association between the FMO3 E308G/E158K haplotype with nicotine-N′-oxidation and smoking behavior (29). However, the mechanism behind the functional role of those polymorphisms in nicotine metabolism and addiction was not addressed in these studies.
The goal of the present study was to examine the variability of NOX formation in smoker's urine, fully characterize FMO involvement in NOX formation, and assess the impact of the most prevalent nonsynonymous polymorphisms in active FMO enzymes on nicotine metabolite production in vitro.
Materials and Methods
Chemicals and materials
Pooled human liver, lung, esophagus, or kidney RNA was purchased from Biochain. Platinum Pfx DNA polymerase, the pcDNA3.1/V5-His-TOPO Mammalian Expression Vector, Lipofectamine 3000, One Shot TOP10 competent E. coli, Superscript II RT Kits, SuperScript VILO Synthesis Kits, anti-V5-HRP antibody, and Luria broth base were obtained from Invitrogen. The anti–calnexin-HRP antibody was purchased from Abcam, while the QuikChange II Site-Directed Mutagenesis Kit used to make FMO-variant vectors was acquired from Agilent. Oligonucleotides used for site-directed mutagenesis were manufactured by Integrated DNA Technologies. RNeasy Kits and QIAquick Gel Extraction Kits were purchased from Qiagen, while the GeneJet Plasmid Mini and Midipreps were purchased from Thermo Fisher Scientific. The Ambion PureLink RNA Mini Kit was purchased from Life Technologies, while the BCA Protein Assays used in protein assessment were purchased from Pierce. DMEM, Dulbecco PBS, FBS, and geneticin (G418) were purchased from Gibco. The NADPH Regeneration System was purchased from Corning. Nicotine tartrate and benzydamine hydrochloride used as substrates for N′-oxidation Activity Assays were purchased from Sigma-Aldrich. Pooled human liver microsomes were purchased from Sekisui XenoTech. TaqMan qPCR probes were purchased from AB Applied Biosystems. High performance liquid chromatography-grade ammonium acetate and acetonitrile were purchased from Thermo Fisher Scientific, while the ACQUITY UPLC BEH-HILIC (1.7 μm 2.1 × 100 mm) column was purchased from Waters. The internal standard d3-NOX was purchased from Toronto Research Chemicals.
Nicotine metabolites in the urine of smokers
Urine samples from 106 Chinese smokers who participated in two cancer epidemiology cohorts were analyzed. All 106 subjects were free of cancer at the latest follow-up in 2015. Among them, 63 were from the Shanghai Cohort Study, a prospective cohort of 18,244 Chinese men in Shanghai, China, when they were 45–64 years of age at enrollment during 1986–1989 (30, 31). The mean (SD) age at the time of urine sample collection was 56.6 (5.5) years. The remaining 43 subjects were from the Singapore Chinese Health Study, a prospective cohort of 63,257 Chinese men and women in Singapore at 45–74 years of age when they were recruited into the study in 1993–1998 (32, 33). Of the 43 subjects, 7 (16%) were women and the average age (SD) at urine collection was 67.0 (6.4) years. The Shanghai cohort subjects smoked an average of 15.3 cigarettes/day for 46.5 years as compared with 17.2 cigarettes/day and 31.7 years for the Singapore cohort subjects. In addition to nicotine, eight of its major metabolites (NOX, cotinine-N-oxide, trans-3-hydroxycotinine, nicotine-N-glucuronide, cotinine-N-glucuronide, trans-3′-hydroxycotinine-O-glucuronide, and 4-hydroxy-4-(3-pyridyl)-butanoid acid) were quantified in urine samples of the 106 smokers using a previously reported LC/MS method (17).
Generation of FMO-overexpressing cell lines
Pooled RNA (2 μg) from lung, liver, esophagus, or kidney was used to generate corresponding cDNAs using the SuperScript II RT-PCR kit, and cDNA corresponding to 10 ng RNA was used as template for PCR amplification of wild-type FMOs 1–5 using 1 U of Pfx polymerase and FMO gene-specific primers (see Supplementary Table S1). PCR was performed with an initial denaturation temperature of 94°C for 2 minutes, followed by 40 cycles of 94°C for 30 seconds, 45 seconds at the specific Tm for each FMO (see Supplementary Table S1), and 68°C for 105 seconds, and a final cycle of 10 minutes at 68°C. Sequences of each PCR-amplified wild-type FMO product were verified by Sanger sequencing, and each were cloned into the pcDNA 3.1/V5-His-TOPO vector in-frame with the plasmid V5 epitope tag at the C-terminus 3′-end prior to transformation into One Shot TOP10 competent E. coli using standard protocols. Insert orientation and sequences of plasmid DNA prepared from individual clones were confirmed by a second round of Sanger sequencing. Lipofectamine 3000 reagent was used to transfect HEK293 cells with 2.5 μg of each FMO-containing pcDNA 3.1/V5-His-TOPO vector, with cells grown in DMEM supplemented with 10% FBS and 1.5% geneticin (G418). The parent cell line used in this study, HEK293, was purchased from ATCC in 2015 and was authenticated by ATCC in 2017 using short-tandem repeat polymorphisms analysis. No test was performed to detect Mycoplasma in these cells.
FMO variants were created by site-directed mutagenesis of wild-type FMO-containing pcDNA3.1/V5-His-TOPO plasmids using QuikChange II XL kits and specific mutagenic FMO primers (see Supplementary Table S2). All FMO variants were confirmed by Sanger sequencing and transfected into HEK293 cells as described above. For cloning of the FMO2 full-length isoform, we used the high-prevalence FMO2*Q472Stop variant as template for site-directed mutagenesis to change the codon from the premature stop to a glutamine residue.
Enzyme kinetic assays and ultra performance liquid chromatography - tandem mass spectrometry analysis
Microsomal membrane fractions of FMO-overexpressing cell lines were prepared by differential centrifugation as described previously (34, 35) and used as the enzyme source for all activity assays. For relative FMO quantification, equal amounts of microsomal protein (20 μg) were loaded on 10% SDS-polyacrylamide gels and FMO protein quantity was determined by Western blot analysis using an anti–V5-HRP antibody in a 1:2,500 dilution. As a loading control for microsomal fractions, an anti–calnexin-HRP antibody was used in a 1:5,000 dilution for all Western blot analyses. ImageJ software (http://rsb.info.nih.gov/ij; NIH) was used to perform densitometry analysis, and the relative expression of each FMO-containing microsomal preparation was used for normalization in N′-oxidation activity assays.
N′-oxidation reactions (final volume = 25 μL) contained 25 μg of total microsomal protein, 50 mmol/L potassium phosphate, a NADPH regenerating system (1.55 mmol/L NADP+, 3.3 mmol/L glucose-6-phosphate, 3.3 mmol/L MgCl2, and 0.5 of 40 U/mL glucose-6-phosphate dehydrogenase), and varying concentrations of nicotine tartrate (0.05–5,000 μmol/L). For reactions with benzydamine, assays were performed as described above for nicotine using benzydamine hydrochloride at concentrations equivalent to the reported Km's for each FMO for this substrate (36). All incubations were performed for 60 minutes at 37°C with shaking at 300 rpm, with reactions terminated by the addition of 25 μL ice-cold acetonitrile containing internal standard (d3-NOX). Supernatants were collected after centrifugation at 16,100 × g for 10 minutes at 4°C prior to subsequent ultra performance liquid chromatography - tandem mass spectrometry (UPLC-MS/MS) analysis. Pooled human liver microsomes and untransfected parent HEK293 cell microsomes were used as the protein source for positive and negative assay controls, respectively.
UPLC-MS/MS was performed with mobile phases that consisted of 5 mmol/L of NH4AC (pH 5.7) in either 50% acetonitrile (volume for volume; buffer A) or 90% acetonitrile (buffer B). An UPLC-BEH-HILIC column (2.1 × 100 mm, 1.7 μm) was used for separation of NOX from nicotine as follows: 20% buffer A for 1.5 minutes, a linear gradient to 100% buffer A from 1.5 to 2.5 minutes, maintenance of 100% buffer A for 3 minutes, and a re-equilibrium step to the initial 20% buffer A conditions from 5.5 to 7 minutes (flow rate 0.4 mL/min). The injection volume was 2 μL using a column temperature of 30°C. MS/MS detection was performed in a Waters ACQUITY XEVO TQD instrument in MRM ESI+ mode. The MS/MS scan was performed using the following mass transitions: m/z 179.1→117.1, 182.1→117.1, and 163.1→105, to monitor NOX, d3-NOX, and nicotine, respectively. The desolvation temperature was 500°C, with 800 L/hour of nitrogen gas. The collision energy was optimized at 28 V, 27 V, and 28 V for NOX, nicotine, and d3-NOX, respectively. A cone voltage of 25 V and 0.5 second dwell time resulted in high-sensitivity detection of NOX and d3-NOX. Metabolite retention times observed in the enzymatic incubations were compared with the retention times of d3 internal standard metabolites. UPLC-MS/MS detection for benzydamine and benzydamine-N-oxide was performed using a previously reported separation method (37).
mRNA expression analysis
Total RNA was isolated from wild-type FMO3- (FMO3158E/308E), FMO3158K/308E-, FMO3158E/308G-, and FMO3158K/308G-expressing HEK293 cells using the Ambion PureLink RNA Mini kit according to the manufacturer's protocol. Total RNA (1 μg) was used to synthesize cDNA using the SuperScript VILO synthesis kit following the manufacturer's protocol. qPCR analysis was performed using a TaqMan probe specific for the FMO3 human gene (Hs00199368_m1) while a TaqMan probe specific for the RPLP0 gene (Hs00420895_gH) was used as a housekeeping gene. Expression of the FMO3 gene was normalized to the internal control RPLP0 and the fold-change was calculated using the delta-delta Ct method (38). All experiments were performed in triplicate.
Statistical analysis
To compare the levels of NOX formation in subjects stratified by the nicotine metabolic ratio (NMR), statistical stratifications by median and quartiles were performed. Once the samples were stratified, relative FMO activity was determined by the levels of NOX/total nicotine equivalents (TNE, equivalent to nicotine plus its eight measured metabolites) in the stratified groups. The distributions of the urinary biomarkers measured or their ratios were markedly skewed toward high values, which were normalized to a large extent by transformation to logarithmic values. Therefore, formal statistical testing was performed on logarithmically transformed values, and geometric (as opposed to arithmetic) means are presented. Statistical analyses were conducted using the statistics program R v.3.4. All analyses were two-sided and considered significant if P < 0.05 for all tests. Continuous outcome variables were analyzed using Student t test or two-way ANOVA using a Tukey post hoc test for multiple corrections for all possible pairwise comparisons (39).
Kinetic parameters were determined from the Michaelis–Menten equation using GraphPad Prism version 6.01 (GraphPad Software). Relative maximum reaction rates (Vmax) were calculated as pmol NOX/L/min, with values normalized to wild-type FMO-V5 microsomal protein as determined by Western blot analysis using ImageJ software as described above. All reported values represent the results (e.g., mean ± SD) of three independent experiments, with triplicates run for each experiment. The activity of each variant enzyme was compared with its corresponding wild-type isoform using the Student t test. A two-tailed P value of less than 0.05 was considered the threshold for statistical significance.
Results
Nicotine metabolites in smokers
Table 1 shows the levels of nicotine and its eight major metabolites in the urine of Chinese smokers. NOX accounted for a mean of 7.4% (95% confidence interval, 6.7%–8.2%) of the urinary TNE in the 106 smokers examined. To examine whether NOX formation differed based on CYP2A6 activity, subjects were stratified based on their NMR, which is the ratio of urinary 3′-hydroxycotinine/cotinine in smokers and which has been shown to be an effective marker of CYP2A6 activity (40). A strong negative correlation was observed between the urinary NMR ratio and both urinary NOX levels alone (r = −0.383; P < 0.001), and urinary NOX measured as a ratio of TNE (r = −0.414; P < 0.001; Fig. 2A). When stratifying urine samples at the median value of NMR, the mean percentage of the urinary NOX/TNE ratio for each stratified group was 9.2% for the low CYP2A6 activity group and 6.0% for the high activity group, a significant (P < 0.001) 1.53-fold increase in NOX levels between groups (Fig. 2B). When stratifying samples by quartiles for NMR, the mean % urinary NOX/TNE for each group (low to high NMR) was 11.8%, 7.3%, 6.7%, and 5.4%, respectively, corresponding to a significant trend (P < 0.001) with a 2.2-fold difference between the lowest versus highest NMR quartiles (Fig. 2C). Furthermore, when subjects with the lowest 5% urinary NMR (n = 5) were compared with the highest 5% urinary NMR samples (n = 5), the percent urinary NOX/TNE were 18.8 in the former compared with 5.5 in the latter, respectively (Supplementary Fig. S1).
FMO activities against nicotine
HEK293 cell lines expressing V5-tagged FMO enzymes were created and FMO expression in microsomal fractions was verified by Western blot analysis using an anti-V5 antibody (Supplementary Fig. S2). Wild-type FMOs 1–5 were successfully expressed in HEK293 cells, with the relative level of FMO expression calculated based on the relative band intensities of the V5 signal and were used to calculate relative Vmax's obtained in N′-oxidation activity assays. Catalytic activities for each of the FMO enzymes were verified using benzydamine as a probe substrate, with all wild-type and variant FMO enzymes demonstrating benzydamine N′-oxidation activity. The Km values found for each FMO against benzydamine (60 ± 8 μmol/L, 440 ± 65 μmol/L, 80 ± 6 μmol/L, >3 mmol/L, and >2 mmol/L for FMOs 1, 2, 3, 4, and 5, respectively) were similar to those reported previously (36).
Using microsomes isolated from FMO-overexpressing HEK293 cell lines, in vitro assays demonstrated that all of the wild-type FMOs (FMO1, FMO2, FMO3, FMO4, and FMO5) exhibited activity against nicotine. As shown by UPLC-MS/MS, a NOX peak was observed at a retention time of 3.84 minutes that was easily discernable from the nicotine substrate peak at 1.40 minutes (Fig. 3A). Representative kinetics plots of NOX formation for individual FMO enzymes are shown in Fig. 3B–F. Enzyme kinetic analysis suggested that both FMOs 1 and 3 exhibited the highest overall activity for nicotine (Vmax/Km = 16 nL/min/mg for both FMO1 and FMO3; Table 2). The Km values of 0.75 and 0.70 mmol/L for FMO1 and FMO3, respectively, were similar to that observed for human liver microsomes (Km = 0.54; Table 2; Fig. 3G). While full-length FMO2 exhibited less overall activity (Vmax/Km = 1.8 nL/min/mg) than either FMO1 or FMO3, the activities of FMO4 and FMO5 were very low, with Km values > 5 mmol/L (Table 2; Fig. 3E and F, respectively).
Effects of FMO1 and FMO3 nonsynonymous polymorphisms on N′-oxidation activity
Prevalent missense polymorphisms for the three most active FMO enzymes (FMO1, FMO2, and FMO3) were identified using the NCBI SNP database. For this study, only nonsynonymous polymorphisms with a MAF > 1% in the 1000 Genomes Project were used (41). Nonsynonymous SNPs were identified for FMO1 (I303V, rs16864314) and FMO3 (N61S, rs72549322; D132H, rs12072582; E158K, rs2266782; V257M, rs1736557; and V277A, rs2066530; Supplementary Table S2). An E158K/E308G haplotype allele with a MAF ≥ 1% was also identified for FMO3. Other than the truncated FMO2472Stop variant, no independent missense polymorphisms with a MAF ≥ 1% were identified for this gene. Variants were cloned and overexpressed in the HEK293 cell line as described for wild-type FMOs.
The FMO1303V-variant exhibited a significantly (P = 0.023) higher Km (1.8 mmol/L) and a significant (P = 0.023) 3.5-fold lower overall catalytic activity (Vmax/Km = 4.5 nL/min/mg) as compared with wild-type FMO1 (Km = 0.8 mmol/L; Vmax/Km = 16 nL/min/mg; Table 2). The FMO2 protein encoded by the common FMO2Q472stop-variant exhibited no activity against nicotine or benzydamine in vitro using up to 5 mmol/L benzydamine. The FMO361S, FMO3132H, FMO3158K, and FMO3257M variants all exhibited significantly (P = 0.006, 0.003, 0.002, and 0.003, respectively) lower overall catalytic efficiencies as determined by Vmax/Km (ranging from 3.9 to 8.9 nL/min/mg) as compared with its wild-type FMO3 counterpart (Vmax/Km = 16 nL/min/mg). While the Km was similar for the wild-type FMO3 and FMO3277A isoforms (0.70 and 0.74 mmol/L, respectively), the Vmax for the FMO3277A variant (6.7 pmol/min/mg) was 40% (P = 0.029) lower than that observed for wild-type FMO3 (11 pmol/min/mg). The Km values for all of the other FMO3-variants were significantly (P < 0.05) higher than that observed for wild-type FMO3, with the FMO361S variant exhibiting a Km that was 6.7-fold higher (4.7 mmol/L). While microsomes containing the FMO3158K/308G haplotype variant exhibited a Km that was 3.1-fold higher than the wild-type isoform (2.2 mmol/L vs. 0.70 mmol/L, respectively), there was no significant (P = 0.565) difference in overall catalytic activity as determined by Vmax/Km (15 vs. 16 nL/min/mg, respectively).
Previous studies suggested that the FMO3158K/308G allele was associated with decreased FMO3 protein expression (42). In the present analysis, FMO3 mRNA expression was determined in the HEK293 cell lines overexpressing wild-type FMO3 (FMO3158E/308E), or FMO3158K/308G variants. Significant (P < 0.05) decreases in both mRNA (2.9- and 1.4-fold) and protein (3.6- and 2.1-fold) expression were observed for the FMO3158K/308E- and FMO3158E/308G-overexpressing HEK293 cell lines, respectively, when compared with cells overexpressing wild-type FMO3 (Fig. 4). In contrast, while a 4.5-fold decrease in mRNA expression was observed for the FMO3158K/308G variant versus wild-type FMO3, a >21-fold lower level of FMO3 protein was observed for cells overexpressing this variant as compared with cells overexpressing wild-type FMO3.
Discussion
The levels of urinary NOX as a percentage of total nicotine metabolites for 106 Chinese smokers showed high variability between subjects. When samples were stratified by CYP2A6 activity (determined by the urinary NMR), the levels of urinary NOX were significantly different between high versus low CYP2A6 activity groups. NOX formation comprised ≥19% of total nicotine metabolites in subjects exhibiting the lowest CYP2A6 activities. These data suggest that nicotine N′-oxidation is a key metabolism pathway for smokers, particularly those individuals whose CYP2A6 activity is diminished or deleted. For example, the NOX/TNE ratio was 18.8% in smokers with the lowest 5% NMR group whereas the corresponding NOX/TNE ratio was 5.5% in smokers with the highest 5% NMR. These results are consistent with a previous report where NOX levels were as high as 31% of total nicotine metabolites in two subjects whose CYP2A6 gene was deleted (18). This is particularly important for Asian populations where the MAF of the CYP2A6 deletion allele is high (15%–20%; ref. 43). While the CYP2A6 deletion allele (CYP2A6 *4) is much less prevalent in other populations including Caucasians and African Americans (44), the prevalence of other alleles that code for CYP2A6 variants with decreased enzyme function is still high in these populations (45, 46). Therefore, the impact of FMO variants on nicotine metabolism needs to be examined in other racial groups.
In this study, the five human FMO enzymes were tested for N′-oxidation activity against nicotine in vitro. While all of the recombinant FMO isoforms mediated NOX formation, FMO1 and FMO3 exhibited the highest activity, with the Km for FMO1 (0.8 mmol/L) similar to that reported previously for this enzyme (1.2 mmol/L; ref. 19). While FMO3 exhibited a Km (0.8 mmol/L) that was similar to that observed for FMO1 in this study, it was significantly lower than that observed in previous studies where the rate of NOX formation was still linear at 2 mmol/L nicotine (19). These differences in kinetics for FMO3 is likely due to the fact that the insect baculosomal system was used as the FMO-overexpression system in previous studies, while the FMOs (including FMO3) were overexpressed in the human HEK293 cell line in this study.
FMOs exhibit tissue-specific expression. FMO1 is an extra-hepatic drug metabolizing enzyme expressed in kidney, intestine, and possibly the brain (47), whereas FMO3 is the FMO isoform most highly expressed in liver (12). Given that the hepatic enzymes FMO4 and FMO5 exhibited minimal levels of activity against nicotine in this study, the data presented in this study suggest that FMO3 may be the primary hepatic enzyme involved in NOX formation, with FMO1 potentially more important for extra-hepatic nicotine metabolism.
The FMO2 variant that expresses a functional full-length protein is observed primarily in people of African descendent ancestry (MAF = 0.133), with most of the remaining human population expressing the FMO2472stop allele which codes for a truncated nonfunctional enzyme (14, 48). This study is the first to demonstrate that the functional full-length FMO2 variant has activity against nicotine. While its Km (4.5 mmol/L) was about 5-to 6-fold higher than that observed for either FMO1 or FMO3, this enzyme has been reported to be primarily expressed in human lung, with FMO2 mRNA detected in human lung specimens from people of African descent (14). Given the low prevalence of the full-length functional FMO2 allele in the Asian population (0%; ref. 41), FMO2 is, however, unlikely to play a role in variability in NOX formation in the Asian populations examined in this study.
Given that FMOs 1, 2, and 3 were the main catalysts for nicotine-N′-oxidation in this study, the potential role of prevalent genetic variants (MAF ≥ 1%) in these enzymes on nicotine metabolism was investigated. The FMO1303V variant exhibited a significant (P = 0.02) decrease in activity when compared with its wild-type counterpart, showing an approximately 3-fold decrease in overall NOX formation as determined by Vmax/Km. In previous studies, Furnes and colleagues (49) tested this variant against several substrates and did not observe differences in sulfoxidation activity against methyl-p-tolyl sulfide, fenthion, and methamizole when compared with the wild-type FMO1 isoform. However, in the same study, N-oxygenation activity was decreased against imipramine, decreasing both its affinity and maximum velocity, supporting the theory that differences observed in catalytic activity could be substrate dependent. The FMO1303V variant is present in African Americans and in the Chinese population (MAF = 10% and 4%, respectively) but not in Caucasians (MAF < 0.0.01%). Because smoking-related GWAS studies have been performed in primarily Caucasian populations, the effect of this variant on nicotine dependence and/or cigarette consumption had not been well studied (19, 26, 28). Interestingly, several noncoding region variants in FMO1 have also been shown to be correlated with nicotine dependence in Caucasian populations. This includes several FMO1 promoter variants which were associated with differences in FMO1 expression and protein levels, subsequently resulting in differences in overall activity between individuals (19).
The FMO2472stop variant was tested using nicotine as a substrate, and as expected, exhibited no activity against nicotine. This finding agrees with previously reported data where the truncated protein is catalytically inactive against a variety of substrates (13, 14, 50).
The FMO361S, FMO3132H, FMO3158K, and FMO3257M variants showed statistically significant (P < 0.05) decreases in overall catalytic activity as compared with the wild-type FMO3 isoform, suggesting that individuals with any of these genetic variants could impair nicotine transformation to NOX in vivo. This decrease was as high as 4-fold for both the FMO361S and FMO3132H variants. The FMO361S variant results in an amidic to hydroxylic amino acid change, and, consistent with this study, previous reports on this variant showed abolished activity towards trymethylamine N-oxidation (51). The FMO3132H variant has one of the most drastic changes in terms of amino acid structure, changing from an acidic to a basic amino acid. Similar to the decreases in nicotine N′-oxidation activity observed in this study, the FMO3132H variant was associated with 30% and 60% decreases in activity against methimazole and trimethylamine, respectively, in previous studies (52).
In this study, the FMO3158K and FMO3257M variants exhibited 50%–60% decreases in overall activity when compared with wild-type FMO3. The FMO3158K variant was shown in different studies to exhibit similar decreases in catalytic activity towards a variety of other substrates (14, 52–54). This variant is widely distributed among different ethnicities (MAF: CEU, 41%; ASW, 33%; MXL, 43%; CDX, 13%; ref. 41), and consistent with the lower activity observed for this variant against nicotine in this study, this polymorphism has been associated with nicotine dependence and/or cigarette consumption in previous studies (26, 28, 29).
For the FMO3257M variant, an aliphatic amino acid is changed to a sulfur-containing amino acid. This polymorphism is widely distributed among different populations (MAF: CEU 7%, ASW, 9%; MXL, 5%; and CDX, 21%; ref. 41). Similar to that observed for the FMO3158K variant, a significant 2.6-fold increase in Km (but no significant alteration in Vmax) was observed as compared with wild-type FMO3, resulting in significant decreases in overall catalytic activity against nicotine. Interestingly, previous studies reported no significant difference in activity between the FMO3257M variant and wild-type FMO3 on trimethylamine-N-oxidation, suggesting substrate specificity for this SNP (51, 53).
The FMO3277A variant, which retains the aliphatic nature of the amino acid at position 277, exhibited a minimal change in overall activity against nicotine (Vmax = 6.7 vs. 11.2 pmol/min/mg, for FMO3277A and wild-type FMO3, respectively). This is similar to that observed in previous studies against other compounds including trimethylamine-N-oxidation (55).
The FMO3158K/308G haplotype variant occurs with high frequency among different populations (MAF: CEU, 22%; ASW, 9%; MXL, 5%; and CDX, 11%). The linkage disequilibrium for the two SNPs at codons 158 and 308 is 98% and the frequency of this haplotype is defined by the MAF of the E308G variant (41). Results from this study indicated no significant difference in overall activity of this haplotype variant against nicotine. However, in vivo studies have shown that individuals with this FMO3-variant exhibit impaired activity toward a number of FMO3 substrates including nicotine (28, 29, 56). A recent study showed that this haplotype was correlated with a 50% reduction in FMO3 protein levels in a panel of human adult liver samples, a pattern not observed for mRNA levels measured in the same samples (42). It was suggested that overall protein instability, resulting in an early ubiquitin–proteasome pathway, was potentially responsible for the in vivo reduced enzymatic activity of this haplotype. Consistent with these results, the HEK293 cell line overexpressing the FMO3158K/308G variant exhibited decreases in FMO3 protein levels that were >4.5-fold less than the decreases observed for FMO3 mRNA as compared with HEK293 cells that overexpressed wild-type FMO3. This differential pattern of expression was not observed when examining the effects on FMO3 protein and mRNA expression of cells overexpressing each of the corresponding missense SNPs at codons 158 and 308 individually and suggest a direct effect on protein stability by the FMO3158K/308G haplotype variant.
While not observed as an independent SNP in the population, the FMO3308G variant has been previously tested against a variety of substrates in vitro, showing a trend of decreased activity (53, 57). This is consistent with the results obtained in this study where this SNP is associated with decreased in vitro NOX formation. In studies implicating the FMO3158K/308G haplotype, due to the high linkage disequilibrium between the FMO3 codons 158 and 308 SNPs, the FMO3308G variant has been significantly correlated with nicotine dependence (time to first cigarette in the morning; ref. 29). In the same study, this variant was associated with aberrant mRNA splicing in the liver and brain, suggesting that this SNP could also potentially be affecting overall FMO protein levels by multiple mechanisms.
This study is, to our knowledge, the most comprehensive study examining nicotine-N′-oxidation by the FMO family of enzymes. Nicotine-N′-oxidation is one of the direct detoxification pathways for nicotine, and while it accounts for a mean of 4%–7% of total urinary nicotine metabolites across multiple studies (16, 17), data from this study demonstrate that NOX formation may be a critical pathway for nicotine elimination in people with diminished or deficient CYP2A6 activity. Furthermore, results from this study suggest that functional SNPs in FMOs 1, 2, and 3 may play an important role in nicotine metabolism variability in people with diminished or deficient CYP2A6 activity. Studies in larger populations will be necessary to better establish the role of FMO genetic variants in overall nicotine metabolism.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: Y.X. Perez-Paramo, G. Chen, P. Lazarus
Development of methodology: Y.X. Perez-Paramo, G. Chen, C.J.W. Watson, P. Lazarus
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): Y.X. Perez-Paramo, Y.-T. Gao, W.-P. Koh, J.-M. Yuan, P. Lazarus
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): Y.X. Perez-Paramo, J.H. Ashmore, S. Nasrin, R. Wang, J.-M. Yuan, P. Lazarus
Writing, review, and/or revision of the manuscript: Y.X. Perez-Paramo, G. Chen, C.J.W. Watson, R. Wang, Y.-T. Gao, W.-P. Koh, J.-M. Yuan, P. Lazarus
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): Y.X. Perez-Paramo, J. Adams-Haduch, Y.-T. Gao, P. Lazarus
Study supervision: G. Chen, Y.-T. Gao, P. Lazarus
Acknowledgments
The authors thank the Mass Spectrometry Core facility at Washington State University, Spokane, for their help with LC/MS. This work is supported by grants from NIH, National Institutes of Environmental Health Sciences (grant R01-ES025460; to P. Lazarus), the Fulbright-Garcia Robles Program (to Y.X. Perez-Paramo), Health Sciences and Services Authority of Spokane, Washington (grant WSU002292 to the College of Pharmacy and Pharmaceutical Sciences, Washington State University), and NIH, NCI (grants R01-CA1440324, UM1-CA182876; to J.-M. Yuan).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.