Tuning Thermodynamics and Electronic Properties of 1-Hexyl-3-Methylimidazolium Based Organic Ionic Liquids with Iron Porphyrin

Ionic liquid-iron porphyrin (IL-FeP) complexes are systems with diverse applications, particularly in oxygen reduction reactions and the biodegradation of ILs by cytochrome P450 enzymes, where FeP serves as the active site. Despite the importance of such systems, there is currently a lack of information on conformational preference for ILs binding to FeP and the influence such binding conformations exert on the electronic properties of FeP. Additionally, the role of binding interactions in facilitating biodegradation of ionic liquids is not yet fully understood. In this article, we address the knowledge gap by employing density functional theory calculations to identify the most stable structures of iIL-FeP complexes. We considered four ILs containing the cation 1-hexyl-3-methylimidazolium paired with acetate, methanesulfonate, ethylsulfate, and butylsulfate anions. Thermodynamic analysis of binding was supplemented with energy decomposition analysis (EDA), non-covalent interactions (NCI) analysis, and natural bond orbital (NBO) analysis to glean insight into the energetic factors driving binding stability. Our results reveal conformations in which one of the oxygen atoms in anions is positioned directly above Fe and the imidazolium ring aligns approximately parallel to the FeP plane exhibit strong binding energies, irrespective of the anion. Electrostatic interactions between the IL and FeP play a dominant role in determining the binding conformations. The propensity to acquire an electron by an IL-FeP complex is significantly reduced when oxygen in the anion is in the close proximity of Fe. These findings enhance our understanding of the fundamental interactions governing IL-FeP complexes, providing a foundation for tailoring ILs for applications in catalysis and biodegradation.