# Redox proteomic study of Bacillus cereus thiol proteome during fermentative anaerobic growth

Fella Hamitouche, Jean-Charles Gaillard, Philippe Schmitt, J. Armengaud,

Catherine Duport, Luc Dedieu

► To cite this version:

Fella Hamitouche, Jean-Charles Gaillard, Philippe Schmitt, J. Armengaud, Catherine Duport, et al. Redox proteomic study of Bacillus cereus thiol proteome during fermentative anaerobic growth. BMC Genomics, 2021, 22 (1), pp.648. 10.1186/s12864-021-07962-y. http-03499925

HAL Id: hal-03499925

https://hal.inrae.fr/hal-03499925

Submitted on 21 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Redox proteomic study of *Bacillus cereus* thiol proteome during fermentative anaerobic growth

Fella Hamitouche<fnref n="1" />, Jean-Charles Gaillard², Philippe Schmitt<fnref n="1" />, Jean Armengaud², Catherine Duport<fnref n="1" />† and Luc Dedieu<fnref n="1" />*†

## Abstract

**Background:** *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.

**Results:** 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.

**Conclusions:** 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.

**Keywords:** *Bacillus cereus*, Thiol proteome, Anaerobiosis, Fermentative growth, Differential labeling strategy

## Background

*Bacillus cereus* is a Gram positive bacterium, which is recognized as a major foodborne pathogen responsible for two types of syndrome: emetic and diarrheal [1]. 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 [4], pH, oxygen level, and oxidoreduction potential (ORP)

[2]. *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 [5]. *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 [6]. 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 [7]. Consequently, both the NAD and NADP systems

<footnote>* Correspondence: luc.dedieu@univ-avignon.fr</footnote>

<footnote>†Catherine Duport and Luc Dedieu should be considered as co-last authors.</footnote>

<footnote n="1">¹Avignon Université, INRAE, UMR SQPOV, Site Agroparc, F-84914 Avignon Cedex 9, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
play major roles in redox homeostasis [8]. Low-molecular weight compounds, such as coenzyme A, free cysteines, and bacillithiol (BSH) are also major contributors to redox homeostasis [9], 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 [10]. 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 [10].

The redox status of Cys SH in proteins contributes significantly to protein folding [11], metal binding [12], and regulating protein function [10]. 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-SO₂H) and sulfonic acid (R-SO₃H) species [10]. The set of proteins containing reversibly modified thiols is referred to as the thiol proteome [13], or redoxome [14].

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 [15], and monitoring protein cysteine oxidation by loss of reactivity with thiol reagents [16], to accurately identifying sites of thiol redox modifications in numerous proteins [17]. 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 [17]. 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 [14]. At the proteome level, reduced and oxidized Cys can be distinguished thanks to the use of different alkylating agents [13]. They can thus be quantified and the abundance of the proteins bearing these residues can be estimated [20].

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* [23], *Saccharomyces cerevisiae* [24], *Escherichia coli* [25, 26], and Firmicutes including

*Bacillus subtilis* [27], 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 [28]; alcohol dehydrogenase (AdhA), which is involved in energy metabolism [23]; and alkyl hydroperoxide reductase (AhpC), a key player in defensive mechanisms [21].

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.

## Results

### 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 [31].

