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

A review of fructosyl-transferases from catalytic characteristics and structural features to reaction mechanisms and product specificity

Carbohydrate-active enzymes are accountable for the synthesis and degradation of glycosidic bonds among diverse carbohydrates. Fructosyl-transferases represent a subclass of these enzymes, employing sucrose as a substrate to generate fructooligosaccharides (FOS) and fructan polymers. This category primarily includes levansucrase (LS, EC 2.4.1.10), inulosucrase (IS, EC 2.4.1.9), and β-fructofuranosidase (Ffase, EC 3.2.1.26). These three enzymes possess a similar five-bladed β-propeller fold and employ an anomer-retaining reaction mechanism mediated by nucleophiles, transition state stabilizers, and general acids/bases. However, they exhibit distinct product profiles, characterized by variations in linkage specificity and molecular mass distribution. Consequently, this article comprehensively explores recent advancements in the catalytic characteristics, structural features, reaction mechanisms, and product specificity of levansucrase, inulosucrase, and β-fructofuranosidase (abbreviated as LS, IS, and Ffase, respectively). Furthermore, it discusses the potential for modifying catalytic properties and product specificity through structure-based design, which enables the rational production of custom fructan and FOS.

Production, effects, and applications of fructans with various molecular weights

Fructan, a widespread functional polysaccharide, has been used in the food, pharmaceutical, cosmetic, and material production fields because of its versatile physicochemical properties and biological activities. Inulin from plants and levan from microorganisms are two of the most extensively studied fructans. Fructans from different plants or microorganisms have inconsistent molecular weights, and the molecular weight of fructan affects its properties, functions, and applications. Recently, increasing attention has been paid to the production and application of fructans having various molecular weights, and biotechnological processes have been explored to produce tailor-made fructans from sucrose. This review encompasses the introduction of extraction, enzymatic transformation, and fermentation production processes for fructans with diverse molecular weights. Notably, it highlights the enzymes involved in fructan biosynthesis and underscores their physiological effects, with a special emphasis on their prebiotic properties. Moreover, the applications of fructans with varying molecular weights are also emphasized.

Unraveling the role of flexible coil near calcium binding site of levansucrase on thermostability and product profile via proline substitution and molecular dynamics simulations

Due to its bioactivity and versatile applications, levan has appeared as a promising biomaterial. Levansucrase is responsible for the conversion of sucrose into levan. With the goal of enhancing levan production, the strategy for enhancing the stability of levansucrase is being intensively studied. To make proteins more stable under high temperatures, proline, the most rigid residue, can be introduced into previously flexible regions. Herein, G249, D250, N251, and H252 on the flexible coil close to the calcium binding site of Bacillus licheniformis levansucrase were replaced with proline. Mutations at G249P greatly enhance both the enzyme's thermodynamic and kinetic stability, while those at H252P improve solely the enzyme's kinetic stability. GPC analysis revealed that G249P synthesize more levan, but H252P generate primarily oligosaccharides. Molecular dynamics simulations (MD) and MM/GBSA analysis revealed that G249P mutation increased not only the stability of levansucrase, but also affinity toward fructan.

Computational redesign of Beta-27 Fab with substantially better predicted binding affinity to the SARS-CoV-2 Omicron variant than human ACE2 receptor

During the COVID-19 pandemic, SARS-CoV-2 has caused large numbers of morbidity and mortality, and the Omicron variant (B.1.1.529) was an important variant of concern. To enter human cells, the receptor-binding domain (RBD) of the S1 subunit of SARS-CoV-2 (SARS-CoV-2-RBD) binds to the peptidase domain (PD) of Angiotensin-converting enzyme 2 (ACE2) receptor. Disrupting the binding interactions between SARS-CoV-2-RBD and ACE2-PD using neutralizing antibodies is an effective COVID-19 therapeutic solution. Previous study found that Beta-27 Fab, which was obtained by digesting the full IgG antibodies that were isolated from a patient infected with SARS-CoV-2 Beta variant, can neutralize Victoria, Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), and Delta (B.1.617.2) variants. This study employed computational protein design and molecular dynamics (MD) to investigate and enhance the binding affinity of Beta-27 Fab to SARS-CoV-2-RBD Omicron variant. MD results show that five best designed Beta-27 Fabs (Beta-27-D01 Fab, Beta-27-D03 Fab, Beta-27-D06 Fab, Beta-27-D09 Fab and Beta-27-D10 Fab) were predicted to bind to Omicron RBD in the area, where ACE2 binds, with significantly better binding affinities than Beta-27 Fab and ACE2. Their enhanced binding affinities are mostly caused by increased binding interactions of CDR L2 and L3. They are promising candidates that could potentially be employed to disrupt the binding between ACE2 and Omicron RBD.

