Computational modeling of factors that modulate the unique FeNO bonding in {FeNO}(6) heme-thiolate model complexes

A density functional theory account of the changes in FeNO bonding that occur in response to both bonded and nonbonded structural perturbations is reported for a series of {FeNO}(6) heme-thiolate model complexes. Using [Fe(porphine)(SCH(3))NO] as the reference complex, we constructed models to mimic equatorial (cis), distal, and proximal influences of protein environments. Overall, the results from these calculations reveal that the Fe-NO and N-O bond strengths change in the same direction upon variations in structure and environment. These bonding changes are manifested in unique direct correlations between the Fe-NO and N-O vibrational frequencies and bond lengths, as evidenced by their positive slopes (slopes of the familiar inverse or backbonding correlations are negative). The electronic origin of the direct correlations appears to derive from the electron density distribution in high-energy molecular orbitals. This variability modulates the FeNO antibonding character throughout the triatomic FeNO moiety. The results of this study suggest that the stabilities and reactivities of {FeNO}(6) centers in heme-thiolate enzymes can be modulated over a significant range through a variety of bonded and nonbonded means.