Computational design of Bacillus licheniformis RN-01 levansucrase for control of the chain length of levan-type fructooligosaccharides

Levansucrase (LS) from Gram-positive bacteria generally produces a large quantity of levan polymer, a polyfructose with glucose at the end (GF_n) but a small quantity of levan-type fructooligosaccharides (LFOs). The properties of levan and LFOs depend on their chain lengths, thereby determining their potential applications in food and pharmaceutical industries such as prebiotics and anti-tumor agents. Therefore, an ability to redesign and engineer the active site of levansucrase for synthesis of products with desired degree of polymerization (DP) is very beneficial. We employed computational protein design, docking and molecular dynamics to redesign and engineer the active site of Bacillus licheniformis RN-01 levansucrase for production of LFOs with DP up to five (GF₄), using two approaches: 1) blocking oligosaccharide binding track of GF₃-LS complex with large aromatic residues and 2) eliminating hydrogen bond interactions between terminal glucose of GF₄ and side chains of binding residues of GF₄-LS complex. The designed enzymes and their product patterns from these two approaches were experimentally characterized. The experimental results show that the first approach was successful in creating N251W and N251W/K372Y mutants that synthesized LFOs with DP up to five. This work illustrates how computer-aided approaches can offer novel opportunities to engineer enzymes for desired products.

Replica exchange molecular dynamics simulations reveal the structural and molecular properties of levan-type fructo-oligosaccharides of various chain lengths

Background

Levan and levan-type fructo-oligosaccharides (LFOs) have various potential applications in pharmaceutical and food industries due to their beneficial properties such as their low intrinsic viscosity and high water solubility. Previous studies showed that they exhibited prebiotic effects, anti-inflammatory and anti-tumor activities against Sarcoma-180 tumor cells of human. Despite their various potential applications, the structural and molecular properties of LFOs of various chain lengths are not well understood.

Results

We employed the replica-exchange molecular dynamics simulations method (REMD) in AMBER14 to elucidate structural and molecular properties of LFOs with chain lengths of 5 (LFO₅), 10 (LFO₁₀) and 15 (LFO₁₅) residues in two models of generalized Born implicit solvent (GB_HCT and GB_OBC1). For LFO₁₀ and LFO₁₅, four distinct conformations (helix-like, partial helix, zig-zag and random structures) were characterized by their upper-middle and lower-middle torsions. For LFO₅, two distinct conformations (partial helix and random structures) were characterized by their middle torsion and molecular angle of residues 1, 3 and 5. To determine hydrogen bonds important for the formation of helix-like structures of LFO₁₀ and LFO₁₅, occurrence frequencies of hydrogen bonds were analyzed, and the O6_(i)--H3O_(i+1) hydrogen bond was found with the highest frequency, suggesting its importance in helix formation. Among three dihedral angles between two fructosyl units [ϕ (O5’-C2’-O6-C6), ψ (C2’-O6-C6-C5) and ω (O6-C6-C5-C4)], dihedral angle distributions showed that ω was the most flexible dihedral angle and probably responsible for conformational differences of LFOs.

Conclusions

Our study provides important insights into the structural and molecular properties of LFOs, which tend to form helical structures as the chain length increases from 5 to 15 residues. This information could be beneficial for the selection of LFOs with appropriate lengths and properties for pharmaceutical and biological applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1182-7) contains supplementary material, which is available to authorized users.