QM/MM Simulations of Enzymatic Hydrolysis of Cellulose: Probing the Viability of an Endocyclic Mechanism for an Inverting Cellulase

Molecular Mechanism of Double-Displacement Retaining β-Kdo Glycosyltransferase WbbB

Glycosyltransferases (GTs) are pivotal enzymes involved in glycosidic bond synthesis, which can lead to either retention or inversion of the glycosyl moiety's anomeric configuration. However, the catalytic mechanism for retaining GTs remains a subject of controversy. In this study, we employ MD and QM/MM metadynamics to investigate the double-displacement catalytic mechanism of the retaining β-Kdo transferase WbbB. Our findings demonstrate that the nucleophile Asp232 initiates the reaction by attacking the sugar ring containing a carboxylate at the anomeric position, forming a covalent adduct. Subsequently, the adduct undergoes a rotational rearrangement, ensuring proper orientation of the anomeric carbon for the acceptor substrate. In the second step, Glu158 acts as the catalytic base to abstract the proton of the acceptor substrate to complete the transglycosylation reaction. Notably, His265 does not function as the anticipated catalytic acid; instead, it stabilizes the phosphate group through H-bonding interactions. Our simulations support the double-displacement mechanism implicated from the crystallographic studies of WbbB. This mechanism deviates from the common SNi-type and retaining glycoside hydrolase mechanisms, thereby expanding our understanding of GT catalytic mechanisms.

DES-ultrasonication treatment of cellulose nanocrystals and the reinforcement in carrageenan biocomposite

CNCs are intensively studied to reinforce biocomposites. However, it remains a challenge to homogeneously disperse the CNC in biocomposites for a smooth film surface. Mechanochemical treatment via ultrasonication in deep eutectic solvent (DES) generated a stable dispersion of CNC before incorporation into carrageenan biocomposite. Shifted peaks of choline chloride (ChCl) methylene groups to 3.95-3.98 ppm in 1H NMR indicated a formation of eutectic mixture between the hydrogen bond acceptor (HBA) and hydrogen bond donor (HBD) at the functional group of CH3···OH. The swelling of CNC in the DES was proven by the formation of intermolecular H-bond at a length of 2.46 Å. The use of DES contributed to a good dispersion of CNC in the solution which increased zeta potential by 43.2 % compared to CNC in deionized water. The ultrasonication amplitude and feed concentration were varied for the best parameters of a stable dispersion of CNC. The crystallinity of 1 wt% of CNC at 20 % sonication amplitude improved from 76 to 81 %. The high crystallinity of CNCDES resulted in an increase in film tensile and capsule loop strength of Carra-CNCDES by 20.7 and 19.4 %, respectively. Improved dispersion of CNCDES reduced the surface roughness of the biocomposite by 21.8 %. H-bond network in CNCDES improved the biocomposite properties for an ingenious reinforcement material.

Catalytic Reaction Mechanism of Bacterial GH92 α-1,2-Mannosidase: A QM/MM Metadynamics Study

The catalytic mechanism of aC a + 2 ${C{a}^{+2}}$ -dependent family 92 α ${{\rm \alpha }}$ -mannosidase, which is abundantly present in human gut flora and malfunctions leading to the lysosomal storage disease α-mannosidosis, has been investigated using quantum mechanics/molecular mechanics and metadynamics methods. Computational efforts show that the enzyme follows a conformational itinerary of and theC a + 2 ${C{a}^{+2}}$ ion serves a dual purpose, as it not only distorts the sugar ring but also plays a crucial role in orchestrating the arrangement of catalytic residues. This orchestration, in turn, contributes to the facilitation of O S 2 ${{{\rm \ }}^{{\rm O}}{{\rm S}}_{2}}$ conformers for the ensuing reaction. This mechanistic insight is well-aligned with the experimental predictions of the catalytic pathway, and the computed energies are of the same order of magnitude as the experimental estimations. Hence, our results extend the mechanistic understanding of glycosidases.

Glycosidase mechanisms: Sugar conformations and reactivity in endo- and exo-acting enzymes

