- Open Access
Redox proteomic study of Bacillus cereus thiol proteome during fermentative anaerobic growth
BMC Genomics volume 22, Article number: 648 (2021)
Bacillus cereus is a notorious foodborne pathogen, which can grow under anoxic conditions. Anoxic growth is supported by endogenous redox metabolism, for which the thiol redox proteome serves as an interface. Here, we studied the cysteine (Cys) proteome dynamics of B. cereus ATCC 14579 cells grown under fermentative anoxic conditions. We used a quantitative thiol trapping method combined with proteomics profiling.
In total, we identified 153 reactive Cys residues in 117 proteins participating in various cellular processes and metabolic pathways, including translation, carbohydrate metabolism, and stress response. Of these reactive Cys, 72 were detected as reduced Cys. The B. cereus Cys proteome evolved during growth both in terms of the number of reduced Cys and the Cys-containing proteins identified, reflecting its growth-phase-dependence. Interestingly, the reduced status of the B. cereus thiol proteome increased during growth, concomitantly to the decrease of extracellular oxidoreduction potential.
Taken together, our data show that the B. cereus Cys proteome during unstressed fermentative anaerobic growth is a dynamic entity and provide an important foundation for future redox proteomic studies in B. cereus and other organisms.
Bacillus cereus is a Gram positive bacterium, which is recognized as a major foodborne pathogen responsible for two types of syndrome: emetic and diarrheal . As well as being a facultative anaerobe, B. cereus can adapt to a wide range of environmental conditions allowing it to multiply in a number of food products and in the human intestine [2, 3]. In these environments, B. cereus adapts its metabolism to variations in temperature , pH, oxygen level, and oxidoreduction potential (ORP) . B. cereus uses aerobic respiratory pathways to grow under aerobiosis, and mixed acid fermentation pathways to grow under anaerobiosis in the absence of an external electron acceptor. The main product of fermentation is lactate, which is synthesized alongside acetate, ethanol, formate, succinate, and small amounts of 2, 3-butanediol . B. cereus re-oxidizes the reducing equivalent NADH, generated by glucose catabolism from NAD+, through the respiratory chain under aerobiosis, and during the formation of end-products under anaerobiosis . Thus, NAD+/ NADH is central to catabolism and energy supply, whereas the NADP+/ NADPH couple plays an important role in biosynthesis and detoxification of cells . Consequently, both the NAD and NADP systems play major roles in redox homeostasis . Low-molecular weight compounds, such as coenzyme A, free cysteines, and bacillithiol (BSH) are also major contributors to redox homeostasis , in particular through their regulation of the redox status of the two amino-acids containing thiol (SH) groups, methionine (Met) and cysteine (Cys). Together, Cys and Met fulfill important roles in cells thanks to their redox chemistry . Although the two residues are subject to redox regulation, the thioether form of the Met sulfur atom is less reactive than the sulfhydryl (thiol, SH) form of the Cys sulfur atom .
The redox status of Cys SH in proteins contributes significantly to protein folding , metal binding , and regulating protein function . Thus, oxidation of Cys SH groups can result in the formation of reversible modifications such as disulfide bridges (S-S), sulfenic acid (S-OH) and nitrosylation (−SNO) groups, and irreversible oxidations such as sulfinic (R-SO2H) and sulfonic acid (R-SO3H) species . The set of proteins containing reversibly modified thiols is referred to as the thiol proteome , or redoxome .
Attempts to study the thiol proteome have taken several approaches. Due to their high lability, cysteine modifications are challenging to analyze. Methods to identify Cys redox modifications have progressed over the years, ranging from the estimation of SH groups using colorimetric tests , and monitoring protein cysteine oxidation by loss of reactivity with thiol reagents , to accurately identifying sites of thiol redox modifications in numerous proteins . The differential thiol labeling method constituted a major advance in the depth of analysis possible. The first step involves irreversibly blocking reduced-Cys thiols with an alkylating agent, subsequently the oxidized thiols are reduced, and finally labeled with a different alkylating agent . This method has been combined with 2D gel-based methods [18, 19] to identify several redox-sensitive proteins. However, this combination presented limitations, mainly in the identification of cysteines. Today, the combination of differential thiol labeling with tandem mass spectrometry is a powerful means to overcome the limitations of previous methods and to provide broad coverage of the thiol proteome . At the proteome level, reduced and oxidized Cys can be distinguished thanks to the use of different alkylating agents . They can thus be quantified and the abundance of the proteins bearing these residues can be estimated .
Large-scale redox proteomic analysis has been used to investigate the relationships between the thiol proteome status in (i) cyanobacteria, focusing on photosynthesis and the response to nutrient limitation [21, 22], and (ii) Clostridium difficile , Saccharomyces cerevisiae , Escherichia coli [25, 26], and Firmicutes including Bacillus subtilis , to investigate responses to oxidative stress. These studies revealed that the thiol proteome modifies key biological processes through alterations to redox-sensitive proteins such as the elongation factor EF-Tu (Tuf), a key component of translation ; alcohol dehydrogenase (AdhA), which is involved in energy metabolism ; and alkyl hydroperoxide reductase (AhpC), a key player in defensive mechanisms .
Previous studies investigated the redox status of Met residues in unstressed B. cereus cells at both the cellular proteome and exoproteome levels [29, 30]. However, the thiol proteome of unstressed B. cereus cells has not yet been investigated. In this study, we used a thiol trapping method combined with shotgun proteomics analyses to decipher the B. cereus thiol proteome during fermentative anaerobic growth. Our results show that the B. cereus thiol proteome is growth-phase-dependent and contains higher number of reduced Cys in the later stages of growth.
Labeling strategy for global analysis of the B. cereus thiol proteome
Fermentative growth of B. cereus ATCC 14579 is illustrated in Fig. 1. It was accompanied by a decrease in extracellular ORP (ΔORP = 217 ± 20 mV). Culture samples were collected anaerobically, in TCA-containing vials at early- (EEP), mid- (MEP) and late-exponential (LEP) growth phases. Mixing with TCA resulted in immediate lysis of cells and precipitation of proteins, thus preserving the in vivo redox status of proteins containing Cysteine (Cys) residues .
Reactive Cys residues were labeled using a three-step sequential strategy as previously described , and named IDN in our study (Fig. 2): First, iodoacetamide (IAM) was used to label and block free Cys thiol groups (SH). Then, reversibly-oxidized thiol (S-OX) were reduced by exposure to DTT, and finally the corresponding thiol groups were labeled using N-ethylmaleimide (NEM). The IN control omitted the DTT reduction step, and was used to assess how efficiently IAM blocked Cys SH. All proteins were digested with trypsin, and the resulting peptides were submitted to extensive LC-MS/MS analysis. A total of 354,684 peptides were identified from the 18 samples (3 cultures × 3 time-points × 2 labels), including 9531 Cys-containing peptides. A total of 20,952 ± 5265, 20,274 ± 4541, and 16,760 ± 2880 peptides were identified at EEP, MEP, and LEP, respectively. These values are not significantly different according to Student’s T-test.
The number of IAM-labeled Cys residues (IAMCys) was not significantly different in IDN-treated samples and IN controls, whatever the growth phase (Fig. 3). In contrast, as expected, the number of NEM- labeled Cys residues (NEMCys) detected was higher in IDN samples than in the IN control, for samples harvested during the EEP and MEP growth phases. Taken together, these results indicate that IAM alkylation stabilizes free sulfhydryl groups, but that a small number of Cys residues escape IAM alkylation and remain available for alkylation with NEM.
B. cereus thiol proteome dynamics
NEMCys in IDN samples correspond to either in vivo oxidized Cys or insensitive IAM residues. To retain only in vivo oxidized Cys residues, we corrected the IDN data by subtracting NEMCys detected in the IN dataset from the list of NEMCys detected in the IDN dataset. The corrected IDN dataset (IDNc) is shown in Table S1.
Analysis of the IDNc dataset shows that the numbers of RedCys (labeled with IAM) and OxCys (labeled with NEM after reduction) were significantly different at both MEP and LEP (Fig. 4A). In addition, the number of OxCys residues decreased during growth and the number of RedCys tend to increase as growth progresses. These data suggest that the thiol proteome evolved toward a highly reduced status as growth progressed. Due to trypsin missed cleavage and differences in abundance, some redox-sensitive Cys residues specified the same cysteine site. Removal of this redundancy resulted in a total of 110, 99 and 93 non-redundant (unique) redox-sensitive Cys residues at EEP, MEP and LEP respectively (Table 1). The number of the non-redundant redox-sensitive Cys residues decreased during growth due to an overall decrease in non-redundant peptide detection. This probably reflects the overall regulation of protein synthesis during growth . However, the labeling efficiency remained stable (3 ± 0.1%, Table 1).
We compared the numbers of RedCys and OxCys of each non-redundant Cys residue. Based on this analysis, we identified 74 Cys residues for which the mean of RedCys number was significantly different from the mean of OxCys number (p < 0.05). Figure 5 shows their growth stage distribution (Fig. 5A), their RedCys vs OxCys distribution and the peptide sequences in which they were identified (Fig. 5B). The data indicate that (i) most of the 74 Cys residues were identified in only one growth phase and only two Cys residues had a higher number of OxCys compared to RedCys: one at EEP and the other at MEP. (ii) LEP samples supported the largest number of Cys residues, and (iii) the 74 Cys residues were contained in 72 non-redundant peptides.
If we consider all the 74 Cys residues, we do not observe any change in the total number of OxCys number during growth. In contrast, we observe a significant increase in the total number of RedCys at MEP and LEP, compared to EEP (Fig. 4B). This overall increase results from an increase in the number of RedCys residues identified in samples from at least two growth stages (hatched bars in Fig. 4B), and RedCys residues identified specifically in samples from only one growth stage (plain bars in Fig. 4B).
Components of the B. cereus thiol proteome
According to our results, the B. cereus thiol proteome includes 153 reactive Cys-containing peptides, of which 70 were detected mainly in the reduced form, and 2 in the oxidized form. We assigned these 153 peptides to their corresponding proteins. The resulting set consisted of 117 proteins (Table S2), of which 64 contained the 74 Cys residues identified above (Table 2). These proteins were mainly involved in translation, amino acid, nucleotide and carbohydrate metabolism, as well as defense mechanisms, whatever the growth stage (Fig. 6). Proteins related to carbohydrate metabolism include the glycolytic enzymes phosphofructokinase PfkA (Q817F3) and phosphoenolpyruvate PtsI (Q819D7), the fermentative enzyme pyruvate formate-lyase PflA (Q81IA2), and alcohol dehydrogenase (AdhA, Q81DX6), which also plays an important role in fermentative metabolism. Three distinct modified cysteines were identified in AdhA: Cys255 was mainly detected as RedCys at MEP, whereas the two others (Cys97 at EEP and Cys91 at MEP) were mainly detected as OxCys.
Several proteins were categorized in functional groups related to defense mechanisms. Superoxide dismutase SodA2 (Q814I6) contains a Cys residue (Cys60) identified as redox-sensitive at MEP, this protein was more abundant during LEP than during EEP (Table S2). The abundance of catalase (Kat, Q81GN6) increased during the LEP compared to both EEP and MEP, it harbors two redox-sensitive residues, Cys357, which are mainly detected as RedCys both at MEP and LEP; and Cys458, which is mainly detected as RedCys both at EEP and LEP. Thiol peroxidase (Tpx, Q817B8) contains two redox-sensitive Cys residues (Cys60 and Cys94), indifferently detected as RedCys and OxCys throughout growth. Finally, alkyl hydroperoxide reductase C (AhpC, Q81IK9) contains one Cys residue (Cys166), which was identified as redox-sensitive in the three growth stages, and mainly detected as RedCys at LEP.
Among the proteins categorized in the translation-related functional group, three are known to bind tRNA. These proteins harbored RedCys residues at EEP and MEP: Cys170 in PheT1 (Q812N7), Cys124 in IleS1 (Q819R4) and Cys153 in ArgS (Q814Q8). The elongation factor FusA (Q814C5) contains a Cys residue (Cys258) that was mainly detected as RedCys at EEP. The elongation factor Tuf (Q814C4) was associated with two redox-sensitive Cys residues (Cys82 and Cys138), which were mainly detected as RedCys at MEP and LEP, respectively.
Several proteins involved in nucleotide metabolism were identified in our study, these included PurD (Q81IP8), Adk (Q81J22) and GuaB (Q81JC8). Interestingly, the three Cys residues present in Inosine-5′-monophosphate dehydrogenase (GuaB) were found to be redox-sensitive at all three growth stages: Cys308, which contributes to the active site, indifferently detected as RedCys and OxCys throughout growth, whereas Cys446 was detected as RedCys at all three growth phases, and Cys327 was detected mainly as RedCys specifically at LEP.
Five proteins involved in amino acid metabolism were also identified. TrpA (Q81GG4), TrpB (Q81GG5), TrpC (Q81GG7), and PabA (Q81GG9). These proteins were more abundant in samples harvested at LEP compared to both EEP and MEP. Anthranilate synthase component I TrpE (Q81GH0) was more abundant at LEP than at MEP (Table 2).
In addition to these classes of proteins, we identified RedCys residues in iron-sulfur cluster proteins such as NifU (Q816B6), which is involved in the biogenesis of bacterial Fe/S proteins, and is thus responsible for several functions, including redox catalysis . The abundance of NifU increased significantly during LEP compared to EEP and MEP.
Interestingly, based on Cys-peptides and spectral count analysis, we detected three proteins involved in the assembly of the bacillibactin siderophore during the LEP. These proteins, the isochorismate synthase DhbC (Q81DQ3), the isochorismatase DhbB (Q81DQ1) and the 2,3-dihydroxybenzoate-AMP ligase DhbE (Q81DQ2) were only detected in samples harvested during the LEP. In DhbE, Cys473 was indifferently detected as RedCys and OxCys. DhbB harbors two Cys residues, Cys67 and Cys197, indifferently detected as RedCys and OxCys. In contrast, Cys14 in DhbC was mainly detected as RedCys.
Thiol redox modifications play important roles in regulating proteins that sustain bacterial growth. However, despite its importance, and mainly due to technical challenges, few studies have attempted to quantify changes in the bacterial thiol redox proteome during growth. The first hurdle relates to the extreme sensitivity of Cys residues to artefactual oxidation. For our study, we circumvented this problem by performing rapid quenching with TCA . A second challenge is to use an efficient strategy to profile thiol redox modifications in the whole organism. We adopted a strategy without peptide enrichment steps to avoid introducing quantitative bias which could lead to deviations from physiological conditions, as we wished to conserve the situation in the cell at the time of harvesting as far as possible. To achieve our aims, we implemented a differential thiol trapping technique, the IDN strategy, which relies on irreversible modifications of reduced Cys residues by IAM, and oxidized Cys residues by NEM after DTT treatment [32, 36]. This strategy does not make it possible to characterize the OxCys forms, and to establish if they are part of an inter- or intra-polypeptide disulfide.
In any differential labeling strategy, the first alkylation – to block free thiols – must be highly efficient. Indeed, incomplete alkylation would increase the number of Cys available for NEM-labeling, leading to incorrect identification of oxidized residues . In a previous study, Shakir et al.  used a sulfhydryl reactive dye (DyLight 550 maleimide) to estimate the proportion of free thiols after IAM labeling as a means to assess the efficiency of IAM alkylation. Here, we included a control sample (IN) to demonstrate and precisely quantify how efficiently IAM blocked free thiols. We were thus able to correct the data and obtain accurate quantifications of the NEMCys target. Through our approach, we identified 153 reactive Cys residues in B. cereus thiol proteome. This relatively low number of redox-target Cys may be explained by the overall low Cys content of proteins in B. cereus (less than 1%) . Like the other facultative anaerobic Gram positive Firmicutes, B. cereus probably promotes cysteine exclusion from both its cytoplasmic and exported proteins due to their high sensitivity to reactive oxygen and nitrogen species (ROS and RNS) . Indeed, any such sensitivity could pose problems for an organism that can survive in a wide variety of, sometimes hostile, environments. Another possibility is that by restricting our analysis to Cys residues identified in all three replicates at each stage of growth to ensure the data generated were reliable and reproducible, we may have excluded some redox-sensitive Cys residues. However, our results are of the same order of magnitude as those reported for Caenorhabditis elegans, for which fewer than 200 redox-sensitive Cys residues were identified during its lifespan .
Our study revealed that the number of detected reduced Cys residues increased as growth progressed. Extracellular ORP decreased at the same time, thus B. cereus growth under anaerobiosis is accompanied by a significant decrease in extracellular ORP, especially between the beginning and the end of exponential growth . The ability to decrease extracellular ORP is common to many bacteria, and may be related to the consumption of oxidizing compounds or the production of reducing end-products, as reported for E. coli . In Lactococcus lactis and Listeria monocytogenes the decrease in extracellular ORP has been shown to depend on the presence of reduced thiol groups present on proteins located at the bacterial cell surface [40, 41]. If the number of reduced thiol groups of these proteins increases as growth progresses, as in the case of the thiol groups in the proteins identified in this study, it could drive the reduction in extracellular ORP. However, this hypothesis will need to be proven through future research.
Overall, redox-sensitive Cys residues were mainly detected in their reduced forms, according to the reducing environment of cytoplasm , and the extracellular milieu in our condition. Only two Cys residues were detected more frequently as OxCys than as RedCys. In addition to being dependent of cellular conditions, the reactivity of Cys residues is partly controlled by the intrinsic properties of the sulfur atom and by the position of the Cys residue in the three-dimensional protein structure [42, 43]. Reactivity can also be modulated by proximity to a metal ion . Interestingly, the two Cys residues, which were mainly detected as OxCys in our samples both belonged to alcohol dehydrogenase AdhA (Q81DX6). The predicted three-dimensional structure of the protein (data not shown) localizes these two residues, Cys91 and Cys97, near a Zn2+ binding site, which is assumed to promote oxidation .
Several redox-sensitive Cys-containing proteins identified in this study (Adk, PpaC, PurL, GuaB) were previously reported to undergo thiol modifications under oxidizing conditions in Bacillus subtilis, Staphylococcus aureus, Corynebacterium glutamicum, and E. coli [27, 28, 46, 47]. These reports suggest that these proteins could undergo specific redox regulation in response to variations in environmental conditions.
Our results showed that the high detection of RedCys residues at LEP was due to differences in abundance for 8 Cys-containing proteins (including 3 Cys-containing proteins that were detected exclusively at LEP) and the detection of 21 Cys-containing proteins with a higher number of RedCys at LEP without abundance variation. Therefore, the growth-phase-specificity of B. cereus thiol proteome may be explained both by a change of proteome composition and a distinct redox sensitivity of Cys residues in individual proteins . The results of the Panther Gene Ontology analysis (http://geneontology.org/) illustrated the composition changes of B. cereus proteome (Figure S3). This analysis showed an enrichment of proteome in proteins involved in siderophore biosynthesis, and tryptophan biosynthesis at LEP, at the expense of proteins involved in deoxyribonucleotide biosynthesis. No proteome change characterizing the onset of sporulation has been detected, according to previous work .
Many abundant Cys-proteins, such as those involved in carbohydrate metabolism (PtsI and PfkA) and translation (FusA and Tuf), harbored cysteines displaying an increasing RedCys residues as growth progressed. These results suggest that these proteins may be redox-regulated to ensure efficient glucose uptake and protein biosynthesis over the course of growth. Moreover, redox regulation of PtsI and PfkA could potentially divert the metabolic flux toward post-translational modifications . Interestingly, the modifications observed on Cys82 and Cys138 in the elongation factor Tuf at MEP, and LEP, respectively, could reflect two different sites of potential redox regulation during anaerobic growth. A growth-phase-dependent modification of Cys residues in Tuf could be the result of a growth-phase-dependent regulation network, as previously described [51, 52]. The redox-sensitive Cys residues in the Tuf protein were previously described to be S-bacillithiolated and S-mycothiolated in response to hypochlorite stress [28, 53]. In addition, Cys82 in Tuf was identified as a conserved S-bacillithiolated site among Bacillus spp. and Staphylococcus carnosus .
Our data showed that some redox-sensitive Cys-containing proteins were detected during LEP, due to their high abundance at this stage of growth, such as Kat, NifU, DhbC, and DhbB, which had at least one Cys detected mainly as RedCys during LEP. From a physiological point of view, these data support the entry of bacteria into the stationary phase [55, 56] in response to stress and nutrient starvation. The Fe-S-containing protein NifU is involved in sensing iron and superoxide levels , catalase is produced to deal with oxidative stress , and the enzymes involved in bacillibactin biosynthesis circumvent iron starvation . To our knowledge, this is the first study to demonstrate the presence of reactive cysteine residues in DhbC and DhbB. Redox regulation of DhbB (Q81DQ1) remains uncharacterized, but might be important since this enzyme uses a thiolation mechanism to activate siderophore assembly .
In summary, we presented an overview of the B. cereus ATCC 14579 thiol proteome, which is mainly defined by a subset of redox-sensitive Cys residues differing from one growth stage to another, and functionally related to growth-stage-dependent cellular events. This works represents an essential step for researchers in the fields of bacterial physiology toward understanding and predicting the effect of changing redox conditions.
Trichloroacetic acid (TCA), dithiothreitol (DTT), Trifluoroacetic acid (TFA), and iodoacetamide (IAM) were purchased from Sigma (St-Quentin-Fallavier, France). N-ethylmaleimide (NEM) and bicinchoninic acid assay (BCA) were purchased from Fisher Scientific (France). Trypsin and ProteaseMax reagent used for proteolysis were purchased from Promega (France).
Regulated batch culture
B. cereus was grown in regulated batch conditions (pH 7, temperature = 37 °C, pO2 = 0%) in MOD medium supplemented with 30 mM glucose . Briefly, B. cereus ATCC 14579 strain  was grown in batch cultures under anoxic conditions in a 3-L-capacity bioreactor (My-control, Applikon technology) equipped with a pH electrode (405-DPAS-SC-K8S/225, Mettler Toledo), a pO2 polarographic oxygen electrode (Z010023525, Applisens), and an AgCl ORP electrode (pT4805-SC-DPAS-K8S/225, Mettler Toledo). All batch culture experiments were conducted in 2 L MOD medium supplemented with 30 mM glucose, as described previously . Continuous sparging of pure N2 (10 L/h) was maintained to generate anoxic conditions (pO2 = 0%). Cultures were stirred at 300 rpm, and temperature was set to 37 °C. The pH of the cell culture medium was maintained at 7 by adding 1 M HCl and 3 M KOH.
B. cereus cultures were performed in biological triplicates and the optical density (OD) at 600 nm was measured every 30 min. ORP values were corrected based on the value recorded by the reference electrode (+ 207 mV at 37 °C). The specific maximum growth rate (μmax) was determined using the Zwietering growth model .
Protein extraction: quenching cellular thiol-disulfide exchange
Culture extracts were harvested anaerobically during the early-exponential growth phase (EEP) (at μmax), mid-exponential growth phase (MEP), and late-exponential growth phase (LEP) in vials containing 25% trichloroacetic acid (TCA). Extracts were incubated with TCA overnight at − 20 °C to maximize protein precipitation. Proteins were pelleted by centrifugation (13,800 x g, 15 min, 4 °C) and washed twice with ice-cold acetone to remove excess TCA. Supernatants were discarded, and protein pellets were air-dried.
Differential thiol trapping method
Differential thiol labeling was performed as described previously [32, 62], with minor modifications. Briefly, protein extracts resuspended in denaturing buffer were divided between two tubes, labeled IDN, and IN (Fig. 2). Free cysteine residues were alkylated with 50 mM iodoacetamide (IAM), at room temperature in the dark for 1 h. Proteins were precipitated once again with 25% (w/v) TCA and washed twice with ice-cold acetone. Protein pellets were resuspended in the same buffer, and the IN samples were incubated with 150 mM N-ethylmaleimide (NEM). IDN samples were first reduced by adding 25 mM dithiothreitol (DTT), and incubating at room temperature for 1 h. Reduced proteins were TCA precipitated (25%), and resuspended before alkylation with 150 mM NEM.
The protein concentration of each sample was estimated using the bicinchoninic acid assay (BCA), according to the manufacturer’s protocol (Thermo-Fisher). A 90-μg aliquot of proteins was loaded onto a 12% SDS-PAGE for a short migration at 90 V (15–20 min). Gels were stained with imperial™ protein stain (Thermo-Fisher) for 1 h and destained with water. Protein bands were excised from gels. IDN and IN samples were digested with trypsin using the ProteaseMax surfactant (Promega) as described elsewhere  without further reduction/alkylation. For all samples, the resulting peptide mixtures were dissolved in 0.1% trifluoroacetic acid (TFA) prior to nano LC-MS/MS analyses.
High-resolution tandem mass spectrometry
Peptides were identified using a Q-Exactive HF mass spectrometer (Thermo Scientific) coupled to an ultimate 3000 nano LC system (Thermo Scientific). The system was essentially operated as previously described  with minor modifications. The peptide mixtures (10 μL) were loaded and desalted online on a reverse-phase precolumn (Acclaim Pepmap 100 C18 5 μm bead size, 100-Å pore size, 5 mm × 300 μm, Thermo). They were then resolved on a nanoscale Acclaim Pepmap 100 C18 column (3 μm bead size, 100-Å pore size, 15 cm × 75 μm) at a flow rate of 200 nl/min using a 120-min gradient combining buffer B (0.1% HCOOH, 80% CH3CN) and buffer A (0.1% HCOOH, 100% H2O): 4–25% B in 100 min, followed by 25–40% B in 20 min. The mass spectrometer was operated in the Top20 data-dependent acquisition mode with full MS scans acquired from 350 to 1800 m/z at 60,000 resolution, and after each scan, selection of the 20 most abundant precursor ions for fragmentation and MS/MS acquisition at 15,000 resolution. An intensity threshold of 9 × 104 was applied. A 10-s dynamic exclusion was used to increase the detection of low-abundance peptides. Only double- and triple-charged ions were selected for MS/MS analysis.
MS/MS spectra interpretation
MS/MS spectra were queried against the theoretical B. cereus ATCC 14579 proteome (5216 sequences) with the Mascot Daemon software, version 2.6.1 (Matrix Science). The parameters were: full trypsin specificity, only 2+ and 3+ peptide charges, a mass tolerance of 5 ppm on the parent ion, a mass tolerance of 0.02 Da on the MS/MS ions, a maximum of two missed cleavages, no static modifications, and dynamic modifications were N-ethylmaleimide (C), Carbamidomethyl (C), Carbamyl (K), Oxidation (M), Deamidation (NQ). All peptide matches with a peptide score associated with a p-value of less than 0.05 were parsed using IRMa 1.31.1c software . Proteins were considered valid when at least two different peptides were detected in the same sample. The proteomics datasets were submitted to the ProteomeXchange Consortium via the PRIDE partner repository , under dataset identifiers PXD022049 and https://doi.org/10.6019/PXD022049. [The dataset can be accessed with the Username: firstname.lastname@example.org and Password: NFX3eY3n].
All experiments were replicated three times. For each Cys residues, we compared the mean value of the number of reduced Cys residues with the mean value of the number of oxidized Cys residues. Statistical analysis was performed using the paired Student t-test, setting the threshold for statistical significance at p less than 0.05. For proteome profiling, data were normalized relative to the total spectral count. Protein abundances were compared between two growth stages, as previously described . Proteins for which abundances changed by more than 1.5 were considered significant when the p-value was less than 0.05. Protein functions were assigned based on information from the Clusters of Orthologous Groups (COG) database .
Availability of data and materials
All data and materials are fully available without restriction. The proteomics datasets were submitted to the ProteomeXchange Consortium via the PRIDE partner repository [Project accession: PXD022049; Project DOI: https://doi.org/10.6019/PXD022049] and can be accessed at the URL [https://www.ebi.ac.uk/pride/login] with the Username: email@example.com and Password: NFX3eY3n.
Ehling-Schulz M, Fricker M, Scherer S. Identification of emetic toxin producing Bacillus cereus strains by a novel molecular assay. FEMS Microbiol Lett. 2004;232(2):189–95.
Duport C, Jobin M, Schmitt P. Adaptation in Bacillus cereus: from stress to disease. Front Microbiol. 2016;7:1550.
Jeßberger N, Rademacher C, Krey VM, Dietrich R, Mohr A-K, Böhm M-E, et al. Simulating intestinal growth conditions enhances toxin production of enteropathogenic Bacillus cereus. Front Microbiol. 2017;8:627.
Diomande SE, Guinebretière M-H, Broussolle V, Brillard J. Role of fatty acids in Bacillus environmental adaptation. Front Microbiol. 2015;6:813.
Rosenfeld E, Duport C, Zigha A, Schmitt P. Characterization of aerobic and anaerobic vegetative growth of the food-borne pathogen Bacillus cereus F4430/73 strain. Can J Microbiol. 2005;51(2):149–58.
Zigha A, Rosenfeld E, Schmitt P, Duport C. Anaerobic cells of Bacillus cereus F4430/73 respond to low oxidoreduction potential by metabolic readjustments and activation of enterotoxin expression. Arch Microbiol. 2006;185(3):222–33.
Agledal L, Niere M, Ziegler M. The phosphate makes a difference: cellular functions of NADP. Redox Rep. 2010;15(1):2–10.
Xiao W, Wang RS, Handy DE, Loscalzo J. NAD(H) and NADP(H) redox couples and cellular energy metabolism. Antioxid Redox Signal. 2018;28(3):251–72.
Reniere ML. Reduce, Induce, Thrive: Bacterial Redox Sensing during Pathogenesis. J Bacteriol. 2018;200(17):e00128–18.
Poole LB. The basics of thiols and cysteines in redox biology and chemistry. Free Radic Biol Med. 2015;80:148–57.
Noiva R. Enzymatic catalysis of disulfide formation. Protein Expr Purif. 1994;5(1):1–13.
Giles NM, Watts AB, Giles GI, Fry FH, Littlechild JA, Jacob C. Metal and redox modulation of cysteine protein function. Chem Biol. 2003;10(8):677–93.
Yang J, Carroll KS, Liebler DC. The expanding landscape of the thiol redox proteome. Mol Cell Proteomics. 2016;15(1):1–11.
Chiappetta G, Ndiaye S, Igbaria A, Kumar C, Vinh J, Toledano MB. Proteome screens for Cys residues oxidation: the redoxome. Methods Enzymol. : Elsevier. 2010;473:199–216.
Ellman GL. Tissue sulfhydryl groups. Arch Biochem Biophys. 1959;82(1):70–7.
Lee S-R, Kwon K-S, Kim S-R, Rhee SG. Reversible inactivation of protein-tyrosine phosphatase 1B in A431 cells stimulated with epidermal growth factor. J Biol Chem. 1998;273(25):15366–72.
Wojdyla K, Rogowska-Wrzesinska A. Differential alkylation-based redox proteomics--lessons learnt. Redox Biol. 2015;6:240–52.
Le Moan N, Clement G, Le Maout S, Tacnet F, Toledano MB. The Saccharomyces cerevisiae proteome of oxidized protein thiols contrasted functions for the thioredoxin and glutathione pathways. J Biol Chem. 2006;281(15):10420–30.
Baty JW, Hampton MB, Winterbourn CC. Detection of oxidant sensitive thiol proteins by fluorescence labeling and two-dimensional electrophoresis. Proteomics: International Edition. 2002;2(9):1261–6.
Wojdyla K, Williamson J, Roepstorff P, Rogowska-Wrzesinska A. The SNO/SOH TMT strategy for combinatorial analysis of reversible cysteine oxidations. J Proteome. 2015;113:415–34.
Guo J, Nguyen AY, Dai Z, Su D, Gaffrey MJ, Moore RJ, et al. Proteome-wide light/dark modulation of thiol oxidation in cyanobacteria revealed by quantitative site-specific redox proteomics. Mol Cell Proteomics. 2014;13(12):3270–85.
Ansong C, Sadler NC, Hill EA, Lewis MP, Zink EM, Smith RD, et al. Characterization of protein redox dynamics induced during light-to-dark transitions and nutrient limitation in cyanobacteria. Front Microbiol. 2014;5:325.
Sievers S, Dittmann S, Jordt T, Otto A, Hochgräfe F, Riedel K. Comprehensive redox profiling of the thiol proteome of Clostridium difficile. Mol Cell Proteomics. 2018;17(5):1035–46.
Brandes N, Reichmann D, Tienson H, Leichert LI, Jakob U. Using quantitative redox proteomics to dissect the yeast redoxome. J Biol Chem. 2011;286(48):41893–903.
Shakir S, Vinh J, Chiappetta G. Quantitative analysis of the cysteine redoxome by iodoacetyl tandem mass tags. Anal Bioanal Chem. 2017;409(15):3821–30.
Leichert LI, Gehrke F, Gudiseva HV, Blackwell T, Ilbert M, Walker AK, et al. Quantifying changes in the thiol redox proteome upon oxidative stress in vivo. Proc Natl Acad Sci. 2008;105(24):8197–202.
Chi BK, Gronau K, Mäder U, Hessling B, Becher D, Antelmann H. S-bacillithiolation protects against hypochlorite stress in Bacillus subtilis as revealed by transcriptomics and redox proteomics. Mol Cell Proteomics. 2011;10(11).
Imber M, Pietrzyk-Brzezinska AJ, Antelmann H. Redox regulation by reversible protein S-thiolation in gram-positive bacteria. Redox Biol. 2019;20:130–45.
Madeira JP, Alpha-Bazin B, Armengaud J, Duport C. Time dynamics of the Bacillus cereus exoproteome are shaped by cellular oxidation. Front Microbiol. 2015;6:342.
Madeira J-P, Armengaud J, Duport C. Methionine residues in exoproteins and their recycling by methionine sulfoxide reductase AB serve as an antioxidant strategy in Bacillus cereus. Front Microbiol. 2017;8:1342.
Chiappetta G, Ndiaye S, Igbaria A, Kumar C, Vinh J, Toledano MB. Proteome Screens for Cys Residues Oxidation. Methods Enzymol. 2010;473:199–216.
Kuo Y-H, Konopko AM, Borotto NB, Majmudar JD, Haynes SE, Martin BR. Profiling protein S-sulfination with maleimide-linked probes. Chembiochem. 2017;18(20):2028.
Bernhardt J, Weibezahn J, Scharf C, Hecker M. Bacillus subtilis during feast and famine: visualization of the overall regulation of protein synthesis during glucose starvation by proteome analysis. Genome Res. 2003;13(2):224–237.
Blanc B, Gerez C, Ollagnier de Choudens S. Assembly of Fe/S proteins in bacterial systems: biochemistry of the bacterial ISC system. Biochim Biophys Acta. 2015;1853(6):1436–47.
Leichert LI, Jakob U. Protein thiol modifications visualized in vivo. PLoS Biol. 2004;2(11):e333.
Hansen RE, Winther JR. An introduction to methods for analyzing thiols and disulfides: reactions, reagents, and practical considerations. Anal Biochem. 2009;394(2):147–58.
Daniels R, Mellroth P, Bernsel A, Neiers F, Normark S, von Heijne G, et al. Disulfide bond formation and cysteine exclusion in gram-positive bacteria. J Biol Chem. 2010;285(5):3300–9.
Knoefler D, Thamsen M, Koniczek M, Niemuth NJ, Diederich AK, Jakob U. Quantitative in vivo redox sensors uncover oxidative stress as an early event in life. Mol Cell. 2012;47(5):767–76.
Jacob H. Redox potential. Methods Microbiol. 1970;2:91–123.
Michelon D, Abraham S, Ebel B, De Coninck J, Husson F, Feron G, et al. Contribution of exofacial thiol groups in the reducing activity of Lactococcus lactis. FEBS J. 2010;277(10):2282–90.
Pillot G, Brillet-Viel A, Prévost H. Implication of exofacial thiol groups in the reducing activity of Listeria monocytogenes. bioRxiv. 2018:353409. https://doi.org/10.1101/353409.
Fowler NJ, Blanford CF, de Visser SP, Warwicker J. Features of reactive cysteines discovered through computation: from kinase inhibition to enrichment around protein degrons. Sci Rep. 2017;7(1):16338.
Netto LE, de Oliveira MA, Monteiro G, Demasi AP, Cussiol JR, Discola KF, et al. Reactive cysteine in proteins: protein folding, antioxidant defense, redox signaling and more. Comp Biochem Physiol Toxicol Pharmacol. 2007;146(1–2):180–93.
Stewart MD, Igumenova TI. Reactive cysteine in the structural Zn2+ site of the C1B domain from PKCα. Biochemistry. 2012;51(37):7263–77.
Tamarit J, Cabiscol E, Ros J. Identification of the major oxidatively damaged proteins in Escherichia coli cells exposed to oxidative stress. J Biol Chem. 1998;273(5):3027–32.
Loi VV, Rossius M, Antelmann H. Redox regulation by reversible protein S-thiolation in bacteria. Front Microbiol. 2015;6:187.
Xie K, Bunse C, Marcus K, Leichert LI. Quantifying changes in the bacterial thiol redox proteome during host-pathogen interaction. Redox Biol. 2019;21:101087.
Go Y-M, Chandler JD, Jones DP. The cysteine proteome. Free Radic Biol Med. 2015;84:227–45.
Abbas AA, Planchon S, Jobin M, Schmitt P. Absence of oxygen affects the capacity to sporulate and the spore properties of Bacillus cereus. Food Microbiol. 2014;42:122–31.
Ralser M, Wamelink MM, Latkolik S, Jansen EE, Lehrach H, Jakobs C. Metabolic reconfiguration precedes transcriptional regulation in the antioxidant response. Nat Biotechnol. 2009;27(7):604–5.
Talavera A, Hendrix J, Versées W, Jurėnas D, Van Nerom K, Vandenberk N, et al. Phosphorylation decelerates conformational dynamics in bacterial translation elongation factors. Sci Adv. 2018;4(3):eaap9714.
Kraal B, Lippmann C, Kleanthous C. Translational regulation by modifications of the elongation factor Tu. Folia Microbiol. 1999;44(2):131–41.
Chi BK, Busche T, Van Laer K, Bäsell K, Becher D, Clermont L, et al. Protein S-mycothiolation functions as redox-switch and thiol protection mechanism in Corynebacterium glutamicum under hypochlorite stress. Antioxid Redox Signal. 2014;20(4):589–605.
Chi BK, Roberts AA, Huyen TTT, Bäsell K, Becher D, Albrecht D, et al. S-bacillithiolation protects conserved and essential proteins against hypochlorite stress in firmicutes bacteria. Antioxid Redox Signal. 2013;18(11):1273–95.
Jaishankar J, Srivastava P. Molecular basis of stationary phase survival and applications. Front Microbiol. 2017;8:2000.
Lee JY, Passalacqua KD, Hanna PC, Sherman DH. Regulation of petrobactin and bacillibactin biosynthesis in Bacillus anthracis under iron and oxygen variation. PLoS One. 2011;6(6):e20777.
Fontecave M. Iron-sulfur clusters: ever-expanding roles. Nat Chem Biol. 2006;2(4):171–4.
Heo YJ, Chung IY, Cho WJ, Lee BY, Kim JH, Choi KH, et al. The major catalase gene (katA) of Pseudomonas aeruginosa PA14 is under both positive and negative control of the global transactivator OxyR in response to hydrogen peroxide. J Bacteriol. 2010;192(2):381–90.
May JJ, Wendrich TM, Marahiel MA. The dhb operon of Bacillus subtilis encodes the biosynthetic template for the catecholic siderophore 2,3-dihydroxybenzoate-glycine-threonine trimeric ester bacillibactin. J Biol Chem. 2001;276(10):7209–17.
Ivanova N, Sorokin A, Anderson I, Galleron N, Candelon B, Kapatral V, et al. Genome sequence of Bacillus cereus and comparative analysis with Bacillus anthracis. Nature. 2003;423(6935):87–91.
Zwietering M, Jongenburger I, Rombouts F, Van't Riet K. Modeling of the bacterial growth curve. Appl Environ Microbiol. 1990;56(6):1875–81.
Wholey W-Y, Jakob U. Hsp33 confers bleach resistance by protecting elongation factor Tu against oxidative degradation in vibrio cholerae. Mol Microbiol. 2012;83(5):981–91.
de Groot A, Dulermo R, Ortet P, Blanchard L, Guérin P, Fernandez B, et al. Alliance of proteomics and genomics to unravel the specificities of Sahara bacterium Deinococcus deserti. PLoS Genet. 2009;5(3).
Klein G, Mathé C, Biola-Clier M, Devineau S, Drouineau E, Hatem E, et al. RNA-binding proteins are a major target of silica nanoparticles in cell extracts. Nanotoxicology. 2016;10(10):1555–64.
Dupierris V, Masselon C, Court M, Kieffer-Jaquinod S, Bruley C. A toolbox for validation of mass spectrometry peptides identification and generation of database: IRMa. Bioinformatics. 2009;25(15):1980–1.
Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47(D1):D442–D50.
Rousset L, Alpha-Bazin B, Château A, Armengaud J, Clavel T, Berge O, et al. Groundwater promotes emergence of asporogenic mutants of emetic Bacillus cereus. Environ Microbiol. 2020;22(12):5248–64.
Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015;43(D1):D261–D9.
The authors would like to thank the Algerian government for supporting Fella Hamitouche’s thesis fellowship. We are also grateful to Bénédicte Doublet for technical assistance and to Joëlle Vinh and Vincent Méjean for helpful discussions.
This work was supported by Avignon-Université and the Algerian government.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Non-redundant Cys-peptides identified in samples harvested at EEP, MEP, and LEP, in the IDNc dataset. Table S2. List of 117 Cys-proteins identified. From column Q to V: abundance levels of Cys-proteins (log2 values of spectral counts for each protein with the adjusted p-value).
Changes in B. cereus proteome composition between EEP, MEP, and LEP obtained by Gene Ontology analysis of pathway enrichment.
About this article
Cite this article
Hamitouche, F., Gaillard, JC., Schmitt, P. et al. Redox proteomic study of Bacillus cereus thiol proteome during fermentative anaerobic growth. BMC Genomics 22, 648 (2021). https://doi.org/10.1186/s12864-021-07962-y
- Bacillus cereus
- Thiol proteome
- Fermentative growth
- Differential labeling strategy