Switching the proton-coupled electron transfer mechanism for non-canonical tyrosine residues in a de novo protein.

The proton-coupled electron transfer (PCET) reactions of tyrosine (Y) are instrumental to many redox reactions in nature. This study investigates how the local environment and the thermodynamic properties of Y influence its PCET characteristics. Herein, 2- and 4-mercaptophenol (MP) are placed in the well-folded α3C protein (forming 2MP-α3C and 4MP-α3C) and oxidized by external light-generated [Ru(L)3]3+ complexes. The resulting neutral radicals are long-lived (>100 s) with distinct optical and EPR spectra. Calculated spin-density distributions are similar to canonical Y˙ and display very little spin on the S-S bridge that ligates the MPs to C32 inside the protein. With 2MP-α3C and 4MP-α3C we probe how proton transfer (PT) affects the PCET rate constants and mechanisms by varying the degree of solvent exposure or the potential to form an internal hydrogen bond. Solution NMR ensemble structures confirmed our intended design by displaying a major difference in the phenol OH solvent accessible surface area (≤∼2% for 2MP and 30-40% for 4MP). Additionally, 2MP-C32 is within hydrogen bonding distance to a nearby glutamate (average O-O distance is 3.2 ± 0.5 Å), which is suggested also by quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations. Neither increased exposure of the phenol OH to solvent (buffered water), nor the internal hydrogen bond, was found to significantly affect the PCET rates. However, the lower phenol pKa values associated with the MP-α3C proteins compared to α3Y provided a sufficient change in PT driving force to alter the PCET mechanism. The PCET mechanism for 2MP-α3C and 4MP-α3C with moderately strong oxidants was predominantly step-wise PTET for pH values, but changed to concerted PCET at neutral pH values and below when a stronger oxidant was used, as found previously for α3Y. This shows how the balance of ET and PT driving forces is critical for controlling PCET mechanisms. The presented results improve our general understanding of amino-acid based PCET in enzymes.

Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify.

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.

Direct Proton-Coupled Electron Transfer between Interfacial Tyrosines in Ribonucleotide Reductase.

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.

"Student-led workshop strengthens perceived discussion skills and community in an interdisciplinary graduate program".

The Integrated Graduate Program in Physical and Engineering Biology (IGPPEB) at Yale University brings together Ph.D. students from the physical, engineering, and biological sciences. The main goals of this program are for students to become comfortable working in an interdisciplinary and collaborative research environment and adept at communicating with scientists and nonscientists. To fill a student-identified learning gap in engaging in inclusive discussions, IGPPEB students developed a communication workshop to improve skills in visual engagement, citing specific content, constructive conversation entrances, and encouragement of peers. Based on short- and long-term assessment of the workshop, 100% of students reported that it should be offered to future cohorts and 63% of students perceived it to be personally helpful. Additionally, 92% of participants reported using one or more of the core skills beyond the course, with skills in "Encouraging peers" and "Constructive conversation entrances" rated the highest in perceived improvement. Based on the highest average rating of 76 ± 24 (on a scale of 0-100), students agreed that the workshop made them feel more welcome in the IGPPEB community. With a rating of 68 ± 13, they also agreed that the workshop had a positive impact on their graduate school experience. Participants provided suggestions for future improvements, such as increasing student involvement in leading discussions of course material. This study demonstrates that a student-led workshop can improve perceived discussion skills and build community across an interdisciplinary program in the sciences.

Kinetic model for reversible radical transfer in ribonucleotide reductase.

The enzyme ribonucleotide reductase (RNR), which catalyzes the reduction of ribonucleotides to deoxynucleotides, is vital for DNA synthesis, replication, and repair in all living organisms. Its mechanism requires long-range radical translocation over ∼32 Šthrough two protein subunits and the intervening aqueous interface. Herein, a kinetic model is designed to describe reversible radical transfer in Escherichia coli RNR. This model is based on experimentally studied photoRNR systems that allow the photochemical injection of a radical at a specific tyrosine residue, Y356, using a photosensitizer. The radical then transfers across the interface to another tyrosine residue, Y731, and continues until it reaches a cysteine residue, C439, which is primed for catalysis. This kinetic model includes radical injection, an off-pathway sink, radical transfer between pairs of residues along the pathway, and the conformational flipping motion of Y731 at the interface. Most of the input rate constants for this kinetic model are obtained from previous experimental measurements and quantum mechanical/molecular mechanical free-energy simulations. Ranges for the rate constants corresponding to radical transfer across the interface are determined by fitting to the experimentally measured Y356 radical decay times in photoRNR systems. This kinetic model illuminates the time evolution of radical transport along the tyrosine and cysteine residues following radical injection. Further analysis identifies the individual rate constants that may be tuned to alter the timescale and probability of the injected radical reaching C439. The insights gained from this kinetic model are relevant to biochemical understanding and protein-engineering efforts with potential pharmacological implications.