The enzymatic breakdown of carbohydrates plays a critical role in several biological events and enables the development of sustainable processes to obtain bioproducts and biofuels. In this scenario, the design of efficient inhibitors for glycosidases that can act as drug targets and the engineering of carbohydrate-active enzymes with tailored catalytic properties is of remarkable importance. To guide rational approaches, it is necessary to elucidate enzyme molecular mechanisms, in particular understanding how the microenvironment modulates the conformational space explored by the substrate. Computer simulations, especially those based on ab initio methods, have provided a suitable atomic description of carbohydrate conformations and catalytic reactions in several glycosidase families. In this review, we will focus on how the active-site topology (pocket or cleft) and mode of cleavage (endo or exo) can affect the catalytic mechanisms adopted by glycosidases, in particular the substrate conformations along the reaction coordinate.

Mechanism of Chiral-Selective Aminoacylation of an RNA Minihelix Explored by QM/MM Free-Energy Simulations

Aminoacylation of a primordial RNA minihelix composed of D-ribose shows L-amino acid preference over D-amino acid without any ribozymes or enzymes. This preference in the amino acylation reaction likely plays an important role in the establishment of homochirality in L-amino acid in modern proteins. However, molecular mechanisms of the chiral selective reaction remain unsolved mainly because of difficulty in direct observation of the reaction at the molecular scale by experiments. For seeking a possible mechanism of the chiral selectivity, quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations of the aminoacylation reactions in a modeled RNA were performed to investigate differences in their free-energy profiles along the reactions for L- and D-alanine and its physicochemical origin. The reaction is initiated by approaching a 3′-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. The QM/MM umbrella sampling MD calculations showed that the height of the free-energy barrier for L-alanine aminoacylation reaction was 17 kcal/mol, which was 9 kcal/mol lower than that for the D-alanine system. At the transition state, the distance between the negatively charged 3′-oxygen and the positively charged amino group of L-alanine was shorter than that of D-alanine, which was caused by the chirality difference of the amino acid. These results indicate that the transition state for L-alanine is more electrostatically stabilized than that for D-alanine, which would be a plausible mechanism previously unexplained for chiral selectivity in the RNA minihelix aminoacylation.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

Purification and enzymatic properties of a new thermostable endoglucanase from Aspergillus oryzae HML366

Aspergillus oryzae HML366 is a newly screened cellulase-producing strain. The endoglucanase HML ED1 from A. oryzae HML366 was quickly purified by a two-step method that combines ammonium sulfate precipitation and strong anion exchange column. SDS-PAGE electrophoresis indicated that the molecular weight of the enzyme was 68 kDa. The optimum temperature of the purified endoglucanase was 60 ℃ and the enzyme activity was stable below 70 ℃. The optimum pH was 6.5, and the enzyme activity was stable at pH between 4.5 and 9.0. The analysis indicated that additional Na⁺, K⁺, Ca²⁺, and Zn²⁺ reduced the catalytic ability of enzyme to the substrate, but Mn²⁺ enhanced its catalytic ability to the substrate.The Km and Vmax of the purified endoglucanase were 8.75 mg/mL and 60.24 μmol/min·mg, respectively. In this study, we report for the first time that A. oryzae HML366 can produce a heat-resistant and wide pH tolerant endoglucanase HML ED1, which has potential industrial application value in bioethanol, paper, food, textile, detergent, and pharmaceutical industries.

Reconstructing reactivity in dynamic host-guest systems at atomistic resolution: amide hydrolysis under confinement in the cavity of a coordination cage

Spatial confinement is widely employed by nature to attain unique efficiency in controlling chemical reactions. Notable examples are enzymes, which selectively bind reactants and exquisitely regulate their conversion into products. In an attempt to mimic natural catalytic systems, supramolecular metal-organic cages capable of encapsulating guests in their cavity and of controlling/accelerating chemical reactions under confinement are attracting increasing interest. However, the complex nature of these systems, where reactants/products continuously exchange in-and-out of the host, makes it often difficult to elucidate the factors controlling the reactivity in dynamic regimes. As a case study, here we focus on a coordination cage that can encapsulate amide guests and enhance their hydrolysis by favoring their mechanical twisting towards reactive molecular configurations under confinement. We designed an advanced multiscale simulation approach that allows us to reconstruct the reactivity in such host-guest systems in dynamic regimes. In this way, we can characterize amide encapsulation/expulsion in/out of the cage cavity (thermodynamics and kinetics), coupling such host-guest dynamic equilibrium with characteristic hydrolysis reaction constants. All computed kinetic/thermodynamic data are then combined, obtaining a statistical estimation of reaction acceleration in the host-guest system that is found in optimal agreement with the available experimental trends. This shows how, to understand the key factors controlling accelerations/variations in the reaction under confinement, it is necessary to take into account all dynamic processes that occur as intimately entangled in such host-guest systems. This also provides us with a flexible computational framework, useful to build structure-dynamics-property relationships for a variety of reactive host-guest systems.

