Antioxidant Mechanism of Lactiplantibacillus plantarum KM1 Under H2O2 Stress by Proteomics Analysis

Lactiplantibacillus plantarum KM1 was screened from natural fermented products, which had probiotic properties and antioxidant function. The survival rate of L. plantarum KM1 was 78.26% at 5 mM H2O2. In this study, the antioxidant mechanism of L. plantarum KM1 was deeply analyzed by using the proteomics method. The results demonstrated that a total of 112 differentially expressed proteins (DEPs) were screened, of which, 31 DEPs were upregulated and 81 were downregulated. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that DEPs participated in various metabolic pathways such as pyruvate metabolism, carbon metabolism, trichloroacetic acid cycle, amino acid metabolism, and microbial metabolism in diverse environments. These metabolic pathways were related to oxidative stress caused by H2O2 in L. plantarum KM1. Therefore, the antioxidant mechanism of L. plantarum KM1 under H2O2 stress provided a theoretical basis for its use as a potential natural antioxidant.


INTRODUCTION
Oxidative stress is an imbalance between the production of reactive oxygen species (ROS) and the antioxidant defense of cells (Jacob et al., 2013). Oxidative stress plays an important role in human health, which is related to many diseases (Preiser, 2012), such as chronic diseases, cardiovascular diseases, atherogenesis, and neurodegenerative diseases (Zhou et al., 2008;Sharifpanah and Sauer, 2016). Antioxidants improve oxidative stress as a target to treat these diseases, which are divided into natural antioxidants and synthetic antioxidants according to their sources. However, most of synthetic antioxidants were proved to have side effects (Kulawik et al., 2013). Therefore, it is necessary to find safe and reliable natural antioxidants. Studies showed that lactic acid bacteria (LAB) have antioxidant effects (Lin and Yen, 1999;Wang et al., 2012Wang et al., , 2021Teng et al., 2020).
Lactic acid bacteria are generally recognized as safe (GRAS) microorganisms (Azat et al., 2016) and also important food resources with certain probiotic activities, including regulating gastrointestinal, lowering cholesterol, alleviating diseases, and relieving lactose intolerance and anti-oxidation (Park et al., 2002;Mishra et al., 2015). As an inexpensive natural antioxidant, Lactiplantibacillus plantarum can scavenge free radicals and exert its antioxidant functions to regulate oxidative stress (Teng et al., 2020). Additionally, the antioxidant properties of L. plantarum were also attributed to the intracellular antioxidant enzyme system (Liu et al., 2021). Under oxidative stress, antioxidant enzymes such as catalase, NADH peroxidase, and thioredoxin reductase were activated to resist oxidative stimulation. L. plantarum MA2 with the antioxidant capacity in vitro was identified, and the antioxidant-related genes including cat, gshR, and npx were upregulated under hydrogen peroxide (H 2 O 2 ) challenge (Tang et al., 2017). The gshR and npx were upregulated in L. plantarum AR113 with probiotic activities after H 2 O 2 treatment (Lin et al., 2020). In addition, many genes in various pathways changed from the global transcriptomics perspective of L. plantarum CAUH2 in response to H 2 O 2 (Zhai et al., 2020). It is a particularly important method to use H 2 O 2 as an oxidant to screen LAB with antioxidant function. In this way, LAB can exert its antioxidant function and use its antioxidant mechanism to alleviate the negative effects of oxidative stress. However, the antioxidant mechanism of LAB is not yet clear, and few studies have reported protein changes associated with oxidative stress.
Moreover, proteomics becomes a hot topic in current genomics research. Proteomics not only characterizes and identifies large-scale proteins but also determines protein activities, functions, and interactions (Vittaladevaram, 2021). More studies have carried out comparative proteomic analysis based on stress conditions. To study the changes and functional roles of the differentially expressed proteins (DEPs), the DEPs of bacteria were quantitatively analyzed through the TMT labeling quantitative proteomics technology (O'Connor et al., 2000). The potential mechanism of nitrite degradation by Lactobacillus fermentum RC4 under nitrite stress was analyzed using proteomics (Zeng et al., 2019). Proteomics technology was used to analyze the DEP changes of Phaffia rhodozyma under TiO 2 stress and to clarify the metabolic mechanism of the strain promoting astaxanthin synthesis (Zhang et al., 2020). When faced with various stress conditions, strains effectively cope with unfavorable conditions under the action of DEPs. Therefore, the mechanism of strains under oxidative stress can be more clearly reflected based on the proteomic analysis.
L. plantarum KM1 isolated from natural fermented food tolerated high concentrations of H 2 O 2 and satisfied the basic criteria for probiotic screening. The antioxidant ability of this strain was evaluated, and its antioxidant mechanism was explored from the proteomics perspective, which provided a theoretical foundation for its use as a natural antioxidant for functional foods.

Strain Isolation and Identification
A piece of naturally fermented food was directly taken into sterilized Mann, Rogosa and Sharpe (MRS) medium (Qingdao Hi-Tech Park Haibo Biotechnology Co., Ltd., China) and Abbreviations: ABTS, 3-ethylbenzothiazoline-6-sulfonic acid; ACN, acetonitrile; DEPC, diethyl pyrocarbonate; DPPH, 1,1-diphenyl-2-picrylhydrazyl radical; DTT, dithiothreitol; GRAS, generally recognized as safe; SDT, 4% sodium lauryl sulfate (SDS), 100mM Tris/HCl pH 7.6, 0.1M DTT; TMT, tandem mass tag. cultivated at 37 • C for 24 h. The strains were inoculated at 3% inoculum into MRS medium containing 10 mM H 2 O 2 for preliminary screening. At this time, the surviving strains were inoculated into a new MRS liquid medium for 24 h, which were the primary screening strains. The strains were repeatedly isolated and purified on MRS agar plates and cultured at 37 • C for 48 h. By observing the morphology of the colony and using a microscope to observe, a single colony was separated and purified to obtain a pure strain. The purified strains were subjected to bacterial 16S-rDNA identification and sequencing. The genomic DNA of the strain was extracted and then performed PCR amplification with universal primers. Finally, the PCR product was recovered and sequenced. The sequencing results came from the Shanghai Biotechnology Engineering Co., Ltd. Sequencing results were compared on the NCBI website for homology comparison, and a phylogenetic tree was constructed.

Determination of H 2 O 2 Concentration
A pure culture of KM1 was inoculated into the MRS broth contained with different concentrations of H 2 O 2 (i.e., 1.25, 2.5, 5.0, and 10.0 mM). The same amount of physiological saline was added to the MRS broth as the control group. KM1 was cultured at 37 • C in incubators (Shanghai Zhicheng Analytical Instrument Manufacturing Co., Ltd.). The samples were shaken briefly, and then the OD value was measured at 600 nm every 2 h.

Determination of Resistance to Acid and Bile Salts
The determination of acid resistance and bile salt resistance was slightly modified according to the reported method (Ribeiro et al., 2018). L. plantarum KM1 was inoculated into the MRS medium (2%, w/w), in which the pH was adjusted to 4.0, 3.0, 2.5, and 2.0 with HCl. The bile salt resistance was determined by inoculating L. plantarum KM1 into the MRS medium containing 0.1 and 0.2% bile salt (Solarbio, China). The normal medium was used as the control group. The above cultures were then cultured at 37 • C, and samples were taken at 0, 2, and 4 h, respectively, and spread on the MRS solid medium for 48 h to count the number of colonies. Each experiment was performed in triplicate.

Antibiotic Sensitivity
According to previous reports, the micro broth dilution method was used to determine the antibiotic sensitivity of strains . Five antibiotics were accurately weighed 5,120 µg, including kanamycin sulfate, streptomycin sulfate, tetracycline HCl, chloramphenicol, and ampicillin trihydrate (Solarbio, China). The above five antibiotics were serially diluted in the range of 256-1 µg/ml, and the final concentration of the antibiotic solution was 256,128,64,32,16,8,4,2, and 1 µg/ml. Then, 2% bacterial solution was added to the diluted solution. Among them, the blank group without antibiotics and bacterial solution and the control group without antibiotics were set up. The cultures were incubated at 37 • C for 24 h to determine its minimum inhibitory concentration (MIC). Each experiment was conducted in triplicate.

ABTS Radical Scavenging
The determination method of ABTS free radicals was modified according to the previous method (Thaipong et al., 2006). ABTS free radical ion working solution was a mixture of 7.4 mM ABTS (Shanghai McLean Biochemical Technology Co., Ltd., China) stock solution and 2.6 mM potassium persulfate solution and must be placed at room temperature in the dark condition for 12 h. The sample was mixed with the diluted ABTS working solution and reacted at 37 • C for 30 min, and an equal volume of 95% ethanol was used as a control group instead of the sample. The absorbance was measured at 734 nm and marked them as A s and A 0 , respectively. Each experiment was conducted in triplicate. The ABTS radical scavenging rate was calculated with the following formula:

Determination of Reducing Activity
The determination of the reducing activity was improved on the basis of the previous reported method (Li et al., 2008). The sample (0.5 ml) was added to 1.5 ml of phosphate buffer (0.2 M, Shanghai Corning Management Co., Ltd., China) and 1.5 ml of 1% (w/v) potassium ferricyanide (Beijing chemical factory, China). The mixture was placed in a water bath at 50 • C for 20 min and quickly cooled to room temperature. The mixture was added to 1.5 ml of trichloroacetic acid (TCA, 10%, Tianjin Guangfu Fine Chemical Research Institute, China) and after shaking well, centrifuged at 8,000 rpm for 10 min to get the supernatant. The supernatant was mixed with the same amount of physiological saline, and 1 ml of ferric chloride solution was added. After standing for 10 min, the absorbance values were read at 700 nm. Different concentrations of L-ascorbic acid (Beijing Dingguo Biotechnology Co., Ltd., China) were made into a standard curve. Each experiment was conducted in triplicate.

Hydroxyl Radical Scavenging
Hydroxyl radical scavenging was conducted according to the Nanjing Jiancheng kit instructions.

DPPH Radical Scavenging Activity Assay
DPPH radical scavenging activity was carried out using the Wei Tang's method (Luan et al., 2021) and modified. The prepared sample was mixed thoroughly at 1:1 with DPPH solution formulated with absolute ethanol. The mixed sample was reacted for 30 min at room temperature under dark conditions and centrifuged at 6,000 rpm for 10 min, and the supernatant was taken to measure the absorbance values (A s ) at 517 nm. At the same time, an equal volume of absolute ethanol replaced the DPPH solution to measure the absorbance value and recorded it as A b . The control group replaced the sample with an equal of distilled water volume, and the measured absorbance value was A c . Each experiment was conducted in triplicate. The DPPH radical scavenging rate was calculated with the following formula:  (Wiśniewski et al., 2009), and then total protein was quantified using the BCA protein assay kit (Thermo Scientific, USA). Protein was completely dissolved in SDS according to the filter-aided proteome preparation (FASP) method (Wiśniewski et al., 2009) for trypsin digestion. The peptides were desalted by C18 Cartridge after enzymatic hydrolysis and then freezedried and reconstituted with adding 40 µl dissolution buffer. The peptides were quantified at the absorbance of OD 280 .

TMT Labeling and Strong Cation Exchange Chromatographic Classification
Furthermore, 100 µg peptides of each sample were labeled according to the instructions of the TMT labeling kit (Thermo Fisher, USA). The labeled peptides from each group were pooled and purified using an AKTA Purifier 100 (GE Healthcare, USA). The sample was loaded to the column equilibrated with buffer A (10 mM KH 2 PO 4 , 25% ACN, pH 3.0) and B (10 mM KH 2 PO 4 , 500 mM KCl, 25% ACN, pH 3.0). The chromatographic column (5 µm, 200 Å, PolyLCInc, Maryland, U.S.A.) was equilibrated with liquid A, and the sample was loaded from the injector to the chromatographic column for separation at a flow rate of 1 ml/min. The liquid-phase gradient was as follows: 0-25 min, the linear gradient of liquid B was from 0 to 10%; 25-32 min, the linear gradient of liquid B was from 10 to 20%; 32-42 min, the linear gradient of liquid B was from 20 to 45%; 42-47 min, the linear gradient of B solution was from 45 to 100%; 47-52 min, 52-60 min B solution was maintained at 100%; after 60 min, B solution was reset to 0%. During the elution process, the absorbance value at 214 nm was monitored, and the elution fractions were collected every 1 min, respectively, lyophilized and desalted using a C18 Cartridge.

Liquid Chromatography With Tandem Mass Spectrometry Data Acquisition
Each sample fraction was separated by high-performance liquid chromatography (HPLC) liquid system Easy nLC (Thermo Scientific, USA) with a separation column (Thermo scientific EASY column, 75 µm × 10 cm, 3 µm, C18). After the sample was chromatographed, it was analyzed using the Q-Exactive mass spectrometer (Thermo, USA) for 60 min. The peptides were detected by the positive ion mode, the scanning range was set to 300-1,800 m/z, the resolution was 70,000 (at m/z 200), automatic gain control (AGC) target was 3e6, maximum IT was 10 ms, and the dynamic exclusion time was 40 s. The mass-to-charge ratios of peptides and peptide fragments were calculated by collecting 10 fragment maps (MS2 scan) after each full scan, and the MS resolution was 17,500 (at m/z 200).

Protein Identification and Quantitative Analysis
Mass spectrometry (MS) analysis was performed with software Mascot 2.2 and Proteome Discoverer 1.4 for library identification and quantitative analysis with MS/MS ion search. The search parameters for protein identification and relative quantification were as follows: The maximum number of missed cleavage sites allowed was 2. Fixed modifications were carbamidomethyl (C), TMT 10plex (N-term), TMT 10plex (K), and variable modifications were oxidation (M), TMT 10plex (Y). The peptide mass tolerance was ±20 ppm, and fragment mass tolerance was 0.1 Da. The database mode used to calculate false discovery rate (FDR) was set to Decoy, and FDR ≤ 0.01 was used as the screening criterion for credible proteins. The ratios (H 2 O 2 /control) were calculated as the median of only unique peptides of the protein normalizes all peptide ratios by the median protein ratio. The median protein ratio should be 1 after the normalization. The relative expression of the DEP was the relative expression of the sample group relative to the internal reference protein, which was the ratio of the signal intensity value of the reporter ion peak of the sample labeling channel and the internal standard sample in the experimental design after correction.

Bioinformatics Analysis
The target protein set used the Omicsbean software for Gene Ontology (GO) function annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) signal pathway. GO and KEGG pathway enrichment analyses of the DEPs were performed by the Fisher's exact test method. The detailed information of GO function annotation was obtained by consulting the related database (http://www.geneontology.org/). Information about the KEGG pathway could be accessed through the KEGG database (http://www.kegg.jp/). The STRING database (http://string-db. org) was used to analyze the protein-protein interactions (PPIs) of L. plantarum KM1 including direct physical interactions and functional correlations.

