Human homogentisate dioxygenase is an Fell-dependent enzyme responsible for aromatic ring cleavage. The mechanism of its catalytic reaction has been investigated with the hybrid density functional method B3LYP. A relatively big model of the active site was first used to determine the substrate binding mode. It was found that binding of the substrate dianion with a vacant position trans to Glu341 is most favorable. The model was then truncated to include only the most relevant parts of the active-site residues involved in iron coordination and substrate binding. Thus, methylimidazole was used to model His292, His335, His365, and His371, while propionate modeled Glu341. The computational results suggest that the catalytic reaction of homogentisate dioxygenases involves three major chemical steps: formation of the peroxo intermediate, homolytic cleavage of the O-O bond leading to an arene oxide radical, and finally, cleavage of the six-membered ring. Calculated barriers for alternative reaction paths are markedly higher than for the proposed mechanism, and thus the computational results successfully explain the product specificity of the enzyme. Interestingly, the results indicate that the type of ring scission, intra or extra with respect to the substituents coordinating to iron, is controlled by the barrier heights for the decay of the arene oxide radical intermediate.