Structural characterization of transglutaminase-catalyzed cross-linking between glyceraldehyde 3-phosphate dehydrogenase and polyglutamine repeats

The accumulation of abnormal polyglutamine-containing protein aggregates within the cytosol and nuclei of affected neurons is a hallmark of the progressive neurodegenerative disorders caused by an elongated (CAG)(n) repeat in the genome. The polyglutamine domains are excellent substrates for the enzyme transglutaminase type 2 (tissue), resulting in the formation of cross-links with polypeptides containing lysyl groups. Enzymatic activity toward the Q(n) domains increases greatly upon lengthening of such Q(n) stretches (n > 40). Among the possible amine donors, the glycolytic enzyme glyceraldehyde-3-phosphate-dehydrogenase was shown to tightly bind several proteins involved in polyglutamine expansion diseases. Recently, the authors have shown that K191, K268, and K331, out of the 26 lysines present in glyceraldehyde-3-phosphate-dehydrogenase, are the reactive amine-donor sites forming cross-links with substance P, which bears the simplest Q(n) domain (n = 2). The present study reports that synthetic peptides of both pathological and nonpathological length (n = 43 and 17, respectively) form cross-links with the same K residues located in the C-terminal region of glyceraldehyde-3-phosphate-dehydrogenase. In addition, it is shown that extra K residues present in the C termini of glyceraldehyde-3-phosphate-dehydrogenase are susceptible to cross-linking in the presence of transglutaminase. The present results indicate a possible modulating effect of Q(n) stretches on tissue transglutaminase substrate specificity and mechanism of recognition.     

Involvement of ATP synthase residues alphaArg-376, betaArg-182, and betaLys-155 in Pi binding

alphaArg-376, betaLys-155, and betaArg-182 are catalytically important ATP synthase residues that were proposed to be involved in substrate Pi binding and subsequent steps of ATP synthesis [Senior, A.E., Nadanaciva, S. and Weber, J. (2002) Biochim. Biophys. Acta 1553, 188-211]. Here, it was shown using purified Escherichia coli F(1)-ATPase that whereas Pi protected wild-type from reaction with 7-chloro-4-nitrobenzo-2-oxa-1,3-diazole, mutations betaK155Q, betaR182Q, betaR182K, and alphaR376Q abolished protection. Therefore, in ATP synthesis initial binding of substrate Pi in open catalytic site betaE is supported by each of these three residues.     

ubiI,a New Gene in Escherichia coli Coenzyme Q Biosynthesis,Is Involved in Aerobic C5-hydroxylation

Coenzyme Q (ubiquinone or Q) is a redox-active lipid found in organisms ranging from bacteria to mammals in which it plays a crucial role in energy-generating processes. Q biosynthesis is a complex pathway that involves multiple proteins. In this work, we show that the uncharacterized conserved visC gene is involved in Q biosynthesis in Escherichia coli, and we have renamed it ubiI. Based on genetic and biochemical experiments, we establish that the UbiI protein functions in the C5-hydroxylation reaction. A strain deficient in ubiI has a low level of Q and accumulates a compound derived from the Q biosynthetic pathway, which we purified and characterized. We also demonstrate that UbiI is only implicated in aerobic Q biosynthesis and that an alternative enzyme catalyzes the C5-hydroxylation reaction in the absence of oxygen. We have solved the crystal structure of a truncated form of UbiI. This structure shares many features with the canonical FAD-dependent para-hydroxybenzoate hydroxylase and represents the first structural characterization of a monooxygenase involved in Q biosynthesis. Site-directed mutagenesis confirms that residues of the flavin binding pocket of UbiI are important for activity. With our identification of UbiI, the three monooxygenases necessary for aerobic Q biosynthesis in E. coli are known.     

The C-terminal end of the Trypanosoma brucei editing deaminase plays a critical role in tRNA binding

