The Richness and Diversity of Catalases in Bacteria

Catalases play a key role in the defense against oxidative stress in bacteria by catalyzing the decomposition of H₂O₂. In addition, catalases are also involved in multiple cellular processes, such as cell development and differentiation, as well as metabolite production. However, little is known about the abundance, diversity, and distribution of catalases in bacteria. In this study, we systematically surveyed and classified the homologs of three catalase families from 2,634 bacterial genomes. It was found that both of the typical catalase and Mn-catalase families could be divided into distinct groups, while the catalase-peroxidase homologs formed a tight family. The typical catalases are rich in all the analyzed bacterial phyla except Chlorobi, in which the catalase-peroxidases are dominant. Catalase-peroxidases are rich in many phyla, but lacking in Deinococcus-Thermus, Spirochetes, and Firmicutes. Mn-catalases are found mainly in Firmicutes and Deinococcus-Thermus, but are rare in many other phyla. Given the fact that catalases were reported to be involved in secondary metabolite biosynthesis in several Streptomyces strains, the distribution of catalases in the genus Streptomyces was given more attention herein. On average, there are 2.99 typical catalases and 0.99 catalase-peroxidases in each Streptomyces genome, while no Mn-catalases were identified. To understand detailed properties of catalases in Streptomyces, we characterized all the five typical catalases from S. rimosus ATCC 10970, the oxytetracycline-producing strain. The five catalases showed typical catalase activity, but possessed different catalytic properties. Our findings contribute to the more detailed classification of catalases and facilitate further studies about their physiological roles in secondary metabolite biosynthesis and other cellular processes, which might facilitate the yield improvement of valuable secondary metabolites in engineered bacteria.

Introduction

Reactive oxygen species (ROS) are generated in the O₂ reduction process, including superoxide anion (O₂^⋅–), hydrogen peroxide (H₂O₂), and hydroxyl radical (OH^⋅) (Montibus et al., 2015; Johnson and Hug, 2019; Kim et al., 2019). These ROS can cause the oxidative damage of cellular macromolecules and lead to dysfunction of cells (Montibus et al., 2015; Johnson and Hug, 2019). Due to the pervasive threat of ROS, bacteria have evolved various scavenging and defense enzymes to mitigate this threat (Johnson and Hug, 2019; Kim et al., 2019). Catalases are one of the families of well-known ROS defense enzymes, which catalyze the decomposition of H₂O₂ (Nicholls et al., 2000; Chelikani et al., 2004; Zamocky et al., 2008; Nicholls, 2012; Johnson and Hug, 2019). In nature, there are three families of catalases with different protein sequence characteristics and co-factor dependences, including the heme-containing monofunctional typical catalases, heme-containing bifunctional catalase-peroxidases, and manganese-containing catalases (Mn-catalases) (Chelikani et al., 2004; Zamocky et al., 2008; Johnson and Hug, 2019).

The typical catalases catalyze dismutation of H₂O₂ through a two-stage mechanism (Reactions 1 and 2) (Chelikani et al., 2004; Zamocky et al., 2008; Sooch et al., 2014). In the first stage, the heme iron in the active center is oxidized by a molecule of H₂O₂ to form the intermediate oxoferryl porphyrin π-cation radical (^⋅+Por-Fe(IV)=O, compound I), together with the production of water (Reaction 1). In the second stage, compound I is rapidly reduced back to ferrous porphyrin form by H₂O₂, with the formation of molecular oxygen and water (Reaction 2). However, the detailed mechanism of the two-electron reduction of compound I by H₂O₂ was still not very clear (Alfonso-Prieto et al., 2012). Previous studies proposed that the reduction process involved the transfer of a hydride ion from H₂O₂ to compound I (Fita and Rossmann, 1985; Kato et al., 2004; Alfonso-Prieto et al., 2012), while the metadynamics simulation analysis suggested that the transfer of one hydrogen atom from H₂O₂ to compound I occurred first, followed by another reduction reaction (Kato et al., 2004; Alfonso-Prieto et al., 2009, 2012).

$\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{III})+{\mathrm{H}}_2{\mathrm{O}}_2{\to }^{\cdot +}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})=\mathrm{O}+{\mathrm{H}}_2\mathrm{O}$

$(\mathrm{Reaction}\mathrm{\hspace{0.17em}1})$

${}^{\cdot +}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})=\mathrm{O}+{\mathrm{H}}_2{\mathrm{O}}_2\to \mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{III})+{\mathrm{O}}_2+{\mathrm{H}}_2\mathrm{O}$

$(\mathrm{Reaction}\mathrm{\hspace{0.17em}2})$

The bifunctional catalase-peroxidases catalyze dismutation of H₂O₂, and decomposition of H₂O₂ using some unidentified organic electron donors (Chelikani et al., 2004; Zamocky et al., 2008). The catalytic mechanism of catalase-peroxidase is complicated and involves several intermediates (Zamocky et al., 2008; Vlasits et al., 2010). Briefly, the heme iron is first oxidized by H₂O₂ to form compound I (Reaction 1). Compound I can be reduced to compound II (Por-Fe(IV)-OH) via transfer of a hydrogen atom from organic electron donors (Reaction 3), followed by another hydrogen atom transfer with the production of water (Reaction 4). In the absence of exogenous electron donor, the porphyrin radical of compound I can be quenched by an electron transferred from protein residues forming compound I^∗ (^{⋅ +} AA Por-Fe(IV)-OH; Reaction 5) (Zamocky et al., 2008; Vlasits et al., 2010), which is reduced by another molecule of H₂O₂ to form compound III^∗ ([^{⋅ +} AA Por-Fe(II)-O₂ ↔ ^{⋅ +} AA Por-Fe(III)-O₂^⋅–]; Reaction 6) (Vlasits et al., 2010). After the release of molecular oxygen, the enzyme returns to its ferric form (Reaction 7).

${}^{\cdot +}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})=\mathrm{O}+{\mathrm{AH}}_2\to \mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})-\mathrm{OH}+{}^{\cdot }\mathrm{AH}$

$(\mathrm{Reaction}\mathrm{\hspace{0.17em}3})$

