Dispelling the effects of a sorceress in enzyme catalysis

Enzyme dynamics-a brief review

Debate has been simmering for some years regarding the importance of internal thermal motions of enzymes to catalysis. Recent developments in protein design may bring resolution of the more contentious points a little closer.

Mitochondrial Peroxiredoxin 3 Is Rapidly Oxidized and Hyperoxidized by Fatty Acid Hydroperoxides

Human peroxiredoxin 3 (HsPrx3) is a thiol-based peroxidase responsible for the reduction of most hydrogen peroxide and peroxynitrite formed in mitochondria. Mitochondrial disfunction can lead to membrane lipoperoxidation, resulting in the formation of lipid-bound fatty acid hydroperoxides (_LFA-OOHs) which can be released to become free fatty acid hydroperoxides (_fFA-OOHs). Herein, we report that HsPrx3 is oxidized and hyperoxidized by _fFA-OOHs including those derived from arachidonic acid and eicosapentaenoic acid peroxidation at position 15 with remarkably high rate constants of oxidation (>3.5 × 10⁷ M^-1s^-1) and hyperoxidation (~2 × 10⁷ M^-1s^-1). The endoperoxide-hydroperoxide PGG₂, an intermediate in prostanoid synthesis, oxidized HsPrx3 with a similar rate constant, but was less effective in causing hyperoxidation. Biophysical methodologies suggest that HsPrx3 can bind hydrophobic structures. Indeed, molecular dynamic simulations allowed the identification of a hydrophobic patch near the enzyme active site that can allocate the hydroperoxide group of _fFA-OOHs in close proximity to the thiolate in the peroxidatic cysteine. Simulations performed using available and herein reported kinetic data indicate that HsPrx3 should be considered a main target for mitochondrial _fFA-OOHs. Finally, kinetic simulation analysis support that mitochondrial _fFA-OOHs formation fluxes in the range of nM/s are expected to contribute to HsPrx3 hyperoxidation, a modification that has been detected in vivo under physiological and pathological conditions.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

A Theoretical Study of the Benzoylformate Decarboxylase Reaction Mechanism

Density functional theory calculations are used to investigate the detailed reaction mechanism of benzoylformate decarboxylase, a thiamin diphosphate (ThDP)-dependent enzyme that catalyzes the nonoxidative decarboxylation of benzoylformate yielding benzaldehyde and carbon dioxide. A large model of the active site is constructed on the basis of the X-ray structure, and it is used to characterize the involved intermediates and transition states and evaluate their energies. There is generally good agreement between the calculations and available experimental data. The roles of the various active site residues are discussed and the results are compared to mutagenesis experiments. Importantly, the calculations identify off-cycle intermediate species of the ThDP cofactor that can have implications on the kinetics of the reaction.