Computational Insights into Cyclodextrin Inclusion Complexes with the Organophosphorus Flame Retardant DOPO

Cyclodextrins (CDs) were used as green char promoters in the formulation of organophosphorus flame retardants (OPFRs) for polymeric materials, and they could reduce the amount of usage of OPFRs and their release into the environment by forming [host:guest] inclusion complexes with them. Here, we report a systematic study on the inclusion complexes of natural CDs (α-, β-, and γ-CD) with a representative OPFR of DOPO using computational methods of molecular docking, molecular dynamics (MD) simulations, and quantum mechanical (QM) calculations. The binding modes and energetics of [host:guest] inclusion complexes were analyzed in details. α-CD was not able to form a complete inclusion complex with DOPO, and the center of mass distance [host:guest] distance amounted to 4-5 Å. β-CD and γ-CD allowed for a deep insertion of DOPO into their hydrophobic cavities, and DOPO was able to frequently change its orientation within the γ-CD cavity. The energy decomposition analysis based on the dispersion-corrected density functional theory (sobEDAw) indicated that electrostatic, orbital, and dispersion contributions favored [host:guest] complexation, while the exchange-repulsion term showed the opposite. This work provides an in-depth understanding of using CD inclusion complexes in OPFRs formulations.