Molecular dynamics study of host–guest interactions in cyclodextrins: methodology and data analysis for a comparison with solution data and the solid-state structure

摘要

A general method is proposed to model the behavior of cyclodextrins (CDs) and of their inclusion compounds through energy minimizations and molecular dynamics (MD) simulations at a constant temperature. In this way, the formation of a host–guest compound is obtained starting from many trial geometries with the guest outside the CD cavity without any a priori assumption. The MD simulation results are analyzed through two functions: (i) the similarity maps of the root-mean-square distances between instantaneous conformations found in the MD runs to recognize different families of conformers; (ii) the pair distribution function PDF, yielding the probability density of finding appropriate atom pairs as a function of their distance at equilibrium. As an example, the inclusion compound formed by β-CD and (?)-menthol-β-d-glucoside is investigated. The lowest-energy inclusion compound is in good agreement with the results of single-crystal X-ray analysis, while at room temperature the MD runs show a closely similar arrangement with thermal fluctuations. In this case, the PDF between diagnostic hydrogen atoms of β-CD and of the guest molecule are fully consistent with the experimental NOE results obtained from NMR measurements in solution.