Reactive Cys residues were labeled using a three-step sequential strategy as previously described [32], 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
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Growth curve and alterations to extracellular oxidoreduction potential (ORP) of *B. cereus* ATCC 14579 in regulated anoxic batch cultures. Circles represent optical densities at 600 nm (OD₆₀₀) and squares correspond to extracellular ORP. Black arrows indicate cell harvesting, at early-exponential growth phase (EEP), mid-exponential growth phase (MEP), and late-exponential growth phase (LEP). Data correspond to mean values calculated for biological triplicates ± SD</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Schematic workflow used to investigate the *B. cereus* thiol proteome. 1/Samples were harvested from the bioreactor at early- (EEP), mid-MEP) and late-exponential growth phase (LEP); 2/Proteins were extracted and thiol-disulfide exchanges were quenched by addition of TCA; 3/Denatured proteins were obtained; 4/Proteins with reactive Cys residues were labeled using two strategies. IDN) Reduced-Cys thiol (SH) were blocked with IAM; oxidized Cys thiols (S-OX) were then reduced with DTT and alkylated with NEM. IN) Reduced-Cys thiol were alkylated with IAM and remaining non-alkylated Cys thiols were blocked with NEM; 5/Proteins were digested with trypsin; 6/Peptides were analyzed by LC-MS/MS; 7/Modified Cys were identified; 8/ Proteins were identified and the thiol proteome profiled. TCA: trichloroacetic acid; IAM: iodoacetamide; DTT: dithiothreitol; NEM: N-ethylmaleimide</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Labeling efficiency of IDN assay compared to IN control. IDN: reduced Cys were labeled with IAM and reversibly-oxidized Cys were reduced with DTT prior to labeling with NEM. IN: reduced Cys were first labeled with IAM and then with NEM without reduction between steps. The labeled Cys-containing peptides and Cys residues were quantified in samples harvested at early-exponential phase (EEP), mid-exponential phase (MEP) and late-exponential phase (LEP). Data correspond to mean values calculated for triplicates. Error bars show standard deviation. Paired Student t-test was performed to determine statistical significance (*: 0.01 < *p* < 0.05; **: 0.001 < *p* < 0.01; ***: *p* < 0.001)</figcaption></figure>

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 [33]. 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
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**Fig. 4** Proportion of reduced Cys (RedCys) and oxidized Cys (OxCys) residues identified in samples harvested at early-exponential, mid-exponential and late-exponential growth phases (EEP, MEP, and LEP, respectively). A Number of modified cysteines: RedCys and OxCys in samples. B Proportion of non-redundant RedCys and non-redundant OxCys residues. Sp.RedCys/Sp.Oxcys: reduced/oxidized cysteine residues that have been detected in only one growth stage (these residues are specific to a growth stage); S.RedCys/S.OxCys: reduced/oxidized cysteine residues that have been detected in two or three growth stages. Data correspond to mean values calculated for three replicates. Error bars show standard deviation. Paired Student t-test was performed to determine statistical significance (*:0.01 < *p* < 0.05; **: 0.001 < *p* < 0.01; ***: *p* < 0.001)

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: Cys<sup>255</sup> was mainly detected as RedCys at MEP, whereas the two others (Cys<sup>97</sup> at EEP and Cys<sup>91</sup> at MEP) were mainly detected as OxCys.
<table><caption>**Table 1** Numbers of non-redundant peptides, Cys-containing peptides and Cys residues identified in all the three replicates at early-exponential phase (EEP), mid-exponential phase (MEP) and late-exponential phase (LEP)</caption><thead><tr><th></th><th>EEP</th><th>MEP</th><th>LEP</th></tr></thead><tbody><tr><td>Non-redundant peptides</td><td>3673</td><td>3330</td><td>3087</td></tr><tr><td>Non-redundant Cys-peptides</td><td>109</td><td>95</td><td>88</td></tr><tr><td>Non-redundant Cys residues</td><td>110</td><td>99</td><td>93</td></tr></tbody></table>

Several proteins were categorized in functional groups related to defense mechanisms. Superoxide dismutase SodA2 (Q814I6) contains a Cys residue (Cys⁶⁰) 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, Cys³⁵⁷, which are mainly detected as RedCys both at MEP and LEP; and Cys⁴⁵⁸, which is mainly detected as RedCys both at EEP and LEP. Thiol peroxidase (Tpx, Q817B8) contains two redox-sensitive Cys residues (Cys⁶⁰ and Cys⁹⁴), indifferently detected as RedCys and OxCys throughout growth. Finally, alkyl hydroperoxide reductase C (AhpC, Q81IK9) contains one Cys residue (Cys¹⁶⁶), 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: Cys¹⁷⁰ in PheT1 (Q812N7), Cys¹²⁴ in IleS1 (Q819R4) and Cys¹⁵³ in ArgS (Q814Q8). The elongation factor FusA (Q814C5) contains a Cys residue (Cys²⁵⁸) that was mainly detected as RedCys at EEP. The elongation factor Tuf (Q814C4) was associated with two redox-sensitive Cys residues (Cys⁸² and Cys¹³⁸), 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: Cys³⁰⁸, which contributes to the active site, indifferently detected as RedCys and OxCys throughout growth, whereas Cys⁴⁴⁶ was detected as RedCys at all three growth phases, and Cys³²⁷ 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 [34]. 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, Cys⁴⁷³ was indifferently detected as RedCys and OxCys. DhbB harbors two Cys residues, Cys⁶⁷ and Cys¹⁹⁷, indifferently detected as RedCys and OxCys. In contrast, Cys¹⁴ in DhbC was mainly detected as RedCys.

