WNDP (Wilson's disease protein) is a copper-transporting ATPase that plays an essential role in human physiology. Mutations in WNDP result in copper accumulation in tissues and cause a severe hepato-neurological disorder known as Wilson's disease. Several mutations were surmised to affect the nucleotide binding and hydrolysis by WNDP; however, how the nucleotides bind to normal and mutated WNDP remains unknown. To aid such studies, we performed the molecular modelling of the spatial structure and dynamics of the ATP-binding domain of WNDP and its interactions with ATP. The three-dimensional models of this domain in two conformations were built using the X-ray structures of the Ca2+-ATPase in the E1 and E2 states. To study the functional aspects of the models, they were subjected to long-term molecular dynamics simulations in an explicit solvent; similar calculations were performed for the ATP-binding domain of Ca2+-ATPase. In both cases, we found large-scale motions that lead to significant changes of distances between several functionally important residues. The ATP docking revealed two possible modes of ATP binding: via adenosine buried in the cleft near residues H1069, R1151 and D1164, and via phosphate moiety 'anchored' by H-bonds with residues in the vicinity of catalytic D1027. Furthermore, interaction of ATP with both sites occurs if they are spatially close to each other. This may be achieved after relative domain motions of the 'closure' type observed in molecular dynamics simulations. The results provide a framework for analysis of disease mutations and for future mutagenesis studies.