Successful Manipulation of the Product Spectrum of the Erwinia amylovora Levansucrase by Modifying the Residues around loop1, Loop 3, and Loop 4

Levansucrase (LS, EC 2.4.1.10) catalyzes the synthesis of levan by successively transferring the fructosyl moiety from sucrose to an elongated fructan chain. Although the product distribution of LS from Erwinia amylovora (Ea-LS) was studied under different sucrose concentrations, the effect of residues on the product formation is yet unknown. The first levanhexaose-complexed structure of LS from Bacillus subtilis (Bs-SacB) provided information on the oligosaccharide binding sites (OB sites), from +1 to +4 subsites. Since Ea-LS would efficiently produce fructooligosaccharides, a substitution mutation of OB sites in Bs-SacB and the corresponding residues of Ea-LS were conducted to investigate how these mutants would influence the product distribution. As a result, a series of mutants with different product spectrum were obtained. Notably, the mutants of G98E, V151F, and N200T around loop 1, loop 3, and loop 4 all showed a significant increase in both the molecular mass and the yield of high-molecular-mass levan, suggesting that the product profile of Ea-LS was significantly modified.

Analysis of carbohydrates and glycoconjugates by matrix-assisted laser desorption/ionization mass spectrometry: An update for 2019-2020

This review is the tenth update of the original article published in 1999 on the application of matrix-assisted laser desorption/ionization (MALDI) mass spectrometry to the analysis of carbohydrates and glycoconjugates and brings coverage of the literature to the end of 2020. Also included are papers that describe methods appropriate to analysis by MALDI, such as sample preparation techniques, even though the ionization method is not MALDI. The review is basically divided into three sections: (1) general aspects such as theory of the MALDI process, matrices, derivatization, MALDI imaging, fragmentation, quantification and the use of arrays. (2) Applications to various structural types such as oligo- and polysaccharides, glycoproteins, glycolipids, glycosides and biopharmaceuticals, and (3) other areas such as medicine, industrial processes and glycan synthesis where MALDI is extensively used. Much of the material relating to applications is presented in tabular form. The reported work shows increasing use of incorporation of new techniques such as ion mobility and the enormous impact that MALDI imaging is having. MALDI, although invented nearly 40 years ago is still an ideal technique for carbohydrate analysis and advancements in the technique and range of applications show little sign of diminishing.

Biochemical and ligand binding properties of recombinant Xenopus laevis cortical granule lectin-1

Intelectins are putative innate immune lectins that are found throughout chordates. The first intelectin reported was Xenopus laevis cortical granule lectin-1 (XCGL-1 or XL-35). XCGL-1 is critical in fertilization membrane development in Xenopus. Here, we explored the biochemical properties of XCGL-1. The cysteines responsible for forming intermolecular disulfide bonds were identified. XCGL-1 adopted a four-lobed structure as observed by electron microscopy. The full-length XCGL-1 and the carbohydrate recognition domain (CRD) bind galactose-containing carbohydrates at nanomolar to micromolar affinities. Molecular modeling suggested that galactoside ligands coordinated the binding site calcium ion and interacted with residues around the groove made available by the non-conserved substitution compared to human intelectin-1. Folding conditions for production of recombinant XCGL-1 CRD were also investigated. Our results not only provide new biochemical insights into the function of XCGL-1, but may also provide foundation for further applications of XCGL-1 as glycobiology tools.

Computational design of Lactobacillus Acidophilus α-L-rhamnosidase to increase its structural stability

