How Do Electrostatic Perturbations of the Protein Affect the Bifurcation Pathways of Substrate Hydroxylation versus Desaturation in the Nonheme Iron-Dependent Viomycin Biosynthesis Enzyme?
Hafiz Saqib AliRichard H HenchmanJim WarwickerSam P De VisserPublished in: The journal of physical chemistry. A (2021)
The viomycin biosynthesis enzyme VioC is a nonheme iron and α-ketoglutarate-dependent dioxygenase involved in the selective hydroxylation of l-arginine at the C3-position for antibiotics biosynthesis. Interestingly, experimental studies showed that using the substrate analogue, namely, l-homo-arginine, a mixture of products was obtained originating from C3-hydroxylation, C4-hydroxylation, and C3-C4-desaturation. To understand how the addition of one CH2 group to a substrate can lead to such a dramatic change in selectivity and activity, we decided to perform a computational study using quantum mechanical (QM) cluster models. We set up a large active-site cluster model of 245 atoms that includes the oxidant with its first- and second-coordination sphere influences as well as the substrate binding pocket. The model was validated against experimental work from the literature on related enzymes and previous computational studies. Thereafter, possible pathways leading to products and byproducts were investigated for a model containing l-Arg and one for l-homo-Arg as substrate. The calculated free energies of activation predict product distributions that match the experimental observation and give a low-energy C3-hydroxylation pathway for l-Arg, while for l-homo-Arg, several barriers are found to be close in energy leading to a mixture of products. We then analyzed the origins of the differences in product distributions using thermochemical, valence bond, and electrostatic models. Our studies show that the C3-H and C4-H bond strengths of l-Arg and l-homo-Arg are similar; however, external perturbations from an induced electric field of the protein affect the relative C-H bond strengths of l-Arg dramatically and make the C3-H bond the weakest and guide the reaction to a selective C3-hydroxylation channel. Therefore, the charge distribution in the protein and the induced electric dipole field of the active site of VioC guides the l-Arg substrate activation to C3-hydroxylation and disfavors the C4-hydroxylation pathway, while this does not occur for l-homo-Arg. Tight substrate positioning and electrostatic perturbations from the second-coordination sphere residues in VioC also result in a slower overall reaction for l-Arg; however, they enable a high substrate selectivity. Our studies highlight the importance of the second-coordination sphere in proteins that position the substrate and oxidant, perturb charge distributions, and enable substrate selectivity.