Why does β-cyclodextrin prefer to bind nucleotides with an adenine base rather than other 2'-deoxyribonucleoside 5'-monophosphates?

β-Cyclodextrin (β-CD), which resides in the α-hemolysin (αHL) protein pore, can act as a molecular adapter in single-molecule exonuclease DNA sequencing approaches, where the different nucleotide binding behavior of β-CD is crucial for base discrimination. In the present contribution, the inclusion modes of β-CD towards four 2'-deoxyribonucleoside 5'-monophosphates (dNMPs) were investigated using quantum mechanics (QM) calculations. The calculated binding energy suggests that the binding affinity of dAMP to β-CD are highest among all the dNMPs in solution, in agreement with experimental results. Geometry analysis shows that β-CD in the dAMP complex undergoes a small conformational change, and weak interaction analysis indicates that there are small steric repulsion regions in β-CD. These results suggest that β-CD has lower geometric deformation energy in complexation with dAMP. Furthermore, topological analysis and weak interaction analysis suggest that the number and strength of intermolecular hydrogen bonds and van der Waals interactions are critical to dAMP binding, and they both make favorable contributions to the lower interaction energy. This work reveals the reason why β-CD prefers to bind dAMP rather than other dNMPs, while opening exciting perspectives for the design of novel β-CD-based molecular adapters in the single-molecule exonuclease method of sequencing DNA. Graphical Abstract The binding affinity of β-cyclodextrin towards four 2'-deoxyribonucleoside 5'-monophosphates was investigated using quantum mechanics calculations.