{'fraction_non_alphanumeric': 0.03709283755687468, 'fraction_numerical': 0.021876035797149487, 'mean_word_length': 5.0928951716541215, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Metabolomic Signatures in Doxorubicin-Induced Metabolites Characterization, Metabolic Inhibition, and Signaling Pathway Mechanisms in Colon Cancer HCT116 Cells

Doxorubicin (DOX) is a chemotherapeutic agent is used for various cancer cells. To characterize the chemical structural components and metabolic inhibition, we applied a DOX to HCT116 colon cancer cells using an independent metabolites profiling approach. Chemical metabolomics has been involved in the new drug delivery systems. Metabolomics profiling of DOX-applied HCT116 colon cancer cellular metabolisms is rare. We used 1H nuclear magnetic resonance (NMR) spectroscopy in this study to clarify how DOX exposure affected HCT116 colon cancer cells. Metabolomics profiling in HCT116 cells detects 50 metabolites. Tracking metabolites can reveal pathway activities. HCT116 colon cancer cells were evenly treated with different concentrations of DOX for 24 h. The endogenous metabolites were identified by comparison with healthy cells. We found that acetate, glucose, glutamate, glutamine, sn-glycero-3-phosphocholine, valine, methionine, and isoleucine were increased. Metabolic expression of alanine, choline, fumarate, taurine, o-phosphocholine, inosine, lysine, and phenylalanine was decreased in HCT116 cancer cells. The metabolic phenotypic expression is markedly altered during a high dose of DOX. It is the first time that there is a metabolite pool and phenotypic expression in colon cancer cells. Targeting the DOX-metabolite axis may be a novel strategy for improving the curative effect of DOX-based therapy for colon cancer cells. These methods facilitate the routine metabolomic analysis of cancer cells.


Introduction
Colorectal carcinoma is the world's third most commonly diagnosed cancer and the second most common cause of death from cancer [1]. Chemotherapeutic drugs of doxorubicin hydrochloride (DOX) have been utilized to eradicate leukemia cells due to their anticancer activity and inhibition of tumorigenesis by reducing cell division. However, DOX-based drug resistance has limited their utility [2,3]. The development of novel cancer therapeutics relies heavily on in vivo systems. At the same time, cancer cell lines provide a comprehensive dataset of the alterations in the metabolite levels in cancer cells and can produce significant information [4]. DOX is used as a chemotherapeutic drug for many cancers. Clinicians frequently need to reduce the dose, reducing efficacy [5,6].
Metabolomics studies have already shown alterations in colon cancer metabolites [7][8][9][10]. Metabolomics is typically performed using millions of cells, cultured cells, and tumor specimens [11,12]. Metabolomics has been a great study for comprehensive small-molecule analysis. The small molecules (e.g., sugars, amino acids, proteins, nucleic acids, and lipids) are essential for improving the overall quality of disease and are required for palliative care [13]. As a 2 of 13 result, metabolites serve as the strongest predictors of cellular phenotypic expression and reveal substantial changes inside cancer cells [14,15]. Additionally, metabolites can interact with other biomolecules expressed in sick cells. Today, NMR spectroscopy's metabolomic profiling in cancer cells is a rapidly advancing technology.
Herein, we describe an evaluating biomimetic pattern in colorectal carcinoma (HCT116 cancer cells) using metabolomics analysis that was proposed for investigating the DOXapplied metabolic pathways in cells. We also evaluate how the DOX treatment of HCT116 cancer cells affected their metabolic pathways. The experimental analyses in this work utilize cytotoxicity and metabolomics approaches. Herein, we demonstrate a practical approach to DOX-applied cellular viability, metabolic pathways, and metabolic discrimination in colon cancer cells. We provide strong evidence that cancer metabolic phenotypes reflect DOX-associated systemic metabolic changes.

DOX Solution Preparation
DOX was distributed in an aqueous media in accordance with earlier findings [16,17]. Briefly, 0.75, 1.5, and 7.5 µg/mL of DOX was added to 1% (w/v) bovine serum albumin (BSA) in 2 mL phosphate buffer solution (PBS). In order to improve the dispersibility of DOX, BSA was utilized. The obtained suspensions were then dispersed using a homogenizer with 200 W for 10-15 min under sterile conditions. After centrifugation at 2000× g for 15 min, insoluble DOX aggregates in the pellets were removed, and the supernatants were further homogenized by ultrasonic dispersive technology at 200 W for 15 min under sterile conditions. Cells were homogenized using a Precellys24 (PeqLab, Erlangen, Germany). The resulting DOX dispersions were kept at 4 • C for further research. Each solution was made with Milli-Q water with a resistivity (ultrapure water: 18.2 MΩ·cm).

WST-8 Cell Viability Assay
According to previous tests [18], HCT116 cells were seeded in 96-well plates at 5 × 10 3 cells per well, leaving them to adhere for 24 h. Then culture medium was changed by DOX different concentrations for 24 h. After incubation, highly sensitive water-soluble tetrazolium salt (WST-8; 10 µL/mL) was added to each well, and cell cytotoxicity was measured using a Cellrix ® viability assay kit (Seoul, Korea). Cell viability was calculated by normalizing the absorbance at 450 nm. The non-DOX-treated cells were considered to be 100% viable.

Isolation of HCT116 Cells and Quenching Metabolome
Samples were prepared for metabolomics analysis as previously described [16,19,20]. First, 2 mL of medium for each sample was collected after 24 h, frozen in liquid nitrogen, and stored at −80 • C. All samples were allowed to thaw for 5 min. Next, 80 µL of chloroform was added to the samples at 4 • C and left to incubate for 30 min. Samples were then homogenized by vortexing and transferred to a 1.5 mL plastic tube. Next, 125 µL of deionized water and 125 µL of chloroform were added at 4 • C and vortexed again. For phase separation, samples were centrifuged at 13,000× g for 20 min at 4 • C. Subsequently, the two phases found in the aqueous phase (top: polar metabolites) and the organic phase (bottom: nonpolar metabolites) were separated into different tubes. The aqueous phase (polar metabolites) was lyophilized overnight to remove methanol and water.
A lot of research was conducted on the outcomes of cell quenching in liquid nitrogen. This typically serves as the initial step in various extraction techniques (immediately following the removal of the growing medium) right before cell washing to remove the medium residues. However, it was found in this study that washing the DOX-applied HCT116 cells after the quenching stage resulted in a considerable loss of the cell metabolites, most likely because the freezing step causes the cell wall to collapse, resulting in metabolite leakage. Quenching and washing were performed in reverse order to solve this issue.
Samples were allowed to thaw for 5 min. Then, 400 µL of all samples was diluted with 100 µL of buffer solution to a final concentration of 150 mM with phosphate buffer (pH: 7.2 in 100% D 2 O) containing 1 mM TSP-d 4 as an internal standard. This solution was then added to the aqueous phase. After being vortexed and centrifuged at 12,000× g for 5 min, the supernatant from the samples was transferred in the amount of 550 µL to a 5 mm NMR tube.

Instrumentation of NMR and Data Files
NMR data were analyzed as previously described [16,21]. Briefly, 1D proton ( 1 H)-NMR spectra were recorded at 26 • C using a Bruker Advance 600 MHz NMR spectrometer (Billerica, MA, USA). Spin-echo sequence (recycle delay-90 • -(τ-180 • -τ) n -acquisition) of t 2 -edited Carr-Purcell-Meiboom-Gill (CPMG) was used. Water-suppressed 1 H-NMR spectra were acquired. Before measurements, for equilibrium maintenance, samples were kept for 5 min inside the NMR probe at 298 K. A total of 128 scans were requested. The width was found at 9615.4 Hz. 1 H-spectral resonance width contained 15 to −1 ppm, and each sample took 8-10 min to be measured. The resonance of the water domain (δ 4.48-4.68) was removed. The phosphate buffer containing TSP-d4 was used as a reference at 0 ppm. The data settings of 1 H (e.g., baseline, filling, and TSP adjustment at 0 ppm) were handled by using the Chenomx NMR Suite software (Chenomx Inc., Edmonton, AB, Canada).

Computational Tools and Statistical Analysis
Use of the in-house MetaboAnalyst 5.0 software, score plot analysis such as principal component analysis (PCA), and orthogonal projections to latent structures discriminant analysis (OPLS-DA) were applied to the NMR data. The goodness of fit (R 2 ), the goodness of prediction (Q 2 ) of the score plot, metabolite sets enrichment analysis (MSEA), and VIP scores (cutoff, >1) were assessed. The difference between R 2 and Q 2 decreases with the score plots' robustness.
The logarithmic transformation of the acquired concentrations was used to reduce the impact of both noise and the high variance of the variables. Pareto scaling was applied. It can divide each variable by the square root of its standard deviation. The normalized metabolite concentrations (mM) are presented as the mean ± standard deviations (SD). p-values < 0.05 were considered to be statistically significant. Figure 1a shows the chemical structure of DOX [22]. As per Martineau et al. [23], 1 H-NMR-based metabolomics analysis to recover the cell metabolites was achieved by combining different steps. Identifying several crucial portions that needed a significant amount of optimization was possible by examining and assessing the various extraction techniques. Figure 1b shows the flowchart design and the experimental pipeline for the metabolomics analysis. All data were autoscaled before analysis. The spectral features were extracted and used for pathway analysis in extracted metabolites.

Optimization of the Quenching and Extraction Procedures for DOX-Treated Cells
of prediction (Q 2 ) of the score plot, metabolite sets enrichment analysis (MSEA), and VIP scores (cutoff, >1) were assessed. The difference between R 2 and Q 2 decreases with the score plots' robustness.
The logarithmic transformation of the acquired concentrations was used to reduce the impact of both noise and the high variance of the variables. Pareto scaling was applied. It can divide each variable by the square root of its standard deviation. The normalized metabolite concentrations (mM) are presented as the mean ± standard deviations (SD). pvalues < 0.05 were considered to be statistically significant. Figure 1a shows the chemical structure of DOX [22]. As per Martineau et al. [23], 1 H-NMR-based metabolomics analysis to recover the cell metabolites was achieved by combining different steps. Identifying several crucial portions that needed a significant amount of optimization was possible by examining and assessing the various extraction techniques. Figure 1b shows the flowchart design and the experimental pipeline for the metabolomics analysis. All data were autoscaled before analysis. The spectral features were extracted and used for pathway analysis in extracted metabolites.

Cell Viability-Treatment with DOX in HCT116 Cells
The cell viability of HCT116 cells cultured in the hypoxic environment decreased by 81 percent from 6 to 72 h. The effect of treatments with DOX and cellular viability are summarized in Figure 2.

Cell Viability-Treatment with DOX in HCT116 Cells
The cell viability of HCT116 cells cultured in the hypoxic environment 81 percent from 6 to 72 h. The effect of treatments with DOX and cellular summarized in Figure 2.  (6,24, and 72 concentration (0, control; 0.75 µg/mL; 1.5 µg/mL; and 7.5 µg/mL) at three time points At least 6 replicates were used for each tested concentration. The bar graphs show µg/mL, blue; 1.5 µg/mL, green; and 7.5 µg/mL, red. Tested with 2-way ANOVA wit tiple comparison test, significant adjusted p < 0.05. * and ** represent p < 0.05 and p tively.
DOX-applied HCT116 cells develop into a subconfluent monolayer, an covered that it was particularly challenging to remove the cell growth med washing process entirely. Because of its residual signals that were disperse the central area of the NMR spectrum, glucose, one of the most prevalent co the medium, presented special difficulties for the spectral interpretation of metabolome. There are now four washing steps instead of three to preve cleaning, the cells could be mechanically scraped off the plates, and the meta of the cells could then be quickly put to rest by liquid nitrogen.

Integration of Targeted Metabolomics Datasets
This study determined and targeted 50 metabolites with DOX concent 1 H-edited spectroscopic peaks. Most 1H signals fall in the 0.5-5.5 ppm regi cused the analysis on it. Phosphate buffer was employed as a reference. Per gle, the spectral data multiplicity of DOX was calculated and calibrated. The tabolites were profiled and quantified using their 1 H nuclei (e.g., CH, CH2, a proton multiplicity, as shown in Table 1. This technique is utilized to estimate probability density functions, whi in practical applications of probability distributions in DOX-applied HCT metabolites and their normalized concentrations (mM) were evaluated base vent. OPLS-DA score plot-based DOX in HCT116 cellular metabolic discrim follows: 0.75 µg/mL, 58.9%; 1.5 µg/mL, 75.8%; and 7.5 µg/mL, 39.5%. The disperse the fundamental level of metabolic variations (Figure 3a-c). The u technique of PCA score plot was used to visualize basic metabolic clarity in . Three DOX concentration (0, control; 0.75 µg/mL; 1.5 µg/mL; and 7.5 µg/mL) at three time points were recorded. At least 6 replicates were used for each tested concentration. The bar graphs show 0, black; 0.75 µg/mL, blue; 1.5 µg/mL, green; and 7.5 µg/mL, red. Tested with 2-way ANOVA with Tukey's multiple comparison test, significant adjusted p < 0.05. * and ** represent p < 0.05 and p < 0.001, respectively. DOX-applied HCT116 cells develop into a subconfluent monolayer, and it was discovered that it was particularly challenging to remove the cell growth media during the washing process entirely. Because of its residual signals that were dispersed throughout the central area of the NMR spectrum, glucose, one of the most prevalent components of the medium, presented special difficulties for the spectral interpretation of the extracted metabolome. There are now four washing steps instead of three to prevent this. After cleaning, the cells could be mechanically scraped off the plates, and the metabolic activity of the cells could then be quickly put to rest by liquid nitrogen.

Integration of Targeted Metabolomics Datasets
This study determined and targeted 50 metabolites with DOX concentrations using 1 Hedited spectroscopic peaks. Most 1H signals fall in the 0.5-5.5 ppm region, so we focused the analysis on it. Phosphate buffer was employed as a reference. Per Pascal's triangle, the spectral data multiplicity of DOX was calculated and calibrated. The potential metabolites were profiled and quantified using their 1 H nuclei (e.g., CH, CH 2 , and CH 3 ) and proton multiplicity, as shown in Table 1.
This technique is utilized to estimate probability density functions, which may result in practical applications of probability distributions in DOX-applied HCT116 cells. The metabolites and their normalized concentrations (mM) were evaluated based on the solvent. OPLS-DA score plot-based DOX in HCT116 cellular metabolic discrimination are as follows: 0.75 µg/mL, 58.9%; 1.5 µg/mL, 75.8%; and 7.5 µg/mL, 39.5%. These data could disperse the fundamental level of metabolic variations (Figure 3a-c). The unsupervised technique of PCA score plot was used to visualize basic metabolic clarity in DOX-treated HCT116 cell data (data not shown). At 1.5 µg/mL of DOX, the highest percentage of metabolic classification was found compared to others. The OPLS-DA score plot indicated that all samples were within the 95% confidence interval (Hotelling T2 ellipse). OPLS-DA analysis was used to classify and filter the metabolites. Note and Abbreviations; a -Singlet (s); doublet (d); triplet (t); doublet of doublets (dd); quartets (q); multiplets (m); molecular weight (MW); parts-per-million (ppm); adenosine monophosphate (AMP); adenosine diphosphate (ADP); adenosine triphosphate (ATP); Trimethylamine N-oxide (TMAO).
The cumulative R 2 and Q 2 values were greater than 0.6, indicating that the OPLS-DA score plot was suitable for determining the difference between control and HCT116 cells (Figure 3d). The metabolic deference was screened from low DOX to high DOX concentrations. This metabolic difference of DOX treatment was further confirmed with a 100-premutation test, which could verify its stability and reliability. Here, the OPLS-DA score plot did not overfit.  The cumulative R 2 and Q 2 values were greater than 0.6, indicating that the OPLS-DA score plot was suitable for determining the difference between control and HCT116 cells (Figure 3d). The metabolic deference was screened from low DOX to high DOX concentrations. This metabolic difference of DOX treatment was further confirmed with a 100premutation test, which could verify its stability and reliability. Here, the OPLS-DA score plot did not overfit.

Volcano Plots and Heatmap Analysis of HCT116 Cells Metabolites
Figure 4a-c show volcano plots of differentially expressed metabolites in HCT116 cancer cells induced by DOX. False discovery rates (FDRs) were plotted against DOXinduced log2 fold changes in HCT116 cancer cells. The 1.5 µg/mL of DOX-applied up-and downregulated metabolites were found to be more than 0.75 µg/mL of DOX. Interestingly, there are no upregulated metabolites in 7.5 µg/mL of DOX-treated in HCT116 cells ( Figure  4c). The DOX-applied volcano plot has been studied in all drug-applied groups. A volcano plot of HCT116 cellular metabolites in positive and negative modes was summarized.
The DOX-treated top 25 metabolites' phenotypic expression is shown as a heatmap in Figure 4d. The metabolites from heatmap analysis of HCT116 cells were identified as candidate biomarkers. The metabolites were filtered and collated with a VIP score > 1 (data not shown). The top 25 metabolites were identified as candidate biometers. Expression levels were assessed and compared to determine the metabolites in 0.75 µg/mL, 1.5 µg/mL, and 7.5 µg/mL of DOX in HCT116 cells. The DOX-treated top 25 metabolites' phenotypic expression is shown as a heatmap in Figure 4d. The metabolites from heatmap analysis of HCT116 cells were identified as candidate biomarkers. The metabolites were filtered and collated with a VIP score > 1 (data not shown). The top 25 metabolites were identified as candidate biometers. Expression levels were assessed and compared to determine the metabolites in 0.75 µg/mL, 1.5 µg/mL, and 7.5 µg/mL of DOX in HCT116 cells.

Tailored Metabolic Adaptations to DOX and Primary Metabolites with Pathways
Among the amino acid metabolites, glutamate, glutamine, sn-glycero-3-PCho, valine, methionine, and isoleucine were increased in HCT116 cells by DOX treatment. Acetate and glucose were increased in cancer cells (Figure 5a). Acetate and glucose are the by-products of energy metabolism. DOX concurrently increased many significant metabolites in HCT116 cells. Several amino acids were significantly decreased by DOX, when amino acids may be involved in autophagy and mitophagy [24]. DOX treatment also induced metabolic downregulation in specific amino acids such as alanine, choline, taurine, o-phosphocholine, inosine, lysine, and phenylalanine (Figure 5b). These observations suggested a decrease in ATP production through glycolysis and oxidative phosphorylation. This could be related to the inhibition of metabolic synthesis in the cells.

Tailored Metabolic Adaptations to DOX and Primary Metabolites with Pathways
Among the amino acid metabolites, glutamate, glutamine, sn-glycero-3-PCho, valine, methionine, and isoleucine were increased in HCT116 cells by DOX treatment. Acetate and glucose were increased in cancer cells (Figure 5a). Acetate and glucose are the byproducts of energy metabolism. DOX concurrently increased many significant metabolites in HCT116 cells. Several amino acids were significantly decreased by DOX, when amino acids may be involved in autophagy and mitophagy [24]. DOX treatment also induced metabolic downregulation in specific amino acids such as alanine, choline, taurine, o-phosphocholine, inosine, lysine, and phenylalanine (Figure 5b). These observations suggested a decrease in ATP production through glycolysis and oxidative phosphorylation. This could be related to the inhibition of metabolic synthesis in the cells.
The metabolites were quantified as before by targeted metabolomics [25]. The increase in hepatic glucose content suggests an enhancement of glycogenesis. Taurine and glutathione (GSH) were most closely associated with liver injury after DOX exposure. GSH, which is biosynthesized from cysteine, glycine, and glutamate, is the most abundant cellular redox molecule and plays a significant role in antioxidant functions [15,19].
The DOX effect appears specific for neurological metabolites, glutamate, and glutamine passed to induce neurogenesis behaviors. The neurological effects of glutamate and glutamine were observed in adult mice C57BL/6J mice treated with indole and DOX [26]. Several molecules are also associated with neurological dysfunction. The relative concen- The metabolites were quantified as before by targeted metabolomics [25]. The increase in hepatic glucose content suggests an enhancement of glycogenesis. Taurine and glutathione (GSH) were most closely associated with liver injury after DOX exposure. GSH, which is biosynthesized from cysteine, glycine, and glutamate, is the most abundant cellular redox molecule and plays a significant role in antioxidant functions [15,19].
The DOX effect appears specific for neurological metabolites, glutamate, and glutamine passed to induce neurogenesis behaviors. The neurological effects of glutamate and glutamine were observed in adult mice C57BL/6J mice treated with indole and DOX [26]. Several molecules are also associated with neurological dysfunction. The relative concentration of eight upregulated and eight downregulated metabolites from the 1 H-NMR analysis were analyzed ( Figure 5). These metabolites showed valid metabolites for OPLS-DA modal discrimination. There are many significant metabolites among DOX-treated groups.
Phosphocholine is a metabolite that partially comes from processing dietary carnitine into trimethylamine and choline. The evidence demonstrated that valine was associated with brain injury, which our study further confirmed by increased valine in DOC treatment. L-valine deficiency is noticeable in neurological defects in the brain [27].
The top metabolites identified by this analysis highlighted changes in selected metabolic pathways in HCT116 cancer cells. Pathways enrichment analysis was employed using the total quantified metabolites (Figure 6a-c). The pathways are represented as circles. We confirmed that metabolites taken from cancer cells undergoing transformation demonstrated significant changes in metabolites and metabolic pathways. Taurine and purine metabolism were highly altered in cancer cells, which could be used at therapeutic intervals. This greater capacity is further confirmed by HCT116 cells having a higher mitochondrial volume than control or DOX cells.
tration of eight upregulated and eight downregulated metabolites from the 1 H-NMR analysis were analyzed ( Figure 5). These metabolites showed valid metabolites for OPLS-DA modal discrimination. There are many significant metabolites among DOX-treated groups.
Phosphocholine is a metabolite that partially comes from processing dietary carnitine into trimethylamine and choline. The evidence demonstrated that valine was associated with brain injury, which our study further confirmed by increased valine in DOC treatment. L-valine deficiency is noticeable in neurological defects in the brain [27]. Figure 5. Essential metabolites are assigned higher priority scores using metabolomics profiling with DOX treatment. (a) Expression of acetate, glucose, glutamate, glutamine, sn-glycero-3-PCho, valine, methionine, and isoleucine were increased metabolites in HCT116 cancer cells. (b) Expression of alanine, choline, fumarate, taurine, o-phosphocholine, inosine, lysine, and phenylalanine were decreased their regulation in HCT116 cancer cells. The black dot represents the control; the blue dot represents 0.75 µg/mL of DOX; the green dot represents 1.5 µg/mL of DOX; and the red dot represents 7.5 µg/mL of DOX. *, ** and *** represent p < 0.05, p < 0.01, and p < 0.001, respectively. The top metabolites identified by this analysis highlighted changes in selected metabolic pathways in HCT116 cancer cells. Pathways enrichment analysis was employed using the total quantified metabolites (Figure 6a-c). The pathways are represented as circles. We confirmed that metabolites taken from cancer cells undergoing transformation demonstrated significant changes in metabolites and metabolic pathways. Taurine and purine metabolism were highly altered in cancer cells, which could be used at therapeutic intervals. This greater capacity is further confirmed by HCT116 cells having a higher mitochondrial volume than control or DOX cells.
The metabolites regulation successfully distinguished and identified the cancer metabolites from the control metabolites. The gut microbe ecosystem modulates and shapes Figure 5. Essential metabolites are assigned higher priority scores using metabolomics profiling with DOX treatment. (a) Expression of acetate, glucose, glutamate, glutamine, sn-glycero-3-PCho, valine, methionine, and isoleucine were increased metabolites in HCT116 cancer cells. (b) Expression of alanine, choline, fumarate, taurine, o-phosphocholine, inosine, lysine, and phenylalanine were decreased their regulation in HCT116 cancer cells. The black dot represents the control; the blue dot represents 0.75 µg/mL of DOX; the green dot represents 1.5 µg/mL of DOX; and the red dot represents 7.5 µg/mL of DOX. *, ** and *** represent p < 0.05, p < 0.01, and p < 0.001, respectively. The metabolites regulation successfully distinguished and identified the cancer metabolites from the control metabolites. The gut microbe ecosystem modulates and shapes many metabolic, immunological, structural, and neurological functions [28]. Most appropriate metabolisms were marked. Figure 6 shows that 0.75 µg/mL of DOX significantly regulated the citrate cycle (TCA cycle); taurine and hypotaurine metabolism; valine, leucine, and isoleucine degradation; and alanine, aspartate, and glutamate metabolism. At 1.5 µg/mL of DOX in colon cancer cells, the alanine, aspartate, and glutamate metabolism; D-glutamine, and D-glutamate metabolisms; citrate cycle (TCA cycle); and taurine and hypotaurine metabolism have been altered. Similarly, arginine and proline metabolism; alanine, aspartate, glutamate metabolism; D-glutamine and D-glutamate metabolism; and glycolysis/gluconeogenesis most robustly contributed to 7.5 µg/mL of DOX treatment.
Study findings support that changes in colon cancer metabolites could contribute to physiological responses to gut microbiota. Human prostate cancer metabolomics via high-resolution NMR spectroscopy has analyzed the metabolites in prostate tissue, seminal fluid, serum, and urine [29,30]. Figure 7 represents the HCT116 cancer cells after DOX treatment. The schematic diagram of glycolysis and TCA cycle imbalance by DOX in colon cancer cells has been summarized.  This study has some limitations. These findings show that DOX resistance exhibits separate basic metabolic vulnerabilities despite shared chemotherapeutic adaptations. Metabolic therapies can address these vulnerabilities to inhibit drug-resistant tumor growth.

Discussion
The study revealed the metabolite regulation of HCT116 cells with DOX treatment. We reduced the noise ratio in all spectrums. We increased the spectral signal resolution by employing NMR spectroscopy. We reduced contamination by eliminating sample drying and increased NMR performance by extracting metabolites with 80% acetonitrile [31,32]. This technique can be used to evaluate any cellular metabolites, but in practice, it Figure 7. The schematic regulation of TCA cycle-associated activity by DOX treatment in HCT116 cancer cells. Red-and green-colored metabolites were represented as increased and decreased in HCT116 cells, respectively. The red sign shows a possible metabolic chemical reaction inhibition in cancer cells. *, ** and *** represent p < 0.05, p < 0.01, and p < 0.001, respectively. This study has some limitations. These findings show that DOX resistance exhibits separate basic metabolic vulnerabilities despite shared chemotherapeutic adaptations. Metabolic therapies can address these vulnerabilities to inhibit drug-resistant tumor growth.

Discussion
The study revealed the metabolite regulation of HCT116 cells with DOX treatment. We reduced the noise ratio in all spectrums. We increased the spectral signal resolution by employing NMR spectroscopy. We reduced contamination by eliminating sample drying and increased NMR performance by extracting metabolites with 80% acetonitrile [31,32]. This technique can be used to evaluate any cellular metabolites, but in practice, it is most helpful when cell counts are modest and enzymatic dissociation is unnecessary. The R 2 Y and Q 2 Y values in OPLS-DA score plots were observed in the distinct metabolic clusters of DOX-treated HCT116 cells [33,34].
Our observation of increased glucose levels by DOX treatment might suggest the inhibition of glycolytic processes and could be linked with the inhibition of cell proliferation. The intervention in glutamate metabolism by glutaminase inhibitors was shown to impact glutathione levels, which could reduce the antioxidative tumor capacity. Due to the high incidence of DOX-mediated anticancer activity, research usually positively affects cellar proliferation. DOX-responsive metabolic profiling detected early toxicity metabolite signature [35]. Evidence suggests that DOX causes these current harmful consequences by producing free radicals, lipid peroxidation, and oxidative damage to tissues and organs [36,37].
An increased amount of acetate and glucose in cancer cells was observed after repeated exposure to DOX. It was concluded that the DOX might remarkably block some signaling pathways. The alteration in glycolysis, TCA cycle, nucleotide synthesis, and choline metabolisms could be associated with growth suppression, which we found in DOX-treated HCT116 cells. These findings underscore the significance of colon cancer metabolites that regulate metabolic reactions that disturb energy, amino acids during metabolic transcriptional changes in mitochondria, and adult neurotransmitter impairment. Herein, we inferred the mechanism of metabolites inhibition enhancing the efficacy of chemotherapeutics could be clarified through metabolomics analysis. Notably, the changes in metabolites indicated that DOX was likely a source of therapeutic drugs for colon cancer cells. The DOX also contains other bioactive therapeutic drug roles.
Breast cancer cells resistant to both medications demonstrated well-known anthracyclines resistance mechanisms in a cancer cell line, such as enhanced drug efflux, lysosomal activity, and oxidative stress response [38,39]. Our in vitro studies display direct, metabolite-specific disorders, demonstrating phenotypic expression [40,41].
Our data gave evidence that DOX inhibits the TCA cycle-medicated mitochondrial damage. However, the reduced fumarate production by DOX may inhibit the ATP molecule production from TCA cycle transporters. Our study shows TCA cycle arrest, metabolite depletion, and mitochondrial dysfunction induced by DOX. In this context, impaired metabolites in mitochondria have been linked to anticancer activity [40,42,43].
We observed a metabolic difference in the abundance of metabolites in HCT116 cells. Functional studies will be required to assess the biological significance of the difference. Cancer cells are undergone metabolic changes to survive oxidative stress during metastasis [44,45]. A better understanding of the metabolic changes could reveal new therapeutic vulnerabilities to colon cancer progression.

Conclusions
We have implemented a fully automated study in HCT116 cancer cells with optimized peak detection, alignment, and an annotation task using 1 H-NMR data. The raw spectra processing produces a peak intensity table and peak annotation table. We found that the metabolites of glutamate, glutamine, sn-glycero-3-PCho, valine, methionine, and isoleucine were increased. The specific amino acids such as alanine, choline, taurine, ophosphocholine, inosine, lysine, and phenylalanine were downregulated in colon cancer HCT116 cells. In the future, we aim to screen more data in HCT116 cells by DOX treatment to support cancer metabolic regulation.