A reoptimized GROMOS force field for hexopyranose-based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers

Refinement of the Drude Polarizable Force Field for Hexose Monosaccharides: Capturing Ring Conformational Dynamics with Enhanced Accuracy

We present a revised version of the Drude polarizable carbohydrate force field (FF), focusing on refining the ring and exocyclic torsional parameters for hexopyranose monosaccharides. This refinement addresses the previously observed discrepancies between calculated and experimental NMR 3J coupling values, particularly in describing ring dynamics and exocyclic rotamer populations within major hexose monosaccharides and their anomers. Specifically, α-MAN, β-MAN, α-GLC, β-GLC, α-GAL, β-GAL, α-ALT, β-ALT, α-IDO, and β-IDO were targeted for optimization. The optimization process involved potential energy scans (PES) of the ring and exocyclic dihedral angles computed using quantum mechanical (QM) methods. The target data for the reoptimization included PES of the inner ring dihedrals (C1-C2-C3-C4, C2-C3-C4-C5, C5-O5-C1-C2, C4-C5-O5-C1, O5-C1-C2-C3, C3-C4-C5-O5) and the exocyclic torsions, other than the pseudo ring dihedrals (O1-C1-O5-C5, O2-C2-C1-O5, and O4-C4-C5-O5) and hydroxyl torsions used in the previous parametrization efforts. These parameters, in conjunction with previously developed Drude parameters for hexopyranose monosaccharides, were validated against experimental observations, including NMR data and conformational energetics, in aqueous environments. The resulting polarizable model is shown to be in good agreement with a range of QM data, experimental NMR data, and conformational energetics of monosaccharides in aqueous solutions. This offers a significant improvement of the Drude carbohydrate force field, wherein the refinement enhances the accuracy of accessing the conformational dynamics of carbohydrates in biomolecular simulations.

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties

The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

Outstanding Properties of the Hydration Shell around β-d-Glucose: A Computational Study

Ab initio molecular dynamics (AIMD) simulations have been performed on aqueous solutions of four simple sugars, α-d-glucose, β-d-glucose, α-d-mannose, and α-d-galactose. Hydrogen-bonding (HB) properties, such as the number of donor- and acceptor-type HB-s, and the lengths and strengths of hydrogen bonds between sugar and water molecules, have been determined. Related electronic properties, such as the dipole moments of water molecules and partial charges of the sugar O atoms, have also been calculated. The hydrophilic and hydrophobic shells were characterized by means of spatial distribution functions. β-d-Glucose was found to form the highest number of hydrophilic and the smallest number of hydrophobic connections to neighboring water molecules. The average sugar-water H-bond length was the shortest for β-d-glucose, which suggests that these are the strongest such H-bonds. Furthermore, β-d-glucose appears to stand out in terms of the symmetry properties of both its hydrophilic and hydrophobic hydration shells. In summary, in all aspects considered here, there seems to be a correlation between the distinct characteristics of β-d-glucose reported here and its outstanding solubility in water. Admittedly, our findings represent only some of the important factors that influence the solubility.

In sight the behavior of natural Bletilla striata polysaccharide hydrocolloids by molecular dynamics method

Plant polysaccharides, distinguished by diverse glycosidic bonds and various cyclic sugar units, constitute a subclass of primary metabolites ubiquitously found in nature. Contrary to common understanding, plant polysaccharides typically form hydrocolloids upon dissolution in water, even though both excessively high and low temperatures impede this process. Bletilla striata polysaccharides (BSP), chosen for this kinetic study due to their regular repeating units, help elucidate the relationship between polysaccharide gelation and temperature. It is suggested that elevated temperatures enhance the mobility of BSP molecular chains, resulting in a notable acceleration of hydrogen bond breakage between BSP and water molecules and consequently, compromising the conformational stability of BSPs to some extent. This study unveils the unique relationship between polysaccharide dissolution processes and temperature from a kinetics perspective. Consequently, the conclusion provides a dynamical basis for comprehending the extraction and preparation of natural plant polysaccharide hydrocolloids, pharmaceuticals and related fields.

The counteracting influence of 2-hydroxypropyl substitution and the presence of a guest molecule on the shape and size of the β-cyclodextrin cavity

The aqueous solubility of β-cyclodextrin (β-CD), a cyclic carbohydrate comprising seven α-D-glucose molecules, is enhanced by 2-hydroxypropyl (2-HP) substitution of the hydroxyl groups at the CD rims. Our thorough analysis of the structural and solvation properties with different degrees of 2-hydroxypropyl substitution on β-CD using molecular dynamics simulations reveals that the solubility is enhanced at the cost of the structural distortion of the CD cyclic structure. Substitution at the secondary rim predominantly enhances the favourable interactions between CD and water by decreasing CD-CD hydrogen bonding and promoting CD-water hydrogen bonding. However, the effect of substitution at the primary rim on the CD-water interactions is minimal; the hydrogen bonds between water and the primary hydroxyl group in native CD merely get replaced by those between water and 2-HP, since the substitution makes the primary hydroxyl oxygen (O6 atom) inaccessible to water. In contrast, substitution at the primary rim maintains the structural integrity of CD, while substitution at the secondary rim results in structural distortion due to the disruption of the intramolecular hydrogen bond belt, even leading to cavity closure. Certain strategic substitutions of the primary hydroxyl groups can help in the reduction of structural distortion, depending upon the degree of substitution at the secondary hydroxyl rim. A detailed inspection of the simulation trajectory revealed that the tilting of glucose units with the primary hydroxyl oxygen (O6) pointing inward is the primary driver for cavity closure. Even though the dynamics of glucose tilting can influence the kinetics of host-guest complex formation, once the guest is well incorporated into the cavity, glucose tilting is inhibited and the cavity opens up as in native β-CD.

