Conformational analysis of hydroxymatairesinol in aqueous solution with molecular dynamics

Molecular dynamics simulations were performed on the naturally occuring lignan hydroxymatairesinol (HMR) using the GROMACS software. The aim of this study was to explore the conformational behavior of HMR in aqueous solution adopting the TIP4P model. The topology of HMR was constructed by hand and HMR was modeled with the OPLS-AA force field implemented in GROMACS. The five torsional angles in HMR were properly analyzed during the simulations. Correlations through certain patterns were observed between the angles. The determining property for the conformation preferred in aqueous solution was found to be the dipole moment and not the lowest energy in gas phase. The solvation effects on HMR was also studied by quantum chemical calculations applying the COnductorlike Screening MOdel (COSMO), the results of which were compared with results from a previous study using the Polarized Continuum Model (PCM). In the present work, COSMO was found to give more credible relative energies than PCM.