## Discussion

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 [35]. 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 [17]. In a previous study, Shakir et al. [25] 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)
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Growth-phase-distribution of the 74 Cys residues detected predominantly either as reduced Cys (RedCys) or as oxidized Cys (OxCys). **A** Venn diagram showing the growth-phase-distribution of the 74 Red/OxCys residues. **B** Heat map visualization of the 72 peptides containing the 74 Red/OxCys residues. EEP, early-exponential growth phase; MEP, mid-exponential growth phase; LEP, late-exponential growth phase. Colors represent the redox form that predominates. Blue: oxidized form; red: reduced form; gray: no significant difference</figcaption></figure>

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%) [23]. 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) [37]. 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 [38].

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 [29]. 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* [39]. 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.
<table><caption>**Table 2** Cys-proteins harboring specific cysteines detected as RedCys or as OxCys</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Cys position</th><th rowspan="2">Protein</th><th rowspan="2">Description</th><th colspan="2">Log<sub>2</sub>(FC)</th></tr><tr><th>EEP/LEP</th><th>MEP/LEP</th></tr></thead><tbody><tr><td colspan="6"><strong>Cys proteins in EEP,MEP,LEP</strong></td></tr><tr><td colspan="6"><strong>Nucleotide metabolism</strong></td></tr><tr><td>Q81IP8</td><td>381</td><td>PurD</td><td>Phosphoribosylamine-glycine ligase</td><td>NS</td><td>NS</td></tr><tr><td>Q81JC8</td><td>446</td><td>GuaB</td><td>Inosine-5'-monophosphate dehydrogenase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Defense mechanisms</strong></td></tr><tr><td>Q81IE6</td><td>123</td><td>TerD</td><td>Tellurium resistance protein</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Translation</strong></td></tr><tr><td>Q812N7</td><td>170</td><td>PheT_1</td><td>Phenylalanine-tRNA ligase beta subunit</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cys proteins in EEP,MEP</strong></td></tr><tr><td colspan="6"><strong>Cofactors and vitamins metabolism</strong></td></tr><tr><td>Q81J74</td><td>60</td><td></td><td>7,8-dihydroneopterin aldolase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Translation</strong></td></tr><tr><td>Q819R4</td><td>124</td><td>IleS1</td><td>Isoleucine-tRNA ligase 1</td><td>NS</td><td>NS</td></tr><tr><td>Q814Q8</td><td>153</td><td>ArgRS</td><td>Arginine-tRNA ligase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Amino acid metabolism</strong></td></tr><tr><td>Q814Q2</td><td>159</td><td>SpeB</td><td>Agmatinase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cys proteins in EEP,LEP</strong></td></tr><tr><td colspan="6"><strong>Defense mechanisms</strong></td></tr><tr><td>Q81GN6</td><td>458</td><td>Kat</td><td>Catalase</td><td>-2,46</td><td>-1,93</td></tr><tr><td colspan="6"><strong>Cys proteins in MEP,LEP</strong></td></tr><tr><td colspan="6"><strong>Nucleotide metabolism</strong></td></tr><tr><td>Q817X3</td><td>64</td><td>Apt</td><td>Adenine phosphoribosyltransferase</td><td>NS</td><td>NS</td></tr><tr><td>Q81IQ3</td><td><u>464</u></td><td>PurL</td><td>Phosphoribosylformylglycinamidine synthase subunit</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Carbohydrate metabolism</strong></td></tr><tr><td>Q818T4</td><td>111</td><td></td><td>2-oxoisovalerate dehydrogenase beta subunit</td><td>NS</td><td>NS</td></tr><tr><td>Q81IA3</td><td>123</td><td></td><td>Formate acetyltransferase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cofactors and vitamins metabolism</strong></td></tr><tr><td>Q818X5</td><td><u>92</u></td><td>RibH</td><td>6,7-dimethyl-8-ribityllumazine synthase</td><td>NS</td><td>NS</td></tr><tr><td>Q818R9</td><td>405</td><td>Dxs</td><td>1-deoxy-D-xylulose-5-phosphate synthase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Amino acid metabolism</strong></td></tr><tr><td>Q817H8</td><td><u>189</u></td><td></td><td>Deblocking aminopeptidase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Lipid and Fatty acid metabolism</strong></td></tr><tr><td>Q81IN1</td><td>38</td><td></td><td>DAGKc domain-containing protein</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Pathogenesis</strong></td></tr><tr><td>Q81I63</td><td><u>184</u></td><td></td><td>Microbial collagenase</td><td>NS</td><td>NS</td></tr><tr><td>Q81JC7</td><td>266</td><td></td><td>D-alanyl-D-alanine carboxypeptidase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Defense mechanisms</strong></td></tr><tr><td>Q81GN6</td><td>357</td><td>Kat</td><td>Catalase</td><td>-2,46</td><td>-1,93</td></tr><tr><td colspan="6"><strong>Transcription</strong></td></tr><tr><td>Q81J66</td><td><u>312</u></td><td>MecB</td><td>Negative regulator of genetic competence</td><td>NS</td><td>NS</td></tr></tbody></table>
<table><caption>**Table 2** Cys-proteins harboring specific cysteines detected as RedCys or as OxCys (Continued)</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Cys position</th><th rowspan="2">Protein</th><th rowspan="2">Description</th><th colspan="2">Log<sub>2</sub>(FC)</th></tr><tr><th>EEP/LEP</th><th>MEP/LEP</th></tr></thead><tbody><tr><td colspan="6"><strong>Cys proteins in EEP</strong></td></tr><tr><td colspan="6"><strong>Nucleotide metabolism</strong></td></tr><tr><td>Q81IQ6</td><td>169</td><td>PurC</td><td>Phosphoribosylaminoimidazole-succinocarboxamide synthase</td><td>NS</td><td>NS</td></tr><tr><td>Q819S3</td><td>1034</td><td>CarB</td><td>Carbamoyl-phosphate synthase large chain</td><td>NS</td><td>NS</td></tr><tr><td>Q819S0</td><td>264</td><td>PyrB</td><td>Aspartate carbamoyltransferase</td><td>NS</td><td>NS</td></tr><tr><td>Q81FV5</td><td>38</td><td>DeoD</td><td>Purine nucleoside phosphorylase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Translation</strong></td></tr><tr><td>Q81IN3</td><td>185</td><td>GatA</td><td>Glutamyl-tRNA (Gln) amidotransferase subunit A</td><td>NS</td><td>NS</td></tr><tr><td>Q814C5</td><td>258</td><td>FusA</td><td>Elongation factor G</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Transcription</strong></td></tr><tr><td>Q81J47</td><td>892</td><td>RpoC</td><td>DNA-directed RNA polymerase subunit beta</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Amino acid metabolism</strong></td></tr><tr><td>Q816W9</td><td>315</td><td></td><td>Deblocking aminopeptidase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Defense mechanisms</strong></td></tr><tr><td>Q81IE5</td><td>110</td><td>TerD</td><td>Tellurium resistance protein</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Function unknown</strong></td></tr><tr><td>Q81IG4</td><td>65</td><td></td><td>Thiamine_BP domain-containing protein</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Carbohydrate metabolism</strong></td></tr><tr><td>Q81DX6</td><td>97</td><td>AdhA</td><td>Alcohol dehydrogenase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cys proteins in MEP</strong></td></tr><tr><td colspan="6"><strong>Nucleotide metabolism</strong></td></tr><tr><td>Q81J22</td><td>153</td><td>Adk</td><td>Adenylate kinase</td><td>NS</td><td>NS</td></tr><tr><td>Q81IS3</td><td>250</td><td>GuaA</td><td>GMP synthase</td><td>NS</td><td>NS</td></tr><tr><td>Q81IQ3</td><td>527</td><td>PurL</td><td>Phosphoribosyl formylglycinamidine synthase subunit</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Translation</strong></td></tr><tr><td>Q814G7</td><td>15</td><td>RpsR</td><td>30S ribosomal protein S18</td><td>NS</td><td>NS</td></tr><tr><td>Q81J29</td><td>27</td><td>RpsZ</td><td>30S ribosomal protein S14 type Z</td><td>NS</td><td>NS</td></tr><tr><td>Q814C4</td><td>82</td><td>Tuf</td><td>Elongation factor Tu</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Transcription</strong></td></tr><tr><td>Q81J47</td><td>818</td><td>RpoC</td><td>DNA-directed RNA polymerase subunit beta</td><td>NS</td><td>NS</td></tr><tr><td>Q81CE5</td><td>158</td><td>PpaC</td><td>Probable manganese-dependent inorganic pyrophosphatase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Amino acid metabolism</strong></td></tr><tr><td>Q81GG4</td><td>115</td><td>TrpA</td><td>Tryptophan synthase alpha chain</td><td>-3.08</td><td>-2.14</td></tr><tr><td colspan="6"><strong>Carbohydrate metabolism</strong></td></tr><tr><td>Q818P6</td><td>194</td><td></td><td>Phosphoglycolate phosphatase</td><td>NS</td><td>NS</td></tr><tr><td>Q81DX6</td><td>91/255</td><td>AdhA</td><td>Alcohol dehydrogenase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cofactors and vitamins metabolism</strong></td></tr><tr><td>Q818X3</td><td>115</td><td>Biob</td><td>Biotin synthase</td><td>NS</td><td>NS</td></tr></tbody></table>
<table><caption>**Table 2** Cys-proteins harboring specific cysteines detected as RedCys or as OxCys (Continued)</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Cys position</th><th rowspan="2">Protein</th><th rowspan="2">Description</th><th colspan="2">Log<sub>2</sub>(FC)</th></tr><tr><th>EEP/LEP</th><th>MEP/LEP</th></tr></thead><tbody><tr><td colspan="6"><strong>Protein folding</strong></td></tr><tr><td>Q819C5</td><td>73</td><td></td><td>Peptidyl-prolyl cis-trans isomerase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cys proteins in LEP</strong></td></tr><tr><td colspan="6"><strong>Amino acid metabolism</strong></td></tr><tr><td>Q81GG5</td><td>378</td><td>TrpB</td><td>Tryptophan synthase beta chain</td><td>-2,24</td><td>-2,26</td></tr><tr><td>Q81GH0</td><td>168/65</td><td>TrpE</td><td>Anthraniolate synthase component 1</td><td>-2,17</td><td>-2,17</td></tr><tr><td>Q81GG7</td><td>108</td><td>TrpC</td><td>Indole-3-glycerol phosphate synthase</td><td>-2,55</td><td>-2,26</td></tr><tr><td>Q81GG9</td><td>62</td><td>PabA (TrpGD)</td><td>Anthraniolate synthase component II</td><td>-1,93</td><td>-1,39</td></tr><tr><td>Q814M2</td><td>364</td><td></td><td>O-acetylhomoserine sulfhydrylase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Carbohydrate metabolism</strong></td></tr><tr><td>Q819M9</td><td>262</td><td></td><td>Pyruvate carboxylase</td><td>NS</td><td>NS</td></tr><tr><td>Q81FZ7</td><td>124</td><td></td><td>Glycerate dehydrogenase</td><td>NS</td><td>NS</td></tr><tr><td>Q81A22</td><td>50</td><td>PorB</td><td>Pyruvate synthase subunit</td><td>NS</td><td>NS</td></tr><tr><td>Q811A2</td><td>12</td><td>PflA</td><td>Pyruvate formate lyase activating enzyme</td><td>NS</td><td>NS</td></tr><tr><td>Q817F3</td><td>283</td><td>PfkA</td><td>ATP-dependent 6-phosphofructokinase</td><td>NS</td><td>NS</td></tr><tr><td>Q819D7</td><td>365</td><td>Ptsi</td><td>Phosphoenolpyruvate-protein phosphotransferase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Cofactors and vitamins metabolism</strong></td></tr><tr><td>Q81JC6</td><td>128/130</td><td>PdxS</td><td>Pyridoxal 5'-phosphate synthase lyase</td><td>NS</td><td>NS</td></tr><tr><td>Q816H9</td><td>71</td><td>MenB</td><td>1,4-dihydroxy-2-naphthoyl-CoA synthase</td><td>NS</td><td>NS</td></tr><tr><td>Q818W9</td><td>69</td><td>BioD</td><td>ATP-dependent dethiobiotin synthetase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Transcription</strong></td></tr><tr><td>Q81J18</td><td>265</td><td>RpoA</td><td>DNA-directed RNA polymerase subunit alpha</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Translation</strong></td></tr><tr><td>Q814C4</td><td>138</td><td>Tuf</td><td>Elongation factor Tu</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Nucleotide metabolism</strong></td></tr><tr><td>Q81JC8</td><td>327</td><td>GuaB</td><td>Inosine-5'-monophosphate dehydrogenase</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Defense mechanisms</strong></td></tr><tr><td>Q81IK9</td><td>166</td><td>AhpC</td><td>Alkyl hydroperoxide reductase C</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Lipid and Fatty acid metabolism</strong></td></tr><tr><td>Q814Y7</td><td>91</td><td>FabZ</td><td>3-hydroxyacyl-[acyl-carrier-protein] dehydratase</td><td>NS</td><td>NS</td></tr><tr><td>Q816B6</td><td>46/49</td><td>NfuA (YutI)</td><td>NifU protein</td><td>-1,9</td><td>-1,68</td></tr><tr><td>Q816D2</td><td>117</td><td></td><td>HESB protein</td><td>NS</td><td>NS</td></tr><tr><td colspan="6"><strong>Secondary metabolites biosynthesis, transport and catabolism</strong></td></tr><tr><td>Q81DQ3</td><td>14</td><td>DhbC</td><td>Isochorismate synthase</td><td>-2,71</td><td>-3,06</td></tr><tr><td>Q81DQ1</td><td>67</td><td>DhbB</td><td>Isochorismatase</td><td>-2,55</td><td>-2,52</td></tr><tr><td colspan="6"><strong>Energy production and conversion</strong></td></tr><tr><td>Q814U7</td><td>135</td><td></td><td>Nitrilotriacetate monooxygenase component B</td><td>NS</td><td>NS</td></tr></tbody></table>

