Many HNO scavenging pathways exist to regulate its in vivo and in vitro biological and pharmacological activities,
including the involvement of numerous ferric heme proteins. Such reactions also build an important basis for HNO probe
development. However, mechanisms of HNO reactions with ferric heme proteins are largely unknown. A computational investigation was
performed to provide the first detailed pathways, using metmyoglobin and catalase as representative ferric heme proteins with
neutral and negatively charged axial ligands. Results well reproduced experimental barriers with an average error of 0.11
kcal/mol. The rate-limiting step was found to be the dissociation of the resting ligand or HNO coordination where there is no
resting ligand. Unlike the non-heme case, the reductive nitrosylation step for both heme proteins was found to be barrierless
proton-coupled electron transfer, providing the major thermodynamic driving force for the overall reaction. Origin of experimental
reactivity difference between metmyoglobin and catalase was revealed. Results will facilitate studies of other heme-based HNO
scavenging systems and probe development.