RT-qPCR Assay
Totally, 50 ml of the cultured bacterial solution was centrifuged at 5,000 g and transferred to a non-enzyme centrifuge tube, and 1 ml of Trizol was added and allowed to stand for 15 min after repeated pipetting and lysis. Furthermore, 200 µl of ice-cold chloroform was added to the above lysate and let stand on ice for 5 min. After centrifugation at 12,000 g for 15 min at 4 • C, the upper aqueous phase was carefully aspirated, an equal volume of precooled isopropanol was added, and the mixture was evenly mixed and placed on ice for 10 min. After centrifugation at 12,000 g for 10 min at 4 • C, the supernatant was carefully discarded, and 75% ethanol prepared with 1 ml of DEPC water was added to wash the white precipitate and let stand for 5 min. After centrifugation at 12,000 g for 5 min at 4 • C, the supernatant was discarded, and the pellet was allowed to stand at room temperature for 5 min to dry the pellet. After adding 20 µl of DEPC water to dissolve the precipitate, the extracted RNA is reverse-transcribed into cDNA according to the instructions of the cDNA kit. According to the instructions of TB Green TM Premix Ex Taq TM II kit (Japan Takara Biological Co., Ltd.), samples were added and mixed well for on-machine test. The amplification conditions are shown in Supplementary Table S1. 16sRNA was used as the internal reference, and the primer design was shown in Supplementary Table S2.

