Structural changes in the molecular sheets along (hk0) planes derived from cellulose Iβ by molecular dynamics simulations

We studied the stability of molecular sheets with four cellotetraoses in an aqueous environment by molecular dynamics simulation to identify the molecular details of first structure as one of the possibilities in the course of crystallization of cellulose I. After simulation, the molecular sheets formed by van der Waals forces along the (11?0) and (110) crystal plane did not change their structures in an aqueous environment, whereas the other ones formed by hydrogen bonds along the (100) and (200) crystal plane changed into a van der Waals associated molecular sheet, similar to the former. These simulated molecular sheets formed by van der Waals forces were structurally stable in water because of their hydrophilic exterior and hydrophobic interior. Therefore, if the molecular sheet structures are formed in the real system, the sheets formed by van der Waals forces are probably the initial structure of crystallization. A close analysis indicated that these sheets could be classified into two groups in terms of the hydrogen bonding networks, camber angle, and main and side chain conformations. One group was the molecular sheets corresponding to the (110) after simulation. This sheet is probably rigid because intramolecular hydrogen bonds of the chains in the sheet are highly developed. The other group was the molecular sheets corresponding to (200), (100), and (11?0) crystal plane: the chains in these sheets seemed to be rather flexible due to their moderately developed intramolecular hydrogen bonds.