$\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})-\mathrm{OH}+{\mathrm{AH}}_2\to \mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{III})+{}^{\cdot }\mathrm{AH}+{\mathrm{H}}_2\mathrm{O}$

$(\mathrm{Reaction}\mathrm{\hspace{0.17em}4})$

${\mathrm{AA}}^{\cdot +}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})=\mathrm{O}+{\mathrm{H}}^{+}{\to }^{\cdot +}\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})-\mathrm{OH}$

$(\mathrm{Reaction}\mathrm{\hspace{0.17em}5})$

${}^{\cdot +}\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{IV})-\mathrm{OH}+{\mathrm{H}}_2{\mathrm{O}}_2\to {[}^{\cdot +}\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{II})-{\mathrm{O}}_2$

${\leftrightarrow }^{\cdot +}\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{III})-{\mathrm{O}}_2^{\cdot -}]+\mathrm{H}{}_2\mathrm{O}+\mathrm{H}{}^{+}(\mathrm{Reaction6})$

${[}^{\cdot +}\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{II})-{\mathrm{O}}_2{\leftrightarrow }^{\cdot +}\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{III})-{\mathrm{O}}_2^{\cdot -}]\to$

$\mathrm{AA}\mathrm{Por}\text{-}\mathrm{Fe}(\mathrm{III})+{\mathrm{O}}_2(\mathrm{Reaction}\mathrm{\hspace{0.17em}7})$

The Mn-catalases are found mainly in bacteria and archaea, which catalyze dismutation of H₂O₂ but with relatively low specific activities compared to the heme-containing typical catalases (Chelikani et al., 2004; Zamocky et al., 2008; Whittaker, 2012). The Mn-catalase has a unique di-manganese cluster (Mn²⁺-Mn²⁺ or Mn³⁺-Mn³⁺) in the active center, and catalyzes a different two-stage reaction (Chelikani et al., 2004; Zamocky et al., 2008; Whittaker, 2012; Sooch et al., 2014). The Mn²⁺-Mn²⁺ cluster is first oxidized by a molecule of H₂O₂ to form the Mn³⁺-Mn³⁺ cluster and two molecules of water (Reaction 8). And then, the Mn³⁺-Mn³⁺ cluster is reduced by another molecule of H₂O₂ and releases molecular oxygen (Reaction 9) (Chelikani et al., 2004; Zamocky et al., 2008).

${\mathrm{Mn}}^{2+}\text{-}{\mathrm{Mn}}^{2+}+2{\mathrm{H}}^{+}+{\mathrm{H}}_2{\mathrm{O}}_2\to {\mathrm{Mn}}^{3+}\text{-}{\mathrm{Mn}}^{3+}$

$+2{\mathrm{H}}_2\mathrm{O}(\mathrm{Reaction}\mathrm{\hspace{0.17em}8})$

${\mathrm{Mn}}^{3+}\text{-}{\mathrm{Mn}}^{3+}+{\mathrm{H}}_2{\mathrm{O}}_2\to {\mathrm{Mn}}^{2+}\text{-}{\mathrm{Mn}}^{2+}+2{\mathrm{H}}^{+}$

$+{\mathrm{O}}_2(\mathrm{Reaction}\mathrm{\hspace{0.17em}9})$

As the major H₂O₂ scavenger, catalases play an important role in the defense against oxidative stress and the maintenance of intracellular H₂O₂ concentrations in bacteria (Montibus et al., 2015; Johnson and Hug, 2019). It is common that multiple catalases are encoded in one bacterium (Mulvey et al., 1990; Engelmann and Hecker, 1996; Cho et al., 2000; Klotz and Loewen, 2003) and they have been shown to behave differently in response to oxidative stress in some studies (Chelikani et al., 2004; Zamocky et al., 2008). For example, in E. coli, the catalase-peroxidase KatG was mainly expressed upon induction of oxidative stress, whereas the typical catalase KatE (HPII) was induced in the stationary growth phase (Mulvey et al., 1990; Zamocky et al., 2008). Studies have also showed that manipulating the expression of catalases could affect the development and secondary metabolite production in strains like streptomycetes (Cho et al., 2000; Beites et al., 2011). Disruption of catB (encoding a typical catalase) in S. coelicolor caused impairment in the formation of aerial mycelium and resistance against osmotic stress, and interestingly led to the reduced production of undecylprodigiosin and hyperproduction of actinorhodin (Cho et al., 2000). On the contrary, disruption of catA (encoding a typical catalase) showed no observed effect on the sporulation of S. coelicolor, but significantly reduced its tolerance to H₂O₂ stress (Cho et al., 2000). These previous studies highlighted the important physiological roles of catalases in bacteria. However, there is a lack of systematic understanding of the richness, diversity, and distribution of catalases in bacteria. Previous studies briefly surveyed the catalases in bacteria, archaea and eukarya, but only a very limited number of catalase homologs (usually 100-200 homologs for each family) were included for the analysis at that time (von Ossowski et al., 1993; Klotz et al., 1997; Klotz and Loewen, 2003; Chelikani et al., 2004; Zamocky et al., 2008; Zámocký et al., 2012).

In this study, we systematically surveyed and analyzed the homologs of three catalase families in bacterial genomes, which enabled us to reveal the richness, diversity, and taxonomic distribution of catalases in bacteria. We further characterized five typical catalases from one representative Streptomyces strain, all of which could catalyze the decomposition of H₂O₂ but with varied catalytic properties. Our study has provided the basis for further investigation of bacterial catalases to better understand their physiological roles in various cellular processes.

Materials and Methods

Bacterial Strains and Growth Conditions

S. rimosus ATCC 10970 was obtained from China General Microbiological Culture Collection Center (CGMCC 4.1438). E. coli JM109 and BL21(DE3) were used for general cloning and protein expression, respectively. LB medium was used for E. coli cultivation. Kanamycin was used at a final concentration of 50 μg/mL. Restriction enzymes and Q5 DNA polymerase were purchased from New England Biolabs (United States). DNA manipulations, competent cell preparation, and transformation were performed as described previously (Sambrook and Russell, 2001). Hydrogen peroxide (30%) solution was purchased from Sigma-Aldrich. The 96-well UV-transparent microplates were purchased from Corning.

Expression and Purification of the Five Catalases

