Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory
Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory
复制标题
在密度泛函理论中模拟配体对 CYP 450 结构的影响
DOI:
10.1021/acs.jpca.2c01783
复制
发表时间:
2022
期刊:
影响因子:
--
通讯作者:
Garashchuk, Sophya
中科院分区:
文献类型:
--
作者:
Dutra, Matthew;McElhenney, Shannon;Manley, Olivia;Makris, Tom;Rassolov, Vitaly;Garashchuk, Sophya
An improved understanding of the P450 structure is relevant to the development of biomimetic catalysts and inhibitors for controlled CH-bond activation, an outstanding challenge of synthetic chemistry. Motivated by the experimental findings of an unusually short Fe–S bond of 2.18 Å for the wild-type (WT) OleT P450 decarboxylase relative to a cysteine pocket mutant form (A369P), a computational model that captures the effect of the thiolate axial ligand on the iron–sulfur distance is presented. With the computational efficiency and streamlined analysis in mind, this model combines a cluster representation of the enzyme─40–110 atoms, depending on the heme and ligand truncation level─with a density functional theory (DFT) description of the electronic structure (ES) and is calibrated against the experimental data. The optimized Fe–S distances show a difference of 0.25 Å between the low and high spin states, in agreement with the crystallographic structures of the OleT WT and mutant forms. We speculate that this difference is attributable to the packing of the ligand; the mutant is bulkier due to an alanine-to-proline replacement, meaning that it is excluded from the energetically favored low-spin minimum because of steric constraints. The presence of pure spin-state pairs and the intersection of the low/high spin states for the enzyme model is indicative of the limitations of single-reference ES methods in such systems and emphasizes the significance of using the proper state when modeling the hydrogen atom transfer (HAT) reaction catalyzed by OleT. At the same time, the correct characterization of both the short and long Fe–S bonds within a small DFT-based model of 42 atoms paves the way for quantum dynamics modeling of the HAT step, which initiates the OleT decarboxylation reaction.