Designing peptides predicted to bind to the omicron variant better than ACE2 via computational protein design and molecular dynamics

Brought about by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), coronavirus disease (COVID-19) pandemic has resulted in large numbers of worldwide deaths and cases. Several SARS-CoV-2 variants have evolved, and Omicron (B.1.1.529) was one of the important variants of concern. It gets inside human cells by using its S1 subunit's receptor-binding domain (SARS-CoV-2-RBD) to bind to Angiotensin-converting enzyme 2 receptor's peptidase domain (ACE2-PD). Using peptides to inhibit binding interactions (BIs) between ACE2-PD and SARS-CoV-2-RBD is one of promising COVID-19 therapies. Employing computational protein design (CPD) as well as molecular dynamics (MD), this study used ACE2-PD's α1 helix to generate novel 25-mer peptide binders (SPB25) of Omicron RBD that have predicted binding affinities (ΔGbind (MM‑GBSA)) better than ACE2 by increasing favorable BIs between SPB25 and the conserved residues of RBD. Results from MD and the MM-GBSA method identified two best designed peptides (SPB25T7L/K11A and SPB25T7L/K11L with ΔGbind (MM‑GBSA) of -92.4 ± 0.4 and -95.7 ± 0.5 kcal/mol, respectively) that have better ΔGbind (MM‑GBSA) to Omicron RBD than ACE2 (-87.9 ± 0.5 kcal/mol) and SPB25 (-71.6 ± 0.5 kcal/mol). Additionally, they were predicted to have slightly higher stabilities, based on their percent helicities in water, than SBP1 (the experimentally proven inhibitor of SARS-CoV-2-RBD). Our two best designed SPB25s are promising candidates as omicron variant inhibitors.

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.

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.

Molecular properties and ligand specificity of zebrafish intelectin-2

Intelectins are immune lectins expressed in chordates, including several fish species, in which intelectins are known to be upregulated upon infection. However, the basic biochemical properties and bacteria binding specificities of several fish intelectins are not well studied. We focus our investigation on zebrafish intelectin-2 (DrIntL-2) that is predominantly expressed in the gastrointestinal tract. The disulfide-linked oligomeric state and the cysteine responsible for intermolecular disulfide bonds are identified. DrIntL-2 is a globular particle of around 30 nm. In addition to the typical exocyclic 1,2-diol ligands, DrIntL-2 binds β-1,3-glucan and recognizes Salmonella typhimurium and Pseudomonas aeruginosa. This investigation not only shed light on the fish innate immunity that will be essential for the aquaculture industry, but will also provide a foundation for further application of DrIntL-2 in bacteria detection and identification.

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.

Biochemical and Structural Investigation of GnnA in the Lipopolysaccharide Biosynthesis Pathway of Acidithiobacillus ferrooxidans

Lipopolysaccharide (LPS) is a crucial component in the outer membrane of Gram-negative bacteria that contributes to both pathogenicity as well as immunity against pathogenic bacteria. Typical LPS contains GlcN disaccharide as the core of lipid A. However, some bacteria such as Acidithiobacillus ferrooxidans and Leptospira interrogans contain GlcN3N in lipid A instead. This modification has been shown to dampen the host immune response and increase resistance to antimicrobial peptides. Therefore, investigation of the enzymes responsible for the biosynthesis of GlcN3N has promising applications in the development of vaccines, antibiotics, or usage of the enzymes in chemoenzymatic synthesis of modified LPS. Here, we describe biochemical and structural investigation of GnnA from A. ferrooxidans (AfGnnA) that is responsible for oxidation of UDP-GlcNAc, which subsequently undergoes transamination to produce UDP-GlcNAc3N as a precursor for LPS biosynthesis. AfGnnA is specific for NAD⁺ and UDP-GlcNAc. The crystal structures of AfGnnA in combination with molecular dynamics simulation and mutational analysis suggest the substrate recognition mode and the catalytic mechanism. K91 or H164 is a potential catalytic base in the oxidation reaction. The results will not only provide insights into the biosynthesis of unusual LPS but will also lay the foundation for development of more immunogenic vaccines, novel antibiotics, or utilization of GnnA in the synthesis of UDP-sugars or modified LPS.

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.