Adenosine to inosine editing at the wobble position allows decoding of multiple codons by a single tRNA. This reaction is catalyzed by adenosine deaminases acting on tRNA (ADATs) and is essential for viability. In bacteria, the anticodon-specific enzyme is a homodimer that recognizes a single tRNA substrate (tRNA(Arg)(ACG)) and can efficiently deaminate short anticodon stem-loop mimics of this tRNA in vitro. The eukaryal enzyme is composed of two nonidentical subunits, ADAT2 and ADAT3, which upon heterodimerization, recognize seven to eight different tRNAs as substrates, depending on the organism, and require a full-length tRNA for activity. Although crystallographic data have provided clues to why the bacterial deaminase can utilize short substrates, residues that provide substrate binding and recognition with the eukaryotic enzymes are not currently known. In the present study, we have used a combination of mutagenesis, binding studies, and kinetic analysis to explore the contribution of individual residues in Trypanosoma brucei ADAT2 (TbADAT2) to tRNA recognition. We show that deletion of the last 10 amino acids at the C terminus of TbADAT2 abolishes tRNA binding. In addition, single alanine replacements of a string of positively charged amino acids (KRKRK) lead to binding defects that correlate with losses in enzyme activity. This region, which we have termed the KR-domain, provides a first glance at key residues involved in tRNA binding by eukaryotic tRNA editing deaminases.     

C-terminal loop of Streptomyces phospholipase D has multiple functional roles

We have recently shown that two flexible loops of Streptomyces phospholipase D (PLD) affect the catalytic reaction of the enzyme by a comparative study of chimeric PLDs. Gly188 and Asp191 of PLD from Streptomyces septatus TH-2 (TH-2PLD) were identified as the key amino acid residues involved in the recognition of phospholipids. In the present study, we further investigated the relationship between a C-terminal loop of TH-2PLD and PLD activities to elucidate the reaction mechanism and the recognition of the substrate. By analyzing chimeras and mutants in terms of hydrolytic and transphosphatidylation activities, Ala426 and Lys438 of TH-2PLD were identified as the residues associated with the activities. We found that Gly188 and Asp191 recognized substrate forms, whereas residues Ala426 and Lys438 enhanced transphosphatidylation and hydrolysis activities regardless of the substrate form. By substituting Ala426 and Lys438 with Phe and His, respectively, the mutant showed not only higher activities but also higher thermostability and tolerance against organic solvents. Furthermore, the mutant also improved the selectivity of the transphosphatidylation activity. The residues Ala426 and Lys438 were located in the C-terminal flexible loop of Streptomyces PLD separate from the highly conserved catalytic HxKxxxxD motifs. We demonstrated that this C-terminal loop, which formed the entrance of the active well, has multiple functional roles in Streptomyces PLD.     

Role of Ser216 in the mechanism of action of membrane-bound lytic transglycosylase B: further evidence for substrate-assisted catalysis

Lytic transglycosylases cleave the beta-(1-->4)-glycosidic bond in the bacterial cell wall heteropolymer peptidoglycan between the N-acetylmuramic acid (MurNAc) and N-acetylglucosamine (GlcNAc) residues with the concomitant formation of a 1,6-anhydromuramoyl residue. Based on sequence alignments, Ser216 in Pseudomonas aeruginosa membrane-bound lytic transglycosylase B (MltB) was targeted for replacement with alanine to delineate its role in the enzyme's mechanism of action. The specific activity of the Ser216-->Ala MltB derivative was less than 12% of that for the wild-type enzyme, while its substrate binding affinity remained virtually unaltered. These data are in agreement with a role of Ser216 in orienting the N-acetyl group on MurNAc at the -1 subsite of MltB for its participation in a substrate-assisted mechanism of action.     

Involvement of Ser94 in RNase HIII from Chlamydophila pneumoniae in the recognition of a single ribonucleotide misincorporated into double-stranded DNA

We recently provided the first report that RNase HIII can cleave a DNA-rN₁-DNA/DNA substrate (rN₁, one ribonucleotide) in vitro. In the present study, mutagenesis analyses and molecular dynamics (MD) simulations were performed on RNase HIII from Chlamydophila pneumoniae AR39 (CpRNase HIII). Our results elucidate the mechanism of ribonucleotide recognition employed by CpRNase HIII, indicating that the G95/K96/G97 motif of CpRNase HIII represents the main surface interacting with single ribonucleotides, in a manner similar to that of the GR(K)G motif of RNase HIIs. However, CpRNase HIII lacks the specific tyrosine required for RNase HII to recognize single ribonucleotides in double-stranded DNA (dsDNA). Interestingly, MD shows that Ser94 of CpRNase HIII forms a stable hydrogen bond with the deoxyribonucleotide at the (5')RNA–DNA(3') junction, moving this nucleotide away from the chimeric ribonucleotide. This movement appears to deform the nucleic acid backbone at the RNA–DNA junction and allows the ribonucleotide to interact with the GKG motif. Based on the inferences drawn from MD simulations, biochemical results indicated that Ser94 was necessary for catalytic activity on the DNA-rN₁-DNA/DNA substrate; mutant S94V could bind this substrate but exhibited no cleavage. Mismatches opposite the single ribonucleotide misincorporated in dsDNA inhibited cleavage by CpRNase HIII to varying degrees but did not interfere with CpRNase/substrate binding. Further MD results implied that mismatches impair the interaction between Ser94 and the deoxyribonucleotide at the RNA–DNA junction. Consequently, recognition of the misincorporated ribonucleotide was disturbed. Our results may help elucidate the distinct substrate-recognition properties of different RNase Hs.     

