Nuclear Quantum Effects in Quantum Mechanical/Molecular Mechanical Free Energy Simulations of Ribonucleotide Reductase

The enzyme ribonucleotide reductase plays a critical role in DNA synthesis and repair. Its mechanism requires long-range radical transfer through a series of proton-coupled electron transfer (PCET) steps. Nuclear quantum effects such as zero-point energy, proton delocalization, and hydrogen tunneling are known to be important in PCET. We present a strategy for efficiently incorporating nuclear quantum effects into multidimensional free energy surfaces and real-time dynamical simulations for condensed-phase systems such as enzymes. This strategy is based on the nuclear-electronic orbital (NEO) method, which treats specified protons quantum mechanically on the same level as the electrons. NEO density functional theory (NEO-DFT) is combined with the quantum mechanical/molecular mechanical finite temperature string method with umbrella sampling via a simple reweighting procedure. Application of this strategy to PCET between two tyrosines, Y731 and Y730, in ribonucleotide reductase illustrates that nuclear quantum effects could either raise or lower the free energy barrier, leading to a range of possible kinetic isotope effects. Real-time time-dependent DFT (RT-NEO-TDDFT) simulations highlight nuclear-electronic quantum dynamics. These approaches enable the incorporation of nuclear quantum effects into a wide range of chemically and biologically important processes.