Conserved Calcium-Binding Residues at the Ca-I Site Involved in Fructooligosaccharide Synthesis by Lactobacillus reuteri 121 Inulosucrase

Inulosucrase is an enzyme that synthesizes inulin-type β-2,1-linked fructooligosaccharides (IFOS) from sucrose. Previous studies have shown that calcium is important for the activity and stability of Lactobacillus reuteri 121 inulosucrase (LrInu). Here, mutational analyses of four conserved calcium-binding site I (Ca-I) residues of LrInu, Asp418, Gln449, Asn488, and Asp520 were performed. Alanine substitution for these residues not only reduced the stability and activity of LrInu, but also modulated the pattern of the IFOS produced. Circular dichroism spectroscopy and molecular dynamics simulation indicated that these mutations had limited impact on the overall conformation of the enzyme. One of Ca-I residues most critical for controlling LrInu-mediated polymerization of IFOS, Asp418, was also subjected to mutagenesis, generating D418E, D418H, D418L, D418N, D418S, and D418W. The activity of these mutants demonstrated that the IFOS chain length could be controlled by a single mutation at the Ca-I site.

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.

Molecular dynamics reveals insight into how N226P and H227Y mutations affect maltose binding in the active site of α-glucosidase II from European honeybee, Apis mellifera

European honeybee, Apis mellifera, produces α-glucosidase (HBGase) that catalyzes the cleavage of an α-glycosidic bond of the non-reducing end of polysaccharides and has potential applications for malt hydrolysis in brewing industry. Characterized by their substrate specificities, HBGases have three isoforms including HBGase II, which prefers maltose to sucrose as a substrate. Previous study found that the catalytic efficiency of maltose hydrolysis of N226P mutant of HBGase II was higher than that of the wild type (WT), and the catalytic efficiency of maltose hydrolysis of WT was higher than those of H227Y and N226P-H227Y mutants. We hypothesized that N226P mutation probably caused maltose to bind with better affinity and position/orientation for hydrolysis than WT, while H227Y and N226P-H227Y mutations caused maltose to bind with worse affinity and position/orientation for hydrolysis than WT. Using this hypothesis, we performed molecular dynamics on the catalytically competent binding conformations of maltose/WT, maltose/N226P, maltose/H227Y, and maltose/N226P-H227Y complexes to elucidate effects of N226P and H227Y mutations on maltose binding in HBGase II active site. Our results reasonably support this hypothesis because the N226P mutant had better binding affinity, higher number of important binding residues, strong and medium hydrogen bonds as well as shorter distance between atoms necessary for hydrolysis than WT, while the H227Y and N226P-H227Y mutants had worse binding affinities, lower number of important binding residues and strong hydrogen bonds as well as longer distances between atoms necessary for hydrolysis than WT. Moreover, results of binding free energy and hydrogen bond interaction of residue 227 support the role of H227 as a maltose preference residue, as proposed by previous studies. Our study provides important and novel insight into how N226P and H227Y mutations affect maltose binding in HBGase II active site. This knowledge could potentially be used to engineer HBGase II to improve its efficiency.

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.

Rational re-design of Lactobacillus reuteri 121 inulosucrase for product chain length control

Fructooligosaccharides (FOSs) are well-known prebiotics that are widely used in the food, beverage and pharmaceutical industries. Inulosucrase (E.C. 2.4.1.9) can potentially be used to synthesise FOSs from sucrose. In this study, inulosucrase from Lactobacillus reuteri 121 was engineered by site-directed mutagenesis to change the FOS chain length. Three variants (R483F, R483Y and R483W) were designed, and their binding free energies with 1,1,1-kestopentaose (GF4) were calculated with the Rosetta software. R483F and R483Y were predicted to bind with GF4 better than the wild type, suggesting that these engineered enzymes should be able to effectively extend GF4 by one residue and produce a greater quantity of GF5 than the wild type. MALDI-TOF MS analysis showed that R483F, R483Y and R483W variants could synthesise shorter chain FOSs with a degree of polymerization (DP) up to 11, 10, and 10, respectively, while wild type produced longer FOSs and in polymeric form. Although the decrease in catalytic activity and the increase of hydrolysis/transglycosylation activity ratio was observed, the variants could effectively synthesise FOSs with the yield up to 73% of substrate. Quantitative analysis demonstrated that these variants produced a larger quantity of GF5 than wild type, which was in good agreement with the predicted binding free energy results. Our findings demonstrate the success of using aromatic amino acid residues, at position D418, to block the oligosaccharide binding track of inulosucrase in controlling product chain length.