Protein identifiers are underlined for proteins with differences in abundance; Cysteine positions are underlined for Cys-peptides detected with a higher proportion of RedCys at LEP; FC fold change; NS not significant
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Functional categorization of reactive Cys-containing proteins at early (EEP), mid (MEP) and late (LEP) exponential growth stages</figcaption></figure>

Overall, redox-sensitive Cys residues were mainly detected in their reduced forms, according to the reducing environment of cytoplasm [35], 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 [44]. 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, Cys⁹¹ and Cys⁹⁷, near a Zn²⁺ binding site, which is assumed to promote oxidation [45].

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 [48]. 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 [49].

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 [50]. Interestingly, the modifications observed on Cys⁸² and Cys¹³⁸ 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, Cys⁸² in Tuf was identified as a conserved S-bacillithiolated site among *Bacillus spp.* and *Staphylococcus carnosus* [54].

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 [57], catalase is produced to deal with oxidative stress [58], and the enzymes involved in bacillibactin biosynthesis circumvent iron starvation [56]. 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 [59].

## Conclusions

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.

## Methods

### Chemicals

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, pO₂ = 0%) in MOD medium supplemented with 30 mM glucose [6]. Briefly, *B. cereus* ATCC 14579 strain [60] 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 pO₂ 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 [29]. Continuous sparging of pure N₂ (10 L/h) was maintained to generate anoxic conditions (pO₂ = 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 [61].

## 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 [63] 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 [64] 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 × 10⁴ 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 [65]. 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 [66], under dataset identifiers PXD022049 and https://doi.org/10.6019/PXD022049. [The dataset can be accessed with the Username: reviewer_pxd022049@ebi.ac.uk and Password: NFX3eY3n].

## Data analysis

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 [67]. 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 [68].

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s12864-021-07962-y.

**Additional file 1: Table S1.** 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).

