A variety of theoretical methods including classical molecular interaction potentials, classical molecular dynamics, and activated molecular dynamics have been used to analyze the substrate recognition mechanisms of peroxisomal catalase from Saccharomyces cerevisiae. Special attention is paid to the existence of channels connecting the heme group with the exterior of the protein. On the basis of these calculations a rationale is given for the unique catalytic properties of this enzyme, as well as for the change in enzyme efficiency related to key mutations. According to our calculations the water is expected to be a competitive inhibitor of the enzyme, blocking the access of hydrogen peroxide to the active site. The main channel is the preferred route for substrate access to the enzyme and shows a cooperative binding to hydrogen peroxide. However, the overall affinity of the main channel for H(2)O(2) is only slightly larger than that for H(2)O. Alternative channels connecting the heme group with the monomer interface and the NADP(H) binding site are detected. These secondary channels might be important for product release.