Abstract Background Endometrial cancer (EMC) is the most common female genital tract malignancy with an increasing prevalence in many countries including Japan, a fact that renders early detection and treatment necessary to protect health and fertility. Although early detection and treatment are necessary to further improve the prognosis of women with endometrial cancer, biomarkers that accurately reflect the pathophysiology of EMC patients are still unclear. Therefore, it is clinically critical to identify biomarkers to assess diagnosis and treatment efficacy to facilitate appropriate treatment and development of new therapies for EMC. Methods In this study, wide-targeted plasma metabolome analysis was performed to identify biomarkers for EMC diagnosis and the prediction of treatment responses. The absolute quantification of 628 metabolites in plasma samples from 142 patients with EMC was performed using ultra-high-performance liquid chromatography with tandem mass spectrometry. Results The concentrations of 111 metabolites increased significantly, while the concentrations of 148 metabolites decreased significantly in patients with EMC compared to healthy controls. Specifically, LysoPC and TGs, including unsaturated fatty acids, were reduced in patients with stage IA EMC compared to healthy controls, indicating that these metabolic profiles could be used as early diagnostic markers of EMC. In contrast, blood levels of amino acids such as histidine and tryptophan decreased as the risk of recurrence increased and the stages of EMC advanced. Furthermore, a marked increase in total TG and a decrease in specific TGs and free fatty acids including polyunsaturated fatty acids levels were observed in patients with EMC. These results suggest that the polyunsaturated fatty acids in patients with EMC are crucial for disease progression. Conclusions Our data identified specific metabolite profiles that reflect the pathogenesis of EMC and showed that these metabolites correlate with the risk of recurrence and disease stage. Analysis of changes in plasma metabolite profiles could be applied for the early diagnosis and monitoring of the course of treatment of EMC patients.

AIM: Uterine cervical cancer (UCC) is the fourth most common cancer in women, responsible for more than 300 000 deaths worldwide. Its early detection, by cervical cytology, and prevention, by vaccinating against human papilloma virus, greatly contribute to reducing cervical cancer mortality in women. However, penetration of the effective prevention of UCC in Japan remains low. Plasma metabolome analysis is widely used for biomarker discovery and the identification of cancer-specific metabolic pathways. Here, we aimed to identify predictive biomarkers for the diagnosis and radiation sensitivity of UCC using wide-targeted plasma metabolomics. METHODS: We analyzed 628 metabolites in plasma samples obtained from 45 patients with UCC using ultra-high-performance liquid chromatography with tandem mass spectrometry. RESULTS: The levels of 47 metabolites were significantly increased and those of 75 metabolites were significantly decreased in patients with UCC relative to healthy controls. Increased levels of arginine and ceramides, and decreased levels of tryptophan, ornithine, glycosylceramides, lysophosphatidylcholine, and phosphatidylcholine were characteristic of patients with UCC. Comparison of metabolite profiles in groups susceptible and non-susceptible to radiation therapy, a treatment for UCC, revealed marked variations in polyunsaturated fatty acid, nucleic acid, and arginine metabolism in the group not susceptible to treatment. CONCLUSIONS: Our findings suggest that the metabolite profile of patients with UCC may be an important indicator for distinguishing these patients from healthy cohorts, and may also be useful for predicting sensitivity to radiotherapy.

The drug 5-fluorouracil (5-FU) is the first-choice chemotherapeutic agent against advanced-stage cancers. However, 10-30% of treated patients experience grade 3-4 toxicity. The deficiency of dihydropyrimidinase (DHPase), which catalyzes the second step of the 5-FU degradation pathway, is correlated with the risk of developing toxicity. Thus, genetic polymorphisms within DPYS, the DHPase-encoding gene, could potentially serve as predictors of severe 5-FU-related toxicity. We identified 12 novel DPYS variants in 3,554 Japanese individuals, but the effects of these mutations on function remain unknown. In the current study, we performed in vitro enzymatic analyses of the 12 newly identified DHPase variants. Dihydrouracil or dihydro-5-FU hydrolytic ring-opening kinetic parameters, Km and Vmax , and intrinsic clearance (CLint = Vmax /Km ) of the wild-type DHPase and eight variants were measured. Five of these variants (R118Q, H295R, T418I, Y448H, and T513A) showed significantly reduced CLint compared with that in the wild-type. The parameters for the remaining four variants (V59F, D81H, T136M, and R490H) could not be determined as dihydrouracil and dihydro-5-FU hydrolytic ring-opening activity was undetectable. We also determined DHPase variant protein stability using cycloheximide and bortezomib. The mechanism underlying the observed changes in the kinetic parameters was clarified using blue-native polyacrylamide gel electrophoresis and three-dimensional structural modeling. The results suggested that the decrease or loss of DHPase enzymatic activity was due to reduced stability and oligomerization of DHPase variant proteins. Our findings support the use of DPYS polymorphisms as novel pharmacogenomic markers for predicting severe 5-FU-related toxicity in the Japanese population. Significance Statement DHPase contributes to the degradation of 5-fluorouracil, and genetic polymorphisms that cause decreased activity of DHPase can cause severe toxicity. In this study, we performed functional analysis of 12 DHPase variants in the Japanese population and identified 9 genetic polymorphisms that cause reduced DHPase function. In addition, we found that the ability to oligomerize and the conformation of the active site are important for the enzymatic activity of DHPase.

Dihydropyrimidine dehydrogenase (DPD), encoded by the DPYD gene, is the rate-limiting enzyme in 5-fluorouracil (5-FU) degradation. In Caucasians, four DPYD risk variants are recognized to be responsible for interindividual variations in the development of 5-FU toxicity. However, these risk variants have not been identified in Asian populations. Recently, 41 DPYD allelic variants, including 15 novel single nucleotide variants, were identified in 3,554 Japanese individuals by analyzing their whole-genome sequences; however, the effects of these variants on DPD enzymatic activity remain unknown. In the present study, an in vitro analysis was performed on 41 DPD allelic variants and three DPD risk variants to elucidate the changes in enzymatic activity. Wild-type and 44 DPD-variant proteins were heterologously expressed in 293FT cells. DPD expression levels and dimerization of DPD were determined by immunoblotting after SDS-PAGE and blue native PAGE, respectively. The enzymatic activity of DPD was evaluated by quantification of dihydro-5-FU, a metabolite of 5-FU, using high-performance liquid chromatography-tandem mass spectrometry. Moreover, we used 3D simulation modeling to analyze the effect of amino acid substitutions on the conformation of DPD. Among the 41 DPD variants, seven exhibited drastically decreased intrinsic clearance (CL int ) compared to the wild-type protein. Moreover, R353C and G926V exhibited no enzymatic activity, and the band patterns observed in the immunoblots after blue native PAGE indicated that DPD dimerization is required for its enzymatic activity. Our data suggest that these variants may contribute to the significant inter-individual variability observed in the pharmacokinetics and pharmacodynamics of 5-FU. In our study, nine DPD variants exhibited drastically decreased or no enzymatic activity due to dimerization inhibition or conformational changes in each domain. Especially, the rare DPYD variants, although at very low frequencies, may serve as important pharmacogenomic markers associated with the severe 5-FU toxicity in Japanese population.

Epithelial ovarian cancer (EOC) is a fatal gynecologic cancer, and its poor prognosis is mainly due to delayed diagnosis. Therefore, biomarker identification and prognosis prediction are crucial in EOC. Altered cell metabolism is a characteristic feature of cancers, and metabolomics reflects an individual's current phenotype. In particular, plasma metabolome analyses can be useful for biomarker identification. In this study, we analyzed 624 metabolites, including uremic toxins (UTx) in plasma derived from 80 patients with EOC using ultra-high-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS). Compared with the healthy control, we detected 77 significantly increased metabolites and 114 significantly decreased metabolites in EOC patients. Especially, decreased concentrations of lysophosphatidylcholines and phosphatidylcholines and increased concentrations of triglycerides were observed, indicating a metabolic profile characteristic of EOC patients. After calculating the parameters of each metabolic index, we found that higher ratios of kynurenine to tryptophan correlates with worse prognosis in EOC patients. Kynurenine, one of the UTx, can affect the prognosis of EOC. Our results demonstrated that plasma metabolome analysis is useful not only for the diagnosis of EOC, but also for predicting prognosis with the variation of UTx and evaluating response to chemotherapy.