Optimization of microwave-expanding pretreatment and microwave-assisted extraction of hemicellulose from bagasse cells with the exploration of the extracting mechanism

Hemicellulose is mainly distributed in the tightly packed S2 layer of the plant cell wall and the middle lamella. This rigid microstructure of wood and interactions among hemicellulose, lignin, and cellulose jointly restrict the separation and transformation of hemicellulose in the wood matrix. To address this issue, a method combined with microwave-expanding pretreatment (MEP) and microwave-assisted extraction (MAE) with a NaOH solution was carried out. We found that the MEP could effectively create new pathways for bagasse cells in mass transferring. More specifically, 195 % of the specific surface area (m2/g) with 193 % of the pores (>50 nm) increased after MEP; the SEM images also confirmed that the microstructure of bagasse was modified. MAE could considerably exfoliate hemicellulose from cellulose fiber and accelerate mass transfer. Additionally, we optimized MEP and MAE by using response surface methodology (RSM). The optimal parameters were 370 K, 3.7 min, 1081 W microwave power, and 9.9 wt% NH4HCO3 consumption for the MEP and 1100 W microwave power, 2.5 wt% NaOH concentration, 34.6 min reaction time for MAE, respectively. Moreover, molecular dynamics (MD) simulation suggests that NaOH could significantly lower the work needed to peel off the xylan chain from cellulose nanofibril.

Molecular modeling of the carbohydrate corona formation on a polyvinyl chloride nanoparticle and its impact on the adhesion to lipid bilayers

The pervasive presence of nanoplastics (NPs) in the environment has gained increasing attention due to their accumulation in living organisms. These emerging contaminants inevitably interact with extracellular polymeric substances along respiratory or gastrointestinal tracts, and diverse organic coating on the surface of NPs, known as bio- or eco-corona, is formed. Although its impact on altering the NP properties and potential cell internalization has been extensively examined, studies on its role in NP partitioning in the cell membrane are elusive yet. In this work, molecular dynamics is used to investigate the formation of chitosan (CT) corona centered on a polyvinyl chloride (PVC) nanoparticle and the uptake of the resulting complex onto lipid membranes. Coarse-grained models compatible with the newly developed Martini 3.0 force field are implemented for the two polymers employing the atomistic properties as targets in the parameterization. The reliability of the coarse-grained polymer models is demonstrated by reproducing the structural properties of the PVC melt and of solvated CT strands, as well as by determining the conformation adopted by the latter at the NP surface. Results show that the spontaneous binding of CT chains of high and intermediate protonation degrees led to the formation of soft and hard corona that modulates the interaction of PVC core with model membranes. The structural changes of the corona adsorbed at the lipid-water interface enable a subsequent transfer of the NP to the center of the saturated lipid membranes and a complete or partial transition to a snorkel conformation depending on the hydrophilic/hydrophobic balance in the CT-PVC complex. Overall, the computational investigation of the coarse-grained model system provides implications for understanding how the eco-corona development influences the uptake and implicit toxicology of NPs.

Understanding Self-Assembly and Molecular Packing in Methylcellulose Aqueous Solutions Using Multiscale Modeling and Simulations

We present a multiscale molecular dynamics (MD) simulation study on self-assembly in methylcellulose (MC) aqueous solutions. First, using MD simulations with a new coarse-grained (CG) model of MC chains in implicit water, we establish how the MC chains self-assemble to form fibrils and fibrillar networks and elucidate the MC chains' packing within the assembled fibrils. The CG model for MC is extended from a previously developed model for unsubstituted cellulose and captures the directionality of H-bonding interactions between the -OH groups. The choice and placement of the CG beads within each monomer facilitates explicit modeling of the exact degree and position of methoxy substitutions in the monomers along the MC chain. CG MD simulations show that with increasing hydrophobic effect and/or increasing H-bonding strength, the commercial MC chains (with degree of methoxy substitution, DS, ∼1.8) assemble from a random dispersed configuration into fibrils. The assembled fibrils exhibit consistent fibril diameters regardless of the molecular weight and concentration of MC chains, in agreement with past experiments. Most MC chains' axes are aligned with the fibril axis, and some MC chains exhibit twisted conformations in the fibril. To understand the molecular driving force for the twist, we conduct atomistic simulations of MC chains preassembled in fibrils (without any chain twists) in explicit water at 300 and 348 K. These atomistic simulations also show that at DS = 1.8, MC chains adopt twisted conformations, with these twists being more prominent at higher temperatures, likely as a result of shielding of hydrophobic methyl groups from water. For MC chains with varying DS, at 348 K, atomistic simulations show a nonmonotonic effect of DS on water-monomer contacts. For 0.0 < DS < 0.6, the MC monomers have more water contacts than at DS = 0.0 or DS > 0.6, suggesting that with few methoxy substitutions, the MC chains are effectively hydrophilic, letting the water molecules diffuse into the fibril to participate in H-bonds with the MC chains' remaining -OH groups. At DS > 0.6, the MC monomers become increasingly hydrophobic, as seen by decreasing water contacts around each monomer. We conclude based on the atomistic observations that MC chains with lower degrees of substitutions (DS ≤ 0.6) should exhibit solubility in water over broader temperature ranges than DS ∼ 1.8 chains.

Conserved residues at the family and subfamily levels determine enzyme activity and substrate binding in glycoside hydrolase family 13

Site-directed mutagenesis is a valuable strategy for modifying enzymes, but the lack of understanding of conserved residues regulating glycosidase function hinders enzyme design. We analyzed 1662 enzyme sequences to identify conserved amino acids in maltohexaose-forming amylase at both family and subfamily levels. Several conserved residues at the family level (G37, P45, R52, Y57, D101, V103, H106, G230, R232, D234, E264, H330, D331, and G360) were found, mutations of which resulted in reduced enzyme activity or inactivation. At the subfamily level, several conserved residues (L65, E67, F68, D111, E114, R126, R147, F154, W156, F161, G163, D165, W218H, V342, W345, and F346) were identified, which primarily facilitate substrate binding in the enzyme's active site, as shown by molecular dynamics and kinetic assays. Our findings provide critical insights into conserved residues essential for catalysis and can inform targeted enzyme design in protein engineering.

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling

Molecular dynamics simulations often struggle to obtain sufficient sampling to study complex molecular events due to high energy barriers separating the minima of interest. Multiple enhanced sampling techniques have been developed and improved over the years to tackle this issue. Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that works by adding a biasing potential, lifting the energy landscape up, and decreasing the height of its barriers. GaMD allows one to increase the sampling of events of interest without the need of a priori knowledge of the system or the relevant coordinates. All required acceleration parameters can be obtained from a previous search run. Upon its development, several improvements for the methodology have been proposed, among them selective GaMD in which the boosting potential is selectively applied to the region of interest. There are currently four selective GaMD methods that have shown promising results. However, all of these methods are constrained on the number, location, and scenarios in which this selective boosting potential can be applied to ligands, peptides, or protein-protein interactions. In this work, we showcase a GROMOS implementation of the GaMD methodology with a fully flexible selective GaMD approach that allows the user to define, in a straightforward way, multiple boosting potentials for as many regions as desired. We show and analyze the advantages of this flexible selective approach on two previously used test systems, the alanine dipeptide and the chignolin peptide, and extend these examples to study its applicability and potential to study conformational changes of glycans and glycosylated proteins.

Distribution of Aromatic Amino Acid Residues in Substrate-Binding Regions Modulates Substrate Specificity of Microbial Debranching Enzymes

Debranching enzymes (DBEs) directly hydrolyze α-1,6-glucosidic linkages in glycogen, starch, and related polysaccharides, making them important in the starch processing industry. However, the ambiguous substrate specificity usually restricts synergistic catalysis with other amylases for improving starch utilization. Herein, a glycogen-debranching enzyme from Saccharolobus solfataricus (SsGDE) and two isoamylases from Pseudomonas amyloderamosa (PaISO) and Chlamydomonas reinhardtii (CrISO) were used to investigate the molecular mechanism of substrate specificity. Along with the structure-based computational analysis, the aromatic residues in the substrate-binding region of DBEs played an important role in binding substrates. The aromatic residues in SsGDE appeared clustered, contributing to a small substrate-binding region. In contrast, the aromatic residues in isoamylase were distributed dispersedly, forming a large active site. The distinct characteristics of substrate-binding regions in SsGDE and isoamylase might explain their substrate preferences for maltodextrin and amylopectin, respectively. By modulating the substrate-binding region of SsGDE, variants Y323F and V375F were obtained with significantly enhanced activities, and the activities of Y323F and V375F increased by 30 and 60% for amylopectin, and 20 and 23% for DE4 maltodextrin, respectively. This study revealed the molecular mechanisms underlying the substrate specificity for SsGDE and isoamylases, providing a route for engineering enzymes to achieve higher catalytic performance.

Understanding Nanocellulose-Water Interactions: Turning a Detriment into an Asset

Modern technology has enabled the isolation of nanocellulose from plant-based fibers, and the current trend focuses on utilizing nanocellulose in a broad range of sustainable materials applications. Water is generally seen as a detrimental component when in contact with nanocellulose-based materials, just like it is harmful for traditional cellulosic materials such as paper or cardboard. However, water is an integral component in plants, and many applications of nanocellulose already accept the presence of water or make use of it. This review gives a comprehensive account of nanocellulose-water interactions and their repercussions in all key areas of contemporary research: fundamental physical chemistry, chemical modification of nanocellulose, materials applications, and analytical methods to map the water interactions and the effect of water on a nanocellulose matrix.

The thermodynamics of enhanced dope stability of cellulose solution in NaOH solution by urea

The addition of urea in pre-cooled alkali aqueous solution is known to improve the dope stability of cellulose solution. However, its thermodynamic mechanism at a molecular level is not fully understood yet. By using molecular dynamics simulation of an aqueous NaOH/urea/cellulose system using an empirical force field, we found that urea was concentrated in the first solvation shell of the cellulose chain stabilized mainly by dispersion interaction. When adding a glucan chain into the solution, the total solvent entropy reduction is smaller if urea is present. Each urea molecule expelled an average of 2.3 water molecules away from the cellulose surface, releasing water entropy that over-compensates the entropy loss of urea and thus maximizing the total entropy. Scaling the Lennard-Jones parameter and atomistic partial charge of urea revealed that direct urea/cellulose interaction was also driven by dispersion energy. The mixing of urea solution and cellulose solution in the presence or absence of NaOH are both exothermic even after correcting for the contribution from dilution.

Impact of Polarization on the Ring Puckering Dynamics of Hexose Monosaccharides

