Ribonucleotide reductase (RNR) is an essential enzyme in DNA synthesis for all living organisms. It reduces ribonucleotides to the corresponding deoxyribonucleotides by a reversible radical transfer mechanism. The active form of E. coli Ia RNR is composed of two subunits, α and β, which form an active asymmetric α 2 β 2 complex. The radical transfer pathway involves a series of proton-coupled electron transfer (PCET) reactions spanning α and β over ∼32 Å. Herein, quantum mechanical/molecular mechanical free energy simulations of PCET between tyrosine residues Y730 and Y731 are performed on the recently solved cryo-EM structure of the active α 2 β 2 complex, which includes a pre-turnover α/β pair with an ordered PCET pathway and a post-turnover α′/β′ pair. The free energy surfaces in both the pre-and post-turnover states are computed. According to the simulations, forward radical transfer from Y731 to Y730 is thermodynamically favored in the pre-turnover state, and backward radical transfer is favored in the post-turnover state, consistent with the reversible mechanism. E623, a glutamate residue that is near these tyrosines only in the pre-turnover state, is discovered to play a key role in facilitating forward radical transfer by thermodynamically stabilizing the radical on Y730 through hydrogen-bonding and electrostatic interactions and lowering the free energy barrier via a proton relay mechanism. Introduction of fluorinated Y731 exhibits expected thermodynamic trends without altering the basic mechanism. These simulations suggest that E623 influences the directionality of PCET between Y731 and Y730 and predict that mutation of E623 will impact catalysis.

Ribonucleotide reductase (RNR) catalyzes the reduction of ribonucleotides to deoxyribonucleotides and is critical for DNA synthesis and repair in all organisms. Its mechanism requires radical transfer along a ∼32 Å pathway through a series of proton-coupled electron transfer (PCET) steps. Previous simulations suggested that a glutamate residue (E623) mediates the PCET reaction between two stacked tyrosine residues (Y730 and Y731) through a proton relay mechanism. This work focuses on the adjacent PCET reaction between Y730 and a cysteine residue (C439). Quantum mechanical/molecular mechanical free energy simulations illustrate that when Y730 and Y731 are stacked, E623 stabilizes the radical on C439 through hydrogen bonding with the Y730 hydroxyl group. When Y731 is flipped away from Y730, a water molecule stabilizes the radical on C439 through hydrogen bonding with Y730 and lowers the free energy barrier for radical transfer from Y730 to C439 through electrostatic interactions with the transferring hydrogen but does not directly accept the proton. These simulations indicate that the conformational motions and electrostatic interactions of the tyrosines, cysteine, glutamate, and water strongly impact the thermodynamics and kinetics of these two coupled PCET reactions. Such insights are important for protein engineering efforts aimed at altering radical transfer in RNR.

Ribonucleotide reductase (RNR) regulates DNA synthesis and repair in all organisms. The mechanism of Escherichia coli RNR requires radical transfer over a proton-coupled electron transfer (PCET) pathway spanning ∼32 Å across two protein subunits. A key step along this pathway is the interfacial PCET reaction between Y356 in the β subunit and Y731 in the α subunit. Herein, this PCET reaction between two tyrosines across an aqueous interface is explored with classical molecular dynamics and quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The simulations suggest that the water-mediated mechanism involving double proton transfer through an intervening water molecule is thermodynamically and kinetically unfavorable. The direct PCET mechanism between Y356 and Y731 becomes feasible when Y731 is flipped toward the interface and is predicted to be approximately isoergic with a relatively low free energy barrier. This direct mechanism is facilitated by the hydrogen bonding of water to both Y356 and Y731. These simulations provide fundamental insights into radical transfer across aqueous interfaces.