Optimization of the hydrated Cu(II)(N7-guanine) structures revealed a number of minima on the potential energy surface. For selected structures, energy decompositions together with the determination of electronic properties (partial charges and electron spin densities) were performed. In the complexes of guanine with the bare copper cation and that with the monoaqua ligated cation, an electron transfer from guanine to Cu(II) was observed, resulting in a Cu(I)-guanine(+) type of complex. Conformers with two aqua ligands are borderline systems characterized by a Cu partial charge of +0.7e and a similar value of the spin density (0.6e) localized on guanine. When tetracoordination of copper was achieved, only then the prevailing electron spin density is unambiguously localized on copper. The energetic preference of diaqua-Cu-(N7,O6-guanine) over triaqua-Cu-(N7-guanine) was found for the four-coordinate structures. However, the energy difference between these two conformations decreases with the number of water molecules present in the systems, and in complexes with five water molecules this preference is preserved only at DeltaG level where thermal and entropy terms are included.