Inversion of the allosteric response of Escherichia coli glucosamine-6-P deaminase to N-acetylglucosamine 6-P, by single amino acid replacements

Amino acid replacements in the active site of glucosamine-6-P deaminase from Escherichia coli (GlcN6P deaminase, EC 3.5.99.6) involving the residues D141 and E148 produce atypical allosteric kinetics. These residues are located in the chain segment 139-156 which is part of the active site and which also forms several intersubunit contacts close to the allosteric site. In the D141N and E148Q mutant forms of this deaminase, there is an inversion of the effect of its physiological allosteric effector, N-acetylglucosamine 6-P, which becomes an inhibitor at substrate concentrations above a critical value. For both mutants, this particular point appears at low substrate concentration and the inhibition by the allosteric activator is the dominant effect in velocity versus substrate curves. These effects are analyzed as a particular case of the concerted allosteric model, assuming that the R state, the conformer displaying the higher affinity for the substrate, is the less catalytic state, thus producing an inverted allosteric response.     

Three subunits contribute amino acids to the active site of tetrameric adenylosuccinate lyase: Lys268 and Glu275 are required

Tetrameric adenylosuccinate lyase (ASL) of Bacillus subtilis catalyzes the cleavage of adenylosuccinate to form AMP and fumarate. We previously reported that two distinct subunits contribute residues to each active site, including the His68 and His89 from one and His141 from a second subunit [Brosius, J. L., and Colman, R. F. (2000) Biochemistry 39, 13336-13343]. Glu(275) is 2.8 A from His141 in the ASL crystal structure, and Lys268 is also in the active site region; Glu275 and Lys268 come from a third, distinct subunit. Using site-directed mutagenesis, we have replaced Lys268 by Arg, Gln, Glu, and Ala, with specific activities of the purified mutant enzymes being 0.055, 0.00069, 0.00028, and 0.0, respectively, compared to 1.56 units/mg for wild-type (WT) enzyme. Glu275 was substituted by Gln, Asp, Ala, and Arg; none of these homogeneous mutant enzymes has detectable activity. Circular dichroism and light scattering reveal that neither the secondary structure nor the oligomeric state of the Lys268 mutant enzymes has been perturbed. Native gel electrophoresis and circular dichroism indicate that the Glu275 mutant enzymes are tetramers, but their conformation is altered slightly. For K268R, the K(m)s for all substrates are similar to WT enzyme. Binding studies using [2-3H]-adenylosuccinate reveal that none of the Glu275 mutant enzymes, nor inactive K268A, can bind substrate. We propose that Lys268 participates in binding substrate and that Glu275 is essential for catalysis because of its interaction with His141. Incubation of H89Q with K268Q or E275Q leads to restoration of up to 16% WT activity, while incubation of H141Q with K268Q or E275Q results in 6% WT activity. These complementation studies provide the first functional evidence that a third subunit contributes residues to each intersubunit active site of ASL. Thus, adenylosuccinate lyase has four active sites per enzyme tetramer, each of which is formed from regions of three subunits.     

A more robust version of the Arginine 210-switched mutant in subunit a of the Escherichia coli ATP synthase

Previous work has shown that the essential R210 of subunit a in the Escherichia coli ATP synthase can be switched with a conserved glutamine Q252 with retention of a moderate level of function, that a third mutation P204T enhances this function, and that the arginine Q252R can be replaced by lysine without total loss of activity. In this study, the roles of P204T and R210Q were examined. It was concluded that the threonine in P204T is not directly involved in function since its replacement by alanine did not significantly affect growth properties. Similarly, it was concluded that the glutamine in R210Q is not directly involved with function since replacement by glycine results in significantly enhanced function. Not only did the rate of ATP-driven proton translocation increase, but also the sensitivity of ATP hydrolysis to inhibition by N,N′-dicyclohexylcarbodiimide (DCCD) rose to more than 50%. Finally, mutations at position E219, a residue near the proton pathway, were used to test whether the Arginine-switched mutant uses the normal proton pathway. In a wild type background, the E219K mutant was confirmed to have greater function than the E219Q mutant, as has been shown previously. This same unusual result was observed in the triple mutant background, P204T/R210Q/Q252R, suggesting that the Arginine-switched mutants are using the normal proton pathway from the periplasm.     

