Theoretical investigation of the catalytic mechanism of the protein arginine deiminase 4 enzyme

What occurs by replacing Mn2+ with Co2+ in human arginase I: first-principles computational analysis

The reaction mechanism of the dinuclear cobalt enzyme arginase is investigated using density functional theory. As an arginase-containing binuclear Mn(2)(2+) cluster, it catalyzes the hydrolysis of L-arginine in L-ornithine and urea. The bridging hydroxide is capable of performing nucleophilic attack on the iminium carbon ion from its bridging position, resulting in the formation of a tetrahedral intermediate, as was already obtained in a previous theoretical study on the manganese enzyme. Our theoretical investigation allows us to obtain an accurate potential energy profile and confirms that the coordination mode of the substrate to the dimetallic center is quite similar to that present in the manganese enzyme. In agreement with the experimental observations, our results show that both Mn- and Co-containing enzymes catalyze the same reaction with quite comparable energy barriers.

On the inhibitor effects of bergamot juice flavonoids binding to the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) enzyme

Density functional theory was applied to study the binding mode of new flavonoids as possible inhibitors of the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), an enzyme that catalyzes the four-electron reduction of HMGCoA to mevalonate, the committed step in the biosynthesis of sterols. The investigated flavonoid conjugates brutieridin and melitidin were recently quantified in the bergamot fruit extracts and identified to be structural analogues of statins, lipids concentration lowering drugs that inhibit HMGR. Computations allowed us to perform a detailed analysis of the geometrical and electronic features affecting the binding of these compounds, as well as that of the excellent simvastatin drug, to the active site of the enzyme and to give better insight into the inhibition process.

Human insulin-degrading enzyme working mechanism

The possible mechanism by which the insulin-degrading enzyme (IDE) zinc-binding protease carries out its catalytic function toward two peptides of different length, simulating a portion of B chain of insulin, was investigated on an enzymatic model consisting of 130 /159 atoms, using the density functional theory method and the hybrid exchange-correlation functional B3LYP in gas phase and in the protein environment. Based on the geometry and relative stabilities of minima and transition states on the potential energy profiles, we determined that proteolysis reaction is exothermic and proceeds quickly as the barrier in the rate-limiting step falls widely within the range of values expected for an enzymatic catalysis, both in vacuum and in protein medium.

Active site cysteine is protonated in the PAD4 Michaelis complex: evidence from Born-Oppenheimer ab initio QM/MM molecular dynamics simulations

The protein arginine deiminase 4 (PAD4) catalyzes the citrullination of the peptidylarginine and plays a critical role in rheumatoid arthritis (RA) and gene regulation. Understanding its catalytic mechanism is not only of fundamental importance but also of significant medical interest for the rational design of new inhibitors. By employing on-the-fly Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, we have demonstrated that it is unlikely for the active site cysteine and histidine to exist as a thiolate-imidazolium ion pair in the PAD4 Michaelis reactant complex. Instead, a substrate-assisted proton transfer mechanism for the deimination reaction step has been characterized: both Cys645 and His471 in the PAD4 active site are neutral prior to the reaction; the deprotonation of Cys645 by the substrate arginine occurs in concert with the nucleophilic addition of the Cys thiolate to Czeta of the substrate, and leads to a covalent tetrahedral intermediate; then, the Czeta-Neta1 bond cleaves and the resulted ammonia is displaced by a solvent water molecule. The initial deprotonation and nucleophilic attack step is found to be rate-determining. The computed free energy barrier with B3LYP(6-31G*) QM/MM MD simulations and umbrella sampling is 20.9 kcal.mol(-1), consistent with the experimental kinetic data. During the deimination, His471 plays an important role in stabilizing the transition state through the formation of the hydrogen bond with the guanidinium group. Our current studies further demonstrated the viability and strength of the ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.