Abstract
The proton-coupled electron transfer (PCET) reaction catalyzed by soybean lipoxygenase has served as a prototype for understanding hydrogen tunneling in enzymes. Herein this PCET reaction is studied with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The free energy surfaces are computed as functions of the proton donor-acceptor (C-O) distance and the proton coordinate, and the potential of mean force is computed as a function of the C-O distance, inherently including anharmonicity. The simulation results are used to calculate the kinetic isotope effects for the wild-type enzyme (WT) and the L546A/L754A double mutant (DM), which have been measured experimentally to be ∼80 and ∼700, respectively. The PCET reaction is found to be exoergic for WT and slightly endoergic for the DM, and the equilibrium C-O distance for the reactant is found to be ∼0.2 Å greater for the DM than for WT. The larger equilibrium distance for the DM, which is due mainly to less optimal substrate binding in the expanded binding cavity, is primarily responsible for its higher kinetic isotope effect. The calculated potentials of mean force are anharmonic and relatively soft at shorter C-O distances, allowing efficient thermal sampling of the shorter distances required for effective hydrogen tunneling. The primarily local electrostatic field at the transferring hydrogen is ∼100 MV/cm in the direction to facilitate proton transfer and increases dramatically as the C-O distance decreases. These simulations suggest that the overall protein environment is important for conformational sampling of active substrate configurations aligned for proton transfer, but the PCET reaction is influenced primarily by local electrostatic effects that facilitate conformational sampling of shorter proton donor-acceptor distances required for effective hydrogen tunneling.