Fluoropyrimidine drugs (FPs), including 5-fluorouracil, tegafur, capecitabine, and doxifluridine, are among the most widely used anticancer agents in the treatment of solid tumors. However, severe toxicity occurs in approximately 30% of patients following FP administration, emphasizing the importance of predicting the risk of acute toxicity before treatment. Three metabolic enzymes, dihydropyrimidine dehydrogenase (DPD), dihydropyrimidinase (DHP), and β-ureidopropionase (β-UP), degrade FPs; hence, deficiencies in these enzymes, arising from genetic polymorphisms, are involved in severe FP-related toxicity, although the effect of these polymorphisms on in vivo enzymatic activity has not been clarified. Furthermore, the clinical usefulness of current methods for predicting in vivo activity, such as pyrimidine concentrations in blood or urine, is unknown. In vitro tests have been established as advantageous for predicting the in vivo activity of enzyme variants. This is due to several studies that evaluated FP activities after enzyme metabolism using transient expression systems in Escherichia coli or mammalian cells; however, there are no comparative reports of these results. Thus, in this review, we summarized the results of in vitro analyses involving DPD, DHP, and β-UP in an attempt to encourage further comparative studies using these drug types and to aid in the elucidation of their underlying mechanisms.

Dihydropyrimidine dehydrogenase (DPD, EC 1.3.1.2), encoded by the DPYD gene, is the rate-limiting enzyme in the degradation pathway of endogenous pyrimidine and fluoropyrimidine drugs such as 5-fluorouracil (5-FU). DPD catalyzes the reduction of uracil, thymine, and 5-FU. In Caucasians, DPYD mutations, including DPYD*2A, DPYD*13, c.2846A>T, and c.1129-5923C>G/hapB3, are known to contribute to interindividual variations in the toxicity of 5-FU; however, none of these DPYD polymorphisms has been identified in the Asian population. Recently, 21 DPYD allelic variants, including some novel single-nucleotide variants (SNVs), were identified in 1070 healthy Japanese individuals by analyzing their whole-genome sequences (WGSs), but the functional alterations caused by these variants remain unknown. In this study, in vitro analysis was performed on 22 DPD allelic variants by transiently expressing wild-type DPD and 21 DPD variants in 293FT cells and characterizing their enzymatic activities using 5-FU as a substrate. DPD expression levels and dimeric forms were determined using immunoblotting and blue-native PAGE, respectively. Additionally, the values of three kinetic parameters-the Michaelis constant (Km ), maximum velocity (Vmax ), and intrinsic clearance (CLint = Vmax/Km )-were determined for the reduction of 5-FU. Eleven variants exhibited significantly decreased intrinsic clearance compared with wild-type DPD. Moreover, the band patterns observed in the immunoblots of blue-native gels indicated that DPD dimerization is required for enzymatic activity in DPD. Thus, the detection of rare DPYD variants might facilitate severe adverse effect prediction of 5-FU-based chemotherapy in the Japanese population.

Sorghum, the fifth major global cereal, has potential as a source crop in temperate regions. To completely use sorghum bagasse, the ideal enzyme cocktail aims to identify and select the contributed enzymatic system. This study investigated the enzymatic system of Clostridium cellulovorans cellulosome and noncellulosomal enzymes using sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and liquid chromatography–tandem mass spectrometry LC-MS/MS. Enzyme solutions from treated and untreated sorghum bagasse were prepared and compared based on carboxymethyl cellulase (CMCase) activity. As a result, the enzyme solution derived from untreated sorghum bagasse had the highest activity. Protein bands from each C. cellulovorans culture showed distinct patterns on SDS-PAGE examination: three enzyme fractions, including culture supernatants, crystalline cellulose (Avicel) bound, and unbound fractions. These results suggested that untreated sorghum bagasse induced a variety of cellulosomal and uncellulosomal proteins. On the other hand, 5% or 10% sorghum supernatants could not induce Avicel-bound proteins, including the cellulosome, although even 5% sorghum juice induced three major bands: 180 kilodalton (kDa), 100 kDa, and 70 kDa, respectively. In contrast, cellobiose induced three major bands, while the total number of all isolated proteins from the cellobiose medium was the most limited among all culture media. More intriguingly, our investigation detected one cellulosomal protein, hydrophobic protein A (HbpA) and three noncellulosomal enzymes, indicating that glycosyl hydrolase family 130 (GH130) was identified as a biomass-induced enzyme in good accord with previously published proteomic studies. Therefore, the proteomic dataset generated in this study provides us a foundation for future computational approaches, including machine learning-based prediction of optimal enzyme cocktails for target biomass degradation.

BACKGROUND: No previous studies have explored metabolites associated with both genetic predispositions to type 2 diabetes (T2DM) and T2DM onset. Therefore, we aimed to explore metabolic profiles using genetic risk to identify pathways for tailored T2DM prevention strategies. METHODS: This prospective community-based cohort study in Japan included a total of 12,461 participants aged ≥ 20 years. Genetic predictors were genome-wide and pathway-based polygenic risk scores (PRSs). We quantified 43 metabolites using nuclear magnetic resonance spectroscopy. T2DM was defined as a non-fasting glucose level of ≥ 200 mg/dL, glycated hemoglobin level ≥ 6.5%, or self-reported T2DM treatment. A modified Poisson regression model was used to examine the associations of PRSs and metabolites with T2DM incidence. Linear associations were tested using the restricted cubic spline, and mediation analysis was conducted to assess the mediating effect of metabolites on the association between PRSs and T2DM incidence. RESULTS: During the 4.3-year median follow-up period, 354 T2DM cases were identified. A higher PRS was associated with incident T2DM (relative risk [RR] 2.09, 95% confidence interval [CI], 1.68-2.60, P < 0.001, 1-standard deviation [SD] increment). The nitrogen compound transport pathway PRS was associated with incident T2DM (RR 1.32, 95% CI 1.03-1.70, P < 0.001, 1-SD increment). Fourteen metabolites like glucose, 2-ketoisocaproic acid, glutamic acid, leucine, 2-aminobutyric acid, 2-hydroxybutyric acid, valine, 3-methyl-2-oxobutyric acid, alanine, 3-hydroxybutyric acid, 3-methyl-2-oxiovaleric acid, formic acid, arginine, and tyrosine were positively associated with the risk of T2DM. Only glycine was inversely associated with the risk of T2DM. Among 43 metabolites, 14 metabolites were positively associated with PRS (P for linear trend < 0.05). 3-hydroxyisobutyric-acid, 2-Aminobutyric acid, 2-ketoisocaproic acid, 2-hydroxybutyric acid, leucine, glycine, and glucose mediated the association between PRS and incident T2DM. CONCLUSIONS: Several metabolites were found to mediate the association between PRS and incident T2DM. These findings may contribute to a better understanding of the metabolic pathways involved in genetic susceptibility to T2DM.

Cancer is a systemic disease rather than a localized pathology and is characterized by widespread effects, including whole-body exhaustion and chronic inflammation. A thorough understanding of cancer pathophysiology requires a systemic approach that accounts for the complex interactions between cancer cells and host tissues. To explore these dynamics, we employed a comprehensive metabolomic analysis of plasma samples from patients with either esophageal or head and neck squamous cell carcinoma (SCC). Plasma samples from 149 patients were metabolically profiled and correlated with clinical data. Among the metabolites identified, lysophosphatidylcholine (LPC) emerged as the sole biomarker strongly correlated with prognosis. A significant reduction in plasma LPC levels was linked to poorer overall survival. Plasma LPC levels demonstrated minimal correlation with patient-specific factors, such as tumor size and general condition, but showed significant association with the response to immune checkpoint inhibitor therapy. Proteomic and cytokine analyses revealed that low plasma LPC levels reflected systemic chronic inflammation, characterized by high levels of inflammatory proteins, the cytokines interleukin-6 and tumor necrosis factor-α, and coagulation-related proteins. These findings indicate that plasma LPC levels may be used as reliable biomarkers for predicting prognosis and evaluating the efficacy of immunotherapy in patients with SCC.