The plasmids and primers used in this study are listed in Supplementary Tables 1, 2, respectively. The five catalase genes of S. rimosus ATCC 10970 were amplified through PCR with their corresponding primer pairs (Supplementary Table 2) using the genomic DNA of S. rimosus as a template. The plasmid pET28a was digested by NdeI and XhoI, and then joined with each of the five catalase gene fragments by Gibson assembly (Gibson et al., 2009) to generate the five overexpression plasmids (Supplementary Table 1). All the plasmids were verified by DNA sequencing and transformed into E. coli BL21(DE3) to obtain the five overexpression strains.

The strains were cultivated in LB supplemented with 50 μg/mL of kanamycin at 37°C and 220 rpm until an OD600 of 0.6 was reached, at which time final concentrations of 100 μM isopropyl-β-D-thiogalactoside (IPTG), 0.5 mM 5-aminolevulinic acid, and 0.2 mM FeCl₃ were added. The cell cultures were further cultivated for 8 h at 28°C and 180 rpm. The five N-terminal His₆-tagged catalases (Cat1-5) were purified by Nickel-NTA affinity chromatography according to the previously reported protocol (Fan et al., 2012), concentrated by centrifugation, and exchanged to ultrafiltration buffer (50 mM Na₂HPO₄-NaH₂PO₄, 10% (v:v) glycerol, pH 7.4) using 10-kDa Amicon Ultra tubes (Centriplus, Merck Millipore). The purities of the proteins were evaluated by 10% SDS-PAGE. The protein concentrations were determined by Bradford assay using bovine serum albumin (BSA) as a standard. The proteins were stored at –70°C.

Heme Occupancy and Thermostability of the Five Catalases by Spectroscopic Analysis

The UV-vis spectra of five catalases were collected on a Beckman Coulter DU800 spectrophotometer using a cell with an optical path of 1 cm. The heme concentrations were calculated using absorption at 406 nm (ε₄₀₆ = 102 mM^–1 cm^–1) (Brill and Sandberg, 1968). The heme occupancy of each catalase was the ratio of heme concentration to the corresponding protein concentration (by Bradford assay).

Circular dichroism (CD) spectra of five catalases were monitored in the far-UV range (195–260 nm) in 50 mM Na₂HPO₄-NaH₂PO₄ (pH 7.4), at 30°C on a Chirascan^TM CD Spectrometer (Applied Photophysics, United Kingdom) using a cell with an optical path of 1 mm. The stability of the five catalases were monitored by CD signals at 222 nm (representing the content of α-helices) from 20 to 90°C (raising 1°C at a time, and holding for 5 s before data collection) in 50 mM Na₂HPO₄-NaH₂PO₄ (pH 7.4).

Enzymatic Reactions of the Five Catalases

The activities of the five catalases were measured using the spectrophotometric method described previously (Beers and Sizer, 1952). A typical 200 μL reaction mixture consisted of 16 mM H₂O₂, purified catalase (1.06-10.32 nM), and 50 mM Na₂HPO₄-NaH₂PO₄ (pH 7.4). The absorbance at 240 nm (representing the concentration of H₂O₂) was monitored continuously in 96-well plates at 30°C using an EnSpire^TM multimode Plate Reader (PerkinElmer Inc., United States). The reaction rates were calculated using the absorbance changes at 240 nm within the first 1 min, and averaged for three independent reaction samples. One unit of enzyme activity is defined as the amount of activity required for the conversion of 1 μmol of H₂O₂ into water and oxygen per minute at 30°C (Zeng et al., 2011).

The catalase activities under different pH conditions were measured at 30°C. A series of buffers were used, including 50 mM NaAc-HAc (pH 4.0 – 6.0), 50 mM NaH₂PO₄-Na₂HPO₄ (pH 6.0 – 8.0), and 50 mM Na₂CO₃-NaHCO₃ (pH 8.0 - 10.0). The reaction rates were calculated using the absorbance changes at 240 nm within the first 1 min. The concentrations of the enzymes used in the assay were Cat1: 1.06 nM, Cat2: 1.35 nM, Cat3: 3.83 nM, Cat4: 1.11 nM, and Cat5: 10.32 nM.

The catalase activities at different temperatures (30 to 65°C) were measured in 50 mM NaH₂PO₄-Na₂HPO₄ (pH 7.4). The reaction rates were calculated using the absorbance changes at 240 nm within the first 1 min. The concentrations of the enzymes used in the assay were Cat1: 1.06 nM, Cat2: 1.35 nM, Cat3: 3.83 nM, Cat4: 1.11 nM, and Cat5: 10.32 nM.

The apparent K_m and k_cat values of the five catalases were determined by non-linear least squares fitting using reaction rates with different H₂O₂ concentrations (7.5 - 50 mM) at 30°C in 50 mM NaH₂PO₄-Na₂HPO₄ (pH 7.4). The specific activities of the five catalases were determined using the reaction rates when 16 mM H₂O₂ was used as the substrate. The reaction rates were calculated using the absorbance changes at 240 nm within the first 1 min. The concentrations of the holoenzymes (deduced from the detected heme occupancy of each purified catalase) used in the reactions were Cat1: 1.15 nM, Cat2: 1.17 nM, Cat3: 1.26 nM, Cat4: 1.05 nM, and Cat5: 6.05 nM.

Bioinformatic Analysis of Catalase Homologs

The catalase homologs were searched using standalone BLASTP program (Camacho et al., 2009) with the 21 known catalases reported previously as query sequences (Supplementary Table 3; Zamocky et al., 2008). Pfam search was performed using HMMER 3.3 (Eddy, 2011) with reported Pfam motifs (typical catalase: PF00199, catalase-peroxidase: PF00141, Mn-catalase: PF05067) (Zamocky et al., 2008; El-Gebali et al., 2019). Sequence comparison was performed using EMBOSS 6.5.0 (Rice et al., 2000). Those homologs with Pfam scores lower than the trusted cutoffs of corresponding Pfam motifs were removed, along with those with shorter lengths (typical catalase clade 1: <450 residues, clade 2: <600 residues, clade 3: <450 residues, catalase-peroxidase: <600 residues, Mn-catalase: <180 residues). The phylogenetic trees were constructed by a neighbor-joining method using MEGA X with the default parameters (Kumar et al., 2018). The information on all catalase homologs, including their GenBank accession numbers and sequences, are provided as Supplementary Dataset 1; the phylogenic trees of three families of catalases are provided as Newick files in Supplementary Dataset 2.

