How Do Variants of Residues in the First Coordination Sphere, Second Coordination Sphere, and Remote Areas Influence the Catalytic Mechanism of Non-Heme Fe(II)/2-Oxoglutarate Dependent Ethylene-Forming Enzyme?

The ethylene-forming enzyme (EFE) is a Fe(II)/2-oxoglutarate (2OG) and l-arginine (l-Arg)-dependent oxygenase that primarily decomposes 2OG into ethylene while also catalyzing l-Arg hydroxylation. While the hydroxylation mechanism in EFE is similar to other Fe(II)/2OG-dependent oxygenases, the formation of ethylene is unique. Various redesign strategies have aimed to increase ethylene production in EFE, but success has been limited, highlighting the need for alternate approaches. It is crucial to incorporate an accurate and comprehensive description of the integrative and multidimensional effects of the protein environment to enhance the redesign strategy in metalloenzymes, particularly in EFE. This involves understanding the role of the second coordination sphere (SCS) and long-range (LR) interacting residues, correlated motions, electronic structure, intrinsic electric field (IntEF), as well as the stabilization of transition states and reaction intermediates. In this study, we employ a molecular dynamics-based quantum mechanics/molecular mechanics approach to examine the integrative effects of the protein environment on reactions catalyzed by EFE variants from the first coordination sphere (FCS, D191E), SCS (A198V and R171A) and LR (E215A). The study uncovers how substitutions at different positions in EFE similarly impact the ethylene-forming reaction while posing distinct effects on the hydroxylation reaction. Results predict the effect of the variants in controlling the 2OG coordination mode in the Fe(II) center. Specifically, the study suggests that D191E uniquely prefers transitioning from an off-line to an in-line 2OG coordination mode before dioxygen binding. However, studies on the 2OG flip in the presence of off-line approaching dioxygen and dioxygen binding in the D191E variant indicate that the 2OG flip might not be feasible in the 5C Fe(II) state. Calculations show the possibility of a hydrogen atom transfer (HAT)-assisted oxygen flip in EFE and its variants (other than D191E). MD simulations elucidate the characteristic dynamic change in the α7 region in the D191E variant that might contribute to its increased hydroxylation reaction. Results indicate the possibility of forming an in-line ferryl from the IM2 (Fe(III)-partial bond intermediate) in the D191E variant. This alternative pathway from IM2 may also exist in WT EFE and other variants, which are yet to be explored. The study also delineates the impact of substitutions on the electronic structure and IntEF. Overall, the calculations support the idea that understanding the integrative and multidimensional effects of the protein environment on the reactions catalyzed by EFE variants provides the basics for improved enzyme redesign protocols of EFE to increase ethylene production. The results of this study will also contribute to the development of alternate redesign strategies for other metalloenzymes.