Molecular dynamics study of the recognition of ATP by nucleic acid aptamers

Despite their great success in recognizing small molecules in vitro, nucleic acid aptamers are rarely used in clinical settings. This is partially due to the lack of structure-based mechanistic information. In this work, atomistic molecular dynamics simulations are used to study the static and dynamic supramolecular structures relevant to the process of the wild-type (wt) nucleic acid aptamer recognition and binding of ATP. The effects brought about by mutation of key residues in the recognition site are also explored. The simulations reveal that the aptamer displays a high degree of rigidity and is structurally very little affected by the binding of ATP. Interaction energy decomposition shows that dispersion forces from π-stacking between ATP and the G6 and A23 nucleobases in the aptamer binding site plays a more important role in stabilizing the supramolecular complex, compared to hydrogen-bond interaction between ATP and G22. Moreover, metadynamics simulations show that during the association process, water molecules act as essential bridges connecting ATP with G22, which favors the dynamic stability of the complex. The calculations carried out on three mutated aptamer structures confirm the crucial role of the hydrogen bonds and π-stacking interactions for the binding affinity of the ATP nucleic acid aptamer.