Structure and molecular dynamics simulation of archaeal prefoldin: the molecular mechanism for binding and recognition of nonnative substrate proteins

Prefoldin (PFD) is a heterohexameric molecular chaperone complex in the eukaryotic cytosol and archaea with a jellyfish-like structure containing six long coiled-coil tentacles. PFDs capture protein folding intermediates or unfolded polypeptides and transfer them to group II chaperonins for facilitated folding. Although detailed studies on the mechanisms for interaction with unfolded proteins or cooperation with chaperonins of archaeal PFD have been performed, it is still unclear how PFD captures the unfolded protein. In this study, we determined the X-ray structure of Pyrococcus horikoshii OT3 PFD (PhPFD) at 3.0 Å resolution and examined the molecular mechanism for binding and recognition of nonnative substrate proteins by molecular dynamics (MD) simulation and mutation analyses. PhPFD has a jellyfish-like structure with six long coiled-coil tentacles and a large central cavity. Each subunit has a hydrophobic groove at the distal region where an unfolded substrate protein is bound. During MD simulation at 330 K, each coiled coil was highly flexible, enabling it to widen its central cavity and capture various nonnative proteins. Docking MD simulation of PhPFD with unfolded insulin showed that the β subunit is essentially involved in substrate binding and that the α subunit modulates the shape and width of the central cavity. Analyses of mutant PhPFDs with amino acid replacement of the hydrophobic residues of the β subunit in the hydrophobic groove have shown that βIle107 has a critical role in forming the hydrophobic groove.     

The role of double covalent flavin binding in chito-oligosaccharide oxidase from Fusarium graminearum

ChitO (chito-oligosaccharide oxidase) from Fusarium graminearum catalyses the regioselective oxidation of N-acetylated oligosaccharides. The enzyme harbours an FAD cofactor that is covalently attached to His94 and Cys154. The functional role of this unusual bi-covalent flavin-protein linkage was studied by site-directed mutagenesis. The double mutant (H94A/C154A) was not expressed, which suggests that a covalent flavin-protein bond is needed for protein stability. The single mutants H94A and C154A were expressed as FAD-containing enzymes in which one of the covalent FAD-protein bonds was disrupted relative to the wild-type enzyme. Both mutants were poorly active, as the k(cat) decreased (8.3- and 3-fold respectively) and the K(m) increased drastically (34- and 75-fold respectively) when using GlcNac as the substrate. Pre-steady-state analysis revealed that the rate of reduction in the mutant enzymes is decreased by 3 orders of magnitude when compared with wild-type ChitO (k(red)=750 s(-1)) and thereby limits the turnover rate. Spectroelectrochemical titrations revealed that wild-type ChitO exhibits a relatively high redox potential (+131 mV) and the C154A mutant displays a lower potential (+70 mV), while the H94A mutant displays a relatively high potential of approximately +164 mV. The results show that a high redox potential is not the only prerequisite to ensure efficient catalysis and that removal of either of the covalent bonds may perturb the geometry of the Michaelis complex. Besides tuning the redox properties, the bi-covalent binding of the FAD cofactor in ChitO is essential for a catalytically competent conformation of the active site.     

Helix 8 and helix 10 are involved in substrate recognition in the rat monocarboxylate transporter MCT1.

Transport of lactate, pyruvate, and the ketone bodies, acetoacetate and beta-hydroxybutyrate, is mediated in many mammalian cells by the monocarboxylate transporter MCT1. To be accepted as a substrate, a carboxyl group and an unpolar side chain are necessary. Site-directed mutagenesis of the rat MCT1 was used to identify residues which are involved in substrate recognition. Helices 8 and 10 but not helix 9 were found to contain critical residues for substrate recognition. Mutation of arginine 306 to threonine in helix 8 resulted in strongly reduced transport activity. Concomitantly, saturable transport was transformed into a nonsaturable dependence of transport activity on lactate concentration, suggesting that binding of the substrate was strongly impaired. Furthermore, proton translocation in the mutant was partially uncoupled from monocarboxylate transport. Mutation of phenylalanine 360 to cysteine in helix 10 resulted in an altered substrate side chain recognition. In contrast to the wild-type transporter, monocarboxylates with more bulky and polar side chains were recognized by the mutated MCT1. Mutation of selected residues in helix 9 and helix 11 (C336A, H337Q, and E391Q) did not cause alterations of the transport properties of MCT1. It is suggested that substrate binding occurs in the carboxy-terminal half of MCT1 and that helices 8 and 10 are involved in the recognition of different parts of the substrate.     