UDP-glucuronosyltransferase (UGT) is an important drug-metabolizing enzyme involved in the detoxification of hydrophobic chemicals by conjugating them with hydrophilic glucuronic acid. Previous studies have revealed that UGTs form oligomers between or among the same or other isoforms, but the molecular mechanism underlying this formation remains unclear. In this study, we used optimized electrophoretic techniques to analyze UGT2B7 homo-oligomer formation. UGT2B7 was expressed in COS-1 cells, and lysates were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). By omitting the heating and reducing steps during SDS-PAGE sample preparation, bands corresponding to dimer, tetramer, and higher-order oligomers were detected in addition to monomeric bands. Since these SDS-stable UGT2B7 oligomer bands disappeared with the addition of reducing agents, we hypothesized that intermolecular disulfide bonds are involved in the formation of UGT oligomers. The cysteine residues important for this oligomer formation were investigated. Analyses using alanine substitution and deletion mutants suggested that three cysteine residues of UGT2B7, Cys127, Cys156, and Cys282, are important not only for oligomer formation but also for glucuronidation ability. We further investigated the oligomerization of UGT2B7 in intact living cells using two membrane-permeable cross-linkers, disuccinimidyl suberate and dithiobis(succinimidyl propionate). UGT2B7-expressing cells were treated with these reagents and analyzed by western blot. This cross-linker treatment markedly reduced the UGT2B7 monomer band and increased the formation of higher-molecular-mass species. These results indicated that the majority of UGT2B7 is present within cells as oligomers, maintaining its enzymatic function, rather than as a monomer.

Advancements in large-scale analysis of metabolites in human peripheral blood samples revealed the links between metabolite concentrations and genetic variations. This field is known as metabolome-genome-wide association study (MGWAS). Although MGWAS is a powerful tool, it has some limitations, particularly in terms of the number of metabolites that can be measured. Whether the observed associations are directly due to genetic variation or indirectly due to changes in unmeasured metabolites is unclear. To address this, we used simulations of metabolic pathway models to investigate the influence of genetic variants on metabolite concentrations and enhance the interpretation of MGWAS results. By systematically adjusting the enzyme reaction rates to simulate genetic variants, we observed changes in the metabolite levels. Our simulations accurately represented most of the variant-metabolite pairs identified by MGWAS with significant p-values, thereby demonstrating the potential of our approach. Furthermore, our simulations revealed additional marked fluctuations in metabolite levels that the MGWAS did not detect, suggesting that some variant-metabolite pairs might become more significant with larger sample sizes. We also categorized the enzymes into three types based on their impact on metabolite concentrations, highlighting enzymes with minimal impact. This indicated that genetic variations in these enzymes may have limited biological significance. Our study not only validates key MGWAS findings, but also provides a systematic framework for understanding enzyme-metabolite relationships. This approach offers valuable insights for future experimental studies and potential therapeutic interventions.

Ykt6 is a conserved SNARE protein involved in multiple membrane trafficking pathways, including intra-Golgi transport and autophagic membrane fusion. We previously demonstrated that mammalian Ykt6 is uniquely modified with farnesyl and geranylgeranyl groups at two C-terminal cysteines through the sequential action of farnesyltransferase (FT) and geranylgeranyltransferase type 3 (GGT3). Although these two cysteines are strictly conserved in all eukaryotes, the evolutionary conservation of Ykt6 double prenylation remains unclear, as budding yeast appears to lack the α subunit of GGT3. In this study, we used structural predictions to identify the uncharacterized protein Ecm9 as the functional α subunit of yeast GGT3. Ecm9 forms a complex with Bet2 and transfers a geranylgeranyl group to mono-farnesylated Ykt6. MALDI-TOF/TOF mass spectrometry confirmed that budding yeast Ykt6 is doubly prenylated with farnesyl and geranylgeranyl groups in wild-type cells but not in ecm9Δ cells. Loss of Ecm9 resulted in fragile cell walls, likely due to mislocalization of Golgi mannosyltransferases. Furthermore, ecm9Δ cells exhibited impaired Ykt6 localization to organelle membranes including autophagosomes, leading to reduced autophagic activity. These findings establish that double prenylation is an evolutionarily conserved structural feature of Ykt6 and is essential for its membrane localization and function.

Functionally impaired human flavin-containing monooxygenase 3 (FMO3) variants are associated with metabolic trimethylaminuria. The present study predicted the level of impairment of catalytic function of eight FMO3 variants newly detected in the updated Tohoku Medical Megabank. Catalytic function was modeled using the in vitro intrinsic trimethylamine N-oxygenation activities of 108 FMO3 variants using a logistic linear regression model. FMO3 proteins with single, double, triple, and quadruple amino acid variants underwent binary classification as unaffected (score 0) or impaired (≤70% of wild-type FMO3, score 1) based on their trimethylamine N-oxygenation activities. Using a model with the sum of the coefficients before and after amino acid substitutions, the mean probabilities of 62 FMO3 variants with impaired catalytic activities (0.78) were significantly higher than those of 46 with unaffected catalytic activities (0.30) (p < 0.01). The area under the corresponding receiver operating characteristic curve was 0.91. This in silico evaluation was applied to eight FMO3 variants (p.Met82Thr, p.Gly180Val, p.Ile212Leu, p.Asn275Ile, p.Gly276Glu, p.Ala377Asp, p.Tyr469Cys, and p.Val502Gly) newly found in the Tohoku Medical Megabank and validated in vitro. The present model could prove useful for estimating the activities of new FMO3 variants by applying a logistic linear regression model based on in vitro-generated catalytic data.

BACKGROUND: Lung cancer is the deadliest disease globally, with more than 120,000 diagnosed cases and more than 75,000 deaths annually in Japan. Several treatment options for advanced lung cancer are available, and the discovery of biomarkers will be useful for personalized medicine. Using metabolome analysis, we aimed to identify biomarkers for diagnosis and treatment response by examining the changes in metabolites associated with lung cancer progression. METHODS: Plasma samples from patients with recurrent or metastatic non-small cell lung carcinomas diagnosed at Tohoku University Hospital between 2019 and 2024 were used in this study. Metabolomic analysis was performed using the Biocrates Life Sciences MxP Quant 500 kit. Multivariate, principal component, and orthogonal partial least squares discriminant analyses were performed. RESULTS: The triglyceride and phosphatidylcholine concentrations were higher in the patients with early than in those with advanced lung adenocarcinomas. However, the cholesterol ester concentrations were higher for the patients with advanced lung cancer. The concentrations of hexosylceramide were higher in patients with early lung adenocarcinoma than in those with squamous cell carcinoma. Relative to epidermal growth factor receptor (EGFR)-mutation negative cases, the EGFR-mutation positive cases showed marked differences between the ceramide and triglyceride concentrations. For the best therapeutic effect of EGFR-TKI treatment, the hexosylceramide (HexCer) (d18:1/24:0), ceramide (Cer) (d18:2/22:0), and ceramide (Cer) (d18:2/24:0) concentrations were higher for the stable and progressive disease groups. The concentrations of phosphatidylcholine (PC) ae C42:2, sphingomyelin (SM) C24:1, and lysophosphatidylcholine (lysoPC) a C18:2 were higher in the partial response group treated with immune checkpoint inhibitors and chemotherapy. CONCLUSION: Metabolomic analysis may be useful for the diagnosis and treatment of lung cancer and may provide clues for new therapeutic strategies. PC ae C42:2, SM C24:1, and lysoPC a C18:2 can serve as predictive biomarkers for monitoring the therapeutic effects of the combination of immune checkpoint inhibitors and chemotherapy.

Cytochrome P450 2B6 (CYP2B6) catalyzes the metabolism of many drugs, including efavirenz and propofol. Genetic polymorphisms in CYP2B6 alter its enzymatic activity and substantially affect its pharmacokinetics. High-frequency variants, such as CYP2B6*6, are associated with the risk of developing side effects due to reduced CYP2B6 activity. However, the impact of rare alterations on enzyme function remains unknown, and some of these variants may significantly decrease the CYP2B6 activity. Therefore, in this study, we evaluated in vitro the functional alterations in 29 missense variants of the CYP2B6 gene identified in 8,380 Japanese individuals. Wild-type CYP2B6 and 29 rare CYP2B6 variants were transiently expressed in mammalian cells. The expression levels of variant CYP2B6 proteins in the microsomal fractions extracted from 293FT cells were assessed using western blotting and reduced-carbon monoxide difference spectroscopy, and a specific peak at 450 nm was detected in the wild-type and 19 variants. Furthermore, kinetic parameters were determined by assaying the reactions with efavirenz and propofol and quantifying the metabolite concentrations. We found that 12 variants had significantly lower or abolished enzymatic activity with both the substrates. In silico three-dimensional docking and molecular-dynamics simulations suggested that these functional changes were due to conformational changes in essential regions, such as the heme-binding site and ligand channels involved in transporting substrates to the active site. These findings have implications for predicting the plasma concentrations of CYP2B6 substrates and controlling their side effects.

