Analysis of Solvation and Gelation Behavior of Methylcellulose Using Atomistic Molecular Dynamics Simulations

We adopt atomistic molecular dynamics (MD) simulations to study the solvation and gelation behavior of homogeneous methylcellulose (MC) and model random oligomers that represent the commercial cellulosic polymer product METHOCEL A in water and acetone solvents. We demonstrate that the two carbohydrate-specific GROMOS force fields, GROMOS 45A4 and 56Acarbo, are capable of reproducing characteristic angle distributions and the persistence length of MC chains reported in the literature. We then use the GROMOS 56Acarbo force field in both single-chain and multiple-chain simulations to characterize their solvation behavior through radial distribution functions, hydrogen-bond counts, and contact map analyses. We find that the un-methylated O6 position on the cellulose ring forms the most hydrogen bonds, followed by O2 and O3, implying that methylation at the 6 position reduces hydrogen bonding more than does methylation at other positions. O6–O6 is the most probable intermolecular hydrogen bond between differen...