Simulation studies of the insolubility of cellulose

Molecular dynamics simulations have been used to calculate the potentials of mean force for separating short cellooligomers in aqueous solution as a means of estimating the contributions of hydrophobic stacking and hydrogen bonding to the insolubility of crystalline cellulose. A series of four potential of mean force (pmf) calculations for glucose, cellobiose, cellotriose, and cellotetraose in aqueous solution were performed for situations in which the molecules were initially placed with their hydrophobic faces stacked against one another, and another for the cases where the molecules were initially placed adjacent to one another in a co-planar, hydrogen-bonded arrangement, as they would be in cellulose Iβ. From these calculations, it was found that hydrophobic association does indeed favor a crystal-like structure over solution, as might be expected. Somewhat more surprisingly, hydrogen bonding also favored the crystal packing, possibly in part because of the high entropic cost for hydrating glucose hydroxyl groups, which significantly restricts the configurational freedom of the hydrogen-bonded waters. The crystal was also favored by the observation that there was no increase in chain configurational entropy upon dissolution, because the free chain adopts only one conformation, as previously observed, but against intuitive expectations, apparently due to the persistence of the intramolecular O3-O5 hydrogen bond.