Density functional molecular dynamics simulations using a QM/MM approach are used to get insight into the binding modes of for-mic acid in catalase. Two ligand binding sites are found, named A and B, in agreement with recent high resolution structures of catalase with bound formic acid. In addition, the calculations show that the His56 residue is protonated and the ligand is present as a formate anion. The lowest energy minimum structure (A) corresponds to the ligand interacting with both the heme iron and the catalytic residues (His56 and Asn129). The second minimum energy structure (B) corresponds to the situation in which the ligand interacts solely with the catalytic residues. A mechanism for the process of formic acid binding in catalase is suggested.