Statistical Analysis
The experimental data were expressed as the mean ± SEM. The significance analysis of experimental data was used t-test and one-way ANOVA (and non-parametric or mixed) by Prism 8.0. The PPI network of DEPs was analyzed using the STRING software 11.0.

Screening and Identification of Hydrogen Peroxide-Tolerant Strains
The strain was separated and purified from natural fermented yogurt, and morphology was observed under a microscope ( Figure 1A). The colony of this strain was round, smooth, dense, and white ( Figure 1B). The strain was identified using 16S-rDNA sequencing, and a phylogenetic tree was constructed based on the sequence (Figure 1C). The results showed that this strain had high system homology (99%) with L. plantarum named L. plantarum KM1.
The growth of L. plantarum KM1 were determined at 1.25-10.00 mM H 2 O 2 in MRS medium. The growth of KM1 was gradually inhibited with increasing the H 2 O 2 concentration, and it grew well under 1.25-5.00 mM H 2 O 2 , while it was completely inhibited at 10.00 mM ( Figure 1D). Therefore, L. plantarum KM1 can tolerate a maximum concentration of 5.00 mM H 2 O 2 .

Acid and Bile Salt Resistance Properties of L. plantarum KM1
It is an important probiotic activity to adapt to the human gastric acid and the bile salt environment after ingestion. The acid resistance of L. plantarum KM1 was tested (Figure 2A). L. plantarum KM1 grew well under the condition of pH 2.5 or higher. Under the condition of pH 2.0, the number of viable colonies of the strain dropped significantly from 7.715 ± 0.041 lg CFU/ml to 4.626 ± 0.009 lg CFU/ml in 2 h and significantly reduced to 3.486 ± 0.257 lg CFU/ml in 4 h. L. plantarum KM1 was cultured in the MRS medium at 0.1 and 0.2% bile salt for 2 and 4 h ( Figure 2B). The number of colonies of L. plantarum KM1 was 7.528 ± 0.003 lg CFU/ml after 2 h when bile salt was 0.1%, and increased to 7.629 ± 0.024 lg CFU/ml after 4 h. When the bile salt concentration was 0.2%, the number of L. plantarum KM1 decreased significantly from 7.664 ± 0.009 lg CFU/ml to 6.113 ± 0.019 lg CFU/ml in 2 h and 6.132 ± 0.037 lg CFU/ml in 4 h (p < 0.0001).

