Revision of the GROMOS 56A6(CARBO) force field: Improving the description of ring-conformational equilibria in hexopyranose-based carbohydrates chains

Structure and Dynamics of Water in Polysaccharide (Alginate) Solutions and Gels Explained by the Core-Shell Model

In both biological and engineered systems, polysaccharides offer a means of establishing structural stiffness without altering the availability of water. Notable examples include the extracellular matrix of prokaryotes and eukaryotes, artificial skin grafts, drug delivery materials, and gels for water harvesting. Proper design and modeling of these systems require detailed understanding of the behavior of water confined in pores narrower than about 1 nm. We use molecular dynamics simulations to investigate the properties of water in solutions and gels of the polysaccharide alginate as a function of the water content and polymer cross-linking. We find that a detailed understanding of the nanoscale dynamics of water in alginate solutions and gels requires consideration of the discrete nature of water. However, we also find that the trends in tortuosity, permeability, dielectric constant, and shear viscosity can be adequately represented using the "core-shell" conceptual model that considers the confined fluid as a continuum.

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.

Investigation of the formation mechanism of the pepper starch-piperine complex

In this study, a new carrier for loading piperine was prepared using pepper starch, and its interaction mechanism was investigated. The porous pepper starch-piperine complex (PPS-PIP) showed higher loading efficiency (76.15 %) compared to the porous corn starch-piperine complex (PCS-PIP (52.34 %)). This may be ascribed to the hemispherical shell structure of porous pepper starch (PPS) compared to the porous structure of porous corn starch (PCS) based on the SEM result. PPS-PIP had smaller particle size (10.53 μm), higher relative crystallinity (38.95 %), and better thermal stability (87.45 °C) than PCS-PIP (17.37 μm, 32.17 %, 74.35 °C). Fourier transform infrared spectroscopy (FTIR) results implied that piperine not only forms a complex with amylose but may also be physically present in porous starch. This was demonstrated by the short-range order and X-ray type. Molecular dynamics simulations confirmed that hydrogen bonding is the primary interaction between amylose and piperine. Besides the formation of the amylose-piperine complex, some of the piperine is also present in physical form.

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.

Exploring Ring Conformation in Uronate Monosaccharides: Insights from Ab Initio Calculations and Classical Molecular Dynamics Simulations

The study focuses on the conformational properties of biologically relevant monosaccharides belonging to the group of uronates: α-l-iduronate, O2-sulfated-α-l-iduronate, and O2-sulfated-α-l-guluronate, either unfunctionalized or O1-methylated. We applied the previously proposed two-step methodology, combining classical MD simulations and subsequent ab initio (QM) calculations, performed on a rationally subsampled set of molecular configurations. We found that, regardless of the number of molecular configurations considered, the level of theory, and the weighting scheme applied, none of the QM approaches is capable of predicting the correct conformational equilibrium of sulfated iduronates as long as the tight counterion binding is not considered. Multicenter, ring-shape-specific binding of either Na+ or Ca2+ ions drastically shifts the conformational equilibrium of the pyranose ring in sulfated iduronates toward 1C4 but does not significantly affect the conformation of non-sulfated compounds. A similar shift is observed upon the protonation of carboxyl groups in all iduronates. In addition, we report a set of average J-coupling constant values related to vicinal protons bound to the pyranose ring of iduronates and corresponding to each of the three main groups of ring conformers, i.e., 4C1, B/S (boat/skew boat), and 1C4. In combination with the conformational energies or with the experimental data, these values allowed the relative proportions of the ring conformers to be estimated and the Karplus-type equations linking the 3JHH-coupling constants to the torsion angles within the pyranose ring to be refined.

Structure-Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities─and their impacts on environmental processes─requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca-Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ∼0.6 nm between alginate strands.

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.

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.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

Carbohydrate Force Fields: The Role of Small Partial Atomic Charges in Preventing Conformational Collapse

