Abstract
The ability of chaperonins to buffer mutations that affect protein folding pathways suggests that their abundance should be evolutionarily advantageous. Here, we investigate the effect of chaperonin overproduction on cellular fitness in Escherichia coli. We demonstrate that chaperonin abundance confers 1) an ability to tolerate higher temperatures, 2) improved cellular fitness, and 3) enhanced folding of metabolic enzymes, which is expected to lead to enhanced energy harvesting potential.
Introduction
Chaperonins are found in nearly every organism across all domains of life, and are essential in all cases tested to date, although in some cases non-essential paralogues are found (; ). The GroE chaperonin system of E. coli, consisting of the 60 kDa GroEL and the 10 kDa GroES proteins assembled into ring complexes of 14 and seven sub-units, respectively, is encoded by the groE operon (Tilly and Georgopoulos, 1982; ; ). This operon is expressed principally from two promoters, one utilized in the presence of housekeeping sigma factor σ70, and the other, which is strongly induced due to the accumulation of unfolded proteins, in the presence of the alternative sigma factor, σ32 (RpoH) (; ; ; ). As σ32 levels respond to unfolded protein, this provides a feedback loop to maintain proteostasis (). When cells are shifted to heat shock temperatures between 42 and 46°C, GroEL levels increase by 5–10 fold, reaching up to 12% of the entire cellular proteome (). These increased levels interact more extensively with the proteome and are assumed to prevent misfolding or assist refolding of heat-stressed proteins (; ; ). Cells that cannot mount an unfolded protein response due to rpoH deletion are extremely temperature sensitive, and selection for pseudo-revertants of these strains at elevated temperatures yields up-promoter mutations in the groE promoter (). GroE is thus important even under normal growth conditions, and indeed GroEL and GroES are respectively the 20th and 21st most abundant proteins in E. coli (excluding ribosomal proteins), with sufficient protein being made under non-stressed conditions to produce approximately 2,800 complexes of GroEL and 5,700 complexes of GroES (). Other chaperones that are also abundant include the ribosome bound trigger factor (TF), which is the 19th, and the Hsp70 homologue, DnaK, which is the 27th most abundant. The high levels of all these chaperones indicates their key roles in cell growth. Although combined loss of TF and DnaK is deleterious to cells, groEL and groES are the only chaperone encoding genes in E. coli that are essential under all conditions ().
GroE (GroEL and GroES) assists the folding of 10–15% cellular proteins (), many of which are essential (). GroE’s ability to fold “folding-compromised” proteins (; ; ; Tokuriki and Tawfik, 2009) is consistent with a “genetic capacitance” function. Many studies with different heterologous proteins have shown that GroE can enhance their folding (Tokuriki et al., 2008; Tokuriki and Tawfik, 2009; Wyganowski et al., 2013; ; ). In addition, some deleterious mutations are retained in the genome upon overexpression of groE, probably due to chaperonin-buffered folding of polypeptides whose folding pathway has been perturbed (Van Dyk et al., 1989; ; Williams and Fares, 2010; ). However, since GroE is an active ATPase, its overproduction could be deleterious to the cell, owing to the depletion of cellular energy pools. Here, we have assessed the effect of GroE overproduction on the growth characteristics and thermal tolerance of E. coli and used proteomics and in silico flux balance analysis (FBA) to determine the likely impact of chaperonin overproduction on the metabolic advantage and consequent fitness of the organism.
Materials and Methods
Materials, Plasmids, Bacterial Strains and Growth Conditions
All chemicals were from Sigma, Inc. Bacterial growth media and media supplements were from HiMedia Laboratories, Inc., Mumbai, India. Phusion polymerase for colony PCR was purchased from New England Biolabs Inc., United States. GroE expression plasmids, pBAD-GSL and pTrc-GSL were generated by cloning GroE operon into NcoI and HindIII sites on plasmids pBAD24 () and pTrc99A (), respectively. The groE conditional mutant strain, E. coli LG6, was a kind gift from Arthur Horwich, Yale University, United States (). This strain produces GroE at levels similar to the wildtype at 30°C upon induction (Supplementary Figure 1). Oligonucleotide primers were purchased from Integrated DNA Technologies, Inc., Coralville, IA, United States.
Construction and Validation of Strains Producing High and Low GroE Levels
To enable control of GroE levels independently from the growth temperature, two strains that differentially express groE were generated from the E. coli strain LG6, in which the chromosomal groE promoter has been replaced with a Plac promoter (). A high level GroE expression strain, GL-Ht (for GroEL High pTrc), was obtained by transforming LG6 with pTrc-GSL and a lower level GroE expression strain, GL-Lt (for GroEL Low pTrc) was obtained by transforming with the control plasmid pTrc99A (). The scheme for the generation of these phenotypes is illustrated in Figure 1. To confirm the expression levels, these strains were cultured in the presence of 0.2% D-lactose to induce chromosome and plasmid borne groE operons, for 3 h at 30°C. The resulting cells were suspended in lysis buffer containing 50 mM HEPES:KOH pH 7.5 and 150 mM NaCl, 1 mM EDTA, and 1 mM PMSF, mixed with Lysing Matrix E and lysed by homogenization in FastPrep (M. P. Biomedicals, Irvine, CA, United States). Lysates were centrifuged at 13,000 rpm for 20 min to obtain soluble lysates. The soluble lysates were resolved on 12.5% SDS-PAGE and 12% Tricine gel followed by Coomassie Brilliant Blue staining to detect the levels of GroEL and GroES, respectively. In parallel, these lysates were probed with an anti-GroEL monoclonal antibody (1.10B) at 1:100 dilution and the blots were developed by BCIP/NBT-Purple Liquid Substrate System (Sigma Aldrich Inc., St. Louis, MO, United States). In addition to these strains, two strains that enable independent regulation of the chromosome and plasmid borne copies of groE operon were generated by transforming LG6 with pBAD-GSL and pBAD24 to result in high and low expression strains, GL-Hb and GL-Lb, respectively. These strains were cultured in the presence of 0.2% lactose plus 0.2% arabinose to obtain the high and low expression levels (Supplementary Figure 2A).
FIGURE 1
Temperature Sensitivity Assessment
The extent to which GroE overproduction enables temperature tolerance was assessed using a complementation assay (; ). Actively growing cultures of GL-Ht and GL-Lt were normalized for OD600, serially diluted, and spotted onto eight LB agar plates supplemented with 0.2% D-lactose. The plates were incubated at 17, 20, 22, 25, 30, 37, 40, 42, 45, 46, and 48°C. Wild type MG1655 harboring pTrc-GSL or empty vector (pTrc99A), respectively, were included as controls.
Competition and Estimation of Relative Fitness
GL-Hb and GL-Lb cells were subjected to competitive serial culturing as described previously (Zambrano et al., 1993; Vulic and Kolter, 2001; ). Briefly, equal number of cells from these two cultures were mixed and grown in fresh LB supplemented with 0.2% L-arabinose and 0.2% D-lactose. This mixed culture was grown to stationary phase at 30°C, recovered, labelled Passage-1 and used to generate the second passage (Figure 3A). Serial sub-culturing was repeated for a further 20 passages (∼700 generations). At each passage, a fraction of the cultures was serially diluted up to 10−7 dilution in LB broth and spread on LB agar plates supplemented with 0.2% D-lactose, which supports the growth of the cells derived from either strain. The resulting colonies at each passage, in the range of 23–28 colonies, were screened using colony PCR to identify whether colonies were derived from either GL-Hb or GL-Lb cells. Colony PCR with the PBADF (5′-CTGTTTCTCCATACCCGTT-3′) and PBADR (5′-CTCATCCGCCAAAACAG-3′) primers, which bind upstream and downstream of the MCS on the parental vector pBAD24, results in the amplification of 2.1 and 0.3 kb fragments from the pBAD-GSL and pBAD24 vectors, harbored by the GL-Hb and GL-Lb cells, respectively. Relative competitive index (CI), a measure of relative fitness, was calculated for each phenotype as the ratio of the proportion of a particular cell type at the final and initial generations (; ; van Opijnen and Camilli, 2013).
Proteomic Analysis
Equal number of cells from exponentially growing cultures (OD600 = ∼0.6) of GL-Hb or GL-Lb strains were harvested, suspended in lysis buffer (50 mM HEPES:KOH pH: 7.5 and 150 mM NaCl, 1 mM EDTA, and 1 mM PMSF), lysed by sonication, and the soluble protein fractions were recovered by centrifugation at 12,000 rpm for 20 min 200 µg protein from the soluble fractions of each lysate were resolved through 2D PAGE following the standard protocols. Briefly, the lysates were resolved on the first dimension through a 7 cm Immobilized pH Gradient (IPG) strip of 3–10 pH range, followed by 10% SDS-PAGE on the second dimension. The separated proteins were stained with Coomassie brilliant blue and intensities of the stained protein spots were compared between the two gels using densitometry. This experiment was repeated three times to identify the spots that exhibited consistent differential enrichment between the strains. Differentially enriched spots between the two lysates were picked and identified by tandem mass-spectrometry in an LTQ Orbitrap Mass Spectrometer (Thermo Fisher Scientific Inc., Waltham, MA, United States). The differentially enriched proteins were identified using MASCOT () search against UniProtKB/TrEMBL (UniProt, 2019) and RefSeq () databases. The spot identification was done in collaboration with the Centre for Cellular and Molecular Platforms, Bangalore, India.
Flux Balance Analysis of the GL-Hb and GL-Lb Strains
E. coli genome-scale metabolic network iJO1366 () was used for performing the FBA simulations. The iJO1366 model was first simulated using a standard energy source (equivalent of a glucose-supplemented minimal media) to obtain the steady state fluxes through each of the reactions (). The objective function of this FBA simulation was to maximize the biomass production, while using some “default constraints” (lower- and upper-bounds of fluxes through each reaction) derived from the literature (). Following this preliminary assessment of the E. coli cell’s metabolic potential, two independent FBA simulations were performed, each of which corresponded to the enzyme expression/enrichment profiles of the GL-Hb and GL-Lb strains. During each of these simulations the reaction flux values were appropriately constrained, based on the results from the preliminary assessment and the corresponding enzyme expression/enrichment profiles (Supplementary File 2). Incorporating enzyme expression profiles into FBA simulations was performed with our software tool “TransFlux,” developed in-house, and housed at http://www.nccs.res.in/TrasFlux/index.jsp. Details of the parameters and the principles applied in FBA are presented in the Supplementary Material methods section.
Results
Construction of GroE Overproducing Strains
To investigate the effect of chaperonin overproduction on E. coli, we constructed two chaperonin producing strains, GL-Ht and GL-Lt, which produce high and low levels of GroE (Figure 1). These strains were derived from strain E. coli LG6 (), in which the PgroE promoter is replaced by the Plac promoter, by transforming with pTrc-GSL, which overexpress groE operon upon induction with lactose, or its parental plasmid pTrc99A. SDS-PAGE confirmed significant overproduction of GroEL (Supplementary Figure 1A) in GL-Ht compared to GL-Lt. From Western blotting of the lysates, we estimate that GroEL levels are twenty-fold greater in GL-Ht than in GL-Lt (Supplementary Figure 1B). The expression levels of GroEL in GL-Lt were lower than the MG1655, where wildtype PgroE promoter drives the expression (Supplementary Figure 1) (). Further, GroES was significantly overproduced in GL-Ht compared to GL-Lt (Supplementary Figure 1C).
GroEL-GroES Overproducing Strains Showed Enhanced Temperature Tolerance
As GroE is involved in protection against thermal stress, we analyzed the impact of different GroE levels in GL-Ht and GL-Lt on growth at temperatures ranging from 17 to 48°C (Figure 2) (). E. coli MG1655 and MG1655 hosting pTrc-GSL were included for comparison. As expected, GL-Lt cells exhibited heat and cold sensitive phenotypes and consequently showed poor growth at many temperatures, consistent with previous observations that sufficient levels of GroE are required for growth over a wide temperature range (). Further, MG1655 showed much better temperature tolerance than GL-Lt, showing the importance of the heat-shock regulation of the PgroE promoter. The strains harboring pTrc-GSL tolerated higher temperatures, up to 48°C, than the vector-only MG1655, where groE expression is temperature regulated, suggesting that higher levels of GroE enable higher temperature tolerance.
FIGURE 2
GroEL-GroES Overproducing Strain Exhibited Fitness Advantage in Competition Culture
Since higher levels of chaperonins led to a growth advantage, we examined whether this translated to a fitness advantage even under low stress conditions, by competing two strains with different GroE levels. Since the two strains showed similar growth profiles on the plates (Figure 2) and in independent liquid cultures at 30°C (Supplementary Figure 3), we chose this temperature for the competition culture. To do these experiments, we needed to be able to control the plasmid borne and chromosomal copies of the groE operon independently. Therefore, we constructed two new strains with a PBAD based plasmid expression system, called GL-Hb (high expression) and GL-Lb (low expression) strains. Similar to GL-Ht, GL-Hb showed several folds higher GroE induction levels (Supplementary Figure 2A) and temperature resistance (Supplementary Figure 2B). The cultures of GL-Hb and GL-Lb were competed for 20 passages (∼700 generations) and their relative fitness(s) were estimated (Figure 3A) as described in Materials and Methods (; ; van Opijnen and Camilli, 2013). The high groE expressing GL-Hb outcompeted GL-Lb (Figure 3B), indicating that chaperonin level is an important fitness determinant.
FIGURE 3
Proteomic Analysis Revealed Preferential Enrichment of Metabolic Enzymes in GroEL-GroES Overproducing Strains
Overproduction of a chaperonin is likely to enrich the levels of folded proteins in the cells, while unfolded or misfolded proteins tend to remain insoluble and thereby targeted to either the inclusion bodies or marked for degradation (). Given this context, we investigated the proteomes of GL-Hb and GL-Lb cells, to identify what might account for the differences in fitness. Both strains were grown under identical conditions and their soluble proteome profiles (on 2D PAGE) were compared for relative abundance (Supplementary Figure 4; Table 1). Many of the identified proteins were known chaperonin clients belonging to either classes I and II (), class IV () or the clients identified exclusively in , which here we have denoted as class V. However, several proteins that were identified as being differentially expressed were not known clients (Table 1), suggesting that either chaperonin overexpression can indirectly affect the folding of these non-client proteins or that the chaperonin client base is larger than currently understood. We noted that none of the obligate class III GroEL clients () were relatively enriched in either strain, showing that there is sufficient chaperonin activity for folding these clients in the GL-Lb strain. Notably, the outer membrane proteins, OmpC and OmpF, which are involved in metabolite import and are known GroE clients (), were enriched in the soluble proteome of GL-Lb. The higher level of OmpC and OmpF in the soluble fraction of GL-Lb suggested a lower proportion of these proteins might be reaching the outer membrane in these strains. We therefore quantified the relative levels of OmpC and OmpF in membrane fractions of both pairs of strains, and confirmed that the levels were lower in both GL-Lb and GL-Lt (Supplementary Figure 5). Further, a higher instability index (obtained from Expasy ProtParam tool), which is a reverse measure of protein stability (; ) was observed for the proteins enriched in GL-Hb strain, suggesting that their enrichment in the chaperonin overexpressing condition may be linked to lower stability and hence a greater chaperonin requirement. The enrichment of TF in GL-Lb (Supplementary Figure 1A; Table 1), is consistent with previously reported interactions between TF and GroE (; ) and suggests TF may be able to partially compensate for low levels of chaperonin function in GL-Lb. Furthermore, enrichment of several metabolic enzymes in the GL-Hb strain, suggested a higher rate of metabolism in this strain. To evaluate this hypothesis, we collated publicly available E. coli proteomic data from the paxdb database (Wang et al., 2012), screened for proteins that were co-enriched with GroE across different experiments and identified 404 proteins that showed significant correlation, in expression levels, with GroE (Pearson correlation co-efficient ≥ 0.7, p < 0.05). Interestingly, a GO enrichment analysis of this set of proteins revealed that majority of these proteins were involved in metabolism and energy production, including multiple GO terms related to carbohydrate metabolism (Table 2).
TABLE 1
| Strain | SwissProt entry | Protein description | MW [kDa] | pI | Unique peptides | Coverage | GroEL client classa | Oligomeric state | COG | SCOP fold class | Protein instabilitya | In vivo location | Gene | Ea | mRNA t1/2 (min) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| GL-Hb | ENO_ECOLI (P0A6P9) | Enolase (EC:4.2.1.11) (2-phosphoglycerate dehydratase) (2-phospho-D-glycerate hydro-lyase) | 45.5 | 5.32 | 19 | 51.68 | One | Homodimer | G | c.1.1.1; d.54.1.1 | 25.64 (Stable) | Cytoplasm, cyto-skeleton, secreted, cell surface | eno (b2779) | 1 | 4.7 |
| 6PGD_ECOLI (P00350) | 6-Phosphogluconate dehydrogenase, decarboxylating (EC:1.1.1.44) | 51.5 | 5.04 | 16 | 57.48 | One | Homodimer | G | a.100.1; c.2.1.6 | 35.98 (Stable) | Cytoplasm | gnd (b2029) | 0 | 10.6 | |
| DLDH_ECOLI (P0A9P0) | Dihydrolipoyl dehydrogenase (EC:1.8.1.4), Dihydrolipoamide dehydrogenase, E3 component of pyruvate and 2-oxoglutarate dehydrogenases complexes | 50.6 | 5.79 | 18 | 44.51 | Four | Homodimer | C | d.87.1.1; c.3.1.5 | 18.84 (Stable) | Cytoplasm, cell inner membrane, peripheral membrane | lpdA (b0116) | 1 | 5.8 | |
| MDH_ECOLI (P61889) | Malate dehydrogenase (EC:1.1.1.37) | 32.3 | 5.28 | 19 | 87.5 | Two | Homodimer | C | d.162.1.1; c.2.1.5 | 30.58 (Stable) | Cytoplasm | mdh (b3236) | 1 | 10.5 | |
| TYPH_ECOLI (P07650) | Thymidine phosphorylase (EC:2.4.2.4) | 51.4 | 5.2 | 19 | 51.36 | Four | Homodimer | F | c.27.1.1; d.41.3.1; a.46.2.1 | 20.63 (Stable) | Cytoplasm | deoA (b4382) | 0 | 15.8 | |
| IDH_ECOLI (P08200) | Isocitrate dehydrogenase [NADP] (EC:1.1.1.42) (Oxalosuccinate decarboxylase) | 45.7 | 5.15 | 19 | 65.28 | — | Homodimer | C | c.77.1.1 | 34.73 (Stable) | Cytoplasm | icd (b1136) | 0 | 5.8 | |
| ACEA_ECOLI (P0A9G6) | Isocitrate lyase (EC:4.1.3.1) | 47.4 | 5.16 | 11 | 40.32 | — | Homo-tetramer | C | c.1.12.7 | 36.53 (Stable) | Cytoplasm | aceA (b4015) | 0 | 11.5 | |
| DPO3B_ECOLI (P0A988) | DNA polymerase III beta subunit protein (EC:2.7.7.7) | 40.5 | 5.45 | 12 | 40.71 | — | Hetero-Oligomer | L | d.131.1.1 | 42.49 (Unstable) | Cytoplasm | dnaN (b3701) | 1 | 2.4 | |
| TALB_ECOLI (P0A870) | Transaldolase B (EC:2.2.1.2) | 35 | 5.11 | 21 | 77.29 | Five | Homodimer | G | c.1.10.1 | 31.71 (Stable) | Cytoplasm | talB (b0008) | 0 | 3.4 | |
| POTD_ECOLI (P0AFK9) | Spermidine/putrescine-binding periplasmic protein | 38.8 | 4.86 | 16 | 50.86 | — | Monomer | E | c.94.1.1 | 21.34 (Stable) | Periplasm | potD (b1123) | 1 | — | |
| RIHA_ECOLI (P41409) | Pyrimidine-specific ribonucleoside hydrolase, RihA (EC:3.2.2.-), Cytidine/uridine-specific hydrolase, ribonucleoside hydrolase 1 | 33.8 | 4.84 | 13 | 77.81 | — | Tetramer | F | C.70.1.0 | 30.57 (Stable) | Cytoplasm | rihA (b0651) | 0 | 4.4 | |
| CH60_ECOLI (P0A6F5) | Chaperonin 60, GroEL | 57 | 4.85 | 11 | 28.89 | Five | Homo-tetradecamer | O | a.129.1.1; d.56.1.1; c.8.5.1 | 29.30 (Stable) | Cytoplasm | groL (b4143) | 1 | 3.5 | |
| BGAL_ECOLI (P00722) | Beta-galactosidase (EC:3.2.1.23) | 116.4 | 5.28 | 51 | 67.68 | — | Homo-tetramer | G | b.30.5.1; c.1.8.3; b.18.1.5; b.1.4.1 | 43.27 (Unstable) | Cytoplasm | lacZ (b0344) | 0 | 10.4 | |
| GL-Lb | TIG_ECOLI (P0A850) | Trigger factor (EC:5.2.1.8) (TF) | 48.2 | 4.83 | 28 | 65.05 | One | Homodimer and monomer | O | i.1.1.2; d.241.2.1; d.26.1.1; a.223.1.1 | 37.21 (Stable) | Cytoplasm | tig (b0436) | 0 | 2.3 |
| RPOA_ECOLI (P0A7Z4) | DNA-directed RNA polymerase subunit alpha (EC:2.7.7.6) (RNAP subunit alpha), RNA polymerase subunit alpha, Transcriptase subunit alpha | 36.5 | 4.97 | 15 | 56.53 | One | Homodimer | K | d.181.1.1; i.8.1.1; a.60.3.1; d.74.3.1 | 41.59 (Unstable) | Cytoplasm | rpoA (b3295) | 1 | 4 | |
| PGK_ECOLI (P0A799) | Phosphoglycerate kinase (EC:2.7.2.3) | 41 | 5.08 | 22 | 73.9 | One | Monomer | G | c.86.1.1; c.1.1.1 | 26.37 (Stable) | Cytoplasm | pgk (b2926) | 1 | 2.5 | |
| OMPC_ECOLI (P06996) | Outer membrane protein C, outer membrane protein 1B, porin, OmpC | 40.3 | 4.48 | 20 | 77.38 | One | Homotrimer | M | f.4.3.1 | 12.86 (Stable) | Outer membrane | ompC (b2215) | 0 | 9.7 | |
| OMPF_ECOLI (P02931) | Outer membrane protein F, outer membrane protein 1A, outer membrane protein B, porin, OmpF | 39.3 | 4.64 | 24 | 82.6 | Two | Homotrimer | M | f.4.3.1 | 13.81 (Stable) | Outer membrane | ompF (b0929) | 0 | 8.5 | |
| ALF_ECOLI (P0AB71) | Fructose-bisphosphate aldolase class II (EC 4.1.2.13) (FBP aldolase), Fructose-1,6-bisphosphate aldolase | 39.1 | 5.52 | 11 | 50.42 | Two | Homodimer | G | c.1.10.2 | 34.82 (Stable) | Cytoplasm | fbaA (b2925) | 1 | 7.2 | |
| GLF_ECOLI (P37747) | UDP-galactopyranose mutase (EC:5.4.99.9), UDP-GALP mutase, Uridine 5-diphosphate galactopyranose mutase | 43 | 6.61 | 27 | 79.02 | Five | Homodimer | M | d.16.1.7; c.4.1.3 | 32.48 (Stable) | Cytoplasm | glf (b2036) | 0 | — | |
| SUCC_ECOLI (P0A836) | Succinyl-CoA ligase [ADP-forming] subunit beta (EC:6.2.1.5), succinyl-CoA synthetase subunit beta | 41.3 | 5.37 | 22 | 76.8 | Five | Hetero-tetramer | C | c.23.4.1; d.142.1.4 | 30.24 (Stable) | Cytoplasm | sucC (b0728) | 0 | 6.7 | |
| MALE_ECOLI (P0AEX9) | Maltose-binding periplasmic protein, MBP, MMBP, Maltodextrin-binding protein | 43.3 | 5.22 | 13 | 51.77 | Five | Hetero-pentamer | G | c.94.1.1 | 18.23 (Stable) | Periplasmic | malE (b4034) | 0 | — | |
| LACI_ECOLI (P03023) | Lactose operon repressor (LacI) | 38.5 | 6.39 | 24 | 80.28 | Five | Homo-tetramer | K | c.93.1.1; a.35.1.5 | 37.37 (Stable) | Cytoplasm | lacI (b0345) | 0 | 5.7 | |
| MANA_ECOLI (P00946) | Mannose-6-phosphate isomerase (EC:5.3.1.8), Phosphohexomutase, Phosphomannose isomerase (PMI) | 42.8 | 5.29 | 16 | 62.92 | — | Monomer | G | b.82.1.3 | 40.37 (Unstable) | Cytoplasm | manA (b1613) | 0 | 3.6 | |
| TREC_ECOLI (P28904) | Trehalose-6-phosphate hydrolase (EC:3.2.1.93), Alphaalpha-phosphotrehalase | 63.8 | 5.51 | 31 | 62.61 | — | — | G | c.87.1.6 | 33.96 (Unstable) | Cytoplasm | treC (b4239) | 0 | 4.3 | |
| AAT_ECOLI (P00509) | Aspartate aminotransferase (EC:2.6.1.1), AspAT, Transaminase A | 43.5 | 5.54 | 23 | 61.62 | Five | Homodimer | E | c.67.1.1 | 29.50 (Stable) | Cytoplasm | aspC (b0928) | 0 | 4.3 |
Properties of the differentially enriched proteins in GL-Hb and GL-Lb strains.
GroEL substrate classes 1–3 are from , class 4 is from , and the proteins exclusive to study were denoted as class 5. Protein stability is depicted as instability index obtained from Expasy Protparam (; ). Column E lists the essential (1) and non-essential (0) genes.
TABLE 2
| Gene ontology terms | Protein count | Fold enrichment | p-value | Bonferroni correction |
|---|---|---|---|---|
| GO:0006091: Generation of precursor metabolites and energy | 76 | 3.810 | 7.95e−24 | 1.28e−21 |
| GO:0044249: Cellular biosynthetic process | 193 | 1.541 | 5.72e−12 | 9.21e−10 |
| GO:0042180: Cellular ketone metabolic process | 84 | 2.130 | 4.30e−11 | 6.92e−09 |
| GO:0006082: Organic acid metabolic process | 82 | 2.118 | 1.04e−10 | 1.68e−08 |
| GO:0009308: Amine metabolic process | 73 | 2.003 | 1.55e−08 | 2.50e−06 |
| GO:0016052: Carbohydrate catabolic process | 42 | 2.525 | 9.87e−08 | 1.59e−05 |
| GO:0022900: Electron transport chain | 26 | 2.992 | 2.28e−06 | 3.68e−04 |
| GO:0006519: Cellular amino acid and derivative metabolic process | 58 | 1.908 | 3.06e−06 | 4.92e−04 |
| GO:0046483: Heterocycle metabolic process | 45 | 2.052 | 8.50e−06 | 1.37e−03 |
| GO:0006793: Phosphorus metabolic process | 29 | 2.405 | 3.27e−05 | 5.26e−03 |
| GO:0019538: Protein metabolic process | 60 | 1.672 | 9.90e−05 | 1.58e−02 |
| GO:0009059: Macromolecule biosynthetic process | 114 | 1.355 | 3.68e−04 | 5.76e−02 |
| GO:0006766: Vitamin metabolic process | 20 | 2.397 | 7.47e−04 | 1.13e−01 |
| GO:0006790: Sulphur metabolic process | 18 | 2.540 | 7.89e−04 | 1.19e−01 |
| GO:0016051: Carbohydrate biosynthetic process | 31 | 1.879 | 1.20e−03 | 1.75e−01 |
| GO:0044248: Cellular catabolic process | 29 | 1.864 | 2.06e−03 | 2.83e−01 |
| GO:0006461: Protein complex assembly | 11 | 3.214 | 2.28e−03 | 3.07e−01 |
| GO:0065003: Macromolecular complex assembly | 11 | 3.189 | 2.41e−03 | 3.22e−01 |
| GO:0005975: Carbohydrate metabolic process | 77 | 1.385 | 2.63e−03 | 3.46e−01 |
| GO:0033014: Tetrapyrrole biosynthetic process | 10 | 3.197 | 4.12e−03 | 4.86e−01 |
| GO:0044255: Cellular lipid metabolic process | 27 | 1.787 | 5.08e−03 | 5.60e−01 |
| GO:0009991: Response to extracellular stimulus | 10 | 3.016 | 6.02e−03 | 6.22e−01 |
| GO:0051186: Cofactor metabolic process | 30 | 1.640 | 9.63e−03 | 7.89e−01 |
| GO:0009057: Macromolecule catabolic process | 14 | 2.247 | 9.88e−03 | 7.98e−01 |
Enriched Gene Ontology terms (level 3 - biological process terms), associated with the 404 proteins that were co-enriched/expressed with GroE across different experiments.
Flux Balance Analysis of Oxidative Phosphorylation in High- and Low-GroEL Strains.
Considering the preferential enrichment of metabolic enzymes upon GroE overproduction, we adopted an FBA approach (; ) to assess how the differential enrichment of metabolic enzymes in the GL-Lb and GL-Hb strains would translate into altered metabolic states and cellular fitness. The FBA simulation analyses were carried out using “TransFlux” (available at: http://www.nccs.res.in/TransFlux/index.jsp), an in-house tool with a module to incorporate gene expression/proteomic profiles in the FBA framework. The proteomic profiles (Table 1) and observations from E. coli gene expression microarray studies, derived from the Many Microbe Microarrays database (M3D, www.m3d.mssm.edu) () were utilized to constrain fluxes though respective reactions, while performing two independent FBA simulations, each of which corresponded to the expression/enrichment profiles of the enzymes enriched in GL-Lb and GL-Hb strains. As expected, higher flux was observed through several pathways of carbon metabolism including glycolysis, gluconeogenesis, citric acid cycle (TCA cycle) and its anaplerotic reactions, and alternate carbon metabolism, in the simulated GL-Hb strain (Table 3). These pathways appear to be supported by enhanced import of glucose and glycerol (Supplementary File 2). Pathways corresponding to several glucogenic amino acids metabolism and energy generating oxidative phosphorylation were enriched in this strain. However, the pathways leading to the toxic methylglyoxal synthesis were also enriched in the GL-Hb strain (Table 3). We also noted that pathways leading to the metabolism of membrane lipids, pyruvic acid, pentose sugars, ubiquinone and salvage of nucleotides are enriched in the GL-Lb strain. Overall, FBA simulations indicated that the metabolic enzymes that were enriched in GL-Hb may lead to higher metabolic flux in this strain (Table 3; Supplementary Material).
TABLE 3
| Strain | Metabolic pathway | Metabolic flux through the pathway (mM/gm-DW/hr)a | |||
|---|---|---|---|---|---|
| Flux in GL-Hb | Flux in GL-Lb | Flux difference | Flux ratio | ||
| GL-Hb | Glycolysis/Gluconeogenesis | 289.2 | 59.6 | 229.6 | 2.3 |
| Citric acid cycle | 183.9 | 100.1 | 83.7 | 0.9 | |
| Oxidative phosphorylation | 103.3 | 48.9 | 54.5 | 1.1 | |
| Threonine and lysine metabolism | 41.2 | 0.6 | 40.6 | 6.0 | |
| Anaplerotic reactions | 33.8 | 1.4 | 32.5 | 4.6 | |
| Inorganic ion transport and metabolism | 62.7 | 31.6 | 31.1 | 1.0 | |
| Methylglyoxal metabolism | 27.2 | 0.0 | 27.2 | NA | |
| Transport, inner membrane | 116.1 | 91.9 | 24.3 | 0.3 | |
| Glutamate metabolism | 25.0 | 3.2 | 21.8 | 3.0 | |
| Alanine and aspartate metabolism | 126.3 | 111.6 | 14.8 | 0.2 | |
| Transport, outer membrane porin | 28.1 | 17.0 | 11.1 | 0.7 | |
| Alternate carbon metabolism | 46.0 | 36.0 | 10.0 | 0.4 | |
| Glycine and serine metabolism | 9.8 | 0.8 | 9.0 | 3.6 | |
| GL-Lb | Membrane lipid metabolism | 0.4 | 0.4 | 0.0 | 0.0 |
| Cofactor and prosthetic group biosynthesis | 0.0 | 0.0 | 0.0 | −4.2 | |
| Nucleotide salvage pathway | 19.7 | 21.4 | −1.7 | −0.1 | |
| Pyruvate metabolism | 452.6 | 478.4 | −25.8 | −0.1 | |
| Unassigned | 0.2 | 30.1 | −29.9 | −7.0 | |
| Pentose phosphate pathway | 218.9 | 266.8 | −47.8 | −0.3 | |
Cumulative metabolic flux through major pathways in simulated GL-Lb and GL-Hb strains as obtained through Flux Balance Analysis. Log two fold-change of fluxes of GL-Hb and GL-Lb are indicated in the Flux Ratio column.
mM/gm-DW/hr, Millimolar Metabolite per Gram Dry Weight of the cell mass per hour.
Discussion
Over- or under-production of chaperonins in several organisms has been demonstrated to perturb rates of proteolysis (), influence growth rates, and alter the expression levels of compensatory chaperones like DnaK (). Here we present a simple model system to study the effects of GroE overproduction (Figure 1). We demonstrate that the overexpression of GroE chaperonin results in enhanced thermal tolerance (Figure 2) and competitive advantage (Figure 3). GroEL is known to be required for growth at low () and high () temperatures. Consistent with this, the GL-Lb and GL-Lt strains exhibited both cold and heat sensitive phenotypes (Figure 2). Proteomic studies (Table 1) followed by FBA (Tables 2,3) suggest that the acquired fitness advantage could be attributed to an enriched set of metabolic enzymes. Chaperonin depletion was observed to induce the enrichment of the compensatory chaperone, TF (Supplementary Figure 1A; Table 1), which may act as a holdase for the GroE client proteins (). Interestingly, while GroE is more abundant than TF in E. coli (Zou et al., 2014), TF is observed to be abundant in mycoplasma which lack the groE operon (; Weiner et al., 2003; ; ), suggesting that higher levels of TF might be needed in such bacteria to compensate for the chaperonin deficiency. The TF - GroEL interplay, owing to their overlapping functions and client-base (; ), has been demonstrated both in vitro () and in vivo in E. coli (). Therefore, it seems likely that TF enrichment in GL-Lb is compensating for GroE depletion and that TF may be acting on some clients as a holdase (). Further, the enrichment of the outer-membrane proteins OmpC and OmpF in the soluble proteome of GL-Lb suggests that these known GroEL client proteins failed to reach their normal final cellular destination (the outer membrane) and may have remained soluble, possibly in a TF-bound state. The reduced levels of these porins in the membranes of GL-Lb and GL-Lt strains (Supplementary Figure 5) might be responsible, in part, for the lower metabolite transport and metabolic flux in this strain (Table 3). TF was not upregulated in the wildtype strain (MG1655), despite lower GroE levels (Supplementary Figure 1), as GroE levels in this strain respond directly to levels of unfolded proteins. Furthermore, a different mode of GroE depletion resulted in the enrichment of DnaK (), which exhibits significant functional overlap with TF (Teter et al., 1999; ; ). The higher fitness of the GroES and GroEL over-producing strains under the conditions of our experiments is likely to be associated with fitness costs under other conditions (Figures 2, 3), otherwise it would be expected that higher expression would have evolved.
We demonstrate a direct relation between chaperonin abundance and competitive fitness. However, the evolution has not selected for intracellular chaperonin levels as high as the ones used in our experiments. The predictions from FBA simulations provide some clues that may explain why this has not occurred. Although enhanced glycolysis, TCA cycle and oxidative phosphorylation in the GL-Hb cells increase cellular energy currency, FBA simulations for the GL-Hb strain predicted an enhanced production of a toxic side product, methylglyoxal (Table 3), a very toxic three-carbon aldehyde that can inhibit E. coli growth at millimolar concentrations (; Weber et al., 2005). Therefore, evolution might have selected a balance in metabolic states between energy production and methylglyoxal toxicity, which would have, in turn, selected for an optimal level of chaperonin production. The fact that chaperonins are active ATPases provides another possible answer to this question. Overabundance of chaperonins might be linked to ATP depletion and consequent reduced growth (). Thus, very high levels of chaperonin expression may have been selected against during the course of evolution. These explanations are not exhaustive, and the final level of chaperonin expression selected for is likely to result from a balance of optimizing fitness, due to multiple different factors.
Our analysis showed that GroE over-production results in several pleiotropic consequences that can enhance cellular fitness under the tested conditions. These observations need to be probed further to enhance our understanding of the precise role of the chaperone-client interactions in influencing fitness and, ultimately, evolution. A similar system could be advantageous in studying the effect of chaperonin overproduction in different microbes, especially the pathogenic bacteria with multiple chaperonins (; ).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
CK did principal experiments. KC did proteomics experiments, AD, TB, VM, and SSM did the FBA, SCM and PL monitored the initial and later parts of the experiments.
Funding
The project is currently funded by the Biotechnology and Biological Sciences Research Council Responsive mode grant (BB/S017526/1). CK was a Newton International Fellow (NF161469) sponsored by the Royal Society, British Academy and Academy of Medical Sciences, United Kingdom. Initial part of this work was supported by grants from the Department of Biotechnology, India (BT/PR3260/BRB/10/967/2011).
Acknowledgments
We would like to acknowledge the support from Biotechnology and Biological Sciences Research Council Responsive mode grant (BB/S017526/1). CK was a Newton International Fellow (NF161469) sponsored by the Royal Society, British Academy and Academy of Medical Sciences, United Kingdom. We would like to thank, Arthur Horwich for the GroE depletion strain, E. coli LG6 and cCAMP, Bangalore for assistance in proteomic studies. We thank Abhijit Sardesai and Gaurang Mahajan for helpful discussions, Melanie Swannell, Amanda Rossiter, Ian Henderson, Anna Schager, Chistopher Icke, Shahida Rafique, Nitin Bayal, Sapna Sugandhi and Surbhi Dhingra for support in the initial studies, and Ishita Verma for helping in designing the TransFlux website. SSM, AD, and TB are employees of TCS Research (Tata Consultancy Services Ltd., Pune, India), and would like to acknowledge TCS for its support. Initial part of this work was supported by grants from the Department of Biotechnology, India (BT/PR3260/BRB/10/967/2011).
Conflict of interest
Authors AD, TB, and SSM are employed by the company TCS Research division in the Tata Consultancy Services Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.669996/full#supplementary-material
References
1
AmannE.OchsB.AbelK.-J. (1988). Tightly Regulated Tac Promoter Vectors Useful for the Expression of Unfused and Fused Proteins in Escherichia coli. Gene69 (2), 301–315. 10.1016/0378-1119(88)90440-4
2
AvellanedaM. J.KoersE. J.NaqviM. M.TansS. J. (2017). The Chaperone Toolbox at the Single-Molecule Level: From Clamping to Confining. Protein Sci.26 (7), 1291–1302. 10.1002/pro.3161
3
BalchinD.Hayer-HartlM.HartlF. U. (2016). In Vivo aspects of Protein Folding and Quality Control. Science353 (6294), aac4354. 10.1126/science.aac4354
4
BangH.PechtA.RaddatzG.SciorT.SolbachW.BruneK.et al (2000). Prolyl Isomerases in a Minimal Cell. Eur. J. Biochem.267 (11), 3270–3280. 10.1046/j.1432-1327.2000.01355.x
5
BhandariV.HouryW. A. (2015). Substrate Interaction Networks of the Escherichia coli Chaperones: Trigger Factor, DnaK and GroEL. Adv. Exp. Med. Biol.883, 271–294. 10.1007/978-3-319-23603-2_15
6
BlaisE. M.ChavaliA. K.PapinJ. A. (2013). Linking Genome-Scale Metabolic Modeling and Genome Annotation. Methods Mol. Biol.985, 61–83. 10.1007/978-1-62703-299-5_4
7
BukauB.HorwichA. L. (1998). The Hsp70 and Hsp60 Chaperone Machines. Cell92 (3), 351–366. 10.1016/s0092-8674(00)80928-9
8
CalloniG.ChenT.SchermannS. M.ChangH.-c.GenevauxP.AgostiniF.et al (2012). DnaK Functions as a central Hub in the E. coli Chaperone Network. Cel Rep.1 (3), 251–264. 10.1016/j.celrep.2011.12.007
9
ChapmanE.FarrG. W.UsaiteR.FurtakK.FentonW. A.ChaudhuriT. K.et al (2006). Global Aggregation of Newly Translated Proteins in an Escherichia coli Strain Deficient of the Chaperonin GroEL. Proc. Natl. Acad. Sci.103 (43), 15800–15805. 10.1073/pnas.0607534103
10
ChilukotiN.KumarC. M. S.MandeS. C. (2016). GroEL2 of Mycobacterium tuberculosis Reveals the Importance of Structural Pliability in Chaperonin Function. J. Bacteriol.198 (3), 486–497. 10.1128/JB.00844-15
11
DeuerlingE.PatzeltH.VorderwülbeckeS.RauchT.KramerG.SchaffitzelE.et al (2003). Trigger Factor and DnaK Possess Overlapping Substrate Pools and Binding Specificities. Mol. Microbiol.47 (5), 1317–1328. 10.1046/j.1365-2958.2003.03370.x
12
DurãoP.AignerH.NagyP.Mueller-CajarO.HartlF. U.Hayer-HartlM. (2015). Opposing Effects of Folding and Assembly Chaperones on Evolvability of Rubisco. Nat. Chem. Biol.11 (2), 148–155. 10.1038/nchembio.1715
13
FaithJ. J.DriscollM. E.FusaroV. A.CosgroveE. J.HayeteB.JuhnF. S.et al (2008). Many Microbe Microarrays Database: Uniformly Normalized Affymetrix Compendia with Structured Experimental Metadata. Nucleic Acids Res.36, D866–D870. 10.1093/nar/gkm815
14
FaresM. A.Ruiz-GonzálezM. X.MoyaA.ElenaS. F.BarrioE. (2002). GroEL Buffers against Deleterious Mutations. Nature417 (6887), 398. 10.1038/417398a
15
FayetO.ZiegelhofferT.GeorgopoulosC. (1989). The groES and groEL Heat Shock Gene Products of Escherichia coli Are Essential for Bacterial Growth at All Temperatures. J. Bacteriol.171 (3), 1379–1385. 10.1128/jb.171.3.1379-1385.1989
16
FerrerM.ChernikovaT. N.YakimovM. M.GolyshinP. N.TimmisK. N. (2003). Chaperonins Govern Growth of Escherichia coli at Low Temperatures. Nat. Biotechnol.21 (11), 1266–1267. 10.1038/nbt1103-1266
17
FujiwaraK.IshihamaY.NakahigashiK.SogaT.TaguchiH. (2010). A Systematic Survey of In Vivo Obligate Chaperonin-dependent Substrates. EMBO J.29 (9), 1552–1564. 10.1038/emboj.2010.52
18
GasteigerE.HooglandC.GattikerA.DuvaudS. e.WilkinsM. R.AppelR. D.et al (2005). “Protein Identification and Analysis Tools on the ExPASy Server,” in The Proteomics Protocols Handbook. Editor WalkerJ. M.1 ed (Totowa, NJ: Humana Press), 571–607. 10.1385/1-59259-890-0:571
19
GenevauxP.KeppelF.SchwagerF.Langendijk‐GenevauxP. S.HartlF. U.GeorgopoulosC. (2004). In Vivo analysis of the Overlapping Functions of DnaK and Trigger Factor. EMBO Rep.5 (2), 195–200. 10.1038/sj.embor.7400067
20
GuisbertE.HermanC.LuC. Z.GrossC. A. (2004). A Chaperone Network Controls the Heat Shock Response in E. coli. Genes Dev.18 (22), 2812–2821. 10.1101/gad.1219204
21
GuruprasadK.ReddyB. V. B.PanditM. W. (1990). Correlation between Stability of a Protein and its Dipeptide Composition: a Novel Approach for Predicting In Vivo Stability of a Protein from its Primary Sequence. Protein Eng. Des. Sel4 (2), 155–161. 10.1093/protein/4.2.155
22
GuzmanL. M.BelinD.CarsonM. J.BeckwithJ. (1995). Tight Regulation, Modulation, and High-Level Expression by Vectors Containing the Arabinose PBAD Promoter. J. Bacteriol.177 (14), 4121–4130. 10.1128/jb.177.14.4121-4130.1995
23
HartlF. U.Hayer-HartlM. (2002). Molecular Chaperones in the Cytosol: from Nascent Chain to Folded Protein. Science295 (5561), 1852–1858. 10.1126/science.1068408
24
HirosawaM.HoshidaM.IshikawaM.ToyaT. (1993). MASCOT: Multiple Alignment System for Protein Sequences Based on Three-Way Dynamic Programming. Bioinformatics9 (2), 161–167. 10.1093/bioinformatics/9.2.161
25
HorwichA. L.LowK. B.FentonW. A.HirshfieldI. N.FurtakK. (1993). Folding In Vivo of Bacterial Cytoplasmic Proteins: Role of GroEL. Cell74 (5), 909–917. 10.1016/0092-8674(93)90470-b
26
HouryW. A.FrishmanD.EckerskornC.LottspeichF.HartlF. U. (1999). Identification of In Vivo Substrates of the Chaperonin GroEL. Nature402 (6758), 147–154. 10.1038/45977
27
IshimotoT.FujiwaraK.NiwaT.TaguchiH. (2014). Conversion of a Chaperonin GroEL-independent Protein into an Obligate Substrate. J. Biol. Chem.289 (46), 32073–32080. 10.1074/jbc.M114.610444
28
KandrorO.ShermanM.MoerschellR.GoldbergA. L. (1997). Trigger Factor Associates with GroEL In Vivo and Promotes its Binding to Certain Polypeptides. J. Biol. Chem.272 (3), 1730–1734. 10.1074/jbc.272.3.1730
29
KandrorO.ShermanM.RhodeM.GoldbergA. L. (1995). Trigger Factor Is Involved in GroEL-dependent Protein Degradation in Escherichia coli and Promotes Binding of GroEL to Unfolded Proteins. EMBO J.14 (23), 6021–6027. 10.1002/j.1460-2075.1995.tb00290.x
30
KayserA.WeberJ.HechtV.RinasU. (2005). Metabolic Flux Analysis of Escherichia coli in Glucose-Limited Continuous Culture. I. Growth-rate-dependent Metabolic Efficiency at Steady State. Microbiology151 (Pt 3), 693–706. 10.1099/mic.0.27481-0
31
KernerM. J.NaylorD. J.IshihamaY.MaierT.ChangH.-C.StinesA. P.et al (2005). Proteome-wide Analysis of Chaperonin-dependent Protein Folding in Escherichia coli. Cell122 (2), 209–220. 10.1016/j.cell.2005.05.028
32
KimY. E.HippM. S.BracherA.Hayer-HartlM.Ulrich HartlF. (2013). Molecular Chaperone Functions in Protein Folding and Proteostasis. Annu. Rev. Biochem.82, 323–355. 10.1146/annurev-biochem-060208-092442
33
KumarC. M. S.KhareG.SrikanthC. V.TyagiA. K.SardesaiA. A.MandeS. C. (2009). Facilitated Oligomerization of Mycobacterial GroEL: Evidence for Phosphorylation-Mediated Oligomerization. J. Bacteriol.191 (21), 6525–6538. 10.1128/JB.00652-09
34
KumarC. M. S.MandeS. C.MahajanG. (2015). Multiple Chaperonins in Bacteria-Novel Functions and Non-canonical Behaviors. Cell Stress and Chaperones20 (4), 555–574. 10.1007/s12192-015-0598-8
35
KumarC. M. S. (2017). “Prokaryotic Multiple Chaperonins: The Mediators of Functional and Evolutionary Diversity,” in Prokaryotic Chaperonins: Multiple Copies and Multitude Functions. Editors KumarC. M. S.MandeS. C. (Singapore: Springer Singapore), 39–51. 10.1007/978-981-10-4651-3_3
36
KusukawaN.YuraT. (1988). Heat Shock Protein GroE of Escherichia coli: Key Protective Roles against thermal Stress. Genes Dev.2 (7), 874–882. 10.1101/gad.2.7.874
37
LemosJ. A.LuzardoY.BurneR. A. (2007). Physiologic Effects of Forced Down-Regulation of dnaK and groEL Expression in Streptococcus Mutans. J. Bacteriol.189 (5), 1582–1588. 10.1128/JB.01655-06
38
LiG.-W.BurkhardtD.GrossC.WeissmanJ. S. (2014). Quantifying Absolute Protein Synthesis Rates Reveals Principles Underlying Allocation of Cellular Resources. Cell157 (3), 624–635. 10.1016/j.cell.2014.02.033
39
LlorcaO.GalánA.CarrascosaJ. L.MugaA.ValpuestaJ. M. (1998). GroEL under Heat-Shock. J. Biol. Chem.273 (49), 32587–32594. 10.1074/jbc.273.49.32587
40
LundP. A. (2001). Molecular Chaperones in the Cell. Oxford University Press.
41
LundP. A. (2009). Multiple Chaperonins in Bacteria - Why So many?FEMS Microbiol. Rev.33 (4), 785–800. 10.1111/j.1574-6976.2009.00178.x
42
MachoA. P.GuidotA.BarberisP.BeuzónC. R.GeninS. (2010). A Competitive index Assay Identifies Several Ralstonia Solanacearum Type III Effector Mutant Strains with Reduced Fitness in Host Plants. Mpmi23 (9), 1197–1205. 10.1094/MPMI-23-9-1197
43
MartinJ.HorwichA.HartlF. (1992). Prevention of Protein Denaturation under Heat Stress by the Chaperonin Hsp60. Science258 (5084), 995–998. 10.1126/science.1359644
44
Martínez-AlonsoM.García-FruitósE.Ferrer-MirallesN.RinasU.VillaverdeA. (2010). Side Effects of Chaperone Gene Co-expression in Recombinant Protein Production. Microb. Cel Fact9, 64. 10.1186/1475-2859-9-64
45
MonkI. R.CaseyP. G.CroninM.GahanC. G.HillC. (2008). Development of Multiple Strain Competitive index Assays for Listeria Monocytogenes Using pIMC; a New Site-specific Integrative Vector. BMC Microbiol.8, 96. 10.1186/1471-2180-8-96
46
MusatovovaO.DhandayuthapaniS.BasemanJ. B. (2006). Transcriptional Heat Shock Response in the Smallest Known Self-Replicating Cell, Mycoplasma Genitalium. J. Bacteriol.188 (8), 2845–2855. 10.1128/JB.188.8.2845-2855.2006
47
O'LearyN. A.WrightM. W.BristerJ. R.CiufoS.HaddadD.McVeighR.et al (2016). Reference Sequence (RefSeq) Database at NCBI: Current Status, Taxonomic Expansion, and Functional Annotation. Nucleic Acids Res.44 (D1), D733–D745. 10.1093/nar/gkv1189
48
OrthJ. D.ConradT. M.NaJ.LermanJ. A.NamH.FeistA. M.et al (2011). A Comprehensive Genome‐scale Reconstruction of Escherichia coli Metabolism-2011. Mol. Syst. Biol.7, 535. 10.1038/msb.2011.65
49
Sabater-MuñozB.Prats-EscricheM.Montagud-MartínezR.López-CerdánA.ToftC.Aguilar-RodríguezJ.et al (2015). Fitness Trade-Offs Determine the Role of the Molecular Chaperonin GroEL in Buffering Mutations. Mol. Biol. Evol.32, 2681–2693. 10.1093/molbev/msv144
50
SamuelsonJ. C. (2011). Recent Developments in Difficult Protein Expression: a Guide to E. coli Strains, Promoters, and Relevant Host Mutations. Methods Mol. Biol.705, 195–209. 10.1007/978-1-61737-967-3_11
51
SchumannW. (2017). “Regulation of the Heat Shock Response in Bacteria,” in Prokaryotic Chaperonins: Multiple Copies and Multitude Functions. Editors KumarC. M. S.MandeS. C. (Singapore: Springer Singapore), 21–36. 10.1007/978-981-10-4651-3_2
52
SinghalK.VreedeJ.MashaghiA.TansS. J.BolhuisP. G. (2015). The Trigger Factor Chaperone Encapsulates and Stabilizes Partial Folds of Substrate Proteins. Plos Comput. Biol.11 (10), e1004444. 10.1371/journal.pcbi.1004444
53
SmithH. L. (2011). Bacterial Competition in Serial Transfer Culture. Math. Biosciences229 (2), 149–159. 10.1016/j.mbs.2010.12.001
54
TeterS. A.HouryW. A.AngD.TradlerT.RockabrandD.FischerG.et al (1999). Polypeptide Flux through Bacterial Hsp70. Cell97 (6), 755–765. 10.1016/s0092-8674(00)80787-4
55
TillyK.GeorgopoulosC. (1982). Evidence that the Two Escherichia coli groE Morphogenetic Gene Products Interact In Vivo. J. Bacteriol.149 (3), 1082–1088. 10.1128/jb.149.3.1082-1088.1982
56
TokurikiN.StricherF.SerranoL.TawfikD. S. (2008). How Protein Stability and New Functions Trade off. Plos Comput. Biol.4 (2), e1000002. 10.1371/journal.pcbi.1000002
57
TokurikiN.TawfikD. S. (2009). Chaperonin Overexpression Promotes Genetic Variation and Enzyme Evolution. Nature459 (7247), 668–673. 10.1038/nature08009
58
UniProtC. (2019). UniProt: a Worldwide Hub of Protein Knowledge. Nucleic Acids Res.47 (D1), D506–D515. 10.1093/nar/gky1049
59
Van DykT. K.GatenbyA. A.LaRossaR. A. (1989). Demonstration by Genetic Suppression of Interaction of GroE Products with many Proteins. Nature342 (6248), 451–453. 10.1038/342451a0
60
van OpijnenT.CamilliA. (2013). Transposon Insertion Sequencing: a New Tool for Systems-Level Analysis of Microorganisms. Nat. Rev. Microbiol.11 (7), 435–442. 10.1038/nrmicro3033
61
VulicM.KolterR. (2001). Evolutionary Cheating in Escherichia coli Stationary Phase Cultures. Genetics158 (2), 519–526.
62
WangM.WeissM.SimonovicM.HaertingerG.SchrimpfS. P.HengartnerM. O.et al (2012). PaxDb, a Database of Protein Abundance Averages across All Three Domains of Life. Mol. Cell Proteomics11 (8), 492–500. 10.1074/mcp.O111.014704
63
WeberJ.KayserA.RinasU. (2005). Metabolic Flux Analysis of Escherichia coli in Glucose-Limited Continuous Culture. II. Dynamic Response to Famine and Feast, Activation of the Methylglyoxal Pathway and Oscillatory Behaviour. Microbiology151 (Pt 3), 707–716. 10.1099/mic.0.27482-0
64
WeinerJ.3rdZimmermanC. U.GöhlmannH. W.HerrmannR. (2003). Transcription Profiles of the Bacterium Mycoplasma Pneumoniae Grown at Different Temperatures. Nucleic Acids Res.31 (21), 6306–6320. 10.1093/nar/gkg841
65
WilliamsT. A.FaresM. A. (2010). The Effect of Chaperonin Buffering on Protein Evolution. Genome Biol. Evol.2, 609–619. 10.1093/gbe/evq045
66
WyganowskiK. T.KaltenbachM.TokurikiN. (2013). GroEL/ES Buffering and Compensatory Mutations Promote Protein Evolution by Stabilizing Folding Intermediates. J. Mol. Biol.425 (18), 3403–3414. 10.1016/j.jmb.2013.06.028
67
ZambranoM.SiegeleD.AlmironM.TormoA.KolterR. (1993). Microbial Competition: Escherichia coli Mutants that Take over Stationary Phase Cultures. Science259 (5102), 1757–1760. 10.1126/science.7681219
68
ZouT.WilliamsN.OzkanS. B.GhoshK. (2014). Proteome Folding Kinetics Is Limited by Protein Halflife. PLoS One9 (11), e112701. 10.1371/journal.pone.0112701
Summary
Keywords
metabolic flux, GroEL, evolution, proteomics, metabolism, competitive index
Citation
Kumar CMS, Chugh K, Dutta A, Mahamkali V, Bose T, Mande SS, Mande SC and Lund PA (2021) Chaperonin Abundance Enhances Bacterial Fitness. Front. Mol. Biosci. 8:669996. doi: 10.3389/fmolb.2021.669996
Received
19 February 2021
Accepted
01 July 2021
Published
26 July 2021
Volume
8 - 2021
Edited by
Kürşad Turgay, Max-Planck-Gesellschaft (MPG), Germany
Reviewed by
Hideki Taguchi, Tokyo Institute of Technology, Japan
Marie-Pierre Castanié-Cornet, UMR5100 Laboratoire de Microbiologie et Génétique Moléculaires (LMGM), France
Updates
Copyright
© 2021 Kumar, Chugh, Dutta, Mahamkali, Bose, Mande, Mande and Lund.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: C. M. Santosh Kumar, s.k.cm@bham.ac.uk
This article was submitted to Protein Folding, Misfolding, and Degradation, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.