The nature of the bonding in model complexes of di-copper metalloenzymes has been analyzed by means of the electronic localization function (ELF) and by the quantum theory of atoms in molecules (QTAIM). The constrained space orbital variations (CSOV) approach has also been used. Density functional theory (DFT) and CASSCF calculations have been carried out on several models of tyrosinase such as...