Transient diquat inclusion and carboxylate-rim binding to carboxymethyl-α-cyclodextrin: a molecular dynamics study

Rotaxane-based supramolecular scavengers have shown promise for diquat detoxification, but the relative contributions of cavity inclusion and interactions at the macrocycle rim remain unclear. We used all-atom molecular dynamics simulations in explicit water to examine diquat binding to carboxymethyl-α-cyclodextrin (CM-α-CD) and the behaviour of CM-α-CD on a poly(ethylene glycol) axle. Five systems were simulated: host-guest mixtures containing a neutral host (pH 3 model) or a hexaanionic host (pH 10 model), an explicitly constructed inclusion complex, and two-macrocycle pseudopolyrotaxanes initiated from pre-threaded or random configurations. The total production sampling was approximately 4.5 μs. A relaxed threaded geometry was obtained, with the diquat centroid 0.03 nm from the cavity axis and 9 of its 14 heavy atoms within the defined cavity cylinder. The guest left this geometry within the first 50 ps sampling interval and subsequently remained associated mainly with the carboxylate rim. During 1.2 μs of unbiased sampling of the hexaanionic host-guest system, only 5 of 802 guest-frame pairs contained a guest atom within the cavity cylinder, and all of these events occurred near a portal. Rim contacts were more frequent: each guest formed an average of 3.8-4.4 host heavy-atom contacts per frame, with 12.8-14.0% of guest heavy atoms in contact with the host. These interactions were substantially reduced for the neutral host. Sodium ions remained predominantly solvent-separated from the carboxylates. In the pseudopolyrotaxane simulations, pre-threaded rings satisfied the threading criterion in 1.6-8.4% of frames, and no spontaneous threading was observed from random configurations during 1 μs. Under the simulated conditions, CM-α-CD therefore bound diquat primarily through dynamic rim interactions rather than persistent cavity inclusion.