Understanding the physiological changes associated with aging and the associated disease risks is essential to establish biomarkers as indicators of biological aging. This study used the NMR-measured plasma metabolome to calculate age-specific metabolite indices. In doing so, the scope of the study was deliberately simplified to capture general trends and insights into age-related changes in metabolic patterns. In addition, changes in metabolite concentrations with age were examined in detail, with the period from 55-59 to 60-64 years being a period of significant metabolic change, particularly in men, and from 45-49 to 50-54 years in females. These results illustrate the different variations in metabolite concentrations by sex and provide new insights into the relationship between age and metabolic diseases.

OBJECTIVE: The discovery of objective indicators for recent epileptic seizures will help confirm the diagnosis of epilepsy and evaluate therapeutic effects. Past studies had shortcomings such as the inclusion of patients under treatment and those with various etiologies that could confound the analysis results significantly. We aimed to minimize such confounding effects and to explore the small molecule biomarkers associated with the recent occurrence of epileptic seizures using urine metabolomics. METHODS: This is a multicenter prospective study. Subjects included pediatric patients aged 2 to 12 years old with new-onset, untreated epilepsy, who had had the last seizure within 1 month before urine collection. Controls included healthy children aged 2 to 12 years old. Those with underlying or chronic diseases, acute illnesses, or recent administration of medications or supplements were excluded. Targeted metabolome analysis of spot urine samples was conducted using gas chromatography (GC)- and liquid chromatography (LC)-tandem mass spectrometry (MS/MS). RESULTS: We enrolled 17 patients and 21 controls. Among 172 metabolites measured by GC/MS/MS and 41 metabolites measured by LC/MS/MS, only taurine was consistently reduced in the epilepsy group. This finding was subsequently confirmed by the absolute quantification of amino acids. No other metabolites were consistently altered between the two groups. CONCLUSIONS: Urine metabolome analysis, which covers a larger number of metabolites than conventional biochemistry analyses, found no consistently altered small molecule metabolites except for reduced taurine in epilepsy patients compared to healthy controls. Further studies with larger samples, subjects with different ages, expanded target metabolites, and the investigation of plasma samples are required.

Genetic variants of human flavin-containing monooxygenase 3 (FMO3) were investigated using an updated Japanese population panel containing 54,000 subjects (the previous panel contained 38,000 subjects). One stop codon mutation and six amino acid-substituted FMO3 variants were newly identified in the updated databank. Of these, two substituted variants (p.Thr329Ala and p.Arg492Trp) were previously identified in compound haplotypes with p.[(Glu158Lys; Glu308Gly)] and were associated with the metabolic disorder trimethylaminuria. Three recombinant FMO3 protein variants (p.Ser137Leu, p.Ala334Val, and p.Ile426Val) expressed in bacterial membranes had similar activities toward trimethylamine N-oxygenation (∼75-125 %) as wild-type FMO3 (117 min-1); however, the recombinant novel FMO3 variant Phe313Ile showed moderately decreased FMO3 catalytic activity (∼20 % of wild-type). Because of the known deleterious effects of FMO3 C-terminal stop codons, the novel truncated FMO3 Gly184Ter variant was suspected to be inactive. To easily identify the four impaired FMO3 variants (one stop codon mutation and three amino-acid substitutions) in the clinical setting, simple confirmation methods for these FMO3 variants are proposed using polymerase chain reaction/restriction fragment length polymorphism or allele-specific PCR methods. The updated whole-genome sequence data and kinetic analyses revealed that four of the seven single-nucleotide nonsense or missense FMO3 variants had moderately or severely impaired activity toward trimethylamine N-oxygenation.

Cytochrome P450 4F2 (CYP4F2) is an enzyme that is involved in the metabolism of arachidonic acid (AA), vitamin E and K (VK), and xenobiotics including drugs. CYP4F2*3 polymorphism (rs2108622; c.1297G>A; p.Val433Met) has been associated with hypertension, ischemic stroke, and variation in the effectiveness of the anticoagulant drug warfarin. In this study, we characterized wild-type CYP4F2 and 28 CYP4F2 variants, including a Val433Met substitution, detected in 8,380 Japanese subjects. The CYP4F2 variants were heterologously expressed in 293FT cells to measure the concentrations of CYP4F2 variant holoenzymes using carbon monoxide-reduced difference spectroscopy, where the wild type and 18 holoenzyme variants showed a peak at 450 nm. Kinetic parameters (Vmax , S50 , and CL int as Vmax /S50 ) of AA ω-hydroxylation were determined for the wild type and 21 variants with enzyme activity. Compared to the wild type, two variants showed significantly decreased CL int values for AA ω-hydroxylation. The values for seven variants could not be determined because no enzymatic activity was detected at the highest substrate concentration used. Three-dimensional structural modeling was performed to determine the reason for reduced enzymatic activity of the CYP4F2 variants. Our findings contribute to a better understanding of CYP4F2 variant-associated diseases and possible future therapeutic strategies. Significance Statement CYP4F2 is involved in the metabolism of AA and VK, and CYP4F2*3 polymorphisms have been associated with hypertension and variation in the effectiveness of the anticoagulant drug warfarin. In this study, we present a functional analysis of 28 CYP4F2 variants identified in Japanese subjects, demonstrating that seven gene polymorphisms cause loss of CYP4F2 function, and propose structural changes that lead to altered function.

BACKGROUND AIMS: Cancer cells reprogram their metabolic pathways to support bioenergetic and biosynthetic needs and to maintain their redox balance. In several human tumors the Keap1-Nrf2 system controls proliferation and metabolic reprogramming by regulating the pentose phosphate pathway (PPP). However, whether this metabolic reprogramming also occurs in normal proliferating cells is unclear. APPROACH AND RESULTS: To define the metabolic phenotype in normal proliferating hepatocytes, we induced cell proliferation in the liver by three distinct stimuli: liver regeneration by partial hepatectomy (PH) and hepatic hyperplasia induced by two direct mitogens, lead nitrate (LN) or triiodothyronine (T3). Following LN treatment, well-established features of cancer metabolic reprogramming including enhanced glycolysis, oxidative PPP, nucleic acid synthesis, NAD+/NADH synthesis and altered amino acid content as well as downregulated oxidative phosphorylation (OXPHOS) occurred in normal proliferating hepatocytes displaying Nrf2 activation. Genetic deletion of Nrf2 blunted LN-induced PPP activation and suppressed hepatocyte proliferation. Moreover, Nrf2 activation and following metabolic reprogramming did not occur when hepatocyte proliferation was induced by PH or T3. CONCLUSION: Many metabolic changes in cancer cells are shared by proliferating normal hepatocytes in response to a hostile environment. Nrf2 activation is essential for bridging metabolic changes with crucial components of cancer metabolic reprogramming including the activation of oxidative PPP. Our study demonstrates that matured hepatocytes exposed to LN undergo a cancer-like metabolic reprogramming and offers a rapid and useful in vivo model to study the molecular alterations underpinning the differences/similarities of metabolic changes in normal and neoplastic hepatocytes.

Single-nucleotide substitutions of human flavin-containing monooxygenase 3 (FMO3) identified in the whole-genome sequences of the updated Japanese population reference panel (now containing 38,000 subjects) were investigated. In this study, 2 stop codon mutations, 2 frameshifts, and 43 amino-acid-substituted FMO3 variants were identified. Among these 47 variants, 1 stop codon mutation, 1 frameshift, and 24 substituted variants were already recorded in the National Center for Biotechnology Information database. Functionally impaired FMO3 variants are known to be associated with the metabolic disorder trimethylaminuria; consequently, the enzymatic activities of the 43 substituted FMO3 variants were investigated. Twenty-seven recombinant FMO3 variants expressed in bacterial membranes had similar activities toward trimethylamine N-oxygenation (~75-125%) to that of wild-type FMO3 (98 min-1). However, 6 recombinant FMO3 variants (Arg51Gly, Val283Ala, Asp286His, Val382Ala, Arg387His, and Phe451Leu) had moderately decreased (≈50%) activities toward trimethylamine N-oxygenation, and 10 recombinant FMO3 variants (Gly11Asp, Gly39Val, Met66Lys, Asn80Lys, Val151Glu, Gly193Arg, Arg387Cys, Thr453Pro, Leu457Trp, and Met497Arg) showed severely decreased FMO3 catalytic activity (<10%). Because of the known deleterious effects of FMO3 C-terminal stop codons, the four truncated FMO3 variants (Val187SerfsTer25, Arg238Ter, Lys416SerfsTer72, and Gln427Ter) were suspected to be inactive with respect to trimethylamine N-oxygenation. The FMO3 p.Gly11Asp and p.Gly193Arg variants were located within the conserved sequences of flavin adenine dinucleotide (positions 9-14) and NADPH (positions 191-196) binding sites, which are important for FMO3 catalytic function. Whole-genome sequence data and kinetic analyses revealed that 20 of the 47 nonsense or missense FMO3 variants had moderately or severely impaired activity toward N-oxygenation of trimethylaminuria. Significance Statement The number of single-nucleotide substitutions in human FMO3 recorded in the expanded Japanese population reference panel database was updated. One novel stop mutation, FMO3 p.Gln427Ter; one novel frameshift (p.Lys416SerfsTer72); and 19 novel amino-acid-substituted FMO3 variants were identified, along with p.Arg238Ter, p.Val187SerfsTer25, and 24 amino-acid-substituted variants already recorded with individual rs numbers. Recombinant FMO3 Gly11Asp, Gly39Val, Met66Lys, Asn80Lys, Val151Glu, Gly193Arg, Arg387Cys, Thr453Pro, Leu457Trp, and Met497Arg variants showed severely decreased FMO3 catalytic activity, possibly associated with the metabolic disorder trimethylaminuria.

