Free energy perturbation and MM/PBSA studies on inclusion complexes of some structurally related compounds with β-cyclodextrin

Abstract

Molecular dynamics (MD) simulations have been conducted to explore time-resolved guest–host interactions involving inclusion complex formation between β-cyclodextrin and organic molecules bearing two peripheral benzene rings in aqueous solution. Moreover, free energy perturbation (FEP) and thermodynamic integration (TI) methods at different simulation times have been employed to estimate the relative free energy of complexation. Also, the less computer-time demanding molecular mechanics/Poisson–Boltzmann surface area (MM/PBSA) method was used to estimate the free energy of complexation based on only 1-ns MD simulation. Results showed that both FEP and TI methods were able to reasonably reproduce the experimental thermodynamic quantities. However, long simulation times (e.g. 15 ns) were needed for benzoin mutating to benzanilide (BAN), while moderately shorter times were sufficient for BAN mutating to phenyl benzoate and for benzilic acid mutating to diphenylacetic acid. The results have been discussed in the light of the differences in the chemical structural and conformational features of the guest molecules. In general, it was apparent that the TI method requires less time for convergence of results than the FEP method. However, the less expensive MM/PBSA method proved capable of producing results that are in agreement with those of the more expensive TI and FEP methods.