Although the quality of current additive all-atom force fields for carbohydrates has been demonstrated in many applications, occasional significant differences reported for the hydrodynamic behavior of specific polysaccharides modeled with different force fields is a cause for concern. In particular, irreversible conformational collapse has been noted for some polysaccharide simulations with the GLYCAM06j force field. Here, we investigate the cause of this phenomenon through comparative simulations of a range of saccharides with both the GLYCAM06j and the CHARMM36 carbohydrate force fields. We find that conformational collapse in GLYCAM06j occurs for saccharide chains containing the deoxy sugar α-l-rhamnose after relatively long simulation intervals. Further, we explore the mechanism of conformational collapse and show that this phenomenon arises because of the anomalous low energy in GLYCAM06j (as compared to quantum mechanical calculations) of a specific orientation of α-l-Rha to α-l-Rha glycosidic linkages, which are subsequently sustained by intramolecular interactions in the saccharide chain. We identify the lack of partial charges on aliphatic hydrogens in GLYCAM as the source of this anomaly, demonstrating that addition of small partial atomic charges on the aliphatic protons in rhamnose removes the conformational collapse phenomenon. This work reveals the large cumulative impact that small partial charges may have on the dynamic behavior of polysaccharides and indicates that future reparameterization of the GLYCAM06j force field should investigate the addition of partial charges on all aliphatic hydrogens.

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.

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.

Structural Determinants of Substrate Recognition and Catalysis by Heparan Sulfate Sulfotransferases

Heparan sulfate (HS) and heparin contain imprinted "sulfation codes", which dictate their diverse physiological and pathological functions. A group of orchestrated biosynthetic enzymes cooperate in polymerizing and modifying HS chains. The biotechnological development of enzymes that can recreate this sulfation pattern on synthetic heparin is challenging, primarily due to the paucity of quantitative data for sulfotransferase enzymes. Herein, we identified critical structural characteristics that determine substrate specificity and shed light on the catalytic mechanism of sugar sulfation of two HS sulfotransferases, 2-O-sulfotransferase (HS2ST) and 6-O-sulfotransferase (HS6ST). Two sets of molecular clamps in HS2ST recognize appropriate substrates; these clamps flank the acceptor binding site on opposite sides. The hexuronic epimers, and not their puckers, have a critical influence on HS2ST selectivity. In contrast, HS6ST recognizes a broader range of substrates. This promiscuity is granted by a conserved tryptophan residue, W210, that positions the acceptor within the active site for catalysis by means of strong electrostatic interactions. Lysines K131 and K132 act in concert with a second tryptophan, W153, shedding water molecules from within the active site, thus providing HS6ST with a binding preference toward 2-O-sulfated substrates. QM/MM calculations provided valuable mechanistic insights into the catalytic process, identifying that the sulfation of both HS2ST and HS6ST follows a SN2-like mechanism. When they are taken together, our findings reveal the molecular basis of how these enzymes recognize different substrates and catalyze sugar sulfation, enabling the generation of enzymes that could create specific heparin epitopes.

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.

Molecular simulation, characteristics and mechanism of thermal-responsive acetylated amylose V-type helical complexes

To explore the thermal-responsive characteristics of acetylated amylose-guest V-type helical complexes (AAGHCs) and their potential use as thermal-responsive drug carriers, different types of AAGHCs were built, in which acetylated amylose was used as a host, and iodine, propofol, or hexane was utilized as the guest molecule. Their thermal-responsive characteristics were investigated through molecular dynamic (MD) simulation and corresponding experiments. MD simulation showed that the thermal-responsive helix-unfolding and guest-release behavior in AAGHCs, and the complete unfolding of AAGHC could be divided into brewing, triggering and collapsing periods. Energy analysis revealed that the Lana-Jones potential is an important binding energy that bridges host and guest molecules and enhances the stability of the helix. The various types or number of guests showed different binding energies. The stronger the binding energy, higher is the temperature required to trigger the unfolding of the helix and the releasing of guests. FT-IR and X-ray diffraction analyses confirmed the structures of AAGHCs. The change in hydrated size, and UV-VIS absorption of AAGHCs at high temperatures both confirmed the thermal-responsiveness of AAGHCs. The fluorescence fluctuation of loaded 7-hydroxycoumarin reflected the same thermal-responsive process and mechanism as MD simulation. This study provides meaningful theoretical guidance for the design of thermal-responsive drug carriers based on acetylated amylose-guest V-type helical complexes.

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.

Oxidative damage to hyaluronan-CD44 interactions as an underlying mechanism of action of oxidative stress-inducing cancer therapy

