Molecular docking predicts groove binding and electrostatic interaction well for the studied copper nucleases.

The DNA binding affinity of copper nucleases influenced by the ligand size, length, functional groups, chelate ring size bound to copper center.

Intercalation modes could be reproduced by “distorted DNA” formed by molecular dynamics simulations.

MM-PBSA approach validated the DNA binding affinity of copper nuclease.