Antioxidant Activity of L. plantarum KM1 in Different Component
The scavenging rate of three free radicals (i.e., ABTS, •OH, and DPPH) and reducing activity of different components (i.e., supernatant, intact cell, and cell-free extract) with H 2 O 2 were measured (Figure 3). The ABTS free radical scavenging rate of L. plantarum KM1 supernatant, intact cells, and cell-free extracts were 79.494 ± 0.984%, 51.129 ± 1.400%, and 53.489 ± 1.416%, respectively. There was no significant effect on the ABTS free radical scavenging rate of the three components of L. plantarum KM1 with 5 mM H 2 O 2 treatment.
L-ascorbic acid was used as a standard curve to determine the standard curve of the reducing ability of L. plantarum KM1 ( Figure 3B). The supernatant of L. plantarum KM1 had the strongest reducing activity (4,133.000 ± 23.334 µmol/L) of three components. The reducing activity of intact cells was 271.778 ± 4.006 µmol/L, and the cell-free extracts were 310.67 ± 12.019 µmol/L.
•OH is the most active free radicals in ROS, which cause oxidative damage (Yang et al., 2017). The •OH scavenging rate of the KM1 supernatant was the highest (77.114 ± 0.595%) among two components (Figure 3C). The •OH scavenging rate of cellfree extract of L. plantarum KM1 was 66.963 ± 1.033%. After treatment with 5 mM H 2 O 2 , the •OH scavenging rate of the cellfree extract was significantly lower than the normal group (p < 0.01). L. plantarum KM1 could scavenge •OH and alleviate oxidative damage caused by •OH.
Based on the measurement of DPPH radicals, the scavenging rate of the supernatant of L. plantarum KM1 was 83.173 ± 0.775% (Figure 3D). The DPPH scavenging of the KM1 cell-free extract in control group was 28.173 ± 1.285%. The DPPH free radical scavenging rate of L. plantarum KM1 cell-free extract  Differentially expressed proteins of L. plantarum KM1 were analyzed. The results showed that a total of 112 proteins were differently expressed, of which 31 DEPs such as phosphoglycerate kinase (PGK), homoserine dehydrogenase, and hypothetical protein were upregulated and 81 DEPs were downregulated such as transaldolase (TAL) and pyruvate oxidase ( Figure 4A). The GO function enrichment of DEPs was mainly composed of three functional categories, namely, biological process (BP), cellular component (CC), and molecular function (MF) (Ashburner et al., 2000). DEPs of L. plantarum KM1 were mainly related to the metabolic process, single-organism process, cellular process, and detoxication of BP, and 91 DEPs participated in the metabolic process of BP. At same time, the DEPs exerted antioxidant activity, catalytic activity, binding, and transport activity in oxidation-reduction process of MF, and 82 of them participated in the catalytic activity of MF after H 2 O 2 treatment ( Figure 4B).
Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis was used to further explore the main pathway of L. plantarum KM1 after H 2 O 2 stress ( Figure 4C). The most significant pathway was pyruvate metabolism, and the second one was carbon metabolism and citrate cycle (TCA cycle) from the KEGG pathway enrichment of DEPs. Moreover, the expression changes in proteins in these pathways may be related to the antioxidation of the strain. In addition, the KEGG enrichment pathway also showed that some DEPs were enriched in glycine, serine, and threonine metabolism. Among the microbial metabolism in diverse environments pathways, the DEPs quantity of L. plantarum KM1 was the largest, and the second one was pyruvate metabolism and carbon metabolism. The results indicated that L. plantarum KM1 responded to the environmental stress of H 2 O 2 by pyruvate metabolism, carbon metabolism, TCA cycle, and microbial metabolism in diverse environments based on KEGG analysis from DEPs.

