Potential of mean force calculations have been performed on the wild-type medium chain acyl-CoA dehydrogenase (MCAD) and two of its mutant forms. Initial simulation and analysis of the active site of the enzyme reveals that an arginine residue (Arg256), conserved in the substrate binding domain of this group of enzymes, exists in two alternate conformations, only one of which makes the enzyme active. This active conformation was used in subsequent computations of the enzymatic reactions. It is known that the catalytic α,β-dehydrogenation of fatty acyl-CoAs consists of two C-H bond dissociation processes: a proton abstraction and a hydride transfer. Energy profiles of the two reaction steps in the wild-type MCAD demonstrate that the reaction proceeds by a stepwise mechanism with a transient species. The activation barriers of the two steps differ by only ∼2 kcal/ mol, indicating that both may contribute to the rate-limiting process. Thus this may be a stepwise dissociation mechanism whose relative barriers can be tuned by suitable alterations of the substrate and/or enzyme. Analysis of the structures along the reaction path reveals that Arg256 plays a key role in maintaining the reaction-center hydrogen-bonding network involving the thioester carbonyl group, which stabilizes transition states as well as the intervening transient species. Mutation of this arginine residue to glutamine increases the activation barrier of the hydride transfer reaction by ∼5 kcal/mol, and the present simulations predict a substantial loss of catalytic activity for this mutant. Structural analysis of this mutant reveals that the orientation of the thioester moiety of the substrate has been changed significantly as compared to that in the wild-type enzyme. In contrast, simulation of the active site of the Thr168Ala mutant shows no significant change in the relative orientation of the substrate and the cofactor in the active site; as a result, this mutation has very little effect on the overall reaction barrier, and this is consistent with the experimental data. This study demonstrates that significant insights of the catalytic mechanism can be obtained by these simulated enzyme catalysis studies whose results can pave the way of designing novel mechanistic probes for the enzyme.Mitochondrial β-oxidation has been extensively studied in the past 15 years because a number of mutations in several of the enzymes involved are responsible for diseases related to fatty acid metabolism. Fatty acid metabolism occurs within the mitochondrial matrix and proceeds by a sequence of steps that remove two carbon units in each cycle. The process involves at a This work was partially supported by grant GM 29344 from the National Institutes of Healths to MTS, by grant GM046736 from National Institutes of Healths to JG, by NSF grant CHE03−49122 to DGT.
NIH-PA Author ManuscriptNIH-PA Author Manuscript NIH-PA Author Manuscript least 12 different enzymes, including the acyl-CoA dehydrogenases, which are flavoenzymes that perform the first step of thi...