Results

Phylogenetic Analysis Revealing the Diversity of Catalases in Bacteria

To explore the abundance and diversity of catalases in bacteria, we carried out a virtual survey of 13,360 completed bacterial genomes in the NCBI Genomes database (as of Sep 2020) using the 21 well studied catalases covering all three families (Supplementary Table 3) as query sequences. As the numbers of sequenced genomes varied dramatically among different taxa, no more than 10 genomes from each genus were randomly selected as representatives for further analysis. As a result, a total of 2,634 genomes were selected, in which 4,931 catalase homologs were identified. The homologs were further filtered by removing 71 proteins with Pfam scores lower than the trusted cutoffs or with abnormal protein sequence lengths (Supplementary Figure 1). Pairwise sequence comparisons of the final 4,860 catalase homologs clearly showed the three distinct protein families consisting of the typical catalases, catalase-peroxidases, and Mn-catalases (Figure 1A).

FIGURE 1

Figure 1. Sequence comparison and phylogenetic analysis of three catalase families. (A) Pairwise sequence identities of all 4,860 catalase homologs. Phylogenetic trees of (B) typical catalases, (C) catalase-peroxidases, and (D) Mn-catalases. The inner strip chart indicates the phylum (red: Actinobacteria, orange: Firmicutes, yellow: Bacteroidetes, cyan: Proteobacteria, green: Chloroflexi, blue: Cyanobacteria, magenta: Acidobacteria, black: others). The middle strip chart denotes the protein length with the corresponding color maps. The outer strip chart indicates the groups of each catalase family.

We next subjected the proteins of each catalase family to phylogenetic analysis, which revealed diverse groups within each family (Figure 1). In the phylogenetic tree, 2,615 typical catalase homologs fell into seven groups (group 1-7, Figure 1B). The proteins of group 5 have apparently larger sizes ranging from 600 to 850 residues. They showed relatively lower sequence similarity to the homologs from other groups (Supplementary Figure 2A). The proteins from groups 6 and 7 have 450-550 residues. The group 6 catalase homologs are mainly from Actinobacteria, while a small portion of them are from Proteobacteria and Firmicutes. However, the majority of group 7 proteins are from Proteobacteria, with a small minority of them coming from Bacteroidetes and Firmicutes. The proteins of groups 2, 3, and 4 are predominantly 500-600 residues in size. Among them, the group 4 catalase homologs are from Proteobacteria, while the proteins of groups 2 and 3 come from various phyla (Figure 1B). The proteins from group 1 are shorter (450-500 residues), and display relatively lower sequence identity to the homologs from other groups of this family (Supplementary Figure 2A). The proteins of this group are mainly from Firmicutes.

Phylogenetic analysis of the 1,681 catalase-peroxidase homologs appeared to form four major groups (Figure 1C). However, these groups are closely related, as they are all approximately 700-850 residues in size, and share high sequence identity (Figure 1C and Supplementary Figures 1B, 2B), forming a tight catalase-peroxidase protein family. Among them, the proteins in group 3 are mainly from Proteobacteria, and a minority of them come from Bacteroidetes and Firmicutes. The proteins in group 4 are from Actinobacteria, Proteobacteria, and Bacteroidetes. The proteins in group 1 come from various phyla, while group 2 proteins are mainly from Proteobacteria.

For the Mn-catalase family, the 564 homologs clearly formed ten groups, and their inter-group sequence identities were typically lower than 40% (Figure 1D and Supplementary Figure 2C). The proteins from groups 1 and 2 are significantly shorter in length (group 1: 150-200 residues, group 2: 200-250 residues), and these proteins are from Firmicutes. Most proteins of group 3 and 4 have medium lengths (250-300 residues), and are mainly from Proteobacteria and Firmicutes, respectively. In contrast, the proteins of groups 5 and 6 are larger (300-350 residues) than others, and they are shown to be from Actinobacteria and Proteobacteria, respectively (Figure 1D). The proteins of group 7 differ in length and taxa distribution. The group 8 protein members are mainly from Carnobacterium and Lactobacillus, while those of group 9 are from Proteobacteria, Bacteroidetes, and Deinococcus. The proteins of group 10 vary over a broad range in length, and their host strains include Actinobacteria, Firmicutes, and Cyanobacteria.

The Taxa Distribution of Catalases

We first surveyed and compared the richness and diversity of catalases among different bacteria at the phylum level. The distribution patterns of the three catalase families varied among different phyla (Figure 2). Overall, the typical catalases are widely distributed except in Chlorobi, and the catalase-peroxidases are present in many phyla, while Mn-catalases mainly exist in Deinococcus-Thermus and Firmicutes. The phyla, such as Acidobacteria, Actinobacteria, Bacteroidetes, Planctomycetes, and Proteobacteria, are rich in typical catalases and catalase-peroxidases but lack Mn-catalases. In contrast, the phylum of Deinococcus-Thermus is rich in Mn-catalases and typical catalases, but not catalase-peroxidases. The majority of catalases in Chlorobi are catalase-peroxidases, while the typical catalases are the ones dominant in Spirochetes. Cyanobacteria and Chloroflexi have relatively fewer catalases than other phyla. The shortage of typical catalases in Cyanobacteria has been also observed previously (Zamocky et al., 2008; Johnson and Hug, 2019).

FIGURE 2

Figure 2. Average counts per genome of the three catalase families (left part) and the different groups (right part) for well-sampled phyla (upper part) and genera (lower part). Data from those phyla with less than 5 genomes and genera with less than 10 genomes are not shown.