Analysis of crystal structures of hexose monosaccharides α-d-mannose (α-MAN), β-d-mannose (β-MAN), α-d-glucose (α-GLC), β-d-glucose (β-GLC), α-d-galactose (α-GAL), β-d-galactose (β-GAL), α-d-altrose (α-ALT), β-d-altrose (β-ALT), α-d-idose (α-IDO), and β-d-idose (β-IDO) reveals that the monosaccharide ring adopts multiple ring conformations. These ring conformations can be broadly classified as chair, half-chair, envelope, boat, and skew-boat conformations. The ability of the monosaccharide ring to adopt multiple conformations has been closely tied with their bioactivity. However, it has been difficult to capture the dynamic information of these conformations from experimental studies. Even from simulations, capturing these different conformations is challenging because of the energy barriers involved in the transitions between the stable ⁴C₁ and ¹C₄ chair forms. In this study, we analyze the influence of the polarizable force field on the ring dynamics of five major types of unsubstituted aldohexoses─glucose, mannose, galactose, altrose, and idose─and their anomers. We simulate microsecond trajectories to capture the influence of the CHARMM36 additive and polarizable carbohydrate force fields on the ring dynamics. The microsecond trajectories allow us to comment on the issues associated with equilibrium molecular dynamics simulations. Further, we use the extended system adaptive biasing force (eABF) method to compare the conformational sampling efficiencies of the additive and polarizable force fields. Our studies reveal that inclusion of polarization enhances the sampling of ring conformations and lowers the energy barriers between the ⁴C₁ and ¹C₄ conformations. Overall, the CHARMM36 additive force field is observed to be rigid and favor the ⁴C₁ conformations. Although the inclusion of polarizability results in enhancing ring flexibility, we observe sampling that does not agree with experimental results, warranting a revision of the polarizable Drude parameters.

Atomic-Resolution Experimental Structural Biology and Molecular Dynamics Simulations of Hyaluronan and Its Complexes

This review summarizes the atomic-resolution structural biology of hyaluronan and its complexes available in the Protein Data Bank, as well as published studies of atomic-resolution explicit-solvent molecular dynamics simulations on these and other hyaluronan and hyaluronan-containing systems. Advances in accurate molecular mechanics force fields, simulation methods and software, and computer hardware have supported a recent flourish in such simulations, such that the simulation publications now outnumber the structural biology publications by an order of magnitude. In addition to supplementing the experimental structural biology with computed dynamic and thermodynamic information, the molecular dynamics studies provide a wealth of atomic-resolution information on hyaluronan-containing systems for which there is no atomic-resolution structural biology either available or possible. Examples of these summarized in this review include hyaluronan pairing with other hyaluronan molecules and glycosaminoglycans, with ions, with proteins and peptides, with lipids, and with drugs and drug-like molecules. Despite limitations imposed by present-day computing resources on system size and simulation timescale, atomic-resolution explicit-solvent molecular dynamics simulations have been able to contribute significant insight into hyaluronan's flexibility and capacity for intra- and intermolecular non-covalent interactions.

Automation of AMOEBA polarizable force field for small molecules: Poltype 2

A next-generation protocol (Poltype 2) has been developed which automatically generates AMOEBA polarizable force field parameters for small molecules. Both features and computational efficiency have been drastically improved. Notable advances include improved database transferability using SMILES, robust torsion fitting, non-aromatic ring torsion parameterization, coupled torsion-torsion parameterization, Van der Waals parameter refinement using ab initio dimer data and an intelligent fragmentation scheme that produces parameters with dramatically reduced ab initio computational cost. Additional improvements include better local frame assignment for atomic multipoles, automated formal charge assignment, Zwitterion detection, smart memory resource defaults, parallelized fragment job submission, incorporation of Psi4 quantum package, ab initio error handling, ionization state enumeration, hydration free energy prediction and binding free energy prediction. For validation, we have applied Poltype 2 to ~1000 FDA approved drug molecules from DrugBank. The ab initio molecular dipole moments and electrostatic potential values were compared with Poltype 2 derived AMOEBA counterparts. Parameters were further substantiated by calculating hydration free energy (HFE) on 40 small organic molecules and were compared with experimental data, resulting in an RMSE error of 0.59 kcal/mol. The torsion database has expanded to include 3543 fragments derived from FDA approved drugs. Poltype 2 provides a convenient utility for applications including binding free energy prediction for computational drug discovery. Further improvement will focus on automated parameter refinement by experimental liquid properties, expansion of the Van der Waals parameter database and automated parametrization of modified bio-fragments such as amino and nucleic acids.

Extending the Martini 3 Coarse-Grained Force Field to Carbohydrates

Carbohydrates play an essential role in a large number of chemical and biochemical processes. High structural diversity and conformational heterogeneity make it problematic to link their measurable properties to molecular features. Molecular dynamics simulations carried out at the level of classical force fields are routinely applied to study the complex processes occurring in carbohydrate-containing systems, while the usefulness of such simulations relies on the accuracy of the underlying theoretical model. In this article, we present the coarse-grained force field dedicated to glucopyranose-based carbohydrates and compatible with the recent version of the Martini force field (v. 3.0). The parameterization was based on optimizing bonded and nonbonded parameters with a reference to the all-atom simulation results and the experimental data. Application of the newly developed coarse-grained carbohydrate model to oligosaccharides curdlan and cellulose displays spontaneous formation of aggregates of experimentally identified features. In contact with other biomolecules, the model is capable of recovering the protective effect of glucose monosaccharides on a lipid bilayer and correctly identifying the binding pockets in carbohydrate-binding proteins. The features of the newly proposed model make it an excellent candidate for further extensions, aimed at modeling more complex, functionalized, and biologically relevant carbohydrates.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Conformational preferences of triantennary and tetraantennary hybrid N-glycans in aqueous solution: Insights from 20 μs long atomistic molecular dynamic simulations

In the current study, we have investigated the conformational dynamics of a triantennary (N-glycan1) and tetraantennary (N-glycan2) hybrid N-glycans found on the surface of the HIV glycoprotein using 20 μs long all-atom molecular dynamics (MD) simulations. The main objective of the present study is to elucidate the influence of adding a complex branch on the overall glycan structural dynamics. Our investigation suggests that the average RMSD value increases when a complex branch is added to N-glycan1. However, the RMSD distribution is relatively wider in the case of N-glycan1 compared to N-glycan2, which indicates that multiple complex branches restrict the conformational variability of glycans. A similar observation is obtained from the principal component analysis of both glycans. All the puckering states (⁴C₁ to ¹C₄) of each monosaccharide except mannose are sampled in our simulations, although the ⁴C₁ chair form is energetically more favorable than ¹C₄. In N-glycan1, the 1-6 linkage in the mannose branch [Man(9)-α(1-6)-Man(5)] stays in the gauche-gauche cluster, whereas it moves towards trans-gauche in N-glycan2. For both glycans, mannose branches are more flexible than the complex branches, and adding a complex branch does not influence the dynamics of the mannose branches. We have noticed that the end-to-end distance of the complex branch shortens by ∼ 10 Å in the presence of another complex branch. This suggests that in the presence of an additional complex branch, the other complex branch adopts a close folded structure. All these conformational changes involve the selective formation of inter-residue and water-mediated hydrogen-bond networks.

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins

The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.

O-Methylation in Carbohydrates: An NMR and MD Simulation Study with Application to Methylcellulose

Methylated carbohydrates are important from both biological and technical perspectives. Specifically, methylcellulose is an interesting cellulose derivative that has applications in foods, materials, cosmetics, and many other fields. While the molecular dynamics simulation technique has the potential for both advancing the fundamental understanding of this polymer and aiding in the development of specific applications, a general drawback is the lack of experimentally validated interaction potentials for the methylated moieties. In the present study, simulations using the GROMOS 56 carbohydrate force field are compared to NMR spin–spin coupling constants related to the conformation of the exocyclic torsion angle ω in d-glucopyranose and derivatives containing a 6-O-methyl substituent and a ¹³C-isotopologue thereof. A ³J_CC Karplus-type relationship is proposed for the C5–C6–O6–C_Me torsion angle. Moreover, solvation free energies are compared to experimental data for small model compounds. Alkylation in the form of 6-O-methylation affects exocyclic torsion only marginally. Computed solvation free energies between nonmethylated and methylated molecules were internally consistent, which validates the application of these interaction potentials for more specialized purposes.

Polymorphism of chitosan-based networks stabilized by phytate investigated by molecular dynamics simulations

Chitosan can associate in the presence of polyphosphates into insoluble hydrogels capable of drug encapsulation and safe and efficient release. On the one hand, chitosan hydrogels were synthesized using the phytate anion as a crosslinking agent and were characterized by employing dynamic light scattering (DLS) and Fourier transform infrared spectroscopy (FTIR). On the other hand, an effective chitosan-phytate model with atomistic details was created to examine the underlying physical crosslinking pattern, and the structure and dynamics of the chitosan-phytate complex were systematically investigated by using molecular dynamics (MD) simulations. To harbor the crosslinker potential for obtaining chitosan-based hydrogels, the impact of the phytate concentration and the functional groups of the chitosan on the reticulation process was addressed. The phytate association was determined by the phosphates' capacity for H-bonding to the amine and hydroxyl groups belonging to two consecutive glucosidic units. The physical crosslinking pattern was determined by the number of chitosan chains bound by one phytate anion and the phytate orientation relative to the glucopyranose neighbors. Cross-linking of two up to six chitosan chains mediated by a phytate anion represented favorable states, and the number distribution of cross-linked chains depended on the phytate concentration. The circular distribution of the cross-linkable phosphates regulated the nearly isotropic orientation of the chitosan chains and phytate at the junction, and the variety of topological crosslinking demonstrated the phytate ion's potential for developing chitosan-based hydrogels with improved structural attributes.

Tautomers of N-acetyl-d-allosamine: an NMR and computational chemistry study

d-Allosamine is a rare sugar in Nature but its pyranoid form has been found α-linked in the core region of the lipopolysaccharide from the Gram-negative bacterium Porphyromonas gingivalis and in the chitanase inhibitor allosamidin, then β-linked and N-acetylated. In water solution the monosaccharide N-acetyl-d-allosamine (d-AllNAc) shows a significant presence of four tautomers arising from pyranoid and furanoid ring forms and anomeric configurations. The furanoid ring forms both showed ³J_H1,H2≈ 4.85 Hz and to differentiate the anomeric configurations a series of chemical shift anisotropy/dipole-dipole cross-correlated relaxation NMR experiments was performed in which the α-anomeric form showed notable different relaxation rates for its components of the H1 doublet, thereby making it possible to elucidate the anomeric configuration of each of the furanoses. The conformational preferences of the different forms of d-AllNAc were investigated by ³J_HH, ²J_CH and ³J_CH coupling constants from NMR experiments, molecular dynamics simulations and density functional theory calculations. The pyranose form resides in the ⁴C₁ conformation and the furanose ring form has the majority of its conformers located on the South-East region of the pseudorotation wheel, with a small population in the Northern hemisphere. The tautomeric equilibrium was quite sensitive to changes in temperature, where the β-anomer of the pyranoid ring form decreased upon a temperature increase while the other forms increased.

ATP hydrolysis and nucleotide exit enhance maltose translocation in the MalFGK2E importer

ATP binding cassette (ABC) transporters employ ATP hydrolysis to harness substrate translocation across membranes. The Escherichia coli MalFGK₂E maltose importer is an example of a type I ABC importer and a model system for this class of ABC transporters. The MalFGK₂E importer is responsible for the intake of malto-oligossacharides in E.coli. Despite being extensively studied, little is known about the effect of ATP hydrolysis and nucleotide exit on substrate transport. In this work, we studied this phenomenon using extensive molecular dynamics simulations (MD) along with potential of mean force calculations of maltose transport across the pore, in the pre-hydrolysis, post-hydrolysis and nucleotide-free states. We concluded that ATP hydrolysis and nucleotide exit trigger conformational changes that result in the decrease of energetic barriers to maltose translocation towards the cytoplasm, with a concomitant increase of the energy barrier in the periplasmic side of the pore, contributing for the irreversibility of the process. We also identified key residues that aid in positioning and orientation of maltose, as well as a novel binding pocket for maltose in MalG. Additionally, ATP hydrolysis leads to conformations similar to the nucleotide-free state. This study shows the contribution of ATP hydrolysis and nucleotide exit in the transport cycle, shedding light on ABC type I importer mechanisms.

