Elucidating esterification reaction during deposition of cutin monomers from classical molecular dynamics simulations

The structural behavior of some cutin monomers, when deposited on mica support, was extensively investigated by our research group. However, other events, such as esterification reaction (ER), are still a way to explore. In this paper, we explore possible ER that could occur when these monomers adsorb on support. Although classical molecular dynamics simulations are not able to capture reactive effects, here, we show that they become valuable strategies to analyze the initial structural configurations to predict the most favorable reaction routes. Thus, when depositing aleuritic acid (ALE), it is observed that the loss of capacity to form self-assembled (SA) systems favors different routes to occur ER. In pure ALE bilayers systems, an ER is given exclusively through the -COOH and primary -OH groups. In pure ALE monolayers systems, the ER does not happen when the system is self-assembled. However, for disorganized systems, it is able to occur by two possible routes: -COOH and primary -OH (route 1) and -COOH and secondary -OH (route 2). When palmitic acid (PAL) is added in small quantities, ALE SAMs can now form an ER. In this case, ER occurs mostly through the -COOH and secondary -OH groups. However, when the presence of PAL is dominant, ER can occur with either of both possibilities, that is, routes 1 and 2. Graphical abstract.