Modulation of fructooligosaccharide chain length and insight into the product binding motif of Lactobacillus reuteri 121 inulosucrase

Inulosucrase (E.C. 2.4.1.9) is a bacterial fructosyltransferase that synthesizes inulin-type fructooligosaccharide, using sucrose as a substrate. We modulated the size of fructooligosaccharide synthesized by Lactobacillus reuteri 121 inulosucrase using rational designed mutagenesis. Nine residues: D478, D479, S482, R483, N543, W551, N555, N561 and D689, were changed based on the active site architecture and amino acids potentially interacting with saccharides. The selected residues were substituted with alanine to investigate the contribution of these residues to FOS chain length. Enzymatic activity assays demonstrated that the transglycosylation/hydrolysis ratios of D479A, R483A, N543A, W551A and N555A mutants were significantly different from that of the wild type. Almost all mutants, except D478A, synthesized oligosaccharides with different size distribution compared to that of wild type. Molecular docking further provides insights into the product binding motif of Lactobacillus reuteri 121 inulosucrase and strengthens an important role of amino acid residues at remote locations from the active site on the enzymatic activity and product specificity.

Molecular dynamics provides insight into how N251A and N251Y mutations in the active site of Bacillus licheniformis RN-01 levansucrase disrupt production of long-chain levan

Produced by levansucrase, levan and levan oligosaccharides (GF_n) have potential applications in food and pharmaceutical industries such as prebiotics, anti-tumor and anti-inflammatory agents. Previous study reported that Bacillus licheniformis RN-01 levansucrase could produce levan oligosaccharides and long-chain levan. However, its N251A and N251Y mutants could effectively produce short-chain oligosaccharides upto GF_3, but they could not produce long-chain levan. We hypothesized that these mutations probably reduced GF₃ binding affinity in levansucrase active site that contains fructosyl-Asp93 intermediate and caused GF₃ to be in an unfavorable orientation for transfructosylation; therefore, levansucrase could not effectively extend GF₃ by one fructosyl residue to produce GF₄ and subsequently long-chain levan. However, these mutations probably did not significantly reduce binding affinity or drastically change orientation of GF₂; therefore, levansucrase could still extend GF₂ to produce GF₃. Using this hypothesis, we employed molecular dynamics to investigate effects of these mutations on GF₂/GF₃ binding in levansucrase active site. Our results reasonably support this hypothesis as N251A and N251Y mutations did not significantly reduce GF₂ binding affinity, as calculated by MM-GBSA technique and hydrogen bond occupations, or drastically change orientation of GF₂ in levansucrase active site, as measured by distance between atoms necessary for transfructosylation. However, these mutations drastically decreased GF₃ binding affinity and caused GF₃ to be in an unfavorable orientation for transfructosylation. Furthermore, the free energy decomposition and hydrogen bond occupation results suggest the importance of Arg255 in GF₂/GF₃ binding in levansucrase active site. This study provides important and novel insight into the effects of N251A and N251Y mutations on GF₂/GF₃ binding in levansucrase active site and how they may disrupt production of long-chain levan. This knowledge could be beneficial in designing levansucrase to efficiently produce levan oligosaccharides with desired length.

Insight into the substrate specificity change caused by the Y227H mutation of α-glucosidase III from the European honeybee (Apis mellifera) through molecular dynamics simulations