Role of Water in Proton-Coupled Electron Transfer between Tyrosine and Cysteine in Ribonucleotide Reductase.

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.

Glutamate Mediates Proton-Coupled Electron Transfer Between Tyrosines 730 and 731 in Escherichia coli Ribonucleotide Reductase.

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.

Computing Proton-Coupled Redox Potentials of Fluorotyrosines in a Protein Environment.

The oxidation of tyrosine to form the neutral tyrosine radical via proton-coupled electron transfer is essential for a wide range of biological processes. The precise measurement of the proton-coupled redox potentials of tyrosine (Y) in complex protein environments is challenging mainly because of the highly oxidizing and reactive nature of the radical state. Herein, a computational strategy is presented for predicting proton-coupled redox potentials in a protein environment. In this strategy, both the reduced Y-OH and oxidized Y-O• forms of tyrosine are sampled with molecular dynamics using a molecular mechanical force field. For a large number of conformations, a quantum mechanical/molecular mechanical (QM/MM) electrostatic embedding scheme is used to compute the free-energy differences between the reduced and oxidized forms, including the zero-point energy and entropic contributions as well as the impact of the protein electrostatic environment. This strategy is applied to a series of fluorinated tyrosine derivatives embedded in a de novo α-helical protein denoted as α3Y. The force fields for both the reduced and oxidized forms of these noncanonical fluorinated tyrosine residues are parameterized for general use. The calculated relative proton-coupled redox potentials agree with experimentally measured values with a mean unsigned error of 24 mV. Analysis of the simulations illustrates that hydrogen-bonding interactions between tyrosine and water increase the redox potentials by ∼100-250 mV, with significant variations because of the fluctuating protein environment. This QM/MM approach enables the calculation of proton-coupled redox potentials of tyrosine and other residues such as tryptophan in a variety of protein systems.

Conformational Motions and Water Networks at the α/ß Interface in E. coli Ribonucleotide Reductase.

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 from Tyrosine in the Interior of a de novo Protein: Mechanisms and Primary Proton Acceptor.

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 (kPCET) 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. kPCET 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.

Editing Domain Motions Preorganize the Synthetic Active Site of Prolyl-tRNA Synthetase.

Prolyl-tRNA synthetases (ProRSs) catalyze the covalent attachment of proline onto cognate tRNAs, an indispensable step for protein synthesis in all living organisms. ProRSs are modular enzymes and the "prokaryotic-like" ProRSs are distinguished from "eukaryotic-like" ProRSs by the presence of an editing domain (INS) inserted between motifs 2 and 3 of the main catalytic domain. Earlier studies suggested the presence of coupled-domain dynamics could contribute to catalysis; however, the role that the distal, highly mobile INS domain plays in catalysis at the synthetic active site is not completely understood. In the present study, a combination of theoretical and experimental approaches has been used to elucidate the precise role of INS domain dynamics. Quantum mechanical/molecular mechanical simulations were carried out to model catalytic Pro-AMP formation by Enterococcus faecalis ProRS. The energetics of the adenylate formation by the wild-type enzyme was computed and contrasted with variants containing active site mutations, as well as a deletion mutant lacking the INS domain. The combined results revealed that two distinct types of dynamics contribute to the enzyme's catalytic power. One set of motions is intrinsic to the INS domain and leads to conformational preorganization that is essential for catalysis. A second type of motion, stemming from the electrostatic reorganization of active site residues, impacts the height and width of the energy profile and has a critical role in fine tuning the substrate orientation to facilitate reactive collisions. Thus, motions in a distal domain can preorganize the active site of an enzyme to optimize catalysis.

