Glutarate Hydroxylation by the Carbon Starvation-Induced Protein D: A Computational Study into the Stereo- and Regioselectivities of the Reaction

The carbon starvation-induced protein D (CsiD) is a recently characterized iron(II)/alpha-ketoglutarate-dependent oxygenase that activates a glutarate molecule as substrate at the C-2 position to exclusively produce (S)-2-hydroxyglutarate products. This selective hydroxylation reaction by CsiD is an important component of the lysine biodegradation pathway in Escherichia coli; however, little is known on the details and the origin of the selectivity of the reaction. So far, experimental studies failed to trap and characterize any short-lived catalytic cycle intermediates. As no computational studies have been reported on this enzyme either, we decided to investigate the chemical reaction mechanism of glutarate activation by an iron(IV)-oxo model of the CsiD enzyme. In this work, we present a density functional theory study on a large active site cluster model of CsiD and investigate the glutarate hydroxylation pathways by a high-valent iron(IV)-oxo species leading to (S)-2-hydroxyglutarate, (R)-2-hydroxyglutarate, and 3-hydroxyglutarate. In agreement with experimental observation, the favorable product channel leads to pro-S C-2-H hydrogen atom abstraction to form (S)-2-hydroxyglutarate. The reaction is stepwise with a hydrogen atom abstraction by an iron(IV)-oxo species followed by OH rebound from a radical intermediate. The work presented in this paper shows that despite the fact that the C-H bond strengths at the C-2 and C-3 positions of glutarate are similar in the gas phase, substrate binding and positioning guide the reaction to an enantioselective reaction process by destabilizing the hydrogen atom abstraction pathways for the pro-R C-2-H and C-3-H positions. Our studies predict the chemical properties of the iron(IV)-oxo species and its rate constants with glutarate and deuterated-glutarate. Moreover, the work shows little protein motions during the catalytic process, while the substrate entrance into the substrate binding pocket appears to be guided by three active site arginine residues that position the substrate for pro-S C2-H hydrogen atom abstraction. Finally, the calculations show that irrespective of the position of the substrate and what C-H bond is closest to the metal center, the lowest energy pathway is for a selective pro-S C-2-H hydrogen atom abstraction.