Honey from the European honeybee, Apis mellifera, is produced by α-glucosidases (HBGases) and is widely used in food, pharmaceutical, and cosmetic industries. Categorized by their substrate specificities, HBGases have three isoforms: HBGase I, II and III. Previous experimental investigations showed that wild-type HBGase III from Apis mellifera (WT) preferred sucrose to maltose as a substrate, while the Y227H mutant (MT) preferred maltose to sucrose. This mutant can potentially be used for malt hydrolysis because it can efficiently hydrolyze maltose. In this work, to elucidate important factors contributing to substrate specificity of this enzyme and gain insight into how the Y227H mutation causes substrate specificity change, WT and MT homology models were constructed, and sucrose/maltose was docked into active sites of the WT and MT. AMBER14 was employed to perform three independent molecular dynamics runs for these four complexes. Based on the relative binding free energies calculated by the MM-GBSA method, sucrose is better than maltose for WT binding, while maltose is better than sucrose for MT binding. These rankings support the experimentally observed substrate specificity that WT preferred sucrose to maltose as a substrate, while MT preferred maltose to sucrose, suggesting the importance of binding affinity for substrate specificity. We also found that the Y227H mutation caused changes in the proximities between the atoms necessary for sucrose/maltose hydrolysis that may affect enzyme efficiency in the hydrolysis of sucrose/maltose. Moreover, the per-residue binding free energy decomposition results show that Y227/H227 may be a key residue for preference binding of sucrose/maltose in the WT/MT active site. Our study provides important and novel insight into the binding of sucrose/maltose in the active site of Apis mellifera HBGase III and into how the Y227H mutation leads to the substrate specificity change at the molecular level. This knowledge could be beneficial in the design of this enzyme for increased production of desired products.

Isomahanine induces endoplasmic reticulum stress and simultaneously triggers p38 MAPK-mediated apoptosis and autophagy in multidrug-resistant human oral squamous cell carcinoma cells

Advanced oral squamous cell carcinoma (OSCC) is typically aggressive and closely correlated with disease recurrence and poor survival. Multidrug resistance (MDR) is the most critical problem leading to therapeutic failure. Investigation of novel anticancer candidates targeting multidrug-resistant OSCC cells may provide a basis for developing effective strategies for OSCC treatment. In the present study, we investigated the cytotoxic mechanism of a carbazole alkaloid, namely isomahanine, in a multidrug‑resistant OSCC cell line CLS-354/DX. We demonstrated that CLS-354/DX cells overexpressing multidrug resistance-associated protein 1 (MRP1) were resistant to anticancer drugs cisplatin and camptothecin. Isomahanine effectively induced cytotoxicity against CLS-354/DX cells regardless of resistance. Apoptosis as determined by FITC‑Annexin V/PI staining and western blot analysis of cleaved caspase-3 and cleaved poly(ADP‑ribose) polymerase (PARP) was significantly induced in a time-dependent manner upon isomahanine treatment. Isomahanine-induced caspase‑dependent apoptosis was determined using z-VAD‑fmk. The effects on autophagy in isomahanine-treated cells were investigated via conversion of LC3B and degradation of p62/SQSTM1 (p62). Isomahanine obviously induced autophagic flux as shown by an increase in punctate GFP-LC3B and the LC3B-II/LC3B-I ratio with a concomitant decrease in p62 levels. Autophagy inhibitors 3-methyladenine (3-MA) and chloroquine (CQ) protected isomahanine-induced cell death, indicating the activation of autophagic cell death. Endoplasmic reticulum (ER) stress and MAPK activation were examined to elucidate the mechanism underlying cell death. The expression levels of PERK, CHOP and phosphorylated MAPK (p38, ERK1/2 and JNK1/2) were upregulated following isomahanine treatment. We found that p38 MAPK inhibitor (SB203580) significantly attenuated isomahanine-induced apoptosis and autophagic flux and this prevented cell death. Collectively, the present study demonstrated that isomahanine was able to induce ER stress and trigger p38 MAPK-mediated apoptosis and autophagic cell death in multidrug-resistant OSCC cells. The potential cytotoxic action of isomahanine may provide the development of anticancer candidates for treating multidrug-resistant cancer.

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.

Design of peptides as inhibitors of human papillomavirus 16 transcriptional regulator E1-E2