OBJECTIVE: The ketogenic diet (KD), a high-fat and low-carbohydrate diet, is effective for a subset of patients with drug-resistant epilepsy, although the mechanisms of the KD have not been fully elucidated. The aims of this observational study were to investigate comprehensive short-term metabolic changes induced by the KD and to explore candidate metabolites or pathways for potential new therapeutic targets. METHODS: Subjects included patients with intractable epilepsy who had undergone the KD therapy (the medium-chain triglyceride [MCT] KD or the modified Atkins diet using MCT oil). Plasma and urine samples were obtained before and at 2-4 weeks after initiation of the KD. Targeted metabolome analyses of these samples were performed using gas chromatography-tandem mass spectrometry (GC/MS/MS) and liquid chromatography-tandem mass spectrometry (LC/MS/MS). RESULTS: Samples from 10 and 11 patients were analysed using GC/MS/MS and LC/MS/MS, respectively. The KD increased ketone bodies, various fatty acids, lipids, and their conjugates. In addition, levels of metabolites located upstream of acetyl-CoA and propionyl-CoA, including catabolites of branched-chain amino acids and structural analogues of γ-aminobutyric acid and lactic acid, were elevated. CONCLUSIONS: The metabolites that were significantly changed after the initiation of the KD and related metabolites may be candidates for further studies for neuronal actions to develop new anti-seizure medications.

CYP3A4, which contributes to the metabolism of more than 30% of clinically used drugs, exhibits high variation in its activity; therefore, predicting CYP3A4 activity before drug treatment is vital for determining the optimal dosage for each patient. We aimed to develop and validate an LC-tandem mass spectrometry (LC-MS/MS) method that simultaneously measures the levels of CYP3A4 activity-related predictive biomarkers (6β-hydroxycortisol (6β-OHC), cortisol (C), 1β-hydroxydeoxycholic acid (1β-OHDCA), and deoxycholic acid (DCA)). Chromatographic separation was achieved using a YMC-Triart C18 column and a gradient flow of the mobile phase comprising deionized water/25% ammonia solution (100 : 0.1, v/v) and methanol/acetonitrile/25% ammonia solution (50 : 50 : 0.1, v/v/v). Selective reaction monitoring in the negative-ion mode was used for MS/MS, and run times of 33 min were used. All analytes showed high linearity in the range of 3-3000 ng/mL. Additionally, their concentrations in urine samples derived from volunteers were analyzed via treatment with deconjugation enzymes, ignoring inter-individual differences in the variation of other enzymatic activities. Our method satisfied the analytical validation criteria under clinical conditions. Moreover, the concentrations of each analyte were quantified within the range of calibration curves for all urine samples. The conjugated forms of each analyte were hydrolyzed to accurately examine CYP3A4 activity. Non-invasive urine sampling employed herein is an effective alternative to invasive plasma sampling. The analytically validated simultaneous quantification method developed in this study can be used to predict CYP3A4 activity in precision medicine and investigate the potential clinical applications of CYP3A4 biomarkers (6β-OHC/C and 1β-OHDCA/DCA ratios).

Oxidative stress is an essential component in the progression of diabetic kidney disease (DKD), and the transcription factor NF-E2-related factor-2 (Nrf2) plays critical roles in protecting the body against oxidative stress. To clarify the roles of Nrf2 in protecting against DKD, in this study we prepared compound mutant mice with diabetes and loss of antioxidative defense. Specifically, we prepared compound Ins2Akita/+ (Akita) and Nrf2 knockout (Akita::Nrf2-/-) or Akita and Nrf2 induction (Akita::Keap1FA/FA) mutant mice. Eighteen-week-old Akita::Nrf2-/- mice showed more severe diabetic symptoms than Akita mice. In the Akita::Nrf2-/- mouse kidneys, the glomeruli showed distended capillary loops, suggesting enhanced mesangiolysis. Distal tubules showed dilation and an increase in 8-hydroxydeoxyguanosine-positive staining. In the Akita::Nrf2-/- mouse kidneys, the expression of glutathione (GSH) synthesis-related genes was decreased, and the actual GSH level was decreased in matrix-assisted laser desorption/ionization mass spectrometry imaging analysis. Akita::Nrf2-/- mice exhibited severe inflammation and enhancement of infiltrated macrophages in the kidney. To further examine the progression of DKD, we compared forty-week-old Akita mouse kidney compounds with Nrf2-knockout or Nrf2 mildly induced (Akita::Keap1FA/FA) mice. Nrf2-knockout Akita (Akita::Nrf2-/-) mice displayed severe medullary cast formation, but the formation was ameliorated in Akita::Keap1FA/FA mice. Moreover, in Akita::Keap1FA/FA mice, tubule injury and inflammation-related gene expression were significantly suppressed, which was evident in Akita::Nrf2-/- mouse kidneys. These results demonstrate that Nrf2 contributes to the protection of the kidneys against DKD by suppressing oxidative stress and inflammation.

Dihydropyrimidinase (DHP) is an enzyme that catabolizes the degradation of pyrimidine and fluoropyrimidine drugs such as 5-fluorouracil. DHP deficiency triggers various clinical symptoms and increases the risk of fluoropyrimidine drug toxicity. Various pathogenic variants of DHP cause DHP deficiency, and their catalytic activities have been well studied. However, the three-dimensional structures of DHP variants have not been clarified. In this study, we investigated the effects of mutations on DHP structures using the molecular dynamics simulations. Simulations of the wild type and 10 variants were performed and compared. In the T68R, D81G, G278D, and L337P variants, the flexibilities of structures related to the interaction for oligomer formation increased in comparison with those of the wild type. W117R, T343A, and R412 M mutations affected the structures of stereochemistry gate loops or the substrate-binding pocket. The three-dimensional structures of W360R and G435R variants were suggested to collapse. On the other hand, only slight structural changes were observed in the R181W variant, whose experimentally observed activity was similar to that of the wild type. The computational results are expected to clarify the relationship between clinical mutations and structural effects of drug-metabolizing enzymes.

The number of single-nucleotide substitutions of human flavin-containing monooxygenase 3 (FMO3) recorded in mega-databases is increasing. Moreover, phenotype-gene analyses have revealed impaired FMO3 variants associated with the metabolic disorder trimethylaminuria. In this study, four novel amino-acid substituted FMO3 variants, namely p.(Gly191Asp), p.(Glu414Gln), p.(Phe510Ser), and p.(Val530CysfsTer1), were identified in the whole-genome sequences in the Japanese population reference panel (8.3K JPN) of the Tohoku Medical Megabank Organization. Additionally, four variants, namely p.(Ile369Thr), p.(Phe463Val), p.(Arg500Gln), and p.(Ala526Thr) FMO3, were found in the 8.3K JPN database but were already recorded in the National Center for Biotechnology Information database. Novel FMO3 variants p.[(Met1Leu)] and p.[(Trp231Ter)] were also identified in phenotype-gene analyses of 290 unrelated subjects with self-reported malodor. Among the eight recombinant FMO3 variants tested (except for p.[(Met1Leu)] and p.[(Trp231Ter)]), Arg500Gln and Gly191Asp FMO3, respectively, had lower and much lower capacities for trimethylamine and/or benzydamine N-oxygenation activities than wild-type FMO3. Because another FMO3 mutation p.[(Gly191Cys)] with diminished recombinant protein activity was previously detected in two independent probands, Gly191 would appear to be important for FMO3 catalytic function. Analysis of whole-genome sequence data and trimethylaminuria phenotypes revealed missense FMO3 variants that severely impaired FMO3-mediated N-oxygenations in Japanese subjects that could be susceptible to low drug clearances.