**Additional file 2: Figure S3.** Changes in *B. cereus* proteome composition between EEP, MEP, and LEP obtained by Gene Ontology analysis of pathway enrichment.

## Acknowledgments

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.

## Authors' contributions

CD and LD conceived and designed the experiments. FH and JCG performed the experiments. JA, JCG, and PS contributed (reagents/materials/ analysis tools). FH, JCG, JA, CD, and LD analyzed the data. FH, CD, and LD wrote the paper. All authors reviewed the manuscript. The authors read and approved the final manuscript.

## Funding

This work was supported by Avignon-Université and the Algerian government.

## 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: reviewer_pxd022049@ebi.ac.uk and Password: NFX3eY3n.

## Declarations

### Ethics approval and consent to participate

'Not applicable'.

### Consent for publication

'Not applicable'.

### Competing interests

The authors declare that they have no competing interests.

### Author details

<footnote n="1">¹Avignon Université, INRAE, UMR SQPOV, Site Agroparc, F-84914 Avignon Cedex 9, France. ²Université Paris-Saclay, CEA, INRAE, Département Médicaments et Technologies pour la Santé (DMTS), SPI, 30200 Bagnols-sur-Cèze, France.</footnote>
Received: 23 March 2021 Accepted: 5 August 2021

Published online: 07 September 2021

