Understanding the Effect of Solvent Environment on the Interaction of Hydronium Ion with Biomass Derived Species: A Molecular Dynamics and Metadynamics Investigation

The addition of aprotic solvents results in higher reactivities and selectivities in many key aqueous phase biomass reactions, including the acid-catalyzed conversion of fructose to 5-hydroxyl methyl furfural (HMF). The addition of certain co-solvents inhibits the formation of humins via preferential solvation of key functional groups and can alter reaction kinetics. An important factor in this context is the relative stability of the hydronium ion (the catalyst) in the vicinity of the biomass moiety as compared to that in bulk, as it could determine its efficacy in the protonation step. Hence, in the present work, molecular dynamics (MD) simulations of HMF (the model product) and fructose (the model reactant) in acidic water and water-DMSO mixtures are performed to analyze their interaction with the hydronium ions. We show that the presence of DMSO favors the interaction of the hydronium ion with fructose, whereas it has a detrimental effect on the interaction of hydronium ion with HMF. Well-tempered metadynamics (WT-MTD) simulations are performed to determine the relative stability of the hydronium ion in the immediate vicinity of fructose and HMF, as compared to that in the bulk solvent phase, as a function of solvent composition. We find that DMSO improves the stabilization of the hydronium ions in the first solvation shell of fructose compared to that in the bulk solvent. On the other hand, hydronium ions become less stable in the immediate vicinity of HMF, as the concentration of DMSO increases.