Both Pendred syndrome (PS) and nonsyndromic hearing loss with an enlarged vestibular aqueduct (EVA) are autosomal recessive disorders caused by SLC26A4 pathogenic variants. The spectrum of SLC26A4 pathogenic variants varies with the ethnic background. Among the patients with EVA in Okinawa, 94% had some combination of NM_000441.2(SLC26A4):c.1707+5G>A and NM_000441.2(SLC26A4):c.2168A>G(p.His723Arg), the two SLC26A4 pathogenic variants that are the most common in this population. We identified these two pathogenic variants using a novel genotyping method that employed an allele-specific polymerase chain reaction (PCR) from a gDNA and single-stranded tag hybridization chromatographic printed-array strip (STH-PAS) in DNA samples obtained from 48 samples in Okinawa, including 34 patients with EVA and 14 carriers of c.1707+5G>A or c.2168A>G. In addition, whole blood and saliva samples were used for analysis in this genotyping method with direct PCR. The results of STH-PAS genotyping were consistent with those obtained using standard Sanger sequencing for all samples. The accuracy of the STH-PAS method is 100% under the optimized conditions. STH-PAS genotyping provided a diagnosis in 30 out of 34 patients (88%) in Okinawan patients with EVA in under 3 h. The turn-around time for STH-PAS genotyping used with direct PCR was 2 h as a result of the omission of the DNA extraction and purification steps. Using information about the ethnic distribution of pathogenic variants in the SLC26A4 gene, STH-PAS genotyping performs a rapid genetic diagnosis that is simple and has a considerably improved efficiency.

In mass spectrometry-based metabolomics, the differences in the analytical results from different laboratories/machines are an issue to be considered because various types of machines are used in each laboratory. Moreover, the analytical methods are unique to each laboratory. It is important to understand the reality of inter-laboratory differences in metabolomics. Therefore, we have evaluated whether the differences in analytical methods, with the exception sample pretreatment and including metabolite extraction, are involved in the inter-laboratory differences or not. In this study, nine facilities are evaluated for inter-laboratory comparisons of metabolomic analysis. Identical dried samples prepared from human and mouse plasma are distributed to each laboratory, and the metabolites are measured without the pretreatment that is unique to each laboratory. In these measurements, hydrophilic and hydrophobic metabolites are analyzed using 11 and 7 analytical methods, respectively. The metabolomic data acquired at each laboratory are integrated, and the differences in the metabolomic data from the laboratories are evaluated. No substantial difference in the relative quantitative data (human/mouse) for a little less than 50% of the detected metabolites is observed, and the hydrophilic metabolites have fewer differences between the laboratories compared with hydrophobic metabolites. From evaluating selected quantitatively guaranteed metabolites, the proportion of metabolites without the inter-laboratory differences is observed to be slightly high. It is difficult to resolve the inter-laboratory differences in metabolomics because all laboratories cannot prepare the same analytical environments. However, the results from this study indicate that the inter-laboratory differences in metabolomic data are due to measurement and data analysis rather than sample preparation, which will facilitate the understanding of the problems in metabolomics studies involving multiple laboratories.

Inter-individual variability in pharmacokinetics and drug response is heavily influenced by single-nucleotide variants (SNVs) and copy-number variations (CNVs) in genes with importance for drug disposition. Nowadays, a plethora of studies implement next generation sequencing to capture rare and novel pharmacogenomic (PGx) variants that influence drug response. To address these issues, we present a comprehensive end-to-end analysis workflow, beginning from targeted PGx panel re-sequencing to in silico analysis pipelines and in vitro validation assays. Specifically, we show that novel pharmacogenetic missense variants that are predicted or putatively predicted to be functionally deleterious, significantly alter protein activity levels of CYP2D6 and CYP2C19 proteins. We further demonstrate that variant priorization pipelines tailored with functional in vitro validation assays provide supporting evidence for the deleterious effect of novel PGx variants. The proposed workflow could provide the basis for integrating next-generation sequencing for PGx testing into routine clinical practice.

Introduction: Pharmacogenomic (PGx) testing results provide valuable information on drug selection and appropriate dosing, maximization of efficacy, and minimization of adverse effects. Although the number of large-scale, next-generation-sequencing-based PGx studies has recently increased, little is known about the risks and benefits of returning PGx results to ostensibly healthy individuals in research settings. Methods: Single-nucleotide variants of three actionable PGx genes, namely, MT-RNR1, CYP2C19, and NUDT15, were returned to 161 participants in a population-based Tohoku Medical Megabank project. Informed consent was obtained from the participants after a seminar on the outline of this study. The results were sent by mail alongside sealed information letter intended for clinicians. As an exception, genetic counseling was performed for the MT-RNR1 m.1555A > G variant carriers by a medical geneticist, and consultation with an otolaryngologist was encouraged. Questionnaire surveys (QSs) were conducted five times to evaluate the participants' understanding of the topic, psychological impact, and attitude toward the study. Results: Whereas the majority of participants were unfamiliar with the term PGx, and none had undergone PGx testing before the study, more than 80% of the participants felt that they could acquire basic PGx knowledge sufficient to understand their genomic results and were satisfied with their potential benefit and use in future prescriptions. On the other hand, some felt that the PGx concepts or terminology was difficult to fully understand and suggested that in-person return of the results was desirable. Conclusions: These results collectively suggest possible benefits of returning preemptive PGx information to ostensibly healthy cohort participants in a research setting.

<title>Abstract</title>Space travel induces stresses that contribute to health problems, as well as inducing the expression of Nrf2 (NF-E2-related factor-2) target genes that mediate adaptive responses to oxidative and other stress responses. The volume of epididymal white adipose tissue (eWAT) in mice increases during spaceflight, a change that is attenuated by <italic>Nrf2</italic> knockout. We conducted metabolome analyses of plasma from wild-type and <italic>Nrf2</italic> knockout mice collected at pre-flight, in-flight and post-flight time points, as well as tissues collected post-flight to clarify the metabolic responses during and after spaceflight and the contribution of Nrf2 to these responses. Plasma glycerophospholipid and sphingolipid levels were elevated during spaceflight, whereas triacylglycerol levels were lower after spaceflight. In wild-type mouse eWAT, triacylglycerol levels were increased, but phosphatidylcholine levels were decreased, and these changes were attenuated in <italic>Nrf2</italic> knockout mice. Transcriptome analyses revealed marked changes in the expression of lipid-related genes in the liver and eWAT after spaceflight and the effects of <italic>Nrf2</italic> knockout on these changes. Based on these results, we concluded that space stress provokes significant responses in lipid metabolism during and after spaceflight; Nrf2 plays critical roles in these responses.

Metabolic profiling is an omics approach that can be used to observe phenotypic changes, making it particularly attractive for biomarker discovery. Although several candidate metabolites biomarkers for disease expression have been identified in recent clinical studies, the reference values of healthy subjects have not been established. In particular, the accuracy of concentrations measured by mass spectrometry (MS) is unclear. Therefore, comprehensive metabolic profiling in large-scale cohorts by MS to create a database with reference ranges is essential for evaluating the quality of the discovered biomarkers. In this study, we tested 8700 plasma samples by commercial kit-based metabolomics and separated them into two groups of 6159 and 2541 analyses based on the different ultra-high-performance tandem mass spectrometry (UHPLC-MS/MS) systems. We evaluated the quality of the quantified values of the detected metabolites from the reference materials in the group of 2541 compared with the quantified values from other platforms, such as nuclear magnetic resonance (NMR), supercritical fluid chromatography tandem mass spectrometry (SFC-MS/MS) and UHPLC-Fourier transform mass spectrometry (FTMS). The values of the amino acids were highly correlated with the NMR results, and lipid species such as phosphatidylcholines and ceramides showed good correlation, while the values of triglycerides and cholesterol esters correlated less to the lipidomics analyses performed using SFC-MS/MS and UHPLC-FTMS. The evaluation of the quantified values by MS-based techniques is essential for metabolic profiling in a large-scale cohort.

