Malin Bergenstråhle-Wohlert

Comparison of Cellulose Iβ Simulations with Three Carbohydrate Force Fields

Molecular dynamics simulations of cellulose have recently become more prevalent due to increased interest in renewable energy applications, and many atomistic and coarse-grained force fields exist that can be applied to cellulose. However, to date no systematic comparison between carbohydrate force fields has been conducted for this important system. To that end, we present a molecular dynamics simulation study of hydrated, 36-chain cellulose Iβ microfibrils at room temperature with three carbohydrate force fields (CHARMM35, GLYCAM06, and Gromos 45a4) up to the near-microsecond time scale. Our results indicate that each of these simulated microfibrils diverge from the cellulose Iβ crystal structure to varying degrees under the conditions tested. The CHARMM35 and GLYCAM06 force fields eventually result in structures similar to those observed at 500 K with the same force fields, which are consistent with the experimentally observed high-temperature behavior of cellulose I. The third force field, Gromos 45a4, produces behavior significantly different from experiment, from the other two force fields, and from previously reported simulations with this force field using shorter simulation times and constrained periodic boundary conditions. For the GLYCAM06 force field, initial hydrogen-bond conformations and choice of electrostatic scaling factors significantly affect the rate of structural divergence. Our results suggest dramatically different time scales for convergence of properties of interest, which is important in the design of computational studies and comparisons to experimental data. This study highlights that further experimental and theoretical work is required to understand the structure of small diameter cellulose microfibrils typical of plant cellulose.

Non-exchanging hydroxyl groups on the surface of cellulose fibrils: The role of interaction with water

The interaction of water with cellulose stages many unresolved questions. Here 2H MAS NMR and IR spectra recorded under carefully selected conditions in 1H-2H exchanged, and re-exchanged, cellulose samples are presented. It is shown here, by a quantitative and robust approach, that only two of the three available hydroxyl groups on the surface of cellulose fibrils are exchanging their hydrogen with the surrounding water molecules. This finding is additionally verified and explained by MD simulations which demonstrate that the 1HO(2) and 1HO(6) hydroxyl groups of the constituting glucose units act as hydrogen-bond donors to water, while the 1HO(3) groups behave exclusively as hydrogen-bond acceptors from water and donate hydrogen to their intra-chain neighbors O(5). We conclude that such a behavior makes the latter hydroxyl group unreactive to hydrogen exchange with water.

Alternative hydrogen bond models of cellulose II and IIII based on molecular force-fields and density functional theory

Alternative hydrogen-bond structures were found for cellulose II and IIII based on molecular dynamics simulations using four force fields and energy optimization based on density functional theory. All the modeling results were in support to the new hydrogen-bonding network. The revised structures of cellulose II and IIII differ with the fiber diffraction models mainly in the orientation of two hydroxyl groups, namely, OH2 and OH6 forming hydrogen-bond chains perpendicular to the cellulose molecule. In the alternative structures, the sense of hydrogen bond is inversed but little difference can be seen in hydrogen bond geometries. The preference of these alternative hydrogen bond structures comes from the local stabilization of hydroxyl groups with respect to the β carbon. On the other hand when simulated fiber diffraction patterns were compared with experimental ones, the current structure of cellulose II with higher energy and the alternative structure of cellulose IIII with lower energy were in better agreement.

Molecular adhesion at clay nanocomposite interfaces depends on counterion hydration-molecular dynamics simulation of montmorillonite/xyloglucan.

Nacre-mimetic clay/polymer nanocomposites with clay platelet orientation parallel to the film surface show interesting gas barrier and mechanical properties. In moist conditions, interfacial adhesion is lowered and mechanical properties are reduced. Molecular dynamic simulations (MD) have been performed to investigate the effects of counterions on molecular adhesion at montmorillonite clay (Mnt)-xyloglucan (XG) interfaces. We focus on the role of monovalent cations K(+), Na(+), and Li(+) and the divalent cation Ca(2+) for mediating and stabilizing the Mnt/XG complex formation. The conformation of adsorbed XG is strongly influenced by the choice of counterion and so is the simulated work of adhesion. Free energy profiles that are used to estimate molecular adhesion show stronger interaction between XG and clay in the monovalent cation system than in divalent cation system, following a decreasing order of K-Mnt, Na-Mnt, Li-Mnt, and Ca-Mnt. The Mnt clay hydrates differently in the presence of different counterions, leading to a chemical potential of water that is highest in the case of K-Mnt, followed by Na-Mnt and Li-Mnt, and lowest in the case of Ca-Mnt. This means that water is most easily displaced from the interface in the case of K-Mnt, which contributes to the relatively high work of adhesion. In all systems, the penalty of replacing polymer with water at the interface gives a positive contribution to the work of adhesion of between 19 and 35%. Our work confirms the important role of counterions in mediating the adsorption of biopolymer XG to Mnt clays and predicts potassium or sodium as the best choice of counterions for a Mnt-based biocomposite design.