Multiple cancer therapies nowadays rely on oxidative stress to damage cancer cells. Here we investigated the biological and molecular effect of oxidative stress on the interaction between CD44 and hyaluronan (HA), as interrupting their binding can hinder cancer progression. Our experiments demonstrated that the oxidation of HA decreased its recognition by CD44, which was further enhanced when both CD44 and HA were oxidized. The reduction of CD44-HA binding negatively affected the proliferative state of cancer cells. Our multi-level atomistic simulations revealed that the binding free energy of HA to CD44 decreased upon oxidation. The effect of HA and CD44 oxidation on CD44-HA binding was similar, but when both HA and CD44 were oxidized, the effect was much larger, in agreement with our experiments. Hence, our experiments and computations support our hypothesis on the role of oxidation in the disturbance of CD44-HA interaction, which can lead to the inhibition of proliferative signaling pathways inside the tumor cell to induce cell death.

Structural changes of rice starch and activity inhibition of starch digestive enzymes by anthocyanins retarded starch digestibility

The effects of anthocyanins on in vitro and in vivo digestibility of rice starch were evaluated in this study. Then, the effects of anthocyanins on physicochemical properties of rice starch and on starch digestive enzymes (α-amylase and α-glucosidase) were investigated to understand the mechanism of the effects of anthocyanins on starch digestibility. Characterization of physicochemical properties of rice starch indicates a structural change due to the presence of anthocyanins, hindering its access to starch digestive enzymes. Besides, anthocyanins inhibited the activities of starch digestive enzymes by binding to their active sites, competing with the substrates and changing the secondary structure of the enzymes. The above stated changes of rice starch and starch digestive enzymes due to the presence of anthocyanins both contributed to retarding the digestibility of rice starch. This study could offer some theoretical guidance to the development of new type rice-based food with low glycemic index.

Relation between glycosidic linkage, structure and dynamics of α- and β-glucans in water

In a molecular dynamics simulation study of several oligosaccharides comprising of the very basic building block of carbohydrate, the α- or β-d glucopyranose units, linked by any one of the 1-3/1-4 or 1-6 glycosidic linkages, we compare and contrast their structural and dynamical properties. Results indicate that the litheness of the oligosaccharide chain is noticeably controlled by the composition, anomeric nature and glycosidic linkage type of the units. In mixed β 1-4/1-3 d-glucopyranosides, as those found in oats and barley, the ratio of the β 1-4 and β 1-3 linked residues is crucial in determining the structural and dynamical attributes. Principal component analysis (PCA) using the internal coordinates of torsion angles subtended by glycosidic oxygen atoms and subsequent K-means clustering of the dynamical space spanned by PC1 to PC2 point to the dynamical and structural disparity in the various types of oligosaccharides studied. The properties simulated in this work are meant to provide a systematic yet comparative understanding of the importance of linkage and anomericity on the oligosaccharide chain properties and are in line with some experimental structural attributes.

The hitchhiker's guide to the periplasm: Unexpected molecular interactions of polymyxin B1 in E. coli

The periplasm of Gram-negative bacteria is a complex, highly crowded molecular environment. Little is known about how antibiotics move across the periplasm and the interactions they experience. Here, atomistic molecular dynamics simulations are used to study the antibiotic polymyxin B1 within models of the periplasm, which are crowded to different extents. We show that PMB1 is likely to be able to "hitchhike" within the periplasm by binding to lipoprotein carriers-a previously unreported passive transport route. The simulations reveal that PMB1 forms both transient and long-lived interactions with proteins, osmolytes, lipids of the outer membrane, and the cell wall, and is rarely uncomplexed when in the periplasm. Furthermore, it can interfere in the conformational dynamics of native proteins. These are important considerations for interpreting its mechanism of action and are likely to also hold for other antibiotics that rely on diffusion to cross the periplasm.

ADP-ribose and analogues bound to the deMARylating macrodomain from the bat coronavirus HKU4

Macrodomains are proteins that recognize and hydrolyze ADP ribose (ADPR) modifications of intracellular proteins. Macrodomains are implicated in viral genome replication and interference with host cell immune responses. They are important to the infectious cycle of Coronaviridae and Togaviridae viruses. We describe crystal structures of the conserved macrodomain from the bat coronavirus (CoV) HKU4 in complex with ligands. The structures reveal a binding cavity that accommodates ADPR and analogs via local structural changes within the pocket. Using a radioactive assay, we present evidence of mono-ADPR (MAR) hydrolase activity. In silico analysis presents further evidence on recognition of the ADPR modification for hydrolysis. Mutational analysis of residues within the binding pocket resulted in diminished enzymatic activity and binding affinity. We conclude that the common structural features observed in the macrodomain in a bat CoV contribute to a conserved function that can be extended to other known macrodomains.

