Potential of mean force for separation of the repeating units in cellulose and hemicellulose.

As a first step toward understanding the energetics of removal of cello-oligomers from the cellulose surface, we have performed umbrella sampling calculations to determine the free energy required for separation of repeating units of cellulose and hemicellulose from each other. Molecular dynamics (MD) simulations were performed for both the stacked and non-stacked arrangements of the cellobiose pair system and the xylobiose pair system. These stacked and non-stacked arrangements were taken as representative systems for the crystalline and amorphous domains in cellulose and hemicellulose. In addition, similar calculations were also carried out to determine the energetics involved in the separation of the cellobiose-xylobiose molecule pair in the non-stacked arrangement. The potential of mean force profiles exhibit a single minimum in all cases and are qualitatively similar. Our results show that the location of the minimum as well as the depth of the well can be correlated with the size of the disaccharide molecules.

[1]  J. Brady,et al.  MODEL DEPENDENCE OF THE ANISOTROPIC STRUCTURING OF SOLVENT WATER AROUND SUGARS IN MOLECULAR DYNAMICS SIMULATIONS , 1997 .

[2]  B. J. Hardy,et al.  Molecular dynamics simulation of cellobiose in water , 1993, J. Comput. Chem..

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

[4]  Bernard R. Brooks,et al.  Solvent-Induced Forces between Two Hydrophilic Groups , 1994 .

[5]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[6]  Peter A. Kollman,et al.  Benzene Dimer: A Good Model for π−π Interactions in Proteins? A Comparison between the Benzene and the Toluene Dimers in the Gas Phase and in an Aqueous Solution , 1996 .

[7]  G. P. Johnson,et al.  Cellulose and the twofold screw axis: modeling and experimental arguments , 2009 .

[8]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[9]  Anna-Pitschna E. Kunz,et al.  A comparison of methods to compute the potential of mean force. , 2007, Chemphyschem : a European journal of chemical physics and physical chemistry.

[10]  Alfred D. French,et al.  Quantum mechanics studies of cellobiose conformations , 2006 .

[11]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[12]  Srinivas Karra,et al.  Modeling Intrinsic Kinetics of Enzymatic Cellulose Hydrolysis , 2008, Biotechnology progress.

[13]  M. Himmel,et al.  Simulation studies of the insolubility of cellulose. , 2010, Carbohydrate research.

[14]  Philippe H. Hünenberger,et al.  Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water , 2010, J. Comput. Chem..

[15]  Mark Holtzapple,et al.  Coordinated development of leading biomass pretreatment technologies. , 2005, Bioresource technology.

[16]  G. P. Johnson,et al.  Conformational flexibility of soluble cellulose oligomers: chain length and temperature dependence. , 2009, Journal of the American Chemical Society.

[17]  John W. Brady,et al.  Calculation of the Potential of Mean Force for the Binding of Glucose to Benzene in Aqueous Solution , 2000 .

[18]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[19]  Alfred D. French,et al.  Disaccharide conformational maps: adiabaticity in analogues with variable ring shapes , 2008 .

[20]  Kenneth M. Merz,et al.  A carbohydrate force field for amber and its application to the study of saccharide to surface adsorption , 1997 .

[21]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[22]  T. Teeri,et al.  The Cellulases Endoglucanase I and Cellobiohydrolase II of Trichoderma reesei Act Synergistically To Solubilize Native Cotton Cellulose but Not To Decrease Its Molecular Size , 1996, Applied and environmental microbiology.

[23]  G. P. Johnson,et al.  Comparison of different force fields for the study of disaccharides. , 2009, Carbohydrate research.

[24]  L. Nilsson,et al.  Structure and Dynamics of the TIP3P, SPC, and SPC/E Water Models at 298 K , 2001 .

[25]  G. P. Johnson,et al.  Advanced conformational energy surfaces for cellobiose** , 2004 .

[26]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[27]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[28]  I. S. Pretorius,et al.  Microbial Cellulose Utilization: Fundamentals and Biotechnology , 2002, Microbiology and Molecular Biology Reviews.

[29]  J. Brady.  Molecular dynamics simulations of β-d-glucopyranose , 1987 .

[30]  Robert J. Woods,et al.  NMR Spectroscopy and Computer Modeling of Carbohydrates: Recent Advances , 2006 .

[31]  L. Lynd,et al.  Toward an aggregated understanding of enzymatic hydrolysis of cellulose: Noncomplexed cellulase systems , 2004, Biotechnology and bioengineering.