Review of computer simulations of isotope effects on biochemical reactions: From the Bigeleisen equation to Feynman's path integral

Estimation of frequency factors for the calculation of kinetic isotope effects from classical and path integral free energy simulations

We use the modified Bigeleisen-Mayer equation to compute kinetic isotope effect values for non-enzymatic phosphoryl transfer reactions from classical and path integral molecular dynamics umbrella sampling. The modified form of the Bigeleisen-Mayer equation consists of a ratio of imaginary mode vibrational frequencies and a contribution arising from the isotopic substitution's effect on the activation free energy, which can be computed from path integral simulation. In the present study, we describe a practical method for estimating the frequency ratio correction directly from umbrella sampling in a manner that does not require normal mode analysis of many geometry optimized structures. Instead, the method relates the frequency ratio to the change in the mass weighted coordinate representation of the minimum free energy path at the transition state induced by isotopic substitution. The method is applied to the calculation of 16/18O and 32/34S primary kinetic isotope effect values for six non-enzymatic phosphoryl transfer reactions. We demonstrate that the results are consistent with the analysis of geometry optimized transition state ensembles using the traditional Bigeleisen-Mayer equation. The method thus presents a new practical tool to enable facile calculation of kinetic isotope effect values for complex chemical reactions in the condensed phase.

Altered Mechanisms for Acid-Catalyzed RNA Cleavage and Isomerization Reactions Models

RNA strand cleavage by 2'-O-transphosphorylation is catalyzed not only by numerous nucleolytic RNA enzymes (ribozymes) but also by hydroxide or hydronium ions. In experiments, both cleavage of the 5'-linked nucleoside and isomerization between 3',5'- and 2',5'-phosphodiesters occur under acidic conditions, while only the cleavage reaction is observed under basic conditions. An ab initio path-integral approach for simulating kinetic isotope effects is used to reveal the reaction mechanisms for RNA cleavage and isomerization reactions under acidic conditions. Moreover, the proposed mechanisms can also be combined through the experimental pH-rate profiles.

Platinum on Carbon-Catalysed Site-Selective H-D Exchange Reaction of Allylic Alcohols Using Alkyl Amines as a Hydrogen Source

We developed platinum on carbon-catalysed deuteration reaction of tert-allylic alcohols using deuterium oxide as a deuterium source. Amylamine was dehydrogenated by platinum on carbon to generate an appropriate amount of...

Racemases and epimerases operating through a 1,1-proton transfer mechanism: reactivity, mechanism and inhibition

Racemases and epimerases catalyse changes in the stereochemical configurations of chiral centres and are of interest as model enzymes and as biotechnological tools. They also occupy pivotal positions within metabolic pathways and, hence, many of them are important drug targets. This review summarises the catalytic mechanisms of PLP-dependent, enolase family and cofactor-independent racemases and epimerases operating by a deprotonation/reprotonation (1,1-proton transfer) mechanism and methods for measuring their catalytic activity. Strategies for inhibiting these enzymes are reviewed, as are specific examples of inhibitors. Rational design of inhibitors based on substrates has been extensively explored but there is considerable scope for development of transition-state mimics and covalent inhibitors and for the identification of inhibitors by high-throughput, fragment and virtual screening approaches. The increasing availability of enzyme structures obtained using X-ray crystallography will facilitate development of inhibitors by rational design and fragment screening, whilst protein models will facilitate development of transition-state mimics.

Silver-catalyzed regioselective deuteration of (hetero)arenes and α-deuteration of 2-alkyl azaarenes

A simple silver-catalyzed regioselective deuteration of (hetero)arenes and α-deuteration of 2-alkyl azaarenes has been described. This strategy provides an efficient and practical avenue to access various deuterated electron-rich arenes, azaarenes and α-deuterated 2-alkyl azaarenes with good to excellent deuterium incorporation utilizing D₂O as the source of deuterium atoms.

A practical silver-catalyzed regioselective deuteration of (hetero)arenes and α-deuteration of 2-alkyl azaarenes utilizing D₂O as a deuterium source has been developed.

Transition-Metal-Free, Selective Reductive Deuteration of Terminal Alkynes with Sodium Dispersions and EtOD- d1

A transition-metal-free single electron transfer reaction has been developed for the synthesis of [D₃]-alkenes from terminal alkynes using sodium dispersions as the electron donor and EtOD- d₁ as the deuterium source. Both reagents are cost-effective and bench-stable. This practical method exhibits remarkable terminal alkyne selectivity and exclusive alkene selectivity. Excellent deuterium incorporations and yields were achieved across a broad range of terminal alkynes without olefin isomerization. Of note, this reaction is highly solvent dependent. n-Hexane provides unique enhancement to this reductive deuteration process.

13C kinetic isotope effects on the reaction of a flavin amine oxidase determined from whole molecule isotope effects

A large number of flavoproteins catalyze the oxidation of amines. Because of the importance of these enzymes in metabolism, their mechanisms have previously been studied using deuterium, nitrogen, and solvent isotope effects. While these results have been valuable for computational studies to distinguish among proposed mechanisms, a measure of the change at the reacting carbon has been lacking. We describe here the measurement of a 13C kinetic isotope effect for a representative amine oxidase, polyamine oxidase. The isotope effect was determined by analysis of the isotopic composition of the unlabeled substrate, N, N'-dibenzyl-1,4-diaminopropane, to obtain a pH-independent value of 1.025. The availability of a 13C isotope effect for flavoprotein-catalyzed amine oxidation provides the first measure of the change in bond order at the carbon involved in this carbon-hydrogen bond cleavage and will be of value to understanding the transition state structure for this class of enzymes.

Enzymatic Kinetic Isotope Effects from Path-Integral Free Energy Perturbation Theory

Path-integral free energy perturbation (PI-FEP) theory is presented to directly determine the ratio of quantum mechanical partition functions of different isotopologs in a single simulation. Furthermore, a double averaging strategy is used to carry out the practical simulation, separating the quantum mechanical path integral exactly into two separate calculations, one corresponding to a classical molecular dynamics simulation of the centroid coordinates, and another involving free-particle path-integral sampling over the classical, centroid positions. An integrated centroid path-integral free energy perturbation and umbrella sampling (PI-FEP/UM, or simply, PI-FEP) method along with bisection sampling was summarized, which provides an accurate and fast convergent method for computing kinetic isotope effects for chemical reactions in solution and in enzymes. The PI-FEP method is illustrated by a number of applications, to highlight the computational precision and accuracy, the rule of geometrical mean in kinetic isotope effects, enhanced nuclear quantum effects in enzyme catalysis, and protein dynamics on temperature dependence of kinetic isotope effects.