Hydrogen-Bonding Network and OH Stretch Vibration of Cellulose: Comparison of Computational Modeling with Polarized IR and SFG Spectra.

Hydrogen bonds play critical roles in noncovalent directional interactions determining the crystal structure of cellulose. Although diffraction studies accurately determined the coordinates of carbon and oxygen atoms in crystalline cellulose, the structural information on hydrogen atoms involved in hydrogen-bonding is still elusive. This could be complemented by vibrational spectroscopy; but the assignment of the OH stretch peaks has been controversial. In this study, we performed calculations using density functional theory with dispersion corrections (DFT-D2) for the cellulose Iß crystal lattices with the experimentally determined carbon and oxygen coordinates. DFT-D2 calculations revealed that the OH stretch vibrations of cellulose are highly coupled and delocalized through intra- and interchain hydrogen bonds involving all OH groups in the crystal. Additionally, molecular dynamics (MD) simulations of a single cellulose microfibril showed that the conformations of OH groups exposed at the microfibril surface are not well-defined. Comparison of the computation results with the experimentally determined IR dichroism of uniaxially aligned cellulose microfibrils and the peak positions of various cellulose crystals allowed unambiguous identification of OH stretch modes observed in the vibrational spectra of cellulose.

Hydration control of the mechanical and dynamical properties of cellulose.

The mechanical and dynamical properties of cellulose, the most abundant biomolecule on earth, are essential for its function in plant cell walls and advanced biomaterials. Cellulose is almost always found in a hydrated state, and it is therefore important to understand how hydration influences its dynamics and mechanics. Here, the nanosecond-time scale dynamics of cellulose is characterized using dynamic neutron scattering experiments and molecular dynamics (MD) simulation. The experiments reveal that hydrated samples exhibit a higher average mean-square displacement above ∼240 K. The MD simulation reveals that the fluctuations of the surface hydroxymethyl atoms determine the experimental temperature and hydration dependence. The increase in the conformational disorder of the surface hydroxymethyl groups with temperature follows the cellulose persistence length, suggesting a coupling between structural and mechanical properties of the biopolymer. In the MD simulation, 20% hydrated cellulose is more rigid than the dry form, due to more closely packed cellulose chains and water molecules bridging cellulose monomers with hydrogen bonds. This finding may have implications for understanding the origin of strength and rigidity of secondary plant cell walls. The detailed characterization obtained here describes how hydration-dependent increased fluctuations and hydroxymethyl disorder at the cellulose surface lead to enhancement of the rigidity of this important biomolecule.