Abstract
The cytochromes P450 constitute a ubiquitous family of metalloenzymes, catalyzing manifold reactions of biological and synthetic importance via a thiolate-ligated iron-oxo (IV) porphyrin radical species denoted compound I (Cpd I). Experimental investigations have implicated this intermediate in a broad spectrum of biophysically interesting phenomena, further augmenting the importance of a Cpd I model system. Ab initio molecular dynamics, including Car-Parrinello and path integral methods, conjoin electronic structure theory with finite temperature simulation, affording tools most valuable to approach such enzymes. These methods are typically driven by density functional theory (DFT) in a plane-wave pseudopotential framework; however, existing studies of Cpd I have been restricted to localized Gaussian basis sets. The appropriate choice of density functional and pseudopotential for such simulations is accordingly not obvious. To remedy this situation, a systematic benchmarking of thiolate-ligated Cpd I is performed using several generalized-gradient approximation (GGA) functionals in the Martins-Troullier and Vanderbilt ultrasoft pseudopotential schemes. The resultant electronic and structural parameters are compared to localized-basis DFT calculations using GGA and hybrid density functionals. The merits and demerits of each scheme are presented in the context of reproducing existing experimental and theoretical results for Cpd I.