Cytochrome P450 1A2 (CYP1A2), which accounts for approximately 13% of the total hepatic cytochrome content, catalyzes the metabolic reactions of approximately 9% of frequently used drugs, including theophylline and olanzapine. Substantial inter-individual differences in enzymatic activity have been observed among patients, which could be caused by genetic polymorphisms. Therefore, we functionally characterized 21 novel CYP1A2 variants identified in 4773 Japanese individuals by determining the kinetic parameters of phenacetin O-deethylation. Our results showed that most of the evaluated variants exhibited decreased or no enzymatic activity, which may be attributed to potential structural alterations. Notably, the Leu98Gln, Gly233Arg, Ser380del Gly454Asp, and Arg457Trp variants did not exhibit quantifiable enzymatic activity. Additionally, three-dimensional (3D) docking analyses were performed to further understand the underlying mechanisms behind variant pharmacokinetics. Our data further suggest that despite mutations occurring on the protein surface, accumulating interactions could result in the impairment of protein function through the destabilization of binding regions and changes in protein folding. Therefore, our findings provide additional information regarding rare CYP1A2 genetic variants and how their underlying effects could clarify discrepancies noted in previous phenotypical studies. This would allow the improvement of personalized therapeutics and highlight the importance of identifying and characterizing rare variants.

The oxygenation of food-derived trimethylamine to its N-oxide is a representative reaction mediated by human flavin-containing monooxygenase 3 (FMO3). Impaired FMO3 enzymatic activity is associated with trimethylaminuria (accumulation of substrate), whereas trimethylamine N-oxide (metabolite) is associated with arteriosclerosis. We previously reported FMO3 single-nucleotide and/or haplotype variants with low FMO3 metabolic capacity using urinary phenotyping and the whole-genome sequencing of Japanese populations. Here, we further analyze Japanese volunteers with self-reported malodor and interrogate an updated Japanese database for novel FMO3 single-nucleotide and/or haplotype variants. After 3 years of follow up, seven probands were found to harbor the known impaired FMO3 variant p.(Gly191Cys) identified in the database or novel variants/haplotypes including p.(Met66Val), p.(Arg223Gln), p.(Glu158Lys;Glu308Gly;Arg492Trp), and p.(Glu158Lys;Glu308Gly;Pro496Ser). The known severe mutation p.(Cys197Ter) (a TG deletion) and four variants including p.(Tyr269His) and p.(Pro496Ser) were first detected in the updated genome panel. Among previously unanalyzed FMO3 variants, the trimethylamine/benzydamine N-oxygenation activities of recombinant p.(Met66Val), p.(Arg223Gln), p.(Tyr269His), p.(Glu158Lys;Glu308Gly;Arg492Trp), and p.(Glu158Lys;Glu308Gly;Pro496Ser) FMO3 variant proteins were severely decreased (Vmax/Km <10% of wild-type). Although the present novel mutations or alleles were relatively rare, both in self-reported Japanese trimethylaminuria sufferers and in the genomic database panel, three common FMO3 missense or deletion variants severely impaired FMO3-mediated N-oxygenation of trimethylamine.