α-L-rhamnosidase catalyzes hydrolysis of the terminal α-L-rhamnose from various natural rhamnoglycosides, including naringin and hesperidin, and has various applications such as debittering of citrus juices in the food industry and flavonoid derhamnosylation in the pharmaceutical industry. However, its activity is lost at high temperatures, limiting its usage. To improve Lactobacillus acidophilus α-L-rhamnosidase stability, we employed molecular dynamics (MD) to identify a highly flexible region, as evaluated by its root mean square fluctuation (RMSF) value, and computational protein design (Rosetta) to increase rigidity and favorable interactions of residues in highly flexible regions. MD results show that five regions have the highest flexibilities and were selected for design by Rosetta. Twenty-one designed mutants with the best ΔΔG at each position and ΔΔG < 0 REU were simulated at high temperature. Eight designed mutants with ΔRMSF of highly flexible regions lower than -10.0% were further simulated at the optimum temperature of the wild type. N88Q, N202V, G207D, Q209M, N211T and Y213K mutants were predicted to be more stable and could maintain their native structures better than the wild type due to increased hydrogen bond interactions of designed residues and their neighboring residues. These designed mutants are promising enzymes with high potential for stability improvement.

Computational redesign of Fab CC12.3 with substantially better predicted binding affinity to SARS-CoV-2 than human ACE2 receptor

SARS-CoV-2 is responsible for COVID-19 pandemic, causing large numbers of cases and deaths. It initiates entry into human cells by binding to the peptidase domain of angiotensin-converting enzyme 2 (ACE2) receptor via its receptor binding domain of S1 subunit of spike protein (SARS-CoV-2-RBD). Employing neutralizing antibodies to prevent binding between SARS-CoV-2-RBD and ACE2 is an effective COVID-19 therapeutic solution. Previous studies found that CC12.3 is a highly potent neutralizing antibody that was isolated from a SARS-CoV-2 infected patient, and its Fab fragment (Fab CC12.3) bound to SARS-CoV-2-RBD with comparable binding affinity to ACE2. To enhance its binding affinity, we employed computational protein design to redesign all CDRs of Fab CC12.3 and molecular dynamics (MD) to validate their predicted binding affinities by the MM-GBSA method. MD results show that the predicted binding affinities of the three best designed Fabs CC12.3 (CC12.3-D02, CC12.3-D05, and CC12.3-D08) are better than those of Fab CC12.3 and ACE2. Additionally, our results suggest that enhanced binding affinities of CC12.3-D02, CC12.3-D05, and CC12.3-D08 are caused by increased SARS-CoV-2-RBD binding interactions of CDRs L1 and L3. This study redesigned neutralizing antibodies with better predicted binding affinities to SARS-CoV-2-RBD than Fab CC12.3 and ACE2. They are promising candidates as neutralizing antibodies against SARS-CoV-2.

Global diversity of the gene encoding the Pfs25 protein-a Plasmodium falciparum transmission-blocking vaccine candidate

Background

Vaccines against the sexual stages of the malarial parasite Plasmodium falciparum are indispensable for controlling malaria and abrogating the spread of drug-resistant parasites. Pfs25, a surface antigen of the sexual stage of P. falciparum, is a leading candidate for transmission-blocking vaccine development. While clinical trials have reported that Pfs25-based vaccines are safe and effective in inducing transmission-blocking antibodies, the extent of the genetic diversity of Pfs25 in malaria endemic populations has rarely been studied. Thus, this study aimed to investigate the global diversity of Pfs25 in P. falciparum populations.

Methods

A database of 307 Pfs25 sequences of P. falciparum was established. Population genetic analyses were performed to evaluate haplotype and nucleotide diversity, analyze haplotypic distribution patterns of Pfs25 in different geographical populations, and construct a haplotype network. Neutrality tests were conducted to determine evidence of natural selection. Homology models of the Pfs25 haplotypes were constructed, subjected to molecular dynamics (MD), and analyzed in terms of flexibility and percentages of secondary structures.

Results

The Pfs25 gene of P. falciparum was found to have 11 unique haplotypes. Of these, haplotype 1 (H1) and H2, the major haplotypes, represented 70% and 22% of the population, respectively, and were dominant in Asia, whereas only H1 was dominant in Africa, Central America, and South America. Other haplotypes were rare and region-specific, resulting in unique distribution patterns in different geographical populations. The diversity in Pfs25 originated from ten single-nucleotide polymorphism (SNP) loci located in the epidermal growth factor (EGF)-like domains and anchor domain. Of these, an SNP at position 392 (GGA/GCA), resulting in amino acid substitution 131 (Gly/Ala), defined the two major haplotypes. The MD results showed that the structures of H1 and H2 variants were relatively similar. Limited polymorphism in Pfs25 could likely be due to negative selection.

