Study of orotidine 5′-monophosphate decarboxylase in complex with the top three OMP, BMP, and PMP ligands by molecular dynamics simulation

Catalytic mechanism of orotidine 5′-monophosphate decarboxylase (OMPDC), one of the nature most proficient enzymes which provides large rate enhancement, has not been fully understood yet. A series of 30 ns molecular dynamics (MD) simulations were run on X-ray structure of the OMPDC from Saccharomyces cerevisiae in its free form as well as in complex with different ligands, namely 1-(5′-phospho-D-ribofuranosyl) barbituric acid (BMP), orotidine 5′-monophosphate (OMP), and 6-phosphonouridine 5′-monophosphate (PMP). The importance of this biological system is justified both by its high rate enhancement and its potential use as a target in chemotherapy. This work focuses on comparing two physicochemical states of the enzyme (protonated and deprotonated Asp91) and three ligands (substrate OMP, inhibitor, and transition state analog BMP and substrate analog PMP). Detailed analysis of the active site geometry and its interactions is properly put in context by extensive comparison with relevant experimental works. Our overall results show that in terms of hydrogen bond occupancy, electrostatic interactions, dihedral angles, active site configuration, and movement of loops, notable differences among different complexes are observed. Comparison of the results obtained from these simulations provides some detailed structural data for the complexes, the enzyme, and the ligands, as well as useful insights into the inhibition mechanism of the OMPDC enzyme. Furthermore, these simulations are applied to clarify the ambiguous mechanism of the OMPDC enzyme, and imply that the substrate destabilization and transition state stabilization contribute to the mechanism of action of the most proficient enzyme, OMPDC.