Molecular Structure of Cefuroxime Axetil Complexes with α-, β-, γ-, and 2-Hydroxypropyl-β-Cyclodextrins: Molecular Simulations and Raman Spectroscopic and Imaging Studies

The formation of cefuroxime axetil+cyclodextrin (CA+CD) complexes increases the aqueous solubility of CA, improves its physico-chemical properties, and facilitates a biomembrane-mediated drug delivery process. In CD-based tablet formulations, it is crucial to investigate the molecular details of complexes in final pharmaceutical preparation. In this study, Raman spectroscopy and mapping were applied for the detection and identification of chemical groups involved in α-, β-, γ-, and 2-hydroxypropyl-β-CD (2-HP- β-CD)+CA complexation process. The experimental studies have been complemented by molecular dynamics-based investigations, providing additional molecular details of CA+CD interactions. It has been demonstrated that CA forms the guest–host type inclusion complexes with all studied CDs; however, the nature of the interactions is slightly different. It seems that both α- and β-CD interact with furanyl and methoxy moieties of CA, γ-CD forms a more diverse pattern of interactions with CA, which are not observed in other CDs, whereas 2HP-β-CD binds CA with the contribution of hydrogen bonding. Apart from supporting this interpretation of the experimental data, molecular dynamics simulations allowed for ordering the CA+CD binding affinities. The obtained results proved that the molecular details of the host–guest complexation can be successfully predicted from the combination of Raman spectroscopy and molecular modeling.

Comparison of Carbohydrate Force Fields in Molecular Dynamics Simulations of Protein-Carbohydrate Complexes

In this paper, we present the results of molecular dynamics simulations aimed at critical comparison of classical, biomolecular force fields (FFs) in the context of their capabilities to describe the structural and thermodynamic features of carbohydrate-protein interactions. We have considered the three main families of FFs (CHARMM, GROMOS, and GLYCAM/AMBER) by applying them to investigate the seven different carbohydrate-protein complexes. The results indicate that although the qualitative pattern of several structural descriptors (intermolecular hydrogen bonding, ligand dynamic location, etc.) is conserved among the compared FFs, there also exists a number of significant divergences (mainly the patterns of contacts between particular amino acid residues and bound carbohydrate). The carbohydrate-protein unbinding free energies also vary from one FF to another, displaying diversified trends in deviations from the experimental data. The magnitude of those deviations is not negligible and indicates the need for refinement in the currently existing combinations of carbohydrate- and protein-dedicated biomolecular force fields. In spite of the lack of explicit functional terms responsible for the corresponding intermolecular forces, all tested FFs are capable of adequately reproducing the CH-π interactions, crucial for carbohydrate-protein binding.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Molecular View into the Cyclodextrin Cavity: Structure and Hydration

We find, through atomistic molecular dynamics simulation of native cyclodextrins (CDs) in water, that although the outer surface of a CD appears like a truncated cone, the inner cavity resembles a conical hourglass because of the inward protrusion of the glycosidic oxygens. Furthermore, the conformations of the constituent α-glucose molecules are found to differ significantly from a free monomeric α-glucose molecule. This is the first computational study that maps the conformational change to the preferential hydrogen bond donating capacity of one of the secondary hydroxyl groups of CD, in consensus with an NMR experiment. We have developed a simple and novel geometry-based technique to identify water molecules occupying the nonspherical CD cavity, and the computed water occupancies are in close agreement with the experimental and density functional theory studies. Our analysis reveals that a water molecule in CD cavity loses out about two hydrogen bonds and remains energetically frustrated but possesses higher orientational degree of freedom compared to bulk water. In the context of CD-drug complexation, these imply a nonclassical, that is, enthalpically driven hydrophobic association of a drug in CD cavity.

Dissipative particle dynamics model of homogalacturonan based on molecular dynamics simulations

In this study we present an alternative dissipative particle dynamics (DPD) parametrization strategy based on data extracted from the united-atom molecular simulations. The model of the homogalacturonan was designed to test the ability of the formation of large-scale structures via hydrogen bonding in water. The extraction of coarse-grained parameters from atomistic molecular dynamics was achieved by means of the proposed molecule aggregation algorithm based on an iterative nearest neighbour search. A novel approach to a time-scale calibration scheme based on matching the average velocities of coarse-grained particles enabled the DPD forcefield to reproduce essential structural features of homogalacturonan molecular chains. The successful application of the proposed parametrization method allowed for the reproduction of the shapes of radial distribution functions, particle velocities and diffusivity of the atomistic molecular dynamics model using DPD force field. The structure of polygalacturonic acid molecules was mapped into the DPD force field by means of the distance and angular bond characteristics, which closely matched the MD results. The resulting DPD trajectories showed that randomly dispersed homogalacturonan chains had a tendency to aggregate into highly organized 3D structures. The final structure resembled a three-dimensional network created by tightly associated homogalacturonan chains organized into thick fibres.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

Temperature-Dependent Blue Shifting of O-H Stretching Frequencies in Crystalline Cellulose Explained

Increasing the temperature of a chemical system generally causes covalent bonds to lengthen and weaken, often the first step in initiating chemical reactions. However, for some hydrogen-bonded systems, infrared (IR) spectroscopy measurements reveal that covalent O-H bonds actually strengthen and therefore shorten when heated. In 1957, Finch and Lippincott proposed a simple one-dimensional (1D) model to explain this effect, in which thermal excitation of intermolecular stretching modes leads to lengthening and weakening of intermolecular O-H···O hydrogen bonds, thereby indirectly strengthening the associated covalent O-H bonds. Taking cellulose (an infinitely repeating polymer of d-glucose) as an example, we use molecular dynamics modeling to show that the same mechanism is responsible for temperature-dependent blue shifting of O-H stretching bands in IR spectra of carbohydrate biopolymers, except that interchain hydrogen bonds are weakened by thermal excitation of chain-separation modes, while intrachain hydrogen bonds are weakened by thermally induced changes in ring puckering and orientation of ring substituents but not reorientation of glucose units relative to one another or overall twisting of the cellulose chains.