Human papillomavirus 16 (HPV 16) is a DNA virus that is capable of infecting humans and causing cervical cancer. HPV16 E2 plays an important role in viral gene regulation. This work aims to predict the binding conformations and interactions between the dodecapeptides and HPV16 E2 as well as to design novel peptide inhibitors that are capable of binding to HPV16 E2 and disrupt the transcriptional regulator E1-E2 complex formation, using computational protein design techniques. Based on previously reported peptide4 (TWFWPYPYPHLP), novel peptide inhibitors were designed and five peptides that showed lower binding energy to HPV16 E2 than that of peptide4, were selected for in vitro experiments. Enzyme-linked immunosorbent (ELISA) assay showed that Y6R, W4H_Y6R, and W4H peptides bound to HPV16 E2 with higher affinity than peptide4 did. Moreover, Y6R, W4H_Y6R, and W4H peptides more effectively inhibited E1-E2 complex formation than peptide4. This work revealed important interactions between the peptides and E1-E2 complex, suggesting a strategy for development of more potent peptide inhibitors.

Molecular dynamics of interactions between rigid and flexible antifolates and dihydrofolate reductase from pyrimethamine-sensitive and pyrimethamine-resistant Plasmodium falciparum

Currently, the usefulness of antimalarials such as pyrimethamine (PYR) is drastically reduced due to the emergence of resistant Plasmodium falciparum (Pf) caused by its dihydrofolate reductase (PfDHFR) mutations, especially the quadruple N51I/C59R/S108N/I164L mutations. The resistance was due to the steric conflict of PYR with S108N. WR99210 (WR), a dihydrotriazine antifolate with a flexible side chain that can avoid such conflict, can overcome this resistance through tight binding with the mutant. To understand factors contributing to different binding affinities of PYR/WR to the wild type (WT) and quadruple mutant (QM), we performed simulations on WR-WT, WR-QM, PYR-WT, and PYR-QM complexes and found that Ile14 and Asp54 were crucial for PYR/WR binding to PfDHFR due to strong hydrogen bonds. The quadruple mutations cause PYR to form, on average, fewer hydrogen bonds with Ile14 and Leu164, and to be displaced from its optimal orientation for Asp54 interaction. The predicted binding affinity ranking (WR-QM ≈ WR-WT ≈ PYR-WT >> PYR-QM) reasonably agrees with the inhibition constant (K(i)) ranking. Our results reveal important residues for tight binding of PYR/WR to WT/QM, which may be used to evaluate the inhibition effectiveness of antimalarials and to provide fundamental information for designing new drugs effective against drug-resistant P. falciparum.

Binding affinity of a small molecule to an amorphous polymer in a solvent. Part 1: free energy of binding to a binding site

Crystallization is commonly used in a separation and purification process in the production of a wide range of materials in various industries. In industry, crystallization usually starts with heterogeneous nucleation on a foreign surface. The complicated mechanism of heterogeneous nucleation is not well understood; however, we hypothesize that there might be a possible correlation between binding affinity to a surface and enhancement of nucleation. Recent studies show that amorphous polymers can be used to control crystallization, selectively produce pharmaceutical polymorphs, and discover novel pharmaceutical polymorphs. To investigate the possible correlation between the binding affinity of one molecule to key binding sites (local binding) and heterogeneous nucleation activity as well as the possibility of using this binding affinity to help guide the selection of polymers that promote heterogeneous nucleation, we computed the free energy of binding of aspirin to four nonporous cross-linked polymers in an ethanol-water 38 v% mixture. These cross-linked polymers are poly(4-acryloylmorpholine) (PAM), poly(2-carboxyethyl acrylate) (PCEA), poly(4-hydroxylbutyl acrylate) (PHBA), and polystyrene (PS); all of them were cross-linked with divinylbenzene (DVB). These systems were used because their heterogeneous nucleation activities are available in literature, and the ranking is PAM > PCEA > PHBA ≈ PS. We generated three independent surfaces for each polymer and computed the free energy of binding of aspirin to the best binding site that we found on each surface. The average free energies of binding to the best sites of PAM, PCEA, PHBA, and PS are -20.4 ± 1.0, -16.7 ± 1.0, -14.4 ± 1.1, and -13.6 ± 1.1 kcal/mol, respectively. We found that the trend of the magnitudes of the average free energies of binding to the best sites is PAM > PCEA > PHBA ≈ PS. This trend is very similar to that of heterogeneous nucleation activity. Our results suggest the importance of the free energy of binding to key sites (local binding) and the possibility of using this quantity to help guide the selection of polymers that promote heterogeneous nucleation.