Conclusions

The study successfully established a Pfs25 sequence database that can become an essential tool for monitoring vaccine efficacy, designing assays for detecting malaria carriers, and conducting epidemiological studies of P. falciparum. The discovery of the two major haplotypes, H1 and H2, and their conserved structures suggests that the current Pfs25-based vaccines could be used globally for malaria control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05078-6.

Computational design of SARS-CoV-2 peptide binders with better predicted binding affinities than human ACE2 receptor

SARS-CoV-2 is coronavirus causing COVID-19 pandemic. To enter human cells, receptor binding domain of S1 subunit of SARS-CoV-2 (SARS-CoV-2-RBD) binds to peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptor. Employing peptides to inhibit binding between SARS-CoV-2-RBD and ACE2-PD is a therapeutic solution for COVID-19. Previous experimental study found that 23-mer peptide (SBP1) bound to SARS-CoV-2-RBD with lower affinity than ACE2. To increase SBP1 affinity, our previous study used residues 21–45 of α1 helix of ACE2-PD (SPB25) to design peptides with predicted affinity better than SBP1 and SPB25 by increasing interactions of residues that do not form favorable interactions with SARS-CoV-2-RBD. To design SPB25 with better affinity than ACE2, we employed computational protein design to increase interactions of residues reported to form favorable interactions with SARS-CoV-2-RBD and combine newly designed mutations with the best single mutations from our previous study. Molecular dynamics show that predicted binding affinities of three peptides (SPB25_Q22R, SPB25_{F8R/K11W/L25R} and SPB25_{F8R/K11F/Q22R/L25R}) are better than ACE2. Moreover, their predicted stabilities may be slightly higher than SBP1 as suggested by their helicities. This study developed an approach to design SARS-CoV-2 peptide binders with predicted binding affinities better than ACE2. These designed peptides are promising candidates as SARS-CoV-2 inhibitors.

Computational Design of Oligosaccharide-Producing Levansucrase from Bacillus licheniformis RN-01 to Increase Its Stability at High Temperature

Levan-type fructooligosaccharides (LFOs) and levan can potentially be used as ingredients in prebiotics, skincare products, and antitumor agents. The Y246S mutant of Bacillus licheniformis RN-01 levansucrase (oligosaccharide-producing levansucrase, OPL) was reported to productively synthesize LFOs; however, OPL's thermostability is low at high temperatures. To enhance OPL structural stability, this study employed molecular dynamics (AMBER) to identify a highly flexible region, as measured by its average root-mean-square fluctuation (RMSF) value, on the OPL surface and computational protein design (Rosetta) to rigidify and increase favorable interactions to increase its structural stability. AMBER identified region nine (residues 277-317) as a highly flexible region that was selected for design because it has the highest number of residues and the second-highest average RMSF, and it is farthest from the active site. Rosetta designed 14 mutants with the best ΔΔG value in each position, where three mutants have better ΔG than OPL. To determine whether their flexibilities and stabilities are lower than those of OPL, all 14 designed mutants were simulated at high temperature (500 K), and we found that K296E, G309S, and A310W mutants were predicted to be more stable and could retain their native structures better than OPL. Our results suggest that enhanced structural stabilities of these mutants are caused by increased hydrogen bond strengths of the designed residues and their neighboring residues. This study designed K296E, G309S, and A310W mutants of OPL with high potential for stability improvement, and they could potentially be used for the effective production of LFOs.

Unravelling Regioselectivity of Leuconostoc citreum ABK-1 Alternansucrase by Acceptor Site Engineering