Diversity of the Hydrogen Bond Network and Its Impact on NMR Parameters of Amylose B Polymorph: A Study Using Molecular Dynamics and DFT Calculations Within Periodic Boundary Conditions

Classical molecular dynamics simulations have been combined with quantum (DFT) calculations of ¹³C NMR parameters in order to relate the experimental spectrum of the double-helix form of the amylose B-polymorph in highly crystalline conditions not only to its 3D structure but also to the arrangement of atoms in the crystal lattice. Structures obtained from these simulations or from geometry optimization procedures at the DFT level have shown the presence of hydrogen bond networks between sugars of the same helix or between residues of the two chains of the double helix. ¹³C NMR parameter calculations have revealed the impact of such a network on the chemical shifts of carbon atoms. In addition, DFT calculations using periodic boundary conditions were compulsory to highlight the presence of two types of sugar within the crystal sample. It allows us to confirm, theoretically, the experimental hypothesis that the existence of two distinct sugar types in the NMR spectrum is a consequence of crystal packing.

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.

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.

Antibiotic Uptake Across Gram-Negative Outer Membranes: Better Predictions Towards Better Antibiotics

Crossing the Gram-negative bacterial membrane poses a major barrier to antibiotic development, as many small molecules that can biochemically inhibit key bacterial processes are rendered microbiologically ineffective by their poor cellular uptake. The outer membrane is the major permeability barrier for many drug-like molecules, and the chemical properties that enable efficient uptake into mammalian cells fail to predict bacterial uptake. We have developed a computational method for accurate prospective prediction of outer membrane uptake of drug-like molecules, which we combine with a new medium-throughput experimental assay of outer membrane vesicle swelling. Parallel molecular dynamics simulations of compound uptake through Escherichia coli (E. coli) OmpF are used to successfully and quantitatively predict experimental permeabilities measured via either outer membrane swelling or prior liposome-swelling measurements. These simulations are analyzed using an inhomogeneous solubility-diffusion model to yield predictions of permeability. For most polar molecules we test, outer membrane permeability also correlates well with whole-cell uptake. The ability to accurately predict and measure outer membrane uptake of a wide variety of small molecules will enable simpler determination of which molecular scaffolds and which derivatives are most promising prior to extensive chemical synthesis. It will also assist in formulating a more systematic understanding of the chemical determinants of outer membrane permeability.

Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations

The conformational flexibility of the glycosaminoglycans (GAGs) is known to be key in their binding and biological function, for example in regulating coagulation and cell growth. In this work, we employ enhanced sampling molecular dynamics simulations to probe the ring conformations of GAG-related monosaccharides, including a range of acetylated and sulfated GAG residues. We first perform unbiased MD simulations of glucose anomers and the epimers glucuronate and iduronate. These calculations indicate that in some cases, an excess of 15 μs is required for adequate sampling of ring pucker due to the high energy barriers between states. However, by applying our recently developed msesMD simulation method (multidimensional swarm-enhanced sampling molecular dynamics), we were able to quantitatively and rapidly reproduce these ring pucker landscapes. From msesMD simulations, the puckering free energy profiles were then compared for 15 further monosaccharides related to GAGs; this includes to our knowledge the first simulation study of sulfation effects on β-GalNAc ring puckering. For the force field employed, we find that in general the calculated pucker free energy profiles for sulfated sugars were similar to the corresponding unsulfated profiles. This accords with recent experimental studies suggesting that variation in ring pucker of sulfated GAG residues is primarily dictated by interactions with surrounding residues rather than by intrinsic conformational preference. As an exception to this, however, we predict that 4-O-sulfation of β-GalNAc leads to reduced ring rigidity, with a significant lowering in energy of the ¹C₄ ring conformation; this observation may have implications for understanding the structural basis of the biological function of β-GalNAc-containing glycosaminoglycans such as dermatan sulfate.

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.

Hydration and Dispersion Forces in Hydroxypropylcellulose Phase Behavior

