Solvent Influence on Cellulose 1,4-β-Glycosidic Bond Cleavage: A Molecular Dynamics and Metadynarnics Study

We explore the influence of two solvents, namely water and the ionic liquid 1-ethyl-3-methylimidazolium acetate (EmimAc), on the conformations of two cellulose models (cellobiose and a chain of 40 glucose units) and the solvent impact on glycosidic bond cleavage by acid hydrolysis by using molecular dynamics and metadynamics simulations. We investigate the rotation around the glycosidic bond and ring puckering, as well as the anomeric effect and hydrogen bonds, in order to gauge the effect on the hydrolysis mechanism. We find that EmimAc eases hydrolysis through stronger solvent cellulose interactions, which break structural and electronic barriers to hydrolysis. Our results indicate that hydrolysis in cellulose chains should start from the ends and not in the centre of the chain, which is less accessible to solvent.