The study of inclusion complexes of Chrysin (ChR) with three forms of cyclodextrins (CDs) α-, β-, and γ-CD was accomplished to examine the stability of ChR inside the central cavities of CDs. The aim of study was to identify the most suitable form of CD to improve the hydro-solubility of poorly soluble ChR bioactive molecule. Microsecond timescale molecular dynamics (MD) simulations were performed on four inclusion complexes (α-CD/ChR, β-CD/ChR, and two conformations of γ-CD/ChR) to examine the dynamics of ChR inside the cavity of CDs. The first conformation of γ-CD/ChR inclusion complex (γ-CD1/ChR) was identified to possess the highest affinity between host and guest molecule on the basis of binding energy calculated by employing Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and umbrella sampling simulations. To further strengthen the claims of classical and biased MD studies, Our own N-layered Integrated molecular Orbital and Molecular mechanics (ONIOM) (wB97XD/6-311+g(d,p):pm7) calculations were performed on the selected inclusion complexes. The ONIOM based complexation energy reaffirmed that ChR had highest affinity for the γ-CD1 host molecule. Further, the non-covalent interaction analysis was conducted using Multiwfn software on QM-optimized inclusion complexes with wB97XD/6-311+G(d,p) model chemistry, revealing non-covalent interactions between ChR and CDs. This atomic level information helped us to gain better insights into critical atoms of ChR and CD that participated in intermolecular interactions and identify γ-CD as a suitable host molecule for improving the hydro-solubulity of ChR. The structural insights would help to derive new derivatives of γ-CD with better host capacity.