Introduction
The Deinococcus species are famous for their ability to survive extreme stresses, such as ultraviolet radiation, ionizing radiation, desiccation, and oxidative stress [2,8,9,20]. IrrE is a novel regulator protein (DR0167, inducer of pleiotropic proteins promoting DNA repair, also named pprI) and was discovered in the analysis of a DNA damage-sensitive strain [8,10]. Compared with wild-type D. radiodurans, the ΔirrE mutant displayed significantly altered phenotypes, including increased sensitivity to UV and γ radiations [23]. IrrE, as previously considered, is not only playing an important role in DNA repair probably by activating the expression of recA and pprA, but is also likely involved with transcriptional regulation, signal transduction, and proteolysis [15,16]. Lu et al. [15] revealed that the IrrE HTH domain can bind directly to the recA and pprA promoters, and disruption of irrE led to increased sensitivity to UV radiation.
D. gobiensis I-0 was isolated by our laboratory from the upper sand layers of the Gobi Desert, Xinjiang, China, where the environment rapidly changes and the bacteria are exposed to cycles of high and low temperatures and prolonged desiccation. Phylogenetic analysis of D. gobiensis I-0 (Fig. 1) indicated that it belongs to the genus Deinococcus. Interestingly, this strain exhibited stronger resistance to γ and UV radiations than D. radiodurans [26]. IrrE in D. gobiensis I-0 is highly conserved with the IrrE sequences in other Deinococcus strains and bears 66% identity to D. radiodurans. Moreover, Dgob-irrE exhibits stronger radiation tolerance than Drad-irrE when used to complement the D. radiodurans ΔirrE mutant.
Fig 1.Unrooted neighbor-joining phylogenetic tree from the alignment of the D. gobiensis I-0 IrrE protein with related proteins.
Crystal structure analysis of IrrE from D. deserti demonstrated three domains: an N-terminal peptidase-like domain, a central helix-turn-helix motif, and a GAF-like putative sensor domain [23]. Current research has shown that IrrE in D. deserti is a metalloprotease that cleaves the repressor protein DdrO [16], which is in accord with the N-terminal function. A previous study showed that expression of HTH domain-deficient Drad-IrrE resulted in increased sensitivity to γ radiation, UV radiation, and mitomycin C stress. Therefore, the IrrE HTH domain is essential to program the DNA repair process and cellular survival of D. radiodurans in response to radiation damage [15]. Previous research on the HTH domain focused on the consensus sequence using site-directed mutagenesis to mutate specific residues, such as Y160, but failed to identify any key residue [23]. However, their research methodology for examining critical sites provided us with an idea.
From a multisequence alignment of IrrE in combination with some fundamental experiments, we hypothesized that three different sites in the Dgob-IrrE HTH domain may be critical for its stronger radiation tolerance. Here, we found that the S131A mutation greatly decreased the radiation tolerance, and the corresponding A184S mutant displayedhigher resistance against UV radiation and desiccation. Taken together, these results suggest that the serine residue is the critical site in the Dgob-IrrE HTH domain, which may contribute to improving the affinity of DNA binding and the enhanced UV radiation tolerance of Drad-IrrE. Moreover, the modified A184S strain possesses a series of remarkable traits, which can be applied in the agricultural field and aid in the production of high-quality crops.
Materials and Methods
Bacterial Strains, Plasmids, and Growth Conditions
The bacterial strains and plasmids used in this study are listed in Table 1. The E. coli strain and its variants were maintained in Luria-Bertani (LB) medium at 37℃ in the presence of the appropriate antibiotics as required. The Deinococcus strains were grown in TGY (1% tryptone, 0.5% yeast extract, 0.1% glucose) medium at 30℃ in the presence of the appropriate antibiotics as required.
Table 1.Bacterial strains and plasmids used in this study.
Cloning and Site-Directed Mutagenesis
Plasmid constructions were carried out according to standard procedures. The D. radiodurans and D. gobiensis I-0 irrE genes (Drad-irrE and Dgob-irrE, respectively) based on the Genbank sequences were amplified with pMD18T-Drad-irrE and pMD18-TDgob-irrE, respectively, as the template and ligated into pRADZ3 using the restriction sites NdeI and SpeI. The resulting recombinant plasmids were named as pRADZ3-E and pRADZ3-IE. Mutagenesis was performed to introduce the following designed changes: A122P, S131A, and A184S. Six primers for each mutation were designed to contain the corresponding nucleotide changes (see Table 2). These primers and the pRADZ3-E and pRADZ3-IE templates were used to introduce mutations using the polymerase chain reaction (PCR) under the following conditions: one cycle of denaturation at 95℃ for 5 min, 30 cycles of denaturation at 95℃ for 30 sec, annealing at 60℃ for 30 sec, and extension at 72℃ for 50 sec, followed by an extra extension at 72℃ for 10 min. The PCR products were purified using the TianGen PCR purification kit and were ligated into the pARDZ3 vector. The nucleotide changes were then sequenced by the Beijing Genomics Institute.
Table 2.Bold residues indicate the mutation sites.
Transformation of the Constructed Plasmids into D. radiodurans ΔirrE Mutant Competent Cells
For the transformation, the ΔirrE mutant was grown to mid-exponential phase (optical density at 0.6-0.8) in the presence of spectinomycin (at 350 μg/ml) and then centrifuged at 6,000 ×g for 3 min. The pellet was suspended in 180 μl of fresh TGY liquid, and 20 μl of 0.3 M CaCl2 was added, and the mixture was incubated at 30℃ for approximately 90 min. A total of 50 μl of the resuspended pellet was mixed with 10 μl (0.03 pM) of the constructed plasmids and was incubated under the same conditions (at 30℃ for approximately 90 min). The total 60 μl volume was pipetted into 2 ml of TGY medium and cultivated with no antibiotics at 30℃ with shaking at 220 rpm overnight. A total of 200 μl of the overnight suspension culture was plated on TGY medium in the presence of chloromycetin (at 3.4 μg/ml) and then incubated for 3 days at 30℃ to select the targeted strains.
Abiotic Stress Resistance Experiments
The strains listed in Table 1 (except for E. coli strains grown in LB medium) were grown in TGY medium at 30℃ to an OD600 of 6. Then 20 ml of the bacteria suspension was centrifuged at 6,000 ×g for 5 min and washed twice with PBS, and then after 20 ml of PBS buffer was added, the mix was treated with UV radiation (254 nm, 5 min 380 J/m2, 10 min 760 J/m2, 15 min 1.14 kJ/m2, 20 min 1.52 kJ/m2). In the desiccation stress experiments, we collected cells when the OD600 reached 6-8 and divided them into many equal parts. We removed the same batch from the desiccator every 10 days and counted the surviving cells. In mitomycin C (MMC) shock experiments, the cells were collected when the OD600 reached 0.6-0.8, added to TGY liquid medium containing MMC (20 μg/ml) for 10 min, and then a diluted suspension was spotted onto TGY plates. All experiments were performed in triplicate.
Quantitative RT-PCR (qRT-PCR)
The bacteria were grown to an OD600 of 6, treated with UV radiation for approximately 10 min, and then harvested. Total RNA was isolated using an RNA extraction kit (FastRNA Pro Blue RNA; Qbiogene). The RNA samples were digested and reverse-transcribed using a cDNA synthesis kit (PrimeScriptRT reagent kit with gDNA Eraser; TaKaRa) according to the manufacturer’s protocol. The expression of recA, pprA, uvrA, and ddrB was determined using the ABI PRISM 7500 system. The primers were designed and optimized using primer software (primer length ranged from 18 to 24 nucleotides and the annealing temperature varied from 55℃ to 60℃). PCR amplification was detected using a SYBR fluorescence dye (TaKaRa). The 16S RNA gene served as an internal control to normalize for the differences in total RNA quantity. The qRT-PCR primers are shown in Table 3.
Table 3.List of primers used in the qRT-PCR study.
Sequence Alignment and Molecular Modeling
The sequences of Drad-IrrE, Dgob-IrrE, and other IrrE proteins were aligned using ClustalW ver. 2.0 and ESPript software [10]. A protein model was built using the online Swiss-Model program after input of the aligned sequences [19] based on the amino acid sequence homology. The three-dimensional (3D) structure of Dgob-IrrE was modeled, based on the reported crystal structure of IrrE from D. deserti (PDB ID 3DTI), which was visualized using the PyMOL software [7]. We also used site-directed mutagenesis and interaction bond manipulation to explain the different phenotypes in detail.
Results
Dgob-IrrE Enhanced the Tolerance of the ΔirrE Mutant Against UV Radiation
A previous study revealed that D. gobiensis I-0 exhibited strong resistance against γ and UV radiations [26]. Owing to the significant contribution of IrrE in the protection against γ and UV radiations, we carried out an experiment to investigate whether Dgob-IrrE improves resistance to UV radiation in the ΔirrE mutant compared with the wild-type DR and the complemented ΔirrE-E strain. We constructed two strains in the ΔirrE mutant background containing the shuttle plasmid pRADZ3 carrying the Dgob-irrE gene and Drad-irrE gene. Four strains were analyzed in this assay: wild-type DR, ΔirrE, ΔirrE-E, and ΔirrE-IE. The four strains were treated under high-intensity UV radiation (254 nm) for 5 min and then grown in TGY plates as shown in Fig.2A . Clearly, ΔirrE-IE displayed higher resistance to UV radiation than the others, and the survival curve (Fig.2B) was in agreement with this phenotype. We also examined other abiotic stresses, such as mitomycin C shock and desiccation; there was also stronger phenotype in ΔirrE-IE compared with ΔirrE-E and DR strains under MMC shock, but no differences under desiccation (data not shown).
Fig. 2.Effect of UV radiation on growth and survival curve of wild-type DR and its variant strains.
Multiple Alignment of IrrE and Computerized Analysis of Site-Directed Mutations
Based on the multiple alignment of the IrrE proteins shown in Fig. 3, there are three different sites in the HTH domain between Dgob-IrrE and Drad-IrrE. Both IrrE 3D structures were built on the basis of homology modeling and aligned and visualized using PyMOL software. The three different sites are marked with an orange sphere in Fig. 4A . The HTH domain of Drad-IrrE ranges from residues 163 to 200, and the region spanning 172-200 is conserved. We analyzed the conformational change of the site-directed A184S mutant (Fig. 4B). The results showed that 184S added double interaction bonds compared with A184, causing the conformational change that stably fastened the lysine residue. Lysine is involved the interaction with the DNA minor groove [18,21]. Therefore, we speculated that the conformational change of the lysine residue may increase the DNA-binding activity. The analysis of the other two sites is shown in Fig. S1.
Fig. 3.Multiple alignment of the IrrE proteins.
Fig. 4.3D structure of the IrrE proteins and the computer analysis of the site-directed mutagenesis.
The A184S Mutant Enhanced Resistance to UV Radiation and MMC Shock
The A184S mutant displayed stronger resistance against UV radiation and MMC shock than the wild-type DR strain. The results in time courses demonstrated that the A184S mutant had significantly increased tolerance against UV radiation compared with the ΔirrE-E strain (Fig.5 A). The survival curve of A184S (Fig. S2A) was in accordance with the growth on TGY plate under UV radiation. A184S also performed stronger tolerance with the treatment of MMC (Fig. 6A). The corresponding site of A184 in Dgob-IrrE i s S131 a nd w as m utated i nto S131A. We also generated the A122P mutant, which is another different site in the Dgob-IrrE HTH domain. Both the A122P and S131A mutants lost resistance against UV radiation (Fig.5B) in the first 5 min coinciding with the survival curve (Fig. S2B). The S131A mutant also showed a dramatic loss of MMC shock (Fig. 6B) and at 50 days desiccation (Fig.S4 ). All of these results indicated that the serine residue is critical for the HTH domain and may cause the conformational change leading to the higher DNA-binding efficiency to improve the expression level of IrrE-dependent genes.
Fig. 5.A184S mutation effect on resistance to UV radiation.
Fig. 6.A184S mutation effect on resistance to MMC shock.
Quantitative RT-PCR Revealed that Crucial Genes of Excision and Recombination Repair Pathways were Induced Significantly in the ΔirrE-IE and A184S Mutants
It was previously reported that both the excision and recombination repair pathways were the main pathways employed by D. radiodurans against UV radiation [24]. Some of the crucial genes involved in both pathways are recA (DR_2340), pprA (DR_A0346), and uvrA (DR_1771). RecA is required for extended synthesis-dependent strand annealing (EMSA) and is essential for radioresistance [3,6,20,27]. PprA was implicated in oxidative stress protection and interacted with double-stranded DNA and DNA gyrase to facilitate the radioresistance of D. radiodurans [13, 14]. UvrA is one of the key components of the nucleotide excision repair (NER) pathway responsible for the repair and removal of damaged DNA lesions caused by UV radiation [22]. In addition, DdrB has been reported to play a critical role at early times after radiation [25]. qRT-PCR results showed that all four genes were induced after 10 min UV radiation (Fig. 7). The recA gene in the ΔirrE-IE and A184S mutants were expressed greater than 1.5-fold in DR and ΔirrE-E (Fig. 7A). It is evident that these four genes in the A184S mutant displayed higher induction than in the ΔirrE-E mutant. qRT-PCR results also revealed that recA was also induced in A184S and suppressed in S131A under MMC shock, as shown in Fig. S5, indicating the importance of the serine residue.
Fig. 7.Relative expression of recA (A), pprA (B), uvrA (C), and ddrB (D) after 10 min UV radiation.
Discussion
A previous study by our laboratory revealed that D. gobiensis I-0 displayed a stronger tolerance than D. radiodurans against UV radiation [26], but the cause of this D. gobiensis I-0 phenotype has not been demonstrated. IrrE, as a general switch, plays a significant role in the extreme radioresistance of Deinococcus, including gamma radiation and UV radiation. Both stresses cause the similar impairment of DNA [12]. The IrrE protein sequence and biological function among the Deinococcus genus is highly conserved [23]. Our work revealed that Dgob-IrrE was not only fully restored but also demonstrated stronger radioresistance in the ΔirrE mutant than Drad-IrrE did against UV radiation. The previous work of Lu et al. [15] showed that the HTH domain of Drad-IrrE specifically bound to the promoters of recA and pprA; both are important elements involved in the DNA repair pathway. The HTH domain-deficient mutant only partially restored resistance of the ΔirrE mutant to γ radiation, UV radiation, and mitomycin C. These results indicated that the HTH domain is essential for IrrE in response to radiation damage [15]. Therefore, we speculated that the stronger resistance of Dgob-IrrE was possibly due to different amino acids in the HTH domain. We undertook a series of biological experiments and bioinformatics analyses to identify the crucial sites.
The crystal structure of D. deserti IrrE was previously solved and was shown to contain an N-terminal peptidase-like domain, a central helix-turn-helix motif, and a GAF-like putative sensor domain [23]. We established the 3D structure of Drad-IrrE and Dgob-IrrE using the online Swiss-Model program based on homology modeling, and we improved the model quality using Chiron and SAVEs software to remove the clashes, as shown in Fig. 4A. The three different sites, A122, S131, and T140, were marked with spheres. Compared with the A184 site, the mutated S184 site added two interaction bonds with the K183 residue, as shown Fig.4 B. A recent report showed that 80% of all protein interactions with the DNA minor groove at AT-tracts are mediated by arginine or lysine [18]. The additional interaction bonds of S184 with K183 may cause a conformational change and increase DNA-binding activity. Conversely, serine is a hydrophilic amino acid and alanine is a hydrophobic residue, and this site-directed mutagenesis also increased the DNA-binding affinity. The third site, T140 of Dgob-IrrE and S193 of Drad-IrrE, had similar chemical properties (Fig. S1), which were not processed by site-directed mutagenesis.
Based on the work of Lu et al. [15], we performed qRT-PCR to identify the regulation of IrrE, as shown in Fig. 7. The genes recA, pprA, uvrA, and ddrB are important in the DNA repair pathway and displayed significant downregulation in the S131A mutant and up-regulation in the A184S mutant under UV radiation and MMC shock. These results verified the importance of the serine residue for the HTH domain. Current research showed that the N-terminal region is a metalloprotease, and it can cleave the repressor protein DdrO, leading to the transcriptional induction of various genes required for repair and survival after exposure of Deinococcus to radiation [16]. We found that the G19 site of the N-terminal region is different among the conserved α2-helix region and has yet to be characterized. The C-terminal domain, also named GAF domain, is a small molecule-binding domain present in many signaling and sensory proteins [11]. The GAF domain contains a V184 site, which is different from the conserved sequences. The site-directed mutant V184A displayed no phenotypic changes (Fig.S3 ) against UV radiation and MMC shock, indicating that V184 may not be the key residue in the GAF domain and is not the reason for the stronger resistance of Dgob-IrrE. Identifying the change in the IrrE-DNA interaction by NMR may provide further information about the HTH domain residues.
Drad-IrrE has been investigated for several years. The tolerance of E. coli harboring Drad-IrrE to osmotic stress, heat stress, and oxidative stress was significantly enhanced [4]. Z. Zhou’s early work revealed that the Drad-IrrE protein can enhance E. coli salt tolerance, and his team transformed the irrE gene into Brassica napus, which improved growth of the plant in high-salinity soil [17]. Interestingly, E. coli strains carrying variants of IrrE selected from a library of randomly generated mutants, such as E1, E28, E30, E79, and E80, conferred enhanced tolerance toward ethanol, butanol, and acetate [5]. These results demonstrate the importance of Drad-IrrE application, so it is meaningful and necessary to generate IrrE variants with increased activity against extreme environmental conditions. The A184S mutant displayed stronger resistance, providing new insights into modifications of the IrrE protein. Therefore, we can potentially apply this excellent IrrE variant to engineered crops to enhance growth traits under extreme environmental conditions.
References
- Adachi M, Hirayama H, Shimizu R, Satoh K, Narumi I, Kuroki R. 2014. Interaction of double-stranded DNA with polymerized PprA protein from Deinococcus radiodurans. Protein Sci. 23: 1349-1358. https://doi.org/10.1002/pro.2519
- Battista JR, Earl AM, Park M-J. 1999. Why is Deinococcus radiodurans so resistant to ionizing radiation? Trends Microbiol. 7: 362-365. https://doi.org/10.1016/S0966-842X(99)01566-8
- Bentchikou E, Servant P, Coste G, Sommer S. 2010. A major role of the RecFOR pathway in DNA double-strand-break repair through ESDSA in Deinococcus radiodurans. PLos Genet. 6. https://doi.org/10.1371/journal.pgen.1000774
- Chen T, Wang J, Yang R, Li J, Lin M, Lin Z. 2011. Laboratory-evolved mutants of an exogenous global regulator, IrrE from Deinococcus radiodurans, enhance stress tolerances of Escherichia coli. PLoS One 6: e16228. https://doi.org/10.1371/journal.pone.0016228
- Chen T, Wang J, Zeng L, Li R, Li J, Chen Y, Lin Z. 2012. Significant rewiring of the transcriptome and proteome of an Escherichia coli strain harboring a tailored exogenous global regulator IrrE. PLoS One 7: e37126. https://doi.org/10.1371/journal.pone.0037126
- Daly MJ, Minton KW. 1996. An alternative pathway of recombination of chromosomal fragments precedes recA-dependent recombination in the radioresistant bacterium Deinococcus radiodurans. J. Bacteriol. 178: 4461-4471. https://doi.org/10.1128/jb.178.15.4461-4471.1996
- DeLano WL. 2002. The PyMOL molecular graphics system. DeLano Scientific, CA.
- Earl AM, Mohundro MM, Mian IS, Battista JR. 2002. The IrrE protein of Deinococcus radiodurans R1 is a novel regulator of recA expression. J. Bacteriol. 184: 6216-6224. https://doi.org/10.1128/JB.184.22.6216-6224.2002
- Eisen JA, Hanawalt PC. 1999. A phylogenomic study of DNA repair genes, proteins, and processes. Mutat. Res. DNA Repair 435: 171-213. https://doi.org/10.1016/S0921-8777(99)00050-6
- Gouet P, Robert X, Courcelle E. 2003. ESPript/ENDscript: extracting and rendering sequence and 3D information from atomic structures of proteins. Nucleic Acids Res. 31: 3320-3323. https://doi.org/10.1093/nar/gkg556
- Ho YSJ, Burden LM, Hurley JH. 2000. Structure of the GAF domain, a ubiquitous signaling motif and a new class of cyclic GMP receptor. EMBO J. 19: 5288-5299. https://doi.org/10.1093/emboj/19.20.5288
- Hua Y, Narumi I, Gao G, Tian B, Satoh K, Kitayama S, Shen B. 2003. PprI: a general switch responsible for extreme radioresistance of Deinococcus radiodurans. Biochem. Biophys. Res. Commun. 306: 354-360. https://doi.org/10.1016/S0006-291X(03)00965-3
- Kota S, Charaka VK, Ringgaard S, Waldor MK, Misra HS. 2014. PprA contributes to Deinococcus radiodurans resistance to nalidixic acid, genome maintenance after DNA damage and interacts with deinococcal topoisomerases. PLoS One 9: e85288. https://doi.org/10.1371/journal.pone.0085288
- Kota S, Misra HS. 2006. PprA: a protein implicated in radioresistance of Deinococcus radiodurans stimulates catalase activity in Escherichia coli. Appl. Microbiol. Biotechnol. 72: 790-796. https://doi.org/10.1007/s00253-006-0340-7
- Lu H, Chen H, Xu G, Shah AM-U-H, Hua Y. 2012. DNA binding is essential for PprI function in response to radiation damage in Deinococcus radiodurans. DNA Repair 11: 139-145. https://doi.org/10.1016/j.dnarep.2011.10.013
- Ludanyi M, Blanchard L, Dulermo R, Brandelet G, Bellanger L, Pignol D, et al. 2014. Radiation response in Deinococcus deserti: IrrE is a metalloprotease that cleaves repressor protein DdrO. Mol. Microbiol. 94: 434-449. https://doi.org/10.1111/mmi.12774
- Pan J, Wang J, Zhou Z, Yan Y, Zhang W, Lu W, et al. 2009. IrrE, a global regulator of extreme radiation resistance in Deinococcus radiodurans, enhances salt tolerance in Escherichia coli and Brassica napus. PLoS One 4: e4422. https://doi.org/10.1371/journal.pone.0004422
- Rohs R, West SM, Sosinsky A, Liu P, Mann RS, Honig B. 2009. The role of DNA shape in protein-DNA recognition. Nature 461: 1248-1253. https://doi.org/10.1038/nature08473
- Schwede T, Kopp J, Guex N, Peitsch MC. 2003. SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 31: 3381-3385. https://doi.org/10.1093/nar/gkg520
- Slade D, Radman M. 2011. Oxidative stress resistance in Deinococcus radiodurans. Microbiol. Mol. Biol. Rev. MMBR 75: 133-191. https://doi.org/10.1128/MMBR.00015-10
- Tawfik D. 2002. The Protein Protocols Handbook. Springer, NY.
- Truglio JJ, Croteau DL, Van Houten B, Kisker C. 2006. Prokaryotic nucleotide excision repair: the UvrABC system. Chem. Rev. 106: 233-252. https://doi.org/10.1021/cr040471u
- Vujicic-Zagar A, Dulermo R, Le Gorrec M, Vannier F, Servant P, Sommer S, et al. 2009. Crystal structure of the IrrE protein, a central regulator of DNA damage repair in Deinococcaceae. J. Mol. Biol. 386: 704-716. https://doi.org/10.1016/j.jmb.2008.12.062
- White O, Eisen JA, Heidelberg JF, Hickey EK, Peterson JD, Dodson RJ, et al. 1999. Genome sequence of the radioresistant bacterium Deinococcus radiodurans R1. Science 286: 1571-1577. https://doi.org/10.1126/science.286.5444.1571
- Xu G, Lu H, Wang L, Chen H, Xu Z, Hu Y, et al. 2010. DdrB stimulates single-stranded DNA annealing and facilitates RecA-independent DNA repair in Deinococcus radiodurans. DNA Repair 9: 805-812. https://doi.org/10.1016/j.dnarep.2010.04.006
- Yuan M, Zhang W, Dai S, Wu J, Wang Y, Tao T, et al. 2009. Deinococcus gobiensis sp. nov., an extremely radiation-resistant bacterium. Int. J. Syst. Evol. Microbiol. 59: 1513-1517. https://doi.org/10.1099/ijs.0.004523-0
- Zahradka K, Slade D, Bailone A, Sommer S, Averbeck D, Petranovic M, et al. 2006. Reassembly of shattered chromosomes in Deinococcus radiodurans. Nature 443: 569-573.
Cited by
- A tamB homolog is involved in maintenance of cell envelope integrity and stress resistance of Deinococcus radiodurans vol.7, pp.None, 2015, https://doi.org/10.1038/srep45929
- Improving Biotransformation Efficiency of Arthrobacter simplex by Enhancement of Cell Stress Tolerance and Enzyme Activity vol.69, pp.2, 2015, https://doi.org/10.1021/acs.jafc.0c06592