Binding affinity of a small molecule to an amorphous polymer in a solvent. Part 2: preferential binding to local sites on a surface

Crystallization, a separation and purification process, is commonly used to produce a wide range of materials in various industries, and it usually begins with heterogeneous nucleation on a foreign surface in industrial practice and most other circumstances. Recent studies show that amorphous polymeric substrates are useful in controlling crystallization and selectively producing pharmaceutical polymorphs. In our previous publication, we investigated the possible correlation of the binding affinity of one molecule to key binding sites (local binding), and the possibility of using this binding affinity to guide the selection of polymers promoting heterogeneous nucleation. The studied systems were aspirin binding to four nonporous cross-linked polymers in ethanol-water 38 v% mixture. Cross-linked with divinylbenzene (DVB), these polymers were poly(4-acryloylmorpholine) (PAM), poly(2-carboxyethyl acrylate) (PCEA), poly(4-hydroxylbutyl acrylate) (PHBA), and polystyrene (PS). We discovered that the trend of the magnitudes of the average free energies of binding to the best sites is very similar to that of heterogeneous nucleation activities. This Article aims to investigate whether or not local binding to key sites is the important variable to describe heterogeneous nucleation as opposed to the overall/average binding affinity of molecules to a surface, and to investigate the possibility of using the overall binding affinity to guide the selection of polymers. We used the polymer surfaces generated from our previous study to calculate the overall binding affinity of aspirin molecules to the surface as measured by the preferential interaction coefficients of aspirin (1 m) to these polymers. We discovered that the trend of the average preferential interaction coefficients does not correlate as well to that of heterogeneous nucleation activities as the free energies of binding to the best sites. We also computed the average numbers of aspirin molecules associated with the areas of the surfaces' best binding sites and found that they correlate better to heterogeneous nucleation activities than the average preferential interaction coefficients. These results further support that local binding is indicative of heterogeneous nucleation. Moreover, we found a weak trend of the distance order parameters of the aspirin molecules to be similar that of heterogeneous nucleation activities. Our results from the two-part study suggest the importance of local binding to heterogeneous nucleation as well as the possibility of using the binding affinity to the local area (the free energy of binding to the best site and the number of nucleating molecules associated with the area of the best binding site) and the distance order parameters to guide the selection of polymers.

Free energy of binding of a small molecule to an amorphous polymer in a solvent

Crystallization is a commonly used purification process in industrial practice. It usually begins with heterogeneous nucleation on a foreign surface. The complicated mechanism of heterogeneous nucleation is not well understood, but we hypothesize that a possible correlation between binding affinity to a surface and nucleation enhancement might exist. Amorphous polymers have been used in controlling crystallization. However, to our knowledge, no attempt has been made to calculate the free energy of binding of a small molecule to an amorphous polymer in a solvent, and to characterize the binding sites/conformations of this system at a molecular level. We developed a two-step approach, first using Adsorption Locator to identify probable binding sites and molecular dynamics to screen for the best binding sites and then using the Blue-Moon Ensemble method to compute the free energy of binding. A system of ethylene glycol, polyvinyl alcohol (PVA), and heavy water (D(2)O) was used for validation, since experimental data exists on a related system. Looking at four independently constructed surfaces, we found that ethylene glycol binds to an indentation on the surface or in a hole beneath the surface. We focused on the indentation binding sites because they are easily accessible and do not have large free energy barriers. The closest system for which experimental data on binding energetics exists is ethylene glycol on PVA in aqueous solutions/gels, and the magnitudes of the free energy of binding to the three best indentation binding sites are close to the experimental value, 0.4-3.7 kcal/mol higher. Our approach offers a way to compute the free energy of binding and characterize the binding sites/conformations, and is general enough to apply to other small molecule/amorphous polymer/solvent systems.