Conformational properties of inulin, levan and arabinan studied by molecular dynamics simulations

Inulin, levan and arabinan are the polysaccharides that consist of exclusively furanose units. To date, their conformational features studied at the molecular scale have remained largely unexplained. To tackle this issue, we have performed a series of explicit-solvent molecular dynamics simulations, carried out within the furanose-dedicated force field. None of the polysaccharides exhibits a single, dominating structure type. Instead, they create a large number of separated conformational states originating from the intensive rotation around the φ and ω glycosidic angles. 2₁-helices are the preferential conformational forms for all compounds but they appear only locally, at the length of several consecutive residues. The flexibility of all three furanose-based polysaccharides is much greater in relation to the (1-4)-linked pyranose polysaccharides and is comparable to that of (1-6)-linked pyranoses. The dynamic geometries of both furanose rings and glycosidic linkages are nearly unchanged independently if considering them at the level of mono-, di- or polysaccharides.

Electron microdiffraction reveals the nanoscale twist geometry of cellulose nanocrystals

Nanocellulose consisting of crystalline cellulose nanoparticles has high potential to serve as a building block for bio-based functional materials. The intrinsic chirality of cellulose provides them with high added values such as optical properties and chiral induction ability. At the nanoscale, this chirality is connected to the right-handed longitudinal twisting of these fibrous crystallites. However, this nanoscale fibrillar twist has been a matter of debate due to contradictory data between ultrastructural observations and molecular simulations and so far, the exact twist geometry has not been elucidated. Here, an electron microdiffraction (μED) analysis under cryogenic conditions reveals the continuous twisting of cellulose nanocrystals (CNCs) in aqueous suspension. This intrinsic regular twist is drastically modified to a discontinuous sharp twist when the CNCs are dried on a flat surface. The present μED-based analysis at the single nanoparticle level allows the establishment of the quantitative structure-property relationship of various solid and colloidal nanocellulose systems.

Construction of cellulose-based Pickering stabilizer as a novel interfacial antioxidant: A bioinspired oxygen protection strategy

Oxygen protection/isolation is imperative to prevent the lipid oxidation since oxygen molecule is an ultimate quencher in photon conversion process. Inspired by the structural buildup of seeds from oil crops, a sustainable solid particle stabilizer with novel antioxidant activity was prepared by using cellulose and polyphenol. In this work, bacterial cellulose (BC) nanofibrils modified by tea polyphenols (TPs) was prepared and used as Pickering emulsifier for the O/W emulsion. BC nanofibirls exhibited excellent adsorption capacity up to 55 μg/mg, and the adsorption kinetics between BC and TPs were further investigated. After modification, the interfacial diffusion rate constant of BC was significantly increased to from 0.43 to 1.21 mN m^-1 s^-0.5. Moreover, the obtained O/W interfacial modulus of the dilatational elasticity was increased from 58 to 130 mN/m. Furthermore, the emulsions exhibited excellent free-radical scavenging activity at oil-water interface, suggesting a potential application in usage to extend the lifespan of the food containing polyunsaturated fats.

Drude Polarizable Force Field Parametrization of Carboxylate and N-Acetyl Amine Carbohydrate Derivatives

In this work, we report the development of Drude polarizable force field parameters for the carboxylate and N-acetyl amine derivatives, extending the functionality of the existing Drude polarizable carbohydrate force field. The force field parameters have been developed in a hierarchical manner, reproducing the quantum mechanical gas-phase properties of small model compounds representing the key functional group in the carbohydrate derivatives, including optimization of the electrostatic and bonded parameters. The optimized parameters were then used to generate the models for carboxylate and N-acetyl amine carbohydrate derivatives. The transferred parameters were further tested and optimized to reproduce crystal geometries and J-coupling data from nuclear magnetic resonance experiments. The parameter development resulted in the incorporation of d-glucuronate, l-iduronate, N-acetyl-d-glucosamine (GlcNAc), and N-acetyl-d-galactosamine (GalNAc) sugars into the Drude polarizable force field. The parameters developed in this study were then applied to study the conformational properties of glycosaminoglycan polymer hyaluronan, composed of d-glucuronate and N-acetyl-d-glucosamine, in aqueous solution. Upon comparing the results from the additive and polarizable simulations, it was found that the inclusion of polarization improved the description of the electrostatic interactions observed in hyaluronan, resulting in enhanced conformational flexibility. The developed Drude polarizable force field parameters in conjunction with the remainder of the Drude polarizable force field parameters can be used for future studies involving carbohydrates and their conjugates in complex, heterogeneous systems.

Relation between the NMR data and the pseudorotational free-energy profile for oxolane

The conformation of five-membered furanose rings is a crucial issue for the structural analysis of many biologically-relevant molecules, including DNA and RNA. Oxolane can be treated as a prototypical furanose, composed only of saturated unsubstituted ring. In spite of its structural simplicity, providing the accurate quantitative description of the oxolane conformational features remains a great challenge for both the experimental and theoretical techniques. Here we show the method of recovering the free-energy profiles describing the conformational equilibrium in the oxolane ring (i.e. pseudorotation) based on the experimentally-inferred NMR data ([Formula: see text] coupling constants). The results remain in agreement with the quantum-mechanical-based molecular dynamics simulations and emphasize the large contributions of all ring conformers, even those located at the free-energy barriers. This includes the significant populations of limiting 3T2/2T3 and OE/EO shapes. Our findings provide another example of a poor applicability of the two-state model, which is routinely applied to analyze the NMR data in terms of population of different ring conformers.

