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.

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.

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.

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.

Bioinspired and highly oriented clay nanocomposites with a xyloglucan biopolymer matrix: extending the range of mechanical and barrier properties.

The development of clay bionanocomposites requires processing routes with nanostructural control. Moreover, moisture durability is a concern with water-soluble biopolymers. Here, oriented bionanocomposite coatings with strong in-plane orientation of clay platelets are for the first time prepared by continuous water-based processing. Montmorillonite (MTM) and a "new" unmodified biological polymer (xyloglucan (XG)) are combined. The resulting nanocomposites are characterized by FE-SEM, TEM, and XRD. XG adsorption on MTM is measured by quartz crystal microbalance analysis. Mechanical and gas barrier properties are measured, also at high relative humidity. The reinforcement effects are modeled. XG dimensions in composites are estimated using atomistic simulations. The nanostructure shows highly oriented and intercalated clay platelets. The reinforcement efficiency and effects on barrier properties are remarkable and are likely to be due to highly oriented and well-dispersed MTM and strong XG-MTM interactions. Properties are well preserved in humid conditions and the reasons for this are discussed.

Overview of computer modeling of cellulose.

Although it has a deceptively simple primary structure, the collective organization of bulk cellulose, particularly as it exists in cellulose fibers in the cell walls of living plants and other organisms, is quite diverse and complex. While some experimental techniques, such as vibrational spectroscopy and diffraction from partially crystalline samples, are able to provide insights into the organization of bulk cellulose, its intrinsic complexity has left many questions still unanswered. For this reason, additional probes of cellulose structure would be highly desirable. With the continuing advances in computer power through massive parallelization, and the steady progress in computer codes and force fields for modeling carbohydrate systems, molecular mechanics simulations have become an attractive means of studying cellulosic systems at the atomic and molecular level. The coming decade will almost certainly see remarkable advances in the understanding of cellulose using such simulations.

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.