## References

1. 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.

2. Duport C, Jobin M, Schmitt P. Adaptation in *Bacillus cereus*: from stress to disease. Front Microbiol. 2016;7:1550.

3. 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.

4. Diomande SE, Guinebretière M-H, Broussolle V, Brillard J. Role of fatty acids in *Bacillus* environmental adaptation. Front Microbiol. 2015;6:813.

5. 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.

6. 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.

7. Agledal L, Niere M, Ziegler M. The phosphate makes a difference: cellular functions of NADP. Redox Rep. 2010;15(1):2-10.

8. Xiao W, Wang RS, Handy DE, Localzo J. NAD(H) and NADP(H) redox couples and cellular energy metabolism. Antioxid Redox Signal. 2018;28(3): 251-72.

9. Reniere ML. Reduce, Induce, Thrive: Bacterial Redox Sensing during Pathogenesis. J Bacteriol. 2018;200(17):e00128-18.

10. Poole LB. The basics of thiols and cysteines in redox biology and chemistry. Free Radic Biol Med. 2015;80:148-57.

11. Noiva R. Enzymatic catalysis of disulfide formation. Protein Expr Purif. 1994; 5(1):1-13.

12. 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.

13. Yang J, Carroll KS, Liebler DC. The expanding landscape of the thiol redox proteome. Mol Cell Proteomics. 2016;15(1):1-11.

