The Keggin anion, PW₁₂O₄₀^3-, is one of the most representative polyoxometalates (POMs). In recent years increasing theoretical work focused on this family of compounds has explained or even predicted some of their properties using quantum mechanics methods. In this report we applied for the first time molecular dynamics (MD) to the title compound to analyze its interactions with water. We used three atomic charge definitions (Mulliken, ChelpG, and formal charges) to carry out MD simulations. The results show that the terminal oxygens of the cluster are invariably most effectively solvated by water because of their prominent position within the framework. On the other hand, bridging oxygens, which are confined in more internal positions, concentrate a smaller portion of the whole solvation. We investigated the hydrogen bonds existing between water and the cluster, confirming that the terminal positions form more contacts with H₂O than any other site of the cluster. In the end, the lifetime of such contacts is longer with bridging oxygens, presumably due to their higher atomic charge.