Propensity for Proton Relay and Electrostatic Impact of Protein Reorganization in Slr1694 BLUF Photoreceptor.

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.

Cyclic Changes in Active Site Polarization and Dynamics Drive the 'Ping-pong' Kinetics in NRH:Quinone Oxidoreductase 2: An Insight from QM/MM Simulations.

Quinone reductases belong to the family of flavin-dependent oxidoreductases. With the redox active cofactor, flavin adenine dinucleotide, quinone reductases are known to utilize a 'ping-pong' kinetic mechanism during catalysis in which a hydride is bounced back and forth between flavin and its two substrates. However, the continuation of this catalytic cycle requires product displacement steps, where the product of one redox half-cycle is displaced by the substrate of the next half-cycle. Using improved hybrid quantum mechanical/molecular mechanical simulations, both the catalytic hydride transfer and the product displacement reactions were studied in NRH:quinone oxidoreductase 2. Initially, the self-consistent charge-density functional tight binding theory was used to describe flavin ring and the substrate atoms, while embedded in the molecular mechanically-treated solvated active site. Then, for each step of the catalytic cycle, a further improvement of energetics was made using density functional theory-based corrections. The present study showcases an integrated interplay of solvation, protonation, and protein matrix-induced polarization as the driving force behind the thermodynamic wheel of the 'ping-pong' kinetics. Reported here is the first-principles model of the 'ping-pong' kinetics that portrays how cyclic changes in the active site polarization and dynamics govern the oscillatory hydride transfer and product displacement in this enzyme.

Insight into the kinetics and thermodynamics of the hydride transfer reactions between quinones and lumiflavin: a density functional theory study.

The kinetics and equilibrium of the hydride transfer reaction between lumiflavin and a number of substituted quinones was studied using density functional theory. The impact of electron withdrawing/donating substituents on the redox potentials of quinones was studied. In addition, the role of these substituents on the kinetics of the hydride transfer reaction with lumiflavin was investigated in detail under the transition state (TS) theory assumption. The hydride transfer reactions were found to be more favorable for an electron-withdrawing substituent. The activation barrier exhibited a quadratic relationship with the driving force of these reactions as derived under the formalism of modified Marcus theory. The present study found a significant extent of electron delocalization in the TS that is stabilized by enhanced electrostatic, polarization, and exchange interactions. Analysis of geometry, bond-orders, and energetics revealed a predominant parallel (Leffler-Hammond) effect on the TS. Closer scrutiny reveals that electron-withdrawing substituents, although located on the acceptor ring, reduce the N-H bond order of the donor fragment in the precursor complex. Carried out in the gas-phase, this is the first ever report of a theoretical study of flavin's hydride transfer reactions with quinones, providing an unfiltered view of the electronic effect on the nuclear reorganization of donor-acceptor complexes.

Effect of stacking interactions on the thermodynamics and kinetics of lumiflavin: a study with improved density functionals and density functional tight-binding protocol.

The π-π stacking interaction between lumiflavin and a number of π-electron-rich molecules has been studied by density functional theory using several new-generation density functionals. Six known lumiflavin-aromatic adducts were used and the models were evaluated by comparing the geometry and energetics with experimental results. The study found that dispersion-corrected and hybrid functionals with larger (>50%) Hartree-Fock exchanges produced superior results in modeling thermodynamic characteristics of these complexes. The functional producing the best energetics for these model systems was used to study the stacking interactions of lumiflavin with biologically relevant aromatic groups. Additionally, the reduction of flavin-in the presence of both a hydride donor and a nondonor π-electronic system was also studied. Weak interactions were observed in the stacked lumiflavin complexes of benzene, phenol, and indole, mimicking phenyl alanine, tryptophan, and tyrosine side chains, respectively, of an enzyme. The stacked complex of naphthalene and flavin showed little change in flavin's redox potential indicating insignificant effect on the thermodynamics of the hydride transfer reaction. In contrast, the hydride transfer reaction with the hydride donor N-methyl nicotinamide tells a different story, as the transition state was found to be strongly impacted by the stacking interactions. A comparison of performance between the density functional theory (DFT) and the computationally less expensive dispersion-corrected self-consistent density functional tight-binding (SCC-DFTB-D) theory revealed that the latter produces consistent energetics for this hydride transfer reaction and additional DFT-computed perturbative corrections could significantly improve these results.