The mechanism for the toxic superoxide radical disproportionation to molecular oxygen and hydrogen peroxide by copper-zinc superoxide dismutase (CuZnSOD) has been studied using the B3LYP hybrid density functional. On the basis of the X-ray structure of the enzyme, the molecular system investigated includes the first-shell protein ligands of the two metal centers as well as the second-shell ligand Asp122. The substrates of the model reaction are two superoxide radical anions, approaching the copper center at the beginning of two half-reactions: the first part of the catalytic cycle involving Cu+ oxidation and the second part reducing Cu2+ back to its initial state. The quantitative free energy profile of the reaction is obtained and discussed in connection with the experimental data on the reduction potentials and CuZnSOD kinetics. The optimized structures are analyzed and compared to the experimental ones. The two transition states alternate the protonation state of His61 and correspond to histidine Cu-His61-Zn bridge rupture/reformation. Modifications applied to the initial model allow the importance of Asp122 for catalysis to be estimated.