Solubility of Cellulose in Supercritical Water Studied by Molecular Dynamics Simulations

The insolubility of cellulose in ambient water and most aqueous systems presents a major scientific and practical challenge. Intriguingly though, the dissolution of cellulose has been reported to occur in supercritical water. In this study, cellulose solubility in ambient and supercritical water of varying density (0.2, 0.7, and 1.0 g cm(-3)) was studied by atomistic molecular dynamics simulations using the CHARMM36 force field and TIP3P water. The Gibbs energy of dissolution was determined between a nanocrystal (4 × 4 × 20 anhydroglucose residues) and a fully dissociated state using the two-phase thermodynamics model. The analysis of Gibbs energy suggested that cellulose is soluble in supercritical water at each of the studied densities and that cellulose dissolution is typically driven by the entropy gain upon the chain dissociation while simultaneously hindered by the loss of solvent entropy. Chain dissociation caused density augmentation around the cellulose chains, which improved water-water bonding in low density supercritical water whereas the opposite occurred in ambient and high density supercritical water.

A molecular dynamics study of the effect of glycosidic linkage type in the hemicellulose backbone on the molecular chain flexibility

The macromolecular conformation of the constituent polysaccharides in lignocellulosic biomass influences their supramolecular interactions, and therefore their function in plants and their performance in technical products. The flexibility of glycosidic linkages from the backbone of hemicelluloses was studied by evaluating the conformational freedom of the φ and ψ dihedral angles using molecular dynamic simulations, additionally selected molecules were correlated with experimental data by nuclear magnetic resonance spectroscopy. Three types of β-(1→4) glycosidic linkages involving the monosaccharides (Glcp, Xylp and Manp) present in the backbone of hemicelluloses were defined. Different di- and tetrasaccharides with combinations of such sugar monomers from hemicelluloses were simulated, and free energy maps of the φ - ψ space and hydrogen-bonding patterns were obtained. The glycosidic linkage between Glc-Glc or Glc-Man (C-type) was the stiffest with mainly one probable conformation; the linkage from Man-Man or Man-Glc (M-type) was similar but with an increased probability for an alternative conformation making it more flexible, and the linkage between two Xyl-units (X-type) was the most flexible with two almost equally populated conformations. Glycosidic linkages of the same type showed essentially the same conformational space in both disaccharides and in the central region of tetrasaccharides. Different probabilities of glycosidic linkage conformations in the backbone of hemicelluloses can be directly estimated from the free energy maps, which to a large degree affect the overall macromolecular conformations of these polymers. The information gained contributes to an increased understanding of the function of hemicelluloses both in the cell wall and in technical products.

Temperature Dependence of Hydroxymethyl Group Rotamer Populations in Cellooligomers

Empirical force fields for computer simulations of carbohydrates are often implicitly assumed to be valid also at temperatures different from room temperature for which they were optimized. Herein, the temperature dependence of the hydroxymethyl group rotamer populations in short oligosaccharides is investigated using molecular dynamics simulations and NMR spectroscopy. Two oligosaccharides, viz., methyl β-cellobioside and β-cellotetraose were simulated using three different carbohydrate force fields (CHARMM C35, GLYCAM06, and GROMOS 56Acarbo) in combination with different water models (SPC, SPC/E, and TIP3P) using replica exchange molecular dynamics simulations. For comparison, hydroxymethyl group rotamer populations were investigated for methyl β-cellobioside and cellopentaose based on measured NMR (3)JH5,H6 coupling constants, in the latter case by using a chemical shift selective NMR-filter. Molecular dynamics simulations in combination with NMR spectroscopy show that the temperature dependence of the hydroxymethyl rotamer population in these short cellooligomers, in the range 263-344 K, generally becomes exaggerated in simulations when compared to experimental data, but also that it is dependent on simulation conditions, and most notably properties of the water model.

Molecular mechanisms for the adhesion of chitin and chitosan to montmorillonite clay

Molecular dynamics simulations have been performed to investigate molecular adhesion of chitin and chitosan oligomers to montmorillonite (Mnt) clay at different degrees of acetylation (DA, 0%, 20%, ...

On the anomalous temperature dependence of cellulose aqueous solubility

The solubility of cellulose in water-based media is promoted by low temperature, which may appear counter-intuitive. An explanation to this phenomenon has been proposed that is based on a temperature-dependent orientation of the hydroxymethyl group. In this paper, this hypothesis is investigated using molecular dynamics computer simulations and NMR spectroscopy, and is discussed in conjunction with alternative explanations based on solvent–solute and solvent–solvent hydrogen bond formation respectively. It is shown that neither simulations nor experiments lend support to the proposed mechanism based on the hydroxymethyl orientation, whereas the two alternative explanations give rise to two distinct contributions to the hydration free energy of cellooligomers.