A novel O-linked glycan modulates Campylobacter jejuni major outer membrane protein-mediated adhesion to human histo-blood group antigens and chicken colonization

Campylobacter jejuni is an important cause of human foodborne gastroenteritis; strategies to prevent infection are hampered by a poor understanding of the complex interactions between host and pathogen. Previous work showed that C. jejuni could bind human histo-blood group antigens (BgAgs) in vitro and that BgAgs could inhibit the binding of C. jejuni to human intestinal mucosa ex vivo. Here, the major flagella subunit protein (FlaA) and the major outer membrane protein (MOMP) were identified as BgAg-binding adhesins in C. jejuni NCTC11168. Significantly, the MOMP was shown to be O-glycosylated at Thr²⁶⁸; previously only flagellin proteins were known to be O-glycosylated in C. jejuni. Substitution of MOMP Thr²⁶⁸ led to significantly reduced binding to BgAgs. The O-glycan moiety was characterized as Gal(β1–3)-GalNAc(β1–4)-GalNAc(β1–4)-GalNAcα1-Thr²⁶⁸; modelling suggested that O-glycosylation has a notable effect on the conformation of MOMP and this modulates BgAg-binding capacity. Glycosylation of MOMP at Thr²⁶⁸ promoted cell-to-cell binding, biofilm formation and adhesion to Caco-2 cells, and was required for the optimal colonization of chickens by C. jejuni, confirming the significance of this O-glycosylation in pathogenesis.     

Substrate Specificity and Oligomerization of Human GMP Synthetase

Guanine monophosphate (GMP) synthetase is a bifunctional two-domain enzyme. The N-terminal glutaminase domain generates ammonia from glutamine and the C-terminal synthetase domain aminates xanthine monophosphate (XMP) to form GMP. Mammalian GMP synthetases (GMPSs) contain a 130-residue-long insert in the synthetase domain in comparison to bacterial proteins. We report here the structure of a eukaryotic GMPS. Substrate XMP was bound in the crystal structure of the human GMPS enzyme. XMP is bound to the synthetase domain and covered by a LID motif. The enzyme forms a dimer in the crystal structure with subunit orientations entirely different from the bacterial counterparts. The inserted sub-domain is shown to be involved in substrate binding and dimerization. Furthermore, the structural basis for XMP recognition is revealed as well as a potential allosteric site. Enzymes in the nucleotide metabolism typically display an increased activity in proliferating cells due to the increased need for nucleotides. Many drugs used as immunosuppressants and for treatment of cancer and viral diseases are indeed nucleobase- and nucleoside-based compounds, which are acting on or are activated by enzymes in this pathway. The information obtained from the crystal structure of human GMPS might therefore aid in understanding interactions of nucleoside-based drugs with GMPS and in structure-based design of GMPS-specific inhibitors.     

Structure of tRNA dimethylallyltransferase: RNA modification through a channel

Dimethylallyltransferase (DMATase) transfers a five-carbon isoprenoid moiety from dimethylallyl pyrophosphate (DMAPP) to the amino group of adenosine at position 37 of certain tRNAs. Reported here are the crystal structures of Pseudomonas aeruginosa DMATase alone and in complex with pyrophosphate at 1.9 A resolution. Surprisingly, the enzyme possesses a central channel spanning the entire width of the enzyme. Both the accepting substrate tRNA and the donating substrate DMAPP appear to enter the channel from opposite sides in an ordered sequence, with tRNA first and DMAPP second, and the RNA modification reaction occurs in the middle of the channel once the two substrates have met. The structure of DMATase is homologous to a class of small soluble kinases involved in biosynthesis of nucleotide precursors for nucleic acids, indicating its possibly evolutionary origin. Furthermore, specific recognition of the pyrophosphate by a conserved loop in DMATase, similar to the P-loop commonly seen in diverse nucleotide-binding proteins, demonstrates that DMATase is structurally and mechanistically distinct from farnesyltransferase, another family of prenyltransferases involved in protein modification.     

Insights into the Catalytic Mechanism of 16S rRNA Methyltransferase RsmE (mU1498) from Crystal and Solution Structures