Alternansucrase (ALT, EC 2.4.1.140) is a glucansucrase that can generate α-(1,3/1,6)-linked glucan from sucrose. Previously, the crystal structure of the first alternansucrase from Leuconostoc citreum NRRL B-1355 was successfully elucidated; it showed that alternansucrase might have two acceptor subsites (W675 and W543) responsible for the formation of alternating linked glucan. This work aimed to investigate the primary acceptor subsite (W675) by saturated mutagenesis using Leuconostoc citreum ABK-1 alternansucrase (LcALT). The substitution of other residues led to loss of overall activity, and formation of an alternan polymer with a nanoglucan was maintained when W675 was replaced with other aromatic residues. Conversely, substitution by nonaromatic residues led to the synthesis of oligosaccharides. Mutations at W675 could potentially cause LcALT to lose control of the acceptor molecule binding via maltose-acceptor reaction-as demonstrated by results from molecular dynamics simulations of the W675A variant. The formation of α-(1,2), α-(1,3), α-(1,4), and α-(1,6) linkages were detected from products of the W675A mutant. In contrast, the wild-type enzyme strictly synthesized α-(1,6) linkage on the maltose acceptor. This study examined the importance of W675 for transglycosylation, processivity, and regioselectivity of glucansucrases. Engineering glucansucrase active sites is one of the essential approaches to green tools for carbohydrate modification.

Cloning and Expression of Levansucrase Gene of Bacillus velezensis BM-2 and Enzymatic Synthesis of Levan

Levan is a versatile and valuable fructose homopolymer, and a few bacterial strains have been found to produce levan. Although levan products have numerous specific functions, their application and promotion were limited by the production capacity and production cost. Bacillus velezensis BM-2 is a levan-synthesizing strain, but its levan production is too low to apply. In this study, the levansucrase gene of B. velezensis BM-2 was cloned to plasmid pET-32a-Acma-zz, and the recombinant plasmids were transferred to Escherichia coli BL21. A transformed clone was selected to express and secrete the fusion enzymes with an Acma-tag efficiently. The expressed products were further purified by a self-developed separating material called bacterial enhancer matrix (BEM) particles. The purification efficiency was 93.4%, with a specific activity of 16.589 U/mL protein. The enzymatic reaction results indicated that the optimal reaction temperature is 50 °C, the optimal pH of the acetate buffer is 5.6, and the buffer system greatly influenced the enzyme activity. The enzyme activity was enhanced to 130% in the presence of 5 mM Ca2+, K+, Zn2+, and Mn2+, whereas it was almost abolished in the case of Cu2+ and Fe3+. The values of Km, kcat, and kcat/Km were 17.41 mM, 376.83 s−1, and 21.64 mM−1s−1, respectively. The enzyme amount of 20 U/g sucrose was added to the system containing 400 g/L sucrose, and the levan products with a concentration of 120 g/L reached after an incubation of 18 h, which was 8 times that of the yield before optimization. The results of molecular docking analysis indicated that the Asp86 might act as a nucleophilic catalytic residue for sucrose, Arg246 and Asp247 act as transition state stabilizer of transfructosylation, and Glu340 and Arg306 were recognized as general acid donors. They formed the catalytic-groups triad. The unique properties and catalytic activity of the levansucrase suggest that it deserves further research and might have good industrial application prospects.

The molecular basis of the nonprocessive elongation mechanism in levansucrases

Levansucrases (LSs) synthesize levan, a β2-6-linked fructose polymer, by successively transferring the fructosyl moiety from sucrose to a growing acceptor molecule. Elucidation of the levan polymerization mechanism is important for using LSs in the production of size-defined products for application in the food and pharmaceutical industries. For a deeper understanding of the levan synthesis reaction, we determined the crystallographic structure of Bacillus subtilis LS (SacB) in complex with a levan-type fructooligosaccharide and utilized site-directed mutagenesis to identify residues involved in substrate binding. The presence of a levanhexaose molecule in the central catalytic cavity allowed us to identify five substrate-binding subsites (−1, +1, +2, +3, and +4). Mutants affecting residues belonging to the identified acceptor subsites showed similar substrate affinity (Km) values to the wildtype (WT) Km value but had a lower turnover number and transfructosylation/hydrolysis ratio. Of importance, compared with the WT, the variants progressively yielded smaller-sized low-molecular-weight levans, as the affected subsites that were closer to the catalytic site, but without affecting their ability to synthesized high-molecular-weight levans. Furthermore, an additional oligosaccharide-binding site 20 Å away from the catalytic pocket was identified, and its potential participation in the elongation mechanism is discussed. Our results clarify, for the first time, the interaction of the enzyme with an acceptor/product oligosaccharide and elucidate the molecular basis of the nonprocessive levan elongation mechanism of LSs.

