Abstract

Molecular dynamics simulations were used to explore the distribution of linkages in coumaryl alcohol and coniferyl alcohol systems, including monomeric systems and monomers interacting with β-O4 dimers, respectively. Studying the interactions of two monolignols and the corresponding dimers sheds light on the preferred mechanism of reaction of the growing lignin polymer from the view of kinetic factors. The energy change upon association was quantified, and the distances between the centers of mass of different molecules, and the relative orientations between the phenol groups were calculated for all the systems. Using a cut-off threshold of 4 Å to indicate association leading to bond formation, it is concluded that the presence of the additional methoxy group on coniferyl alcohol assists in promoting interaction of the O4 group with the second moiety. Based on the computed data it is furthermore concluded that in aqueous solution, the most likely model of polymerization is that involving initial dimer formation, followed by dimer-dimer association.