Metabolomics has been widely used for investigating the biological functions of disease expression and has the potential to discover biomarkers in circulating biofluids or tissue extracts that reflect in phenotypic changes. Metabolic profiling has advantages because of the use of unbiased techniques, including multivariate analysis, and has been applied in pharmacological studies to predict therapeutic and adverse reactions of drugs, which is called pharmacometabolomics (PMx). Nuclear magnetic resonance (NMR)- and mass spectrometry (MS)-based metabolomics has contributed to the discovery of recent disease biomarkers; however, the optimal strategy for the study purpose must be selected from many established protocols, methodologies and analytical platforms. Additionally, information on molecular localization in tissue is essential for further functional analyses related to therapeutic and adverse effects of drugs in the process of drug development. MS imaging (MSI) is a promising technology that can visualize molecules on tissue surfaces without labeling and thus provide localized information. This review summarizes recent uses of MS-based global and wide-targeted metabolomics technologies and the advantages of the MSI approach for PMx and highlights the PMx technique for the biomarker discovery of adverse drug effects.

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 (K m , k cat , 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.

Cytochrome P450 2C9 (CYP2C9) is an important drug-metabolizing enzyme that contributes to the metabolism of approximately 15% of clinically used drugs, including warfarin, which is known for its narrow therapeutic window. Interindividual differences in CYP2C9 enzymatic activity caused by CYP2C9 genetic polymorphisms lead to inconsistent treatment responses in patients. Thus, in this study, we characterized the functional differences in CYP2C9 wild-type (CYP2C9.1), CYP2C9.2, CYP2C9.3, and 12 rare novel variants identified in 4773 Japanese individuals. These CYP2C9 variants were heterologously expressed in 293FT cells, and the kinetic parameters (Km, kcat, Vmax, catalytic efficiency, and CLint) of (S)-warfarin 7-hydroxylation and tolbutamide 4-hydroxylation were estimated. From this analysis, almost all novel CYP2C9 variants showed significantly reduced or null enzymatic activity compared with that of the CYP2C9 wild-type. A strong correlation was found in catalytic efficiencies between (S)-warfarin 7-hydroxylation and tolbutamide 4-hydroxylation among all studied CYP2C9 variants. The causes of the observed perturbation in enzyme activity were evaluated by three-dimensional structural modeling. Our findings could clarify a part of discrepancies among genotype-phenotype associations based on the novel CYP2C9 rare allelic variants and could, therefore, improve personalized medicine, including the selection of the appropriate warfarin dose.

Next-generation sequencing (NGS) has identified variations in cytochrome P450 (CYP) 2D6 associated with drug responses. However, determination of novel haplotypes is difficult because of the short reads generated by NGS. We aimed to identify novel CYP2D6 variants in the Japanese population and predict the CYP2D6 phenotype based on in vitro metabolic studies. Using a targeted NGS panel (PKSeq), 990 Japanese genomes were sequenced, and then novel CYP2D6 haplotypes were determined. Km, Vmax, and intrinsic clearance (Vmax/Km) of N-desmethyl-tamoxifen 4-hydroxylation were calculated by in vitro metabolic studies using cDNA-expressed CYP2D6 proteins. After determination of the CYP2D6 diplotypes, phenotypes of the individuals were predicted based on the in vitro metabolic activities. Targeted NGS identified 14 CYP2D6 variants not registered in the Pharmacogene Variation Consortium (PharmVar) database. Ten novel haplotypes were registered as CYP2D6*128 to *137 alleles in the PharmVar database. Based on the Vmax/Km value of each allele, *128, *129, *130, *131, *132, and *133 were predicted to be nonfunctional alleles. According to the results of the present study, six normal metabolizers (NM) and one intermediate (IM) metabolizers were designated as IM and poor metabolizers (PM), respectively. Our findings provide important insights into novel haplotypes and haplotypes of CYP2D6 and the effects on in vitro metabolic activities.

In the Tohoku Medical Megabank project, genome and omics analyses of participants in two cohort studies were performed. A part of the data is available at the Japanese Multi Omics Reference Panel (jMorp; https://jmorp.megabank.tohoku.ac.jp) as a web-based database, as reported in our previous manuscript published in Nucleic Acid Research in 2018. At that time, jMorp mainly consisted of metabolome data; however, now genome, methylome, and transcriptome data have been integrated in addition to the enhancement of the number of samples for the metabolome data. For genomic data, jMorp provides a Japanese reference sequence obtained using de novo assembly of sequences from three Japanese individuals and allele frequencies obtained using whole-genome sequencing of 8,380 Japanese individuals. In addition, the omics data include methylome and transcriptome data from ∼300 samples and distribution of concentrations of more than 755 metabolites obtained using high-throughput nuclear magnetic resonance and high-sensitivity mass spectrometry. In summary, jMorp now provides four different kinds of omics data (genome, methylome, transcriptome, and metabolome), with a user-friendly web interface. This will be a useful scientific data resource on the general population for the discovery of disease biomarkers and personalized disease prevention and early diagnosis.

Aims: Pathogenic variants in mitochondrial DNA are known to be associated with sensorineural hearing loss (SNHL) and aminoglycoside-induced HL. Among them, the m.1555A>G mutation is the most common. Thus, a rapid and easy companion diagnostic method for this mutation would be desirable to prevent HL caused by aminoglycoside therapy. In this study, we report an improved protocol for the single-stranded tag hybridization chromatographic printed-array strip (STH-PAS) method for identifying the m.1555A>G mutation. Methods: To evaluate the accuracy of a novel diagnostic for the m.1555A>G mutation we analyzed 378 DNA samples with or without the m.1555A>G mutation, as determined by Invader assay, and calculated the sensitivity, specificity, and false negative and false positive ratios of this new method. Results: The newly developed protocol was robust; we, obtained the same results using multiple DNA concentrations, differing annealing temperatures, and different polymerase chain reaction thermal cyclers. The diagnostic sensitivity based on the STH-PAS method was 0.99, and the specificity was 1.00. The false negative and false positive ratios were 0 and 0.01, respectively. Conclusion: We improved the genotyping method for m.1555A>G mutations. This assays will be useful as a rapid companion diagnostic before aminoglycoside use.

Activating mutations in KEAP1-NRF2 are frequently found in tumors of the lung, esophagus, and liver, where they are associated with aggressive growth, resistance to cancer therapies, and low overall survival. Despite the fact that NRF2 is a validated driver of tumorigenesis and chemotherapeutic resistance, there are currently no approved drugs which can inhibit its activity. Therefore, there is an urgent clinical need to identify NRF2-selective cancer therapies. To this end, we developed a novel synthetic lethal assay, based on fluorescently labeled isogenic wild-type and Keap1 knockout cell lines, in order to screen for compounds which selectively kill cells in an NRF2-dependent manner. Through this approach, we identified three compounds based on the geldanamycin scaffold which display synthetic lethality with NRF2. Mechanistically, we show that products of NRF2 target genes metabolize the quinone-containing geldanamycin compounds into more potent HSP90 inhibitors, which enhances their cytotoxicity while simultaneously restricting the synthetic lethal effect to cells with aberrant NRF2 activity. As all three of the geldanamycin-derived compounds have been used in clinical trials, they represent ideal candidates for drug repositioning to target the currently untreatable NRF2 activity in cancer.

The evaluation of Cytochrome P450 (CYP) enzymatic activity is essential to estimate drug pharmacokinetics. Numerous CYP allelic variants have been identified; the functional characterisation of these variants is required for their application in precision medicine. Results from heterologous expression systems using mammalian cells can be integrated in in vivo studies; however, other systems such as E. coli, bacteria, yeast, and baculoviruses are generally used owing to the difficulty in expressing high CYP levels in mammalian cells. Here, by optimising transfection and supplementing conditions, we developed a heterologous expression system using 293FT cells to evaluate the enzymatic activities of three CYP isoforms (CYP1A2, CYP2C9, and CYP3A4). Moreover, we established co-expression with cytochrome P450 oxidoreductase and cytochrome b5. This expression system would be a potential complementary or beneficial alternative approach for the pharmacokinetic evaluation of clinically used and developing drugs in vitro.

While CYP2D6 allele and phenotype frequencies have been extensively studied, currently, very little ethnically specific data is available regarding the East African and South Pacific region, including Kenya and Vanuatu. The absence of information regarding gene polymorphisms and their resulting clinical effects in these populations may hinder treatment strategies and patient outcome. Given the scarceness of CYP2D6 related data in these populations, the purpose of this study was to perform a pharmacogenomic analysis of the Kenyan and Ni-Vanuatu population and ultimately characterize the enzymatic properties of eight novel CYP2D6 variant proteins expressed in 293FT cells in vitro using dextromethorphan as a substrate. Our study revealed a prevalence of functional alleles in both populations a low frequency for decreased function defining genotypes in the Ni-Vanuatu population, with approximately 36% of our Kenyan subjects presenting substrate-dependent decreased function alleles. Additionally, 6 variants (P171L, G306R, V402L, K1, K2, and K3) showed significantly reduced intrinsic clearance compared to wild-type CYP2D6.1. Our findings aid in efforts to bridge the gap between pharmacogenomic analysis and clinical application, by providing useful information in the development of ethnic-specific strategies as well as stressing the importance of population-specific genotyping when conducting multi-regional clinical trials and designing therapeutic strategies.

Flavin-containing monooxygenase 3 (FMO3) is a polymorphic xenobiotic- and dietary compound-metabolizing enzyme associated with the genetic disorder trimethylaminuria. We phenotyped 428 Japanese subjects using traditional urinary phenotyping assays and identified two subjects with <20% FMO3 metabolic capacity. Both subjects had novel frameshift mutations. Proband 1 harbored a novel CC deletion resulting in p.[(Pro153Gln fs; Phe166Ter)] FMO3, which was in trans configuration with p.(Cys197Ter). Proband 2 harbored a novel T deletion resulting in p.[(Met211Arg fs; Val220Ter)] FMO3, which was in trans configuration with p.[(Val257Met; Met260Val)]. We also analyzed a new large Japanese database for novel single nucleotide substitutions of FMO3 and identified the following variants with very low frequencies (<∼0.1%): p.(Lys56Glu), p.(Ser112Asn), p.(Asn164Lys), p.(Gly191Cys), p.(Ile199Ser), p.(Pro248Thr), p.(Pro248Leu), p.(Asp286Tyr), and p.(Ala311Pro). Recombinant FMO3 proteins of the above and unanalyzed variants underwent kinetic analysis of their trimethylamine/benzydamine N-oxygenation activities. Gly191Cys, Ile199Ser, Asp286Tyr, and Ala311Pro variant FMO3 proteins exhibited severely decreased activities (Vmax/Km <5% of wild-type). Although these new variants were rare alleles in Japanese self-reported trimethylaminuria sufferers and in the large genomic database, we found that most Japanese individuals compound heterozygous or homozygous for any of these missense FMO3 variants or known severe mutations [e.g., p.(Cys197Ter)] had impaired FMO3-dependent N-oxygenation of malodorous trimethylamine.

Genetic polymorphisms contribute to inter-individual variability in the metabolism of multiple clinical drugs, including warfarin, thiopurines, primaquine, and aminoglycosides. A rapid and sensitive clinical assessment of various genome biomarkers is, therefore, required to predict the individual responsiveness of each patient to these drugs. In this study, we developed a novel genotyping method for the detection of nine pharmacogene variants that are important in the prediction of drug efficiency and toxicity. This genotyping method uses competitive allele-specific PCR and a single-stranded tag hybridization chromatographic printed-array strip (STH-PAS) that can unambiguously determine the presence or absence of the gene variant by displaying visible blue lines on the chromatographic printed-array strip. Notably, the results of our STH-PAS method were in 100% agreement with those obtained using standard Sanger sequencing and KASP assay genotyping methods for CYP4F2 gene deletion. Moreover, the results were obtained within 90 min, including the PCR amplification and signal detection processes. The sensitive and rapid nature of this novel method make it ideal for clinical genetic testing to predict drug efficacy and toxicity, and in doing so will aid in the development of individualized medicine and better patient care.

Genetic variations within cytochrome P450 2B6 (CYP2B6) contribute to inter-individual variation in the metabolism of clinically important drugs, including cyclophosphamide, bupropion, methadone and efavirenz (EFZ). In this study, we performed an in vitro analysis of 40 CYP2B6 allelic variant proteins including seven novel variants identified in 1070 Japanese individuals. Wild-type and 39 variant proteins were heterologously expressed in 293FT cells to estimate the kinetic parameters (Km, Vmax, and CLint) of EFZ 8-hydroxylation and 7-ethoxy-4-trifluoromethylcoumarin (7-ETC) O-deethylation activities. The concentrations of CYP2B6 variant holo-enzymes were measured by using carbon monoxide (CO)-reduced difference spectroscopy, and the wild-type and 28 variants showed a peak at 450 nm. The kinetic parameters were measured for the wild-type and 24 variant proteins. The values for the remaining 15 variants could not be determined because the enzymatic activity was not detected at the highest substrate concentration used. Compared to wild-type, six variants showed significantly decreased EFZ 8-hydroxylation CLint values, while these values were significantly increased in another six variants, including CYP2B6.6. Although 7-ETC O-deethylation CLint values of CYP2B6 variants did not differ significantly from that of CYP2B6.1, the CLint ratios obtained for 7-ETC O-deethylation were highly correlated with EFZ 8-hydroxylation. Furthermore, three-dimensional structural modeling analysis was performed to elucidate the mechanism of changes in the kinetics of CYP2B6 variants. Our findings could provide evidence of the specific metabolic activities of the CYP2B6 proteins encoded by these variant alleles.