Evaluation of coarse-grained mapping schemes for polysaccharide chains in cellulose.

A fundamental understanding of the intermolecular forces that bind polysaccharide chains together in cellulose is crucial for designing efficient methods to overcome the recalcitrance of lignocellulosic biomass to hydrolysis. Because the characteristic time and length scales for the degradation of cellulose by enzymatic hydrolysis or chemical pretreatment span orders of magnitude, it is important to closely integrate the molecular models used at each scale so that, ultimately, one may switch seamlessly between quantum, atomistic, and coarse-grained descriptions of the system. As a step towards that goal, four multiscale coarse-grained models for polysaccharide chains in a cellulose-Iα microfiber are considered. Using the force matching method, effective coarse-grained forces are derived from all-atom trajectories. Performance of the coarse-grained models is evaluated by comparing the intrachain radial distribution functions with those obtained using the all-atom reference data. The all-atom simulation reveals a double peak in the radial distribution function for sites within each glucose residue that arises from the distinct conformations sampled by the primary alcohol group in the glucose residues. The three-site and four-site coarse-grained models have sufficient degrees of freedom to predict this double peak while the one-site and two-site models do not. This is the first time that coarse-grained models have been shown to reproduce such subtle, yet important, molecular features in a polysaccharide chain. The relative orientations between glucose residues along the polysaccharide chain are evaluated and it is found that the four-site coarse-grained model is best at reproducing the glucose-glucose conformations observed in the all-atom simulation. The success of the four-site coarse-grained model underscores the importance of decoupling the pyranose ring from the oxygen atom in the glycosidic bond when developing all-atom to coarse-grained mapping schemes for polysaccharides.

[1]  T. Higuchi,et al.  Transformation of native cellulose crystals induced by saturated steam at high temperatures , 1987 .

[2]  W. Goddard,et al.  M3B: A Coarse Grain Force Field for Molecular Simulations of Malto-Oligosaccharides and Their Water Mixtures , 2004 .

[3]  Alexander D. MacKerell,et al.  CHARMM Additive All-Atom Force Field for Glycosidic Linkages between Hexopyranoses. , 2009, Journal of chemical theory and computation.

[4]  C. Wyman,et al.  Features of promising technologies for pretreatment of lignocellulosic biomass. , 2005, Bioresource technology.

[5]  M. Himmel,et al.  Comparison of Cellulose Iβ Simulations with Three Carbohydrate Force Fields. , 2012, Journal of chemical theory and computation.

[6]  James B. Adams,et al.  Interatomic Potentials from First-Principles Calculations: The Force-Matching Method , 1993, cond-mat/9306054.

[7]  Michael L Klein,et al.  Self-assembly and properties of diblock copolymers by coarse-grain molecular dynamics , 2004, Nature materials.

[8]  Jeremy C. Smith,et al.  REACH coarse-grained simulation of a cellulose fiber. , 2012, Biomacromolecules.

[9]  Harry J. Gilbert,et al.  Cellulosomes: highly efficient nanomachines designed to deconstruct plant cell wall complex carbohydrates. , 2010, Annual review of biochemistry.

[10]  J. Sugiyama,et al.  Combined infrared and electron diffraction study of the polymorphism of native celluloses , 1991 .

[11]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[12]  Mark S. Gordon,et al.  A coarse-grained model for β-d-glucose based on force matching , 2012, Theoretical Chemistry Accounts.

[13]  Gregory A Voth,et al.  Multiscale coarse-graining of monosaccharides. , 2007, The journal of physical chemistry. B.

[14]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[15]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[16]  Gregory A Voth,et al.  The multiscale coarse-graining method. II. Numerical implementation for coarse-grained molecular models. , 2008, The Journal of chemical physics.

[17]  Edward A Bayer,et al.  Applications of computational science for understanding enzymatic deconstruction of cellulose. , 2011, Current opinion in biotechnology.

[18]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[19]  M. Himmel,et al.  Computer simulation studies of microcrystalline cellulose Iβ , 2006 .

[20]  Gregory A. Voth,et al.  The multiscale coarse-graining method. I. A rigorous bridge between atomistic and coarse-grained models. , 2008, The Journal of chemical physics.

[21]  B. Brooks,et al.  Multiscale methods for macromolecular simulations. , 2008, Current opinion in structural biology.

[22]  Xiaolin Cheng,et al.  A Solvent-Free Coarse Grain Model for Crystalline and Amorphous Cellulose Fibrils. , 2011, Journal of chemical theory and computation.

[23]  Ilpo Vattulainen,et al.  Analysis of twisting of cellulose nanofibrils in atomistic molecular dynamics simulations. , 2011, The journal of physical chemistry. B.

[24]  Benjamin Lindner,et al.  Scaling of Multimillion-Atom Biological Molecular Dynamics Simulation on a Petascale Supercomputer. , 2009, Journal of chemical theory and computation.

[25]  E. Bayer,et al.  The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. , 2004, Annual review of microbiology.

[26]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[27]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[28]  Paul Langan,et al.  Crystal structure and hydrogen-bonding system in cellulose Ibeta from synchrotron X-ray and neutron fiber diffraction. , 2002, Journal of the American Chemical Society.

[29]  Gregory A Voth,et al.  Multiscale modeling of biomolecular systems: in serial and in parallel. , 2007, Current opinion in structural biology.

[30]  S. Ramakrishnan,et al.  Enzymatic hydrolysis of cellulose dissolved in N-methyl morpholine oxide/water solutions. , 2010, Bioresource technology.

[31]  Antti-Pekka Hynninen,et al.  Coarse-Grain Model for Glucose, Cellobiose, and Cellotetraose in Water. , 2011, Journal of chemical theory and computation.

[32]  Seema Singh,et al.  Visualization of biomass solubilization and cellulose regeneration during ionic liquid pretreatment of switchgrass , 2009, Biotechnology and bioengineering.

[33]  Christopher H. Chang,et al.  The energy landscape for the interaction of the family 1 carbohydrate-binding module and the cellulose surface is altered by hydrolyzed glycosidic bonds. , 2009, The journal of physical chemistry. B.

[34]  I. Tanaka,et al.  Structure of bacterial cellulose synthase subunit D octamer with four inner passageways , 2019 .

[35]  B. Dale,et al.  Coarse-grained model for the interconversion between native and liquid ammonia-treated crystalline cellulose. , 2012, The journal of physical chemistry. B.

[36]  Lars Berglund,et al.  A Coarse-Grained Model for Molecular Dynamics Simulations of Native Cellulose , 2011 .

[37]  Siewert J Marrink,et al.  Martini Coarse-Grained Force Field: Extension to Carbohydrates. , 2009, Journal of chemical theory and computation.

[38]  Gregory A Voth,et al.  Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching. , 2004, The Journal of chemical physics.