RsmE is the founding member of a new RNA methyltransferase (MTase) family responsible for methylation of U1498 in 16S ribosomal RNA in Escherichia coli. It is well conserved across bacteria and plants and may play an important role in ribosomal intersubunit communication. The crystal structure in monomer showed that it consists of two distinct but structurally related domains: the PUA (pseudouridine synthases and archaeosine‐specific transglycosylases)-like RNA recognition and binding domain and the conserved MTase domain with a deep trefoil knot. Analysis of small-angle X-ray scattering data revealed that RsmE forms a flexible dimeric conformation that may be essential for substrate binding. The S‐adenosyl‐l‐methionine (AdoMet)-binding characteristic determined by isothermal titration calorimetry suggested that there is only one AdoMet molecule bound in the subunit of the homodimer. In vitro methylation assay of the mutants based on the RsmE-AdoMet-uridylic acid complex model showed key residues involved in substrate binding and catalysis. Comprehensive comparisons of RsmE with closely related MTases, combined with the biochemical experiments, indicated that the MTase domain of one subunit in dimeric RsmE is responsible for binding of one AdoMet molecule and catalytic process while the PUA-like domain in the other subunit is mainly responsible for recognition of one substrate molecule (the ribosomal RNA fragment and ribosomal protein complex). The methylation process is required by collaboration of both subunits, and dimerization is functionally critical for catalysis. In general, our study provides new information on the structure-function relationship of RsmE and thereby suggests a novel catalytic mechanism.     

The nucleobase-ascorbate transporter (NAT) signature motif in UapA defines the function of the purine translocation pathway

UapA, a member of the NAT/NCS2 family, is a high affinity, high capacity, uric acid-xanthine/H+ symporter of Aspergillus nidulans. We have previously presented evidence showing that a highly conserved signature motif ([Q/E/P]408-N-X-G-X-X-X-X-T-[R/K/G])417 is involved in UapA function. Here, we present a systematic mutational analysis of conserved residues in or close to the signature motif of UapA. We show that even the most conservative substitutions of residues Q408, N409 and G411 modify the kinetics and specificity of UapA, without affecting targeting in the plasma membrane. Q408 substitutions show that this residue determines both substrate binding and transport catalysis, possibly via interactions with position N9 of the imidazole ring of purines. Residue N409 is an irreplaceable residue necessary for transport catalysis, but is not involved in substrate binding. Residue G411 determines, indirectly, both the kinetics (K(m), V) and specificity of UapA, probably due to its particular property to confer local flexibility in the binding site of UapA. In silico predictions and a search in structural databases strongly suggest that the first part of the NAT signature motif of UapA (Q(408)NNG(411)) should form a loop, the structure of which is mostly affected by mutations in G411. Finally, substitutions of residues T416 and R417, despite being much better tolerated, can also affect the kinetics or the specificity of UapA. Our results show that the NAT signature motif defines the function of the UapA purine translocation pathway and strongly suggest that this might occur by determining the interactions of UapA with the imidazole part of purines.     

A network-based approach to identify substrate classes of bacterial glycosyltransferases

Background

Bacterial interactions with the environment- and/or host largely depend on the bacterial glycome. The specificities of a bacterial glycome are largely determined by glycosyltransferases (GTs), the enzymes involved in transferring sugar moieties from an activated donor to a specific substrate. Of these GTs their coding regions, but mainly also their substrate specificity are still largely unannotated as most sequence-based annotation flows suffer from the lack of characterized sequence motifs that can aid in the prediction of the substrate specificity.

Results

In this work, we developed an analysis flow that uses sequence-based strategies to predict novel GTs, but also exploits a network-based approach to infer the putative substrate classes of these predicted GTs. Our analysis flow was benchmarked with the well-documented GT-repertoire of Campylobacter jejuni NCTC 11168 and applied to the probiotic model Lactobacillus rhamnosus GG to expand our insights in the glycosylation potential of this bacterium. In L. rhamnosus GG we could predict 48 GTs of which eight were not previously reported. For at least 20 of these GTs a substrate relation was inferred.

Conclusions

We confirmed through experimental validation our prediction of WelI acting upstream of WelE in the biosynthesis of exopolysaccharides. We further hypothesize to have identified in L. rhamnosus GG the yet undiscovered genes involved in the biosynthesis of glucose-rich glycans and novel GTs involved in the glycosylation of proteins. Interestingly, we also predict GTs with well-known functions in peptidoglycan synthesis to also play a role in protein glycosylation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-349) contains supplementary material, which is available to authorized users.