Many-body polarization and hydration forces can strongly affect the equilibrium structure and energetics of mixed phases. Accurately reproducing both forces presents a challenge to force field models because it requires balancing hydrogen bonding at short range with many-body orientational order and dispersive attraction at long range. This work reports the first comparison of experimental measurements of the pressure-area isotherm for hydroxypropylcellulose (HPC) against molecular dynamics results with four different force field models-united-atom, all-atom (OPLS and CHARMM), and Drude oscillator models. All force fields exhibit the experimentally determined, exponentially shaped repulsive force at short range. Above a critical temperature of about 40 °C and a lattice spacing of around 14 Å, HPC experiments show a reversible, heat-induced polymer aggregation into an ordered phase driven by loss of water. The nonpolarizable force fields do not display the critical point and instead show biphasic behavior at all temperatures tested. This indicates net attractive forces at intermediate lattice spacings. In contrast, the Drude polarizable force field shows positive osmotic pressure and a single, homogeneous phase over all temperatures and spacings tested. Analysis of structural data from our simulations provides several clues to help interpret these findings. Although all force fields show similar water-water hydrogen bond numbers in the mixed phase, the polarizable model predicts that water-HPC hydrogen bonds are much more favorable than HPC-HPC hydrogen bonds when polymers are dispersed. At high density, water is driven out and replaced by HPC-HPC hydrogen bonds. The polarizable force field shows that both effects have a stronger dependence on polymer density than any of the nonpolarizable models. Our observations support the conclusion that hydration forces are coupled to the polymer coordination number by local, structural waters and that long-range dispersive attraction is overestimated by pairwise additive models.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

Structural Aspects of the O-glycosylation Linkage in Glycopeptides via MD Simulations and Comparison with NMR Experiments

A powerful conformational searching and enhanced sampling simulation method, and unbiased molecular dynamics simulations have been used along with NMR spectroscopic observables to provide a detailed structural view of O-glycosylation. For four model systems, the force-field parameters can accurately predict experimental NMR observables (J couplings and NOE's). This enables us to derive conclusions based on the generated ensembles, in which O-glycosylation affects the peptide backbone conformation by forcing it towards to an extended conformation. An exception is described for β-GalNAc-Thr where the α content is increased and stabilized via hydrogen bonding between the sugar and the peptide backbone, which was not observed in the rest of the studied systems. These observations might offer an explanation for the evolutionary preference of α-linked GalNAc glycosylation instead of a β link.

Describing Molecular Polarizability by a Bond Capacity Model

We propose a bond capacity model for describing molecular polarization in force field energy functions at the charge-only level. Atomic charges are calculated by allowing charge to flow between atom pairs according to a bond capacity and a difference in electrostatic potential. The bond capacity is closely related to the bond order and decays to zero as the bond distance is increased. The electrostatic potential is composed of an intrinsic potential, identified as the electronegativity, and a screened Coulomb potential from all other charges. The bond capacity model leads to integer fragment charges upon bond dissociation and displays linear scaling of the polarizability with system size. Bond capacity parameters can be derived from reference molecular polarizabilities, while electronegativity parameters can be derived from reference atomic charges or a reference molecular electrostatic potential. Out-of-plane polarization for planar systems is modeled by off-nuclei charge sites. The model is shown to be able to reproduce anisotropic reference molecular polarizabilities with an accuracy of ∼10% using a limited set of bond capacity parameters and can describe both inter- and intramolecular polarization.

Perspective on computational simulations of glycosaminoglycans

Glycosaminoglycans (GAGs) represent a formidable frontier for chemists, biochemists, biologists, medicinal chemists and drug delivery specialists because of massive structural complexity. GAGs are arguably the most complex, natural linear biopolymers with theoretical diversity orders of magnitude higher than proteins and nucleic acids. Yet, this diversity remains generally untapped. Computational approaches offer major routes to understand GAG structure and dynamics so as to enable novel applications of these biopolymers. In fact, computational algorithms, softwares, online tools and techniques have reached a level of sophistication that help understand atomistic details of conformational variation and protein recognition of individual GAG sequences. This review describes current approaches and challenges in computational study of GAGs. It presents a history of major findings since the earliest mention of GAGs (the 1960s), the development of parameters and force fields specific for GAGs, and the application of these tools in understanding GAG structure-function relationship. This review also presents a section on how to perform simulation of GAGs, which is directed toward researchers interested in entering this promising field with potential to impact therapy.