Ribonucleotide reductases (RNRs) catalyze the conversion of all four ribonucleotides to deoxyribonucleotides and are essential for DNA synthesis in all organisms. The active form of E. coli Ia RNR is composed of two homodimers that form the active α2β2 complex. Catalysis is initiated by long-range radical translocation over a ∼32 Å proton-coupled electron transfer (PCET) pathway involving Y356β and Y731α at the interface. Resolving the PCET pathway at the α/β interface has been a long-standing challenge due to the lack of structural data. Herein, molecular dynamics simulations based on a recently solved cryogenic-electron microscopy structure of an active α2β2 complex are performed to examine the structure and fluctuations of interfacial water, as well as the hydrogen-bonding interactions and conformational motions of interfacial residues along the PCET pathway. Our free energy simulations reveal that Y731 is able to sample both a flipped-out conformation, where it points toward the interface to facilitate interfacial PCET with Y356, and a stacked conformation with Y730 to enable collinear PCET with this residue. Y356 and Y731 exhibit hydrogen-bonding interactions with interfacial water molecules and, in some conformations, share a bridging water molecule, suggesting that the primary proton acceptor for PCET from Y356 and from Y731 is interfacial water. The conformational flexibility of Y731 and the hydrogen-bonding interactions of both Y731 and Y356 with interfacial water and hydrogen-bonded water chains appear critical for effective radical translocation along the PCET pathway. These simulations are consistent with biochemical and spectroscopic data and provide previously unattainable atomic-level insights into the fundamental mechanism of RNR.

Proton-coupled electron transfer (PCET) from tyrosine produces a neutral tyrosyl radical (Y•) that is vital to many catalytic redox reactions. To better understand how the protein environment influences the PCET properties of tyrosine, we have studied the radical formation behavior of Y32 in the α3Y model protein. The previously solved α3Y solution NMR structure shows that Y32 is sequestered ∼7.7 ± 0.3 Å below the protein surface without any primary proton acceptors nearby. Here we present transient absorption kinetic data and molecular dynamics (MD) simulations to resolve the PCET mechanism associated with Y32 oxidation. Y32 • was generated in a bimolecular reaction with [Ru­(bpy)3]3+ formed by flash photolysis. At pH > 8, the rate constant of Y32 • formation (k PCET) increases by one order of magnitude per pH unit, corresponding to a proton-first mechanism via tyrosinate (PTET). At lower pH < 7.5, the pH dependence is weak and shows a previously measured KIE ≈ 2.5, which best fits a concerted mechanism. k PCET is independent of phosphate buffer concentration at pH 6.5. This provides clear evidence that phosphate buffer is not the primary proton acceptor. MD simulations show that one to two water molecules can enter the hydrophobic cavity of α3Y and hydrogen bond to Y32, as well as the possibility of hydrogen-bonding interactions between Y32 and E13, through structural fluctuations that reorient surrounding side chains. Our results illustrate how protein conformational motions can influence the redox reactivity of a tyrosine residue and how PCET mechanisms can be tuned by changing the pH even when the PCET occurs within the interior of a protein.

Photoreceptor proteins play a vital role in a wide range of light-regulated processes. The formation of the light-adapted state of blue light using flavin (BLUF) photoreceptors is thought to involve rearrangements of hydrogen-bonding networks upon photoexcitation. Free energy simulations with partial charges corresponding to relevant ground and excited states of the Slr1694 BLUF domain characterize conformations prior to and following photoexcitation. The simulations indicate that Trp91 is thermodynamically favored to be in the active site, although it is also able to sample conformations outside the active site. For experimentally observed conformations of Trp91, Gln50 is thermodynamically favored to be oriented for a proton relay bridging Tyr8 and the flavin. When Trp91 is rotated such that it can donate a hydrogen bond to Gln50, as observed in other BLUF domains, the proton relay is not thermodynamically favored in the ground state, providing a possible explanation for the relatively fast photocycle of the Slr1694 BLUF domain. Photoexcitation to the locally excited (LE) state of the flavin induces the formation of the proton relay if it is not already formed. Electrostatically embedded time-dependent density functional theory calculations indicate that the proton relay reduces the energy gap between the LE state and the charge-transfer (CT) state associated with electron transfer from Tyr8 to the flavin. Although the CT state is higher in energy than the LE state prior to photoexcitation, the protein environment can reorganize in a manner that stabilizes the CT state so that it is lower than the LE state, enabling the LE to CT state transition. An electrostatic analysis identifies motions of individual residues, such as Arg65, that stabilize electron transfer from Tyr8 to the flavin. These conformational changes facilitate the critical proton-coupled electron transfer reaction in the BLUF photocycle.

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.