Computational Design of 25-mer Peptide Binders of SARS-CoV-2

SARS-CoV-2 is the novel coronavirus causing the COVID-19 pandemic. To enter human cells, the receptor-binding domain (RBD) of the S1 subunit of SARS-CoV-2 (SARS-CoV-2-RBD) initially binds to the peptidase domain of angiotensin-converting enzyme 2 receptor (ACE2-PD). Using peptides to inhibit SARS-CoV-2-RBD binding to ACE2 is a potential therapeutic solution for COVID-19. A previous study identified a 23-mer peptide (SBP1) that bound to SARS-CoV-2-RBD with comparable KD to ACE2. We employed computational protein design and molecular dynamics (MD) to design SARS-CoV-2-RBD 25-mer peptide binders (SPB25) with better predicted binding affinity than SBP1. Using residues 21-45 of the α1 helix of ACE2-PD as the template, our strategy is employing Rosetta to enhance SPB25 binding affinity to SARS-CoV-2-RBD and avoid disrupting existing favorable interactions by using residues that have not been reported to form favorable interactions with SARS-CoV-2-RBD as designed positions. Designed peptides with better predicted binding affinities, by Rosetta, than SPB25 were subjected to MD validation. The MD results show that five designed peptides (SPB25F8N, SPB25F8R, SPB25L25R, SPB25F8N/L25R, and SPB25F8R/L25R) have better predicted binding affinities, by the MM-GBSA method, than SPB25 and SBP1. This study developed an approach to design SARS-CoV-2-RBD peptide binders, and these peptides may be promising candidates as potential SARS-CoV-2 inhibitors.

Computational design of oligosaccharide producing levansucrase from Bacillus licheniformis RN-01 to improve its thermostability for production of levan-type fructooligosaccharides from sucrose

Levansucrase catalyzes production of levan and levan-type fructooligosaccharides (LFOs) with potential applications in food and pharmaceutical industries such as prebiotics and anti-tumor agents. Previous study found that Y246S mutant of Bacillus licheniformis RN-01 levansucrase (oligosaccharide producing levansucrase, OPL) could effectively produce LFOs but its thermostability is limited at high temperature. In this study, molecular dynamics (MD) and computational protein design were used to create mutants with higher thermostability than OPL by rigidifying highly flexible residues on enzyme surface. MD results show that highly flexible residues suitable for design are K82, N83, D179, and Q308. Two approaches were employed to improve their interactions by allowing them to be amino acids that could potentially form favorable interactions with their neighboring residues or natural amino acids except G, P and C. Flexibilities of designed residues of K82H, N83R, Q308S and K82H/N83R mutants are lower than those of OPL. Experimental results show that characteristics and product patterns of designed mutants are relatively similar to those of OPL. K82H/N83R mutant has higher thermostability than OPL with 1.7-fold increase in t_1/2. Circular dichroism result suggests that designed mutations do not drastically affect secondary structures. This study shows how computational technique can engineer enzyme for thermostability improvement.

Levansucrase from Bacillus amyloliquefaciens KK9 and Its Y237S Variant Producing the High Bioactive Levan-Type Fructooligosaccharides

Levan-typed fructooligosaccharide (LFOS), a β-2,6 linked oligofructose, displays the potential application as a prebiotic and therapeutic dietary supplement. In the present study, LFOS was synthesized using levansucrase from Bacillus amyloliquefaciens KK9 (LsKK9). The wild-type LsKK9 was cloned and expressed in E. coli, and purified by cation exchanger chromatography. Additionally, Y237S variant of LsKK9 was constructed based on sequence alignment and structural analysis to enhance the LFOS production. High-performance anion-exchange chromatography coupled with pulsed amperometric detection (HPAEC-PAD) analysis indicated that Y237S variant efficiently produced a higher amount of short-chain LFOS than wild type. Also, the concentration of enzyme and sucrose in the reactions was optimized. Finally, prebiotic activity assay demonstrated that LFOS produced by Y237S variant had higher prebiotic activity than that of the wild-type enzyme, making the variant enzyme attractive for food biotechnology.