Transition metal doping of ZnO is considered as a promising way to obtain a diluted magnetic semiconducting oxide. In this work we investigate copper doping of ZnO by means of density functional theory, using a hybrid exchange-correlation functional and a periodic approach with localized atomic basis functions. Isolated copper species, such as copper substitutional to zinc, Cu s , and Cu interstitial, Cu i , are analyzed in terms of transition energy levels and hyperfine coupling constants with reference to available spectroscopic data. We also examine the potential magnetic interaction between copper species, their interaction with oxygen vacancies, and the possibility of copper clustering. The relative stability of the various copper impurities considered in this study is finally compared on the basis of their formation energy at different oxygen chemical potentials and Fermi level values.