Computational and experimental analysis of Luteolin-β-cyclodextrin supramolecular complexes: Insights into conformational dynamics and phase solubility

Investigating the structural stability of poorly-soluble luteolin (LuT) after encapsulation within cyclodextrins (CDs) is crucial for unlocking the therapeutic potential of LuT bioactive molecule. Herein, native and modified β-CD were employed to investigate LuT inclusion complex formation. Molecular mechanics (MM) and quantum mechanics (QM) were utilized for structural dynamics analysis. Microsecond timescale MD simulations yielded insights into LuT-CD interactions. The binding affinity between LuT and selected β-CDs was assessed by calculating the binding free energy using MM-PBSA and umbrella sampling simulations. The MM-PBSA results indicated that Heptakis-O-(2-hydroxypropyl)-β-CD (HP-β-CD) (-82.59+/-11.67 kJ/mol) and Di-O-methyl-β-CD (DM-β-CD) (-54.01+/-11.07 kJ/mol) exhibited good binding affinity for LuT. Subsequently, derivative screening of HP-β-CD revealed that only 2-HP-β-CD (HP-β-CD-1)/LuT (-21.38 kJ/mol) displayed a superior binding free energy (obtained from umbrella sampling) than HP-β-CD/LuT (-16.55 kJ/mol) inclusion complex. We conducted QM calculations on the top three inclusion complexes namelly HP-β-CD, DM-β-CD, and HP-β-CD-1 employing wB97X-D/6-311 + G(d,p) model chemistry to strengthen the MM results. The computational analysis aligns with experimental findings (phase solubility analysis), validating HP-β-CD-1 as most effective cavitand molecule for improving the solubility of LuT. This study offers critical structural insights for developing novel HP-β-CD derivatives with enhanced host capacity to encapsulate guest molecules efficiently.