Computer Simulation to Rationalize “Rational” Engineering of Glycoside Hydrolases and Glycosyltransferases

Glycoside hydrolases and glycosyltransferases are the main classes of enzymes that synthesize and degrade carbohydrates, molecules essential to life that are a challenge for classical chemistry. As such, considerable efforts have been made to engineer these enzymes and make them pliable to human needs, ranging from directed evolution to rational design, including mechanism engineering. Such endeavors fall short and are unreported in numerous cases, while even success is a necessary but not sufficient proof that the chemical rationale behind the design is correct. Here we review some of the recent work in CAZyme mechanism engineering, showing that computational simulations are instrumental to rationalize experimental data, providing mechanistic insight into how native and engineered CAZymes catalyze chemical reactions. We illustrate this with two recent studies in which (i) a glycoside hydrolase is converted into a glycoside phosphorylase and (ii) substrate specificity of a glycosyltransferase is engineered toward forming O-, N-, or S-glycosidic bonds.

Fungal cellulases: protein engineering and post-translational modifications

Enzymatic degradation of lignocelluloses into fermentable sugars to produce biofuels and other biomaterials is critical for environmentally sustainable development and energy resource supply. However, there are problems in enzymatic cellulose hydrolysis, such as the complex cellulase composition, low degradation efficiency, high production cost, and post-translational modifications (PTMs), all of which are closely related to specific characteristics of cellulases that remain unclear. These problems hinder the practical application of cellulases. Due to the rapid development of computer technology in recent years, computer-aided protein engineering is being widely used, which also brings new opportunities for the development of cellulases. Especially in recent years, a large number of studies have reported on the application of computer-aided protein engineering in the development of cellulases; however, these articles have not been systematically reviewed. This article focused on the aspect of protein engineering and PTMs of fungal cellulases. In this manuscript, the latest literatures and the distribution of potential sites of cellulases for engineering have been systematically summarized, which provide reference for further improvement of cellulase properties. KEY POINTS: •Rational design based on virtual mutagenesis can improve cellulase properties. •Modifying protein side chains and glycans helps obtain superior cellulases. •N-terminal glutamine-pyroglutamate conversion stabilizes fungal cellulases.

Unlocking the secret of lignin-enzyme interactions: Recent advances in developing state-of-the-art analytical techniques

Bioconversion of renewable lignocellulosics to produce liquid fuels and chemicals is one of the most effective ways to solve the problem of fossil resource shortage, energy security, and environmental challenges. Among the many biorefinery pathways, hydrolysis of lignocellulosics to fermentable monosaccharides by cellulase is arguably the most critical step of lignocellulose bioconversion. In the process of enzymatic hydrolysis, the direct physical contact between enzymes and cellulose is an essential prerequisite for the hydrolysis to occur. However, lignin is considered one of the most recalcitrant factors hindering the accessibility of cellulose by binding to cellulase unproductively, which reduces the saccharification rate and yield of sugars. This results in high costs for the saccharification of carbohydrates. The various interactions between enzymes and lignin have been explored from different perspectives in literature, and a basic lignin inhibition mechanism has been proposed. However, the exact interaction between lignin and enzyme as well as the recently reported promotion of some types of lignin on enzymatic hydrolysis is still unclear at the molecular level. Multiple analytical techniques have been developed, and fully unlocking the secret of lignin-enzyme interactions would require a continuous improvement of the currently available analytical techniques. This review summarizes the current commonly used advanced research analytical techniques for investigating the interaction between lignin and enzyme, including quartz crystal microbalance with dissipation (QCM-D), surface plasmon resonance (SPR), attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy, atomic force microscopy (AFM), nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy (FLS), and molecular dynamics (MD) simulations. Interdisciplinary integration of these analytical methods is pursued to provide new insight into the interactions between lignin and enzymes. This review will serve as a resource for future research seeking to develop new methodologies for a better understanding of the basic mechanism of lignin-enzyme binding during the critical hydrolysis process.