Glycoside hydrolase family 32 is present in Bacillus subtilis phages

Bacteriophages are the most abundant organisms on Earth; estimates are at 10³¹ phage particles in the biosphere [1]. Moreover, they play a major role in the lateral gene transfer (LGT) [2]. For example, the genomes of some phages (cyanophages) have been shown to contain host-like genes known as auxiliary metabolic genes (AMGs) [3]. These genes are presumed to have been acquired from the virus’ host by lateral gene transfer (LGT) and may confer a selective advantage for persistence under certain environmental conditions [4]. The present report focuses on an AMG, a glycoside hydrolase 32 (GH32) family protein, observed for the first time in a viral genome and shown to have been acquired by LGT from the bacterial host of a phage.

Glycoside hydrolases (GH) 32 cleave the glycosidic bond between two monosaccharides or between a carbohydrate and an aglycone moiety [5]. Structurally, in addition to the catalytic five bladed β-propeller fold, GH32 enzymes are characterized by an additional β-sandwich in the C-terminal region [6]. The active site is composed of a WMNDPNG motif as the nucleophile and the EC motif as the acid/base catalyst [7]. The aspartate in the RDP motif may not to be directly implicated in the catalytic mechanism and presumably plays a role in substrate recognition and stabilization of the transition-state [8, 9].

Levan is a β-2,6-linked polymeric fructose that constitutes a carbohydrate reservoir in some plants, bacteria and fungi [10]. In microbes, it participates in the formation of the non-charged extracellular polysaccharide (EPS) matrix and plays a role in microbial biofilm formation [11]. Levan fructotransferase (LFTase; EC 4.2.2.16) is a member of GH32 family of enzymes. LFTase converts β-2,6-linked levan into DFA-IV (di-β-D-fructofuranose-2,6′:6,2′-dianhydride) [12]. DFAs production is of great importance because of their beneficial effects on human health [13]. When Bacillus subtilis is grown in batch cultures in sucrose-rich growth medium, it synthesizes levan [14]. In addition, Dogsa et al. [15] showed that levan, despite the fact that it is not essential for biofilm formation, could play a structural and possibly stabilizing component of B. subtilis floating biofilms. Bacillus species are found in terrestrial environments and are important industrial microorganisms used in traditional fermented foods [16, 17].

During our study of β-fructofuranosidase proteins of budworm (unpublished results), we noticed during the blast search that some bacillus phages also contain homolog to β-fructofuranosidase. The search in the CAZy database (http://www.cazy.org/) showed the presence of one Bacillus phage phiNIT1 sequence (accession AP013029.1) in the section of GH32 family. To determine if other homolog sequences exist in other viruses, we searched for the presence of GH32 enzyme homologs in the complete sequenced genomes of viruses in GenBank using BLASTp, tBLASTn and HMM profiles. Among 4602 (1449 phages) complete viral genomes available at NCBI (as of June 5^th, 2015); only three Bacillus phage (SP10, ϕNIT1 and SPG24) genomes contain GH32 homologs. The amino acid sequence of SPG24 and ϕNIT1 shows 98 % of identity and SP10 presents 76 % of identity with SPG24 and ϕNIT1. We used the amino acid sequence of Bacillus phage SP10 GH32 to search by BlastP for close homologues in bacteria, fungi, plants and animals. The sequences extracted from databases were aligned with Mafft [18] (Fig. 1 and Additional file 1: Figure S1). Analysis of sequences alignment showed that in contrast to other enzymes of the GH32 family, the three enzymes from phages did not possess a signal peptide. However, using Secretome 2.0 Server (http://www.cbs.dtu.dk/services/SecretomeP/) that predicts non-classically secreted proteins in Gram-positive and Gram-negative bacteria, we obtained high scores of secretion, 0.90, 0.86 and 0.84 for enzymes from phages SP10, ϕNIT1 and SPG24 and is in accordance with secreted enzymes. Moreover, the GH32 enzymes from the three phages possess the catalytic triad, D19, D150 and E198 (Figs. 1 and 3a). We noted that in these three phage enzymes, the last three amino acid residues in the catalytic nucleophile WMNDPNG motifs are changed to IQR residues (Fig. 1 and Additional file 1: Figure S1).

To establish the phylogenetic relationships between GH32 of phages and those of prokaryotes and eukaryotes, the sequences were aligned with Mafft (Additional file 1: Figure S1) and a phylogenetic tree was constructed using PhyML [19] and BioNJ [20]. The three phage GH32 enzymes are phylogenetically closer to enzyme GH32 of Sporolactobacillus laevolacticus than B. subtilis (Fig. 2). Interestingly, S. laevolacticus was first isolated and described as Bacillus laevolacticus by Nakayama and Yanoshi [21] and confirmed by Andersch et al. [22]. In 2006, Hatayama et al. [23] phylogenetically reclassified Bacillus laevolacticus as Sporolactobacillus laevolacticus supported by chemotaxonomic and physiological characterizations. In addition to phylogeny, among 67 complete Bacillus phage genomes, only three contain the GH32 enzymes. Therefore, we can speculate that the Bacillus phages SP10, ϕNIT1 and SPG24 acquired the GH32 gene by lateral gene transfer (LGT) from S. laevolacticus. In metazoans, enzymes of the GH32 family are present only in arthropoda that acquired this gene from bacteria by LGT [24]. Bacillus phage GH32 enzymes are also close to the Clostridium genus. The cluster (Clostridium and Sporolactobacillus) closest to GH32 phage enzymes (Fig. 2) has the canonical catalytic nucleophile WMXDI/VQR-motif instead of the classic WMNDPNG motif (Additional file 1: Figure S1).

To begin assessing the function of the GH32 phage proteins, a three-dimensional structure model was built and the carbohydrate substrate was docked into the active site. The 3D model of phage SP10 GH32 was constructed by homology using GH32 of Arthrobacter ureafaciens (PDB accession: 4FFG) as template. The modeled structure of the Bacillus phage SP10 GH32 enzyme is similar to structures of GH32 family. The enzyme structure consists of an N-terminal domain, a five bladed β-propeller, and C-terminal, β-sandwich (Fig. 3a). The active site model contains a pocket with two compartments (subsite −1 and −2) that can accommodate a disaccharide (Fig. 3b). The active site is suitable for the exo-type cleavage of disaccharide from polysaccharide [12]. For docking, the coordinates of a levantriose carbohydrate molecule were extracted from a A. ureafaciens GH32 structure (PDB: 4FFI) and docked into the 3D enzyme model from phage SP10 using the software AutoDock Vina [25]. Results of the docking simulations are in accordance with the catalytic mechanism proposed for levan fructotransferase [12]. In summary, this mechanism consists of the levan substrate binding by its nonreducing terminal fructose in the −2 subsite, and the preceding fructosyl moiety acts at −1 subsite closes at the D19 nucleophile. The β-2,6-glycosidic bond between −1 and +1 subsite (interacts with the terminal fructose moiety of levantriose) and is cleaved by a nucleophilic attack by D19 followed by the E198 acid /base catalyst (Fig. 3b and c).