Extreme tunnelling in methylamine dehydrogenase revealed by hybrid QM/MM calculations: potential energy surface profile for methylamine and ethanolamine substrates and kinetic isotope effect values

Abstract

The rate-determining proton transfer step in the amine reduction reaction catalysed by the enzyme methylamine dehydrogenase has been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) model. Variational transition state theory, combined with multidimensional tunnelling corrections, has been employed to calculate reaction rate constants, and hence deuterium kinetic isotope effects (KIE). To render these calculations computationally feasible, the electronic structure was described using a PM3 method with specific reaction parameters obtained by a fit to energetics obtained at a high level for a small model system. Compared to the use of standard parameters, these revised parameters result in a considerable improvement in the predicted KIE values and activation energy. For both methylamine and ethanolamine substrates, through-barrier, rather than over-barrier, motion is found to dominate with KIE values that are large and close to the experimental values. A major difference between the two substrates is that, for ethanolamine, different hydrogen bonding structures involving the substrate hydroxyl are possible, leading to very different potential energy surfaces with KIE values covering a considerable range. We speculate that this is the origin of the differing temperature behaviour observed for the KIEs of the two substrates.