Carbohydrate Transport and Metabolism of L. plantarum KM1 Under H 2 O 2 Stress
Microorganisms use carbon sources as nutrients to carry out life activities. The transportation and metabolism of carbohydrates are vital to the growth and metabolism of the strain. The phosphotransferase system (PTS) system galactitolspecific EIIC component was upregulated by 2.86-fold after H 2 O 2 stimulation ( Table 2). Beta-phosphoglucomutase (β-PGM) was downregulated by 0.64-fold (Table 2), which catalyzed the mutual conversion of β-D-glucose-1-phosphate and β-Dglucose-6-phosphate, while β-D-glucose-1-phosphate played a role in the biosynthesis and degradation of polysaccharides and was the product of the phosphating effect of maltose and trehalose. Meanwhile, maltodextrin-binding protein MdxE and maltose transport system permease protein MalF were downregulated by 0.35-and 0.63-fold (Table 2), respectively, which participated in the transport and metabolism of maltose. Trehalose import ATP-binding protein SugC involved in the transport and metabolism of trehalose was downregulated by 0.59-fold under the stress of H 2 O 2 .
Phosphoglycerate kinase and Xaa-Pro dipeptidyl-peptidase participated in the metabolic pathways of glycolysis. PGK could catalyze the conversion of glycerate-1, 3 P 2 to glycerate-3 P, which was conducive to produce more ATP in glycolysis  Table 2). Xaa-Pro dipeptidyl-peptidase catalyzed glyceraldehyde-3-phosphate production and provided NADH and H + for the subsequent conversion of pyruvate to lactic acid. It was upregulated by 1.5-fold after stimulation by H 2 O 2 ( Table 2).
Transaldolase and other DEPs indirectly participated in glycolysis and TCA cycle through pentose phosphate pathway (PPP). TAL was downregulated by 0.44-fold in L. plantarum KM1 under H 2 O 2 stress ( The cell should maintain a certain oxidation-reduction homeostasis to resist oxidative stress (Shimizu and Matsuoka, 2019). Homoserine dehydrogenase, anaerobic ribonucleosidetriphosphate reductase, and hypothetical protein all participated in the oxidation-reduction process, which were upregulated by 1.57-, 1.60-, and 5.32-fold under H 2 O 2 stress ( Table 2), respectively. In addition, some DEPs also participated in the oxidation-reduction process and were downregulated by 0.31-0.67 times under H 2 O 2 stress, including HTHtype transcriptional regulator GmuR, D-arabitol-phosphate dehydrogenase, acyl carrier protein, succinate-semialdehyde dehydrogenase [NADP (+)] 1, thioredoxin-like protein YtpP, NAD-dependent malic enzyme, flavodoxin, SsrA-binding protein, and NADH peroxidase ( Table 2). SsrA-binding protein in KM1 slightly decreased under 5 mM H 2 O 2 stimulation. In the process of maintaining redox balance, homoserine dehydrogenase, anaerobic ribonucleoside-triphosphate reductase, and hypothetical protein in L. plantarum KM1 played a major role.

Amino Acid Biosynthesis of L. plantarum KM1 Under H 2 O 2 Stress
Under H 2 O 2 stress, homoserine dehydrogenase and aspartokinase 3 increased by 1.57-and 2.04-fold (Figures 5A,B), respectively. Phosphoserine aminotransferase catalyzed the production of D-aspartate to resist H 2 O 2 , and the enzyme was upregulated by 1.63-fold in L. plantarum KM1 at 5 mM H 2 O 2 ( Figure 5C). Additionally, the TCA cycle was weakened when the pyruvate oxidase and pyruvate dehydrogenase E1 components were downregulated. To improve the adaptability to H 2 O 2 , L. plantarum KM1 upregulated proteins related to amino acid synthesis.

DNA Repair of L. plantarum KM1 Under H 2 O 2 Stress
Oxidative stress causes DNA damage (Preiser, 2012). Anaerobic ribonucleoside-triphosphate reductase [EC:1.1.98.6] in L. plantarum KM1 was upregulated by 1.60-fold under H 2 O 2 stress ( Figure 5D) and catalyzed the production of dGTP for DNA synthesis and repair. Phosphoribosyl-ATP pyrophosphatase for DNA repair and recombination in response to H 2 O 2 -induced oxidative stress were upregulated by 2.27-fold ( Figure 5E). Under H 2 O 2 stress, HTH-type transcriptional regulator AdhR and HTH-type transcriptional repressor GlcR were upregulated by 1.66-and 2.05-fold, respectively (Figures 5F,G). In addition, protein RecA was upregulated by 1.63-fold, which performed homologous recombination DNA repair ( Figure 5H). The Holliday junction resolvase RecU of L. plantarum KM1 increased by 1.57-fold under H 2 O 2 stress, which enabled RecA to catalyze DNA strand exchange ( Figure 5I). Therefore, L. plantarum KM1 could alleviate the DNA damage caused by H 2 O 2 through the DNA repair system.

PPI Networks of DEPs in L. plantarum KM1 Under H 2 O 2 Stress
To observe the interaction of the DEPs of L. plantarum KM1 under H 2 O 2 stress, the PPI network was analyzed (Figure 6). The most obvious PPI network of L. plantarum KM1 was pyruvate metabolism and TCA cycle (pox, pdhA, pdhB, phhC, pdhD). Among them, pyruvate oxidase (pox) played an important role in the detoxification of H 2 O 2 . NADH peroxidase (npr2), NAD-dependent malic enzyme (mae), and SsrA-binding protein (kat) interacted more with other proteins and played an antioxidant role. In addition, the upregulated DEPs related to DNA repair (recA, recU) and amino acid biosynthesis (hom2, thrA2) formed an interaction network that could be involved in reducing the cell damage caused by oxidative stress to L. plantarum KM1, which was consistent with all the above descriptions.

Transcription Level Validation in L. plantarum KM1 Under H 2 O 2 Stress
To verify the reliability of the protein level changes of L. plantarum KM1 under H 2 O 2 stress, the relative expression changes of the genes were measured at the transcriptional level (Figure 7). The results showed that under H 2 O 2 stress, the relative expression levels of 9 genes that played key roles in each pathway were verified to be significantly upregulated. Among them, the gene pepX encoding Xaa-Pro dipeptidyl-peptidase was involved in the glycolytic metabolic pathway and was up-regulated 74.89-fold (p < 0.0001); the genes pox3 (pyruvate oxidase) and pdhA ( (Tang et al., 2017). In addition, L. plantarum needed to have probiotic properties that could tolerate the human gastrointestinal environment before it could exert important antioxidant functions. L. plantarum KM1 grew normally under the condition of pH 2.5, and the survival rate of culturing for 2 h was 96.75%. However, the number of viable cells of L. plantarum L3C1E8 decreased from 8.84 ± 0.23 to 2.25 ± 0.14 log CFU/ml at pH 2.5 for 2 h, and its survival rate was only 25.45% (Ribeiro et al., 2018), which was evidently lower than that of L. plantarum KM1. The survival rate of L. plantarum KM1 cultured at 0.2% bile salt concentration for 2 h was 80.26%, which was significantly higher than L. plantarum RW1 at the same concentration for 3 h (75%)  (Raheem et al., 2021). Antibiotics including kanamycin sulfate, streptomycin sulfate, tetracycline HCl, and ampicillin trihydrate inhibited protein synthesis (Oh and Jung, 2015). Therefore, coadministration with these antibiotics inhibited the growth of L. plantarum KM1 and reduced protein synthesis. Moreover, the MIC of L. plantarum KM1 against the five antibiotics all met the breakpoints provided by EFSA, where kanamycin sulfate was higher than the breakpoint of EFSA (FEEDAP, 2012). The MIC of L. plantarum KM1 on kanamycin sulfate and streptomycin sulfate was consistent with the breakpoints provided by Danielsen and Wind (2003). It was determined by in vitro tests that L. plantarum KM1 had certain probiotic properties and met the screening standards for probiotics. Then, the antioxidant capacity of L. plantarum KM1 was evaluated, and the results showed the antioxidant capacity from strong to weak was fermentation supernatant, cellfree extract, and intact cells. The reducing capacity of L. plantarum KM1 cell-free extracts (310.67 ± 12.019 µmol/L) Therefore, the downregulation of TAL expression increased L. plantarum KM1 resistance to oxidative stress induced by H 2 O 2 . Moreover, the oxidation part of the PPP is the process of converting glucose-6-phosphate into ribulose-5-phosphate, which is also accompanied by NADPH production (Tang, 2019). The putative oxidoreductase acted on NAD + or NADP + acceptor, and its downregulation reduced the use of NAD + or NADP + that increased NADPH production of PPP to reduce cell damage.
The pyruvate produced by glycolysis is oxidatively decarboxylated to acetyl-CoA under the catalysis of the pyruvate dehydrogenase complex and then enters the TCA cycle through acetyl-CoA. Pyruvate metabolism, carbon metabolism and TCA cycle were affected by the oxidation-reduction state (Shetty and Varshney, 2021). Pyruvate metabolism inhibits the production of H 2 O 2 and pyruvate as the final product of glycolysis participated in the TCA cycle (Crestanello et al., 1998;Liu et al., 2009). FAD and TPP were used by pyruvate oxidase as cofactors to catalyze the reaction of pyruvate and phosphate to form H 2 O 2 (Lorquet et al., 2004). The downregulation of it greatly reduced endogenous H 2 O 2 , which was beneficial for the strain to cope with H 2 O 2 stress. The pyruvate dehydrogenase complex mainly converted pyruvate to acetyl-CoA (Grabsztunowicz et al., 2020). The downregulation of DEPs associated with the pyruvate dehydrogenase complex reduced the NADH and ROS production by mitochondrial respiration and maintained cell homeostasis (Anand et al., 2021). However, the TCA cycle captured most of the energy released by glucose oxidation, while downregulation of DEPs on the TCA cycle pathway weakened the TCA cycle and reduced energy production (Fernie et al., 2004). In addition, bacterial growth is inseparable from energy. Furthermore, H 2 O 2 treatment downregulated DEPs associated with the TCA cycle. Therefore, the growth of the strain slowed down significantly with the increase in H 2 O 2 concentration (Figure 1D), which might be due to insufficient energy.
Attenuation of the TCA cycle shifted phosphoenolpyruvate toward amino acid biosynthesis. Studies showed that glycine, serine, and threonine metabolism were closely related to the antioxidant activity (Alhasawi et al., 2015). L. plantarum synthesized amino acids including methionine and D-aspartate to resist oxidative stress (Zhai et al., 2020;Kajitani et al., 2021). Homoserine dehydrogenase and aspartokinase 3 were both involved in the synthesis of L-homoserine, which was the precursor of methionine, and L-homoserine had certain stress resistance (Wei et al., 2019;Liu et al., 2020;Tang et al., 2021). Thus, their marked upregulation increased the ability of the strains to resist oxidative stress. In addition, the weakening of TCA cycle increased the conversion of phosphoenolpyruvate to aspartate and finally improved the homoserine synthesis, which was consistent with previous works .
Finally, L. plantarum KM1 alleviated DNA damage through the DNA repair system. HTH-type transcriptional regulator AdhR and HTH-type transcriptional repressor GlcR belonged to the MerR family. The N-terminal DNA binding region of this family was highly conserved, which played an important role in DNA binding (Brown et al., 2003;Hayashi et al., 2013;Romero-Rodriguez et al., 2015). Particularly, both the protein expression level and the PPI network showed that the protein RecA played a vital role in the DNA repair system, infusing the strain with powerful antioxidant capacity. The expression trend of the gene recA of L. plantarum KM1 at the transcriptional level was consistent with that of its corresponding DEPs (protein RecA). In addition, the expression of gene recA of L. plantarum KM1 in this study was consistent with the corresponding gene transcription level of L. plantarum CAUH2 (Zhai et al., 2020).

CONCLUSION
L. plantarum KM1 with certain intestinal viability was derived from natural fermented food, which exhibited good antioxidant activity in vitro. The proteomics analysis indicated that DEPs participated in various metabolic pathways including pyruvate metabolism, carbon metabolism, TCA cycle, amino acid metabolism, and microbial metabolism in diverse environments in response to oxidative stress caused by H 2 O 2 in L. plantarum KM1. In addition, glycolysis and DNA repair also played an important role in dealing with H 2 O 2 stress.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors. The data presented in the study are deposited in the ProteomeXchange repository, accession number PXD032910. The names of the repository/repositories and accession number(s) can be found at: https://www.ebi.ac.uk/pride/archive, PXD032910.

AUTHOR CONTRIBUTIONS
YT: conceptualization, investigation, formal analysis, writingoriginal draft, and preparation. YuW, JZ, and BN: writingreview and editing. NZ, MX, XX, YZ, and YF: technical support. XL, YuhW, and JL: conceptualization, resources, supervision, funding acquisition, project administration, and writing-review and editing. All authors contributed to the article and approved the submitted version.