Complexation Mechanisms of Aqueous Amylose: Molecular Dynamics Study Using 3-Pentadecylphenol

Molecular dynamics (MD) simulations of linear amylose fragments containing 10 to 40 glucose units were used to study the complexation of the prototypical compound, 3-pentadecylphenol (PDP)─a natural product with surfactant-like properties─in aqueous solution. The amylose-PDP binding leverages mainly hydrophobic interactions together with excluded volume effects. It was found that while the most stable complexes contained PDP inside the helical structure of the amylose in the expected guest-host (inclusion) complexation manner, at higher temperatures, the commonly observed PDP-amylose complexes often involved more nonspecific interactions than inclusion complexation. In the case where a stoichiometric excess of PDP was added to the simulation box, self-aggregation of the small molecule precluded its ability to enter the internal helical part of the oligosaccharide, and as a result, inclusion complexation became ineffective. MD simulation trajectories were analyzed preliminarily using cluster analysis (CA), followed by more rigorous solvent accessible surface area (SASA) determination over the temperature range spanning from 277 to 433 K. It was found that using the SASA of PDP corrected for its intrinsic conformational changes, together with a generic hidden Markov model (HMM), an adequate quantification of the different types of PDP-amylose aggregates was obtained to allow further analysis. The enthalpy change associated with the guest-host binding equilibrium constant (K_gh) in aqueous solution was estimated to be -75 kJ/mol, which is about twice as high as one might expect based on experimentally measured values of similar complexes in the solid state where the (unsolvated) helical structure of amylose remains rigid. On the other hand, the nonspecific binding (K_ns) enthalpy change associated with PDP-amylose interactions in the same solution environment was found to be about half of the inclusion complexation value.