14. 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.

15. Ellman GL. Tissue sulfhydryl groups. Arch Biochem Biophys. 1959;82(1):70-7.

16. 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.

17. Wojdyla K, Rogowska-Wrzesinska A. Differential alkylation-based redox proteomics—lessons learnt. Redox Biol. 2015;6:240-52.

18. 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.

19. 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.

20. 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.

21. 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.

22. 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.

23. 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.

24. 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.

25. 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.

26. 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.

27. 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).

28. Imber M, Pietrzyk-Brzezinska AJ, Antelmann H. Redox regulation by reversible protein S-thiolation in gram-positive bacteria. Redox Biol. 2019;20:130-45.

29. 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.

30. 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.

31. Chiappetta G, Ndiaye S, Igbaria A, Kumar C, Vinh J, Toledano MB. Proteome Screens for Cys Residues Oxidation. Methods Enzymol. 2010;473:199-216.

32. 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.

33. 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.

34. 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.

35. Leichert LI, Jakob U. Protein thiol modifications visualized in vivo. PLoS Biol. 2004;2(11):e333.

36. 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.

37. 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.

38. 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.

39. Jacob H. Redox potential. Methods Microbiol. 1970;2:91-123.

40. 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.

41. 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.

42. Fowler NJ, Blanford CF, de Visser SP, Warwick J. Features of reactive cysteines discovered through computation: from kinase inhibition to enrichment around protein degrons. Sci Rep. 2017;7(1):16338.

43. 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.

44. Stewart MD, Iqumenova TI. Reactive cysteine in the structural Zn2+ site of the C1B domain from PKCa. Biochemistry. 2012;51(37):7263-77.

45. 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.

46. Loi VV, Rossius M, Antelmann H. Redox regulation by reversible protein S-thiolation in bacteria. Front Microbiol. 2015;6:187.

47. 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.

48. Go Y-M, Chandler JD, Jones DP. The cysteine proteome. Free Radic Biol Med. 2015;84:227-45.

49. 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.

50. 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.

51. 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.
52. Kraal B, Lippmann C, Kleanthous C. Translational regulation by modifications of the elongation factor Tu. Folia Microbiol. 1999;44(2):131-41.

53. 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.

54. 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.

55. Jaishankar J, Srivastava P. Molecular basis of stationary phase survival and applications. Front Microbiol. 2017;8:2000.

56. 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.

57. Fontecave M. Iron-sulfur clusters: ever-expanding roles. Nat Chem Biol. 2006; 2(4):171-4.

58. 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.

59. 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.

60. 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.

61. Zwietering M, Jongenburger I, Rombouts F, Van't Riet K. Modeling of the bacterial growth curve. Appl Environ Microbiol. 1990;56(6):1875-81.

62. 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.

63. 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).

64. 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.

65. 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.

66. 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.

67. 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.

68. 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.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_10.png" /></figure>