The systematic influence of solvent on the conformational features of furanosides

The endo- and exo-anomeric effects are the two most recognizable stereoelectronic effects exhibited by carbohydrates. Their presence relies on the interactions between ring substituent(s) and ring oxygen atoms. Here, we report the finding of a new effect that partially controls the conformational properties of furanose rings and can be ascribed to the influence of the solvent on the electronic structure of a molecule. In contrast to anomeric effects, it is not dependent on either presence or absence of ring substituents. Its origins lie in a solvent-induced flux of atomic charges that involves atoms forming the furanose ring. This systematically changes the energy of the whole molecular system and alters the ring-distortion free energies by ∼2.5-6.5 kJ mol-1, favoring the geometries close to the twist 3T2/2T3 conformers and disfavoring the envelope OE/EO-like shapes. This intriguing effect has never been reported before, although it is expected to exist in all furanose rings. Along with more recognized stereoelectronic effects, this phenomenon contributes to a wide applicability of the two-state (north vs. south) model of pseudorotation in furanosides and, in the case of extremely flexible furanose rings, may change the preferred conformation type in comparison with the gas-phase-oriented predictions.

Mechanism and theory of d-glucopyranose homogeneous acid catalysis in the aqueous solution phase

A detailed systematic theoretical study of the mechanism of the homogeneous Brønsted-acid catalysis of d-glucose in aqueous solution phase ("acid hydrolysis") is reported. G4MP2 with the SMD solvation model at B3LYP/6-31G(2df,p) are employed to compute the free energies of the molecular and ionic species pertaining to the isomerization, protonation, hydrogen cation transfer and decomposition processes of d-glucopyranose in aqueous solution phase. This information is used to hypothesise a reaction mechanism that is of improved accuracy and completeness from the existing art. It is found that rotation of the d-glucose alkyl carbon-carbon bond is a facile process and is very important to the subsequent catalytic mechanism. This rotation produces two rotameric isomers which are of notably different thermodynamic stability and reactivity, even with regard to the products of this acid catalysis. As a low energy process (ΔG‡ = ∼3.8-6.7 kcal mol-1), the alkyl carbon-carbon bond may rotate toward the hydroxyl group at the adjacent "4" position reducing the energy required to protonate that position by 3.0-7.2 kcal mol-1 (or 15-30%). The combination of two rotomeric isomers with the six structural isomers owing to the oxygen atoms, means that protonated d-glucose cations embark on a complex competition of interconversion and decomposition that is both thermodynamically and kinetically influenced. The calculations support the hypothesis that the acid-catalysed hydrolysis of d-glucose may yield a number of platform chemicals that have not previously been suggested. These include the prospect of three isomers of 5-hydroxymethylfurfural (HMF); 5-(hydroxymethyl)furan-2-carbaldehyde, 5-(hydroxymethyl)furan-3-carbaldehyde and 5-(hydroxymethyl)furan-4-carbaldehyde. Vibrational spectra of these HMF isomers are also computed and compared to experimentally determined infrared spectra of "humins". On this basis, it is cautiously speculated that the alternative HMF isomers, may be monomeric constituent of the polymeric "humins".

Effects of Tyr555 and Trp678 on the processivity of cellobiohydrolase A from Ruminiclostridium thermocellum: A simulation study

Cellobiohydrolase A from Ruminiclostridium thermocellum (Cbh9A) is a processive exoglucanase from family 9 and is an important cellobiohydrolase that hydrolyzes cello-oligosaccharide into cellobiose. Residues Tyr555 and Trp678 considerably affect catalytic activity, but their mechanisms are still unknown. To investigate how the Tyr555 and Trp678 affect the processivity of Cbh9A, conventional molecular dynamics, steered molecular dynamics, and free energy calculation were performed to simulate the processive process of wild type (WT)-Cbh9A, Y555S mutant, and W678G mutant. Analysis of simulation results suggests that the binding free energies between the substrate and WT-Cbh9A are lower than those of Y555S and W678G mutants. The pull forces and energy barrier in Y555S and W678G mutants also reduced significantly during the steered molecular dynamics (SMD) simulation compared with that of the WT-Cbh9A. And the potential mean force calculations showed that the pulling energy barrier of Y555S and W678G mutants is much lower than that of WT-Cbh9A.

Sugar-based bactericides targeting phosphatidylethanolamine-enriched membranes

Anthrax is an infectious disease caused by Bacillus anthracis, a bioterrorism agent that develops resistance to clinically used antibiotics. Therefore, alternative mechanisms of action remain a challenge. Herein, we disclose deoxy glycosides responsible for specific carbohydrate-phospholipid interactions, causing phosphatidylethanolamine lamellar-to-inverted hexagonal phase transition and acting over B. anthracis and Bacillus cereus as potent and selective bactericides. Biological studies of the synthesized compound series differing in the anomeric atom, glycone configuration and deoxygenation pattern show that the latter is indeed a key modulator of efficacy and selectivity. Biomolecular simulations show no tendency to pore formation, whereas differential metabolomics and genomics rule out proteins as targets. Complete bacteria cell death in 10 min and cellular envelope disruption corroborate an effect over lipid polymorphism. Biophysical approaches show monolayer and bilayer reorganization with fast and high permeabilizing activity toward phosphatidylethanolamine membranes. Absence of bacterial resistance further supports this mechanism, triggering innovation on membrane-targeting antimicrobials.

Bacillus anthracis causes the infectious disease anthrax. Here, the authors synthesized deoxy glycosides that are effective against B. anthracis and related bacteria and found that these amphiphilic compounds kill bacteria via an unusual mechanism of action.