This work presents a systematic theoretical study on Cu(I) and Cu(II) cations in variable hydrogen sulfide-aqua-aminine ligand fields. These ligands model the biologically most common environment for Cu ions. Molecular structures of the complexes were optimized at the density functional theory (DFF) level. Subsequent thorough energy analyses revealed the following trends: (i) The ammine complexes are the most stable, followed by those containing the aqua and hydrogen sulfide ligands, which are characterized by similar stabilization energies. (ii) The most preferred Cu(I) coordination number is 2 in ammine or aqua ligand fields. A qualitatively different binding picture was obtained for complexes with H2S ligands where the 4-coordination is favored. (iii) The 4- and 5-coordinated structures belong to the most stable complexes for Cu(II), regardless of the ligand types. Vertical and adiabatic ionization potentials of Cu(I) complexes were calculated. Charge distribution (using the natural population analysis (NPA) method) and molecular orbital analyses were performed to elucidate the nature of bonding in the examined systems. The results provide in-depth insight into the Cu-binding properties and can be, among others, used for the calibration of bioinorganic force fields.