We then further analyzed the distributions of three catalase families in 81 well-sampled genera (each genus with ten selected genomes). A notable high frequency of typical catalases was observed in soil bacteria Streptomyces, Massilia, Rhodococcus, and the pathogenic bacteria Pseudomonas (Figure 2). Indeed, all the selected genera under the phylum of Actinobacteria are high in typical catalases (mainly groups 5 and 6), and majority of them are also rich in catalase-peroxidase (mainly group 4) (Figure 2). Some genera of Actinobacteria also contain Mn-catalases, such as Amycolatopsis. The Bacillus and Paenibacillus strains are rich in various groups of typical catalases and Mn-catalases, while their close relatives, the thermophilic Geobacillus strains, mainly contain catalase-peroxidases and Mn-catalases. The Staphylococcus strains are rich in typical catalases (groups 5 and 7) but lack of catalase-peroxidases and Mn-catalases. The facultative or strict anaerobes under the Firmicutes phylum, including Enterococcus, Lactobacillus, and Clostridium, are also short of catalase-peroxidases but rich in typical catalases and Mn-catalases. The representative opportunistic pathogenic bacteria, such as Enterobacter, Klebsiella, and Proteus strains, possess all three catalase families, including the typical catalases (mainly group 5), catalase-peroxidases (group 4), and Mn-catalases (group 3). However, the richness in all the three catalase families is not observed in other pathogenic bacteria like Yersinia and Pseudomonas strains.

A Comprehensive Survey of Catalases in Streptomyces

In particular, we were interested to know the abundance and diversity of catalases in Streptomyces strains, which have complex life cycles and are the prolific producers of natural products with therapeutic applications (Bentley et al., 2002; Flärdh and Buttner, 2009; Musiol-Kroll et al., 2019). To get an overarching view of the distribution of catalases in Streptomyces, we extracted all the catalase homologs in 271 annotated Streptomyces genomes available in NCBI Genomes database, resulting in 811 typical catalase homologs and 268 catalase-peroxidase homologs, with no Mn-catalase homologs identified. In other words, there are 2.99 typical catalases and 0.99 catalase-peroxidases on average in each Streptomyces genome. Specifically, the catalase homologs in four well-characterized antibiotic-producing Streptomyces strains, including S. coelicolor A3(2), S. venezuelae ATCC 10712 (ISP5230), S. avermitilis MA-4680, and S. rimosus ATCC 10970, are listed in Supplementary Table 4. It is worth noting that, although an average of 0.99 catalase-peroxidase homolog per genome was found in Streptomyces, it was clear that some Streptomyces strains contain no catalase-peroxidase (such as S. avermitilis MA-4680 and S. rimosus ATCC 10970), whereas some other strains can have more than one homolog (like S. venezuelae ATCC 10712; Supplementary Table 4).

Remarkably, sequence analyses revealed that all the identified 268 catalase-peroxidase homologs belong to group 4 (Figure 3). For the 811 typical catalase homologs, the majority of them belong to group 6 (71%), with the remaining homologs falling into groups 5 (19%), 3 (8%), and 7 (3%). This distribution pattern was consistent with that of the ten randomly selected Streptomyces genomes (Figure 2). As shown in Figure 3, the group 5 and 6 typical catalases in Streptomyces may be further divided into two subgroups based on sequence analysis.

FIGURE 3

Figure 3. Protein sequence identities of catalase homologs from 271 Streptomyces genomes. Groups are shown in different colors. Orange: typical catalase group 3; green: typical catalase group 5; blue: typical catalase group 6; violet: typical catalase group 7; magenta: catalase-peroxidase group 4.

Functional Characterization of Five Catalases From S. rimosus ATCC 10970

The above bioinformatics analysis revealed the richness of typical catalase homologs in Streptomyces. To verify the predicted catalase function of these homologs, we selected and characterized all five typical catalase homologs (Cat1-5) identified in S. rimosus ATCC 10970, the oxytetracycline producing strain (Zhang et al., 2006). The five proteins with N-terminal His₆-tags were overproduced in E. coli BL21(DE3), and purified to homogeneity (Supplementary Figure 3). The purified proteins were reddish brown in solution, and their UV-vis spectra showed characteristic Soret bands at 405 nm (Supplementary Figure 4), supporting that they were heme-binding enzymes. All five enzymes could catalyze the decomposition of H₂O₂, as the absorbance at 240 nm of the reaction mixtures decreased continuously for all five enzymes while such an absorbance change was not observed for the reaction without enzymes (Supplementary Figure 5). These results confirmed that all five proteins are catalases. The specific activities of these catalases were determined, in which the Cat5 (group 5, Clade 2) showed relatively low activity (Table 1). Overall, the specific activities of these typical catalases were comparable to those of reported catalases (Switala and Loewen, 2002; Sooch et al., 2014).

TABLE 1

Table 1. The catalytic activities and kinetic parameters of the five typical catalases from S. rimosus ATCC 10970.

We then further analyzed the properties of these catalases by evaluating the optimal reaction pH and temperature for each of them. The Cat1, Cat3, and Cat4, all belonging to group 6, showed optimal pH at 7.4, and their activities decreased significantly when pH changed (Figure 4). Cat2 (group 3) showed high activities under alkaline pHs, with an optimal pH of 9.0. The Cat1-4 showed highest activities at 35°C, and their activities decreased significantly above 35°C. Surprisingly, Cat5 (group 5) showed nearly constant activity over a broad pH range of 5.0-10.0, and remained active at temperatures from 30 to 65°C (Figure 4). The thermostabilities of five catalases were also characterized by circular dichroism (CD) spectroscopic analysis. All five catalases showed negative peaks at 208 and 222 nm (Supplementary Figure 6), which corresponded to the α-helices of the conserved α + β catalase fold (Nicholls et al., 2000; Chelikani et al., 2004; Díaz et al., 2012). The thermal denaturation processes of five catalases were monitored by the CD signals at 222 nm (Supplementary Figure 6). Among them, Cat1 (Tm: 52°C) and Cat3 (Tm: 53°C) were less resistant to heat, while Cat4 (Tm: 75°C) displayed moderate heat tolerance. Surprisingly, although Cat2 showed complete loss of activity at 65°C, no significant change of the CD signal at 222 nm was observed up to 80°C. These results indicated that, under the above assay conditions, the tertiary structure of Cat2 was likely damaged, but most of the α-helix secondary structures remained unchanged. Cat5 displayed excellent thermostability up to 80°C, which explained why this enzyme remained catalytically active at 65°C as shown in Figure 4.

FIGURE 4

Figure 4. The relative specific activities of five catalases under different pH and temperatures. Specific activity of each catalase at 30°C and pH 7.4 was set as 100%. The catalytic activities at different pH values were determined at 30°C using several buffer systems: 50 mM NaAc-HAc, pH 4.0 – 6.0 (up triangle), 50 mM NaH₂PO₄-Na₂HPO₄, pH 6.0 – 8.0 (circle), and 50 mM Na₂CO₃-NaHCO₃, pH 8.0 - 10.0 (down triangle). Activities at different temperatures were determined at pH 7.4.

Due to the two-stage reaction mechanism, catalases did not follow Michaelis-Menten kinetics except at low H₂O₂ concentrations (Switala and Loewen, 2002; Zamocky et al., 2008). Therefore, the apparent K_m and k_cat values for the five catalase reactions were determined at low substrate concentrations (not more than 50 mM) (Table 1 and Supplementary Figure 7). Among them, Cat5 displayed lower catalytic efficiency. Overall, the kinetic values were comparable to the reported data of known typical catalases (Switala and Loewen, 2002; Jia et al., 2016, 2017).

Discussion

Catalases are widely distributed in most aerobic organisms and also in some anaerobic microorganisms (von Ossowski et al., 1993; Whittaker, 2012; Zámocký et al., 2012). An early phylogenetic analysis of 20 typical catalases suggested that small subunit typical catalases of animal and fungi were derived from one ancestor, while those catalases of plant originated from another ancestor (von Ossowski et al., 1993). The study also revealed that the typical catalases from bacteria did not cluster together (von Ossowski et al., 1993). The following studies with more catalase proteins showed that bacterial typical catalases fell into three clades (Klotz et al., 1997; Klotz and Loewen, 2003), which was further supported by later molecular evolution analyses using more catalase homologs (Zamocky et al., 2008; Zámocký et al., 2012). The Clade 1 (small subunit) typical catalases were mainly from plants and bacteria (Firmicutes and Proteobacteria), the Clade 3 (small subunit) typical catalases were from animal, fungi, archaea, and various bacteria, while the Clade 2 (large subunit) typical catalases were shown to be from bacteria and fungi (Zámocký et al., 2012). Molecular evolution analysis showed that the catalase-peroxidases fell into two clades, the (main) Clade 1 and the (minor) Clade 2, while the Mn-catalases were separated into five clades (Zámocký et al., 2012).

In bacteria, catalases are the key players in the defense against oxidative stress (Zámocký et al., 2012; Johnson and Hug, 2019; Kim et al., 2019). In addition, they also participate in many other cellular processes, such as cell development and differentiation, and production of metabolites (Cho et al., 2000; Aguirre et al., 2005; Matsuura et al., 2014; Montibus et al., 2015; Bibián et al., 2020). In this study, by mining the abundant genomic information available in the GenBank database, we were able to know the overall abundance, diversity, and distribution of catalases in bacteria. Compared to previous studies (von Ossowski et al., 1993; Klotz et al., 1997; Klotz and Loewen, 2003; Chelikani et al., 2004; Zamocky et al., 2008; Zámocký et al., 2012), our study focused on bacteria and included thousands of catalase homologs for analysis. Unlike the catalase-peroxidase family in which the proteins were closely related, both typical catalase and Mn-catalase families could be classified into distinct groups.

Based on sequence similarities, the typical catalase groups 3 and 4, group 5, and groups 6 and 7 in our study, might correspond to the previous Clade 1 (small subunit), Clade 2 (large subunit), and Clade 3 (small subunit), respectively. The catalase-peroxidase groups 3 and 4 corresponded to the previously identified Clade 1 (main clade), and the group 2 seemed to be correlated to the Clade 2 (minor clade). The Mn-catalase groups 4 and 5, group 7, and groups 8-10, appeared to correspond to the previously reported Clade 2, Clade 1, and Clade 3 of Mn-catalases, respectively, while group 3 seemed to correspond to both Clades 4 and 5 (Zámocký et al., 2012). Meanwhile, the homologs in typical catalase group 1 and 2, and Mn-catalase group 1 and 6 showed relatively low sequence similarities to known catalases in RedoxiBase, a database containing the catalases used for analyses in previous studies (Zámocký et al., 2012; Savelli et al., 2019; Supplementary Figure 8).

In particular, we paid great attention to the Streptomyces strains, which are well-known for their complicated cell cycles and capability for producing various valuable secondary metabolites (Bentley et al., 2002; Flärdh and Buttner, 2009; Musiol-Kroll et al., 2019). An extreme richness of typical catalases and catalase-peroxidases was observed in Streptomyces. On average, each Streptomyces strain contains 2.99 typical catalases and 0.99 catalase-peroxidases. But remarkably, there were no Mn-catalases found in any of the Streptomyces genomes. Sequence comparison of these catalase homologs clearly showed several subgroups of typical catalases (Figure 3), which might be the result of gene duplication. To support our bioinformatic analyses, the five typical catalase homologs identified in S. rimosus ATCC 10970 were characterized in vitro. All of them showed the catalytic activities of H₂O₂ decomposition. Among them, Cat5 (group 5, Clade 2) had the lowest specific activity but showed remarkable pH and thermal tolerances. The typical catalases (Clade 2) from Aspergillus niger and E. coli were also reported to show low activities, and in addition, the catalase from E. coli exhibited good thermal stability (Switala and Loewen, 2002). Another typical catalase (Clade 2) was isolated from Aspergillus terreus MTCC 6324, and it was extremely active and stable in broad pH range (4-12) and at temperatures up to 90°C (Vatsyayan and Goswami, 2016).

Taken together, our results set the stage for further studies of bacterial catalases to understand their roles in the adaptation of host strains to internal and external environmental changes. Furthermore, previous studies showed that typical catalases were involved in the regulation of development and differentiation of streptomycetes, and affected the production of several secondary metabolites (Cho et al., 2000; Beites et al., 2011), suggesting the manipulation of these catalases as a potential approach to modulate the secondary metabolite production in strains like streptomycetes.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

GP, KF, and ZL conceived the idea. FY, YX, LX, and SY performed the experiments. FY, KF, and GP analyzed the data and wrote the manuscript with the input of all authors. KF, FY, and HW carried out the bioinformatics analyses. All authors edited the manuscript and approved its final version.

Funding

This work was supported by funding from the Institute of Microbiology, Chinese Academy of Sciences, and Laboratory for Marine Drugs and Bioproducts of Qingdao National Laboratory for Marine Science and Technology (LMDBKF201801).

Conflict of Interest

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

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.645477/full#supplementary-material

References

Aguirre, J., Ríos-Momberg, M., Hewitt, D., and Hansberg, W. (2005). Reactive oxygen species and development in microbial eukaryotes. Trends Microbiol. 13, 111–118. doi: 10.1016/j.tim.2005.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Alfonso-Prieto, M., Biarnés, X., Vidossich, P., and Rovira, C. (2009). The molecular mechanism of the catalase reaction. J. Am. Chem. Soc. 131, 11751–11761. doi: 10.1021/ja9018572

PubMed Abstract | CrossRef Full Text | Google Scholar

Alfonso-Prieto, M., Vidossich, P., and Rovira, C. (2012). The reaction mechanisms of heme catalases: an atomistic view by ab initio molecular dynamics. Arch. Biochem. Biophys. 525, 121–130. doi: 10.1016/j.abb.2012.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Beers, R. F., and Sizer, I. W. (1952). A spectrophotometric method for measuring the breakdown of hydrogen peroxide by catalase. J. Biol. Chem. 195, 133–140.

Google Scholar

Beites, T., Pires, S. D. S., Santos, C. L., Osório, H., Moradas-Ferreira, P., and Mendes, M. V. (2011). Crosstalk between ROS homeostasis and secondary metabolism in S. natalensis ATCC 27448: modulation of pimaricin production by intracellular ROS. PLoS One 6:e27472. doi: 10.1371/journal.pone.0027472

PubMed Abstract | CrossRef Full Text | Google Scholar

Bentley, S. D., Chater, K. F., Cerdeño-Tárraga, A.-M., Challis, G. L., Thomson, N. R., James, K. D., et al. (2002). Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature 417, 141–147. doi: 10.1038/417141a

PubMed Abstract | CrossRef Full Text | Google Scholar

Bibián, M. E., Pérez-Sánchez, A., Mejía, A., and Barrios-González, J. (2020). Penicillin and cephalosporin biosyntheses are also regulated by reactive oxygen species. Appl. Microbiol. Biotechnol. 104, 1773–1783. doi: 10.1007/s00253-019-10330-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Brill, A. S., and Sandberg, H. E. (1968). Spectral studies of iron coordination in hemeprotein complexes: difference spectroscopy below 250 millimicrons. Biophys. J. 8, 669–690. doi: 10.1016/s0006-3495(68)86514-2

CrossRef Full Text | Google Scholar

Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinform. 10:421. doi: 10.1186/1471-2105-10-421

PubMed Abstract | CrossRef Full Text | Google Scholar

Chelikani, P., Fita, I., and Loewen, P. C. (2004). Diversity of structures and properties among catalases. Cell. Mol. Life Sci. 61, 192–208. doi: 10.1007/s00018-003-3206-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, Y. H., Lee, E. J., and Roe, J. H. (2000). A developmentally regulated catalase required for proper differentiation and osmoprotection of Streptomyces coelicolor. Mol. Microbiol. 35, 150–160. doi: 10.1046/j.1365-2958.2000.01685.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Díaz, A., Loewen, P. C., Fita, I., and Carpena, X. (2012). Thirty years of heme catalases structural biology. Arch. Biochem. Biophys. 525, 102–110. doi: 10.1016/j.abb.2011.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Eddy, S. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol. 7:e1002195. doi: 10.1371/journal.pcbi.1002195

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Gebali, S., Mistry, J., Bateman, A., Eddy, S. R., Luciani, A., Potter, S. C., et al. (2019). The Pfam protein families database in 2019. Nucleic Acids Res. 47, D427–D432. doi: 10.1093/nar/gky995

PubMed Abstract | CrossRef Full Text | Google Scholar

Engelmann, S., and Hecker, M. (1996). Impaired oxidative stress resistance of Bacillus subtilis sigB mutants and the role of katA and katE. FEMS Microbiol. Lett. 145, 63–69. doi: 10.1111/j.1574-6968.1996.tb08557.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, K., Pan, G., Peng, X., Zheng, J., Gao, W., Wang, J., et al. (2012). Identification of JadG as the B ring opening oxygenase catalyzing the oxidative C-C bond cleavage reaction in jadomycin biosynthesis. Chem. Biol. 19, 1381–1390. doi: 10.1016/j.chembiol.2012.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Fita, I., and Rossmann, M. G. (1985). The active center of catalase. J. Mol. Biol. 185, 21–37. doi: 10.1016/0022-2836(85)90180-9

CrossRef Full Text | Google Scholar

Flärdh, K., and Buttner, M. J. (2009). Streptomyces morphogenetics: dissecting differentiation in a filamentous bacterium. Nat. Rev. Microbiol. 7, 36–49. doi: 10.1038/nrmicro1968

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibson, D. G., Young, L., Chuang, R.-Y., Venter, J. C., Hutchison, C. A., and Smith, H. O. (2009). Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 6, 343–345. doi: 10.1038/nmeth.1318

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, X., Chen, J., Lin, C., and Lin, X. (2016). Cloning, expression, and characterization of a novel thermophilic monofunctional catalase from Geobacillus sp. CHB1. Biomed. Res. Int. 2016:7535604. doi: 10.1155/2016/7535604

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, X., Lin, X., Tian, Y., Chen, J., and You, M. (2017). High production, purification, biochemical characterization and gene analysis of a novel catalase from the thermophilic bacterium Ureibacillus thermosphaericus FZSF03. Int. J. Biol. Macromol. 103, 89–98. doi: 10.1016/j.ijbiomac.2017.05.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, L. A., and Hug, L. A. (2019). Distribution of reactive oxygen species defense mechanisms across domain bacteria. Free Radic. Biol. Med. 140, 93–102. doi: 10.1016/j.freeradbiomed.2019.03.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, S., Ueno, T., Fukuzumi, S., and Watanabe, Y. (2004). Catalase reaction by myoglobin mutants and native catalase: mechanistic investigation by kinetic isotope effect. J. Biol. Chem. 279, 52376–52381. doi: 10.1074/jbc.M403532200

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. Y., Park, C., Jang, H.-J., Kim, B.-O., Bae, H.-W., Chung, I.-Y., et al. (2019). Antibacterial strategies inspired by the oxidative stress and response networks. J. Microbiol. 57, 203–212. doi: 10.1007/s12275-019-8711-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Klotz, M. G., Klassen, G. R., and Loewen, P. C. (1997). Phylogenetic relationships among prokaryotic and eukaryotic catalases. Mol. Biol. Evol. 14, 951–958. doi: 10.1093/oxfordjournals.molbev.a025838

PubMed Abstract | CrossRef Full Text | Google Scholar

Klotz, M. G., and Loewen, P. C. (2003). The molecular evolution of catalatic hydroperoxidases: evidence for multiple lateral transfer of genes between prokaryota and from bacteria into eukaryota. Mol. Biol. Evol. 20, 1098–1112. doi: 10.1093/molbev/msg129

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsuura, H. N., Rau, M. R., and Fett-Neto, A. G. (2014). Oxidative stress and production of bioactive monoterpene indole alkaloids: biotechnological implications. Biotechnol. Lett. 36, 191–200. doi: 10.1007/s10529-013-1348-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Montibus, M., Pinson-Gadais, L., Richard-Forget, F., Barreau, C., and Ponts, N. (2015). Coupling of transcriptional response to oxidative stress and secondary metabolism regulation in filamentous fungi. Crit. Rev. Microbiol. 41, 295–308. doi: 10.3109/1040841X.2013.829416

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulvey, M. R., Switala, J., Borys, A., and Loewen, P. C. (1990). Regulation of transcription of katE and katF in Escherichia coli. J. Bacteriol. 172, 6713–6720. doi: 10.1128/jb.172.12.6713-6720.1990

PubMed Abstract | CrossRef Full Text | Google Scholar

Musiol-Kroll, E. M., Tocchetti, A., Sosio, M., and Stegmann, E. (2019). Challenges and advances in genetic manipulation of filamentous actinomycetes - the remarkable producers of specialized metabolites. Nat. Prod. Rep. 36, 1351–1369. doi: 10.1039/c9np00029a

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicholls, P. (2012). Classical catalase: ancient and modern. Arch. Biochem. Biophys. 525, 95–101. doi: 10.1016/j.abb.2012.01.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicholls, P., Fita, I., and Loewen, P. C. (2000). “Enzymology and structure of catalases. Adv. Inorganic Chem. 51, 51–106. doi: 10.1016/S0898-8838(00)51001-0

CrossRef Full Text | Google Scholar

Rice, P., Longden, I., and Bleasby, A. (2000). EMBOSS: the european molecular biology open software suite. Trends Genet. 16, 276–277.

Google Scholar

Sambrook, J., and Russell, D. W. (2001). Molecular Cloning: A Laboratory Manual, 3rd Edn. New York, NY: Cold Spring Harbor Laboratory Press.

Google Scholar

Savelli, B., Li, Q., Webber, M., Jemmat, A. M., Robitaille, A., Zamocky, M., et al. (2019). RedoxiBase: a database for ROS homeostasis regulated proteins. Redox Biol. 26:101247. doi: 10.1016/j.redox.2019.101247

PubMed Abstract | CrossRef Full Text | Google Scholar

Sooch, B. S., Kauldhar, B. S., and Puri, M. (2014). Recent insights into microbial catalases: isolation, production and purification. Biotechnol. Adv. 32, 1429–1447. doi: 10.1016/j.biotechadv.2014.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Switala, J., and Loewen, P. C. (2002). Diversity of properties among catalases. Arch. Biochem. Biophys. 401, 145–154. doi: 10.1016/S0003-9861(02)00049-8

CrossRef Full Text | Google Scholar

Vatsyayan, P., and Goswami, P. (2016). Highly active and stable large catalase isolated from a hydrocarbon degrading Aspergillus terreus MTCC 6324. Enzyme Res. 2016:4379403. doi: 10.1155/2016/4379403

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlasits, J., Jakopitsch, C., Bernroitner, M., Zamocky, M., Furtmüller, P. G., and Obinger, C. (2010). Mechanisms of catalase activity of heme peroxidases. Arch. Biochem. Biophys. 500, 74–81. doi: 10.1016/j.abb.2010.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

von Ossowski, I., Hausner, G., and Loewen, P. C. (1993). Molecular evolutionary analysis based on the amino acid sequence of catalase. J. Mol. Evol. 37, 71–76. doi: 10.1007/BF00170464

PubMed Abstract | CrossRef Full Text | Google Scholar

Whittaker, J. W. (2012). Non-heme manganese catalase–the “other” catalase. Arch. Biochem. Biophys. 525, 111–120. doi: 10.1016/j.abb.2011.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Zamocky, M., Furtmüller, P. G., and Obinger, C. (2008). Evolution of catalases from bacteria to humans. Antioxid. Redox Signal. 10, 1527–1548. doi: 10.1089/ars.2008.2046

PubMed Abstract | CrossRef Full Text | Google Scholar

Zámocký, M., Gasselhuber, B., Furtmüller, P. G., and Obinger, C. (2012). Molecular evolution of hydrogen peroxide degrading enzymes. Arch. Biochem. Biophys. 525, 131–144. doi: 10.1016/j.abb.2012.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, H.-W., Cai, Y.-J., Liao, X.-R., Zhang, F., and Zhang, D.-B. (2011). Production, characterization, cloning and sequence analysis of a monofunctional catalase from Serratia marcescens SYBC08. J. Basic Microbiol. 51, 205–214. doi: 10.1002/jobm.201000147

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Ames, B. D., Tsai, S.-C., and Tang, Y. (2006). Engineered biosynthesis of a novel amidated polyketide, using the malonamyl-specific initiation module from the oxytetracycline polyketide synthase. Appl. Environ. Microbiol. 72, 2573–2580. doi: 10.1128/AEM.72.4.2573-2580.2006

PubMed Abstract | CrossRef Full Text | Google Scholar