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

The critical role of a conserved lysine residue in periplasmic nitrate reductase catalyzed reactions

Periplasmic nitrate reductase NapA from Campylobacter jejuni (C. jejuni) contains a molybdenum cofactor (Moco) and a 4Fe-4S cluster and catalyzes the reduction of nitrate to nitrite. The reducing equivalent required for the catalysis is transferred from NapC → NapB → NapA. The electron transfer from NapB to NapA occurs through the 4Fe-4S cluster in NapA. C. jejuni NapA has a conserved lysine (K79) between the Mo-cofactor and the 4Fe-4S cluster. K79 forms H-bonding interactions with the 4Fe-4S cluster and connects the latter with the Moco via an H-bonding network. Thus, it is conceivable that K79 could play an important role in the intramolecular electron transfer and the catalytic activity of NapA. In the present study, we show that the mutation of K79 to Ala leads to an almost complete loss of activity, suggesting its role in catalytic activity. The inhibition of C. jejuni NapA by cyanide, thiocyanate, and azide has also been investigated. The inhibition studies indicate that cyanide inhibits NapA in a non-competitive manner, while thiocyanate and azide inhibit NapA in an uncompetitive manner. Neither inhibition mechanism involves direct binding of the inhibitor to the Mo-center. These results have been discussed in the context of the loss of catalytic activity of NapA K79A variant and a possible anion binding site in NapA has been proposed.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

Probing Nonadiabaticity of Proton-Coupled Electron Transfer in Ribonucleotide Reductase

The enzyme ribonucleotide reductase, which is essential for DNA synthesis, initiates the conversion of ribonucleotides to deoxyribonucleotides via radical transfer over a 32 Å pathway composed of proton-coupled electron transfer (PCET) reactions. Previously, the first three PCET reactions in the α subunit were investigated with hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations. Herein, the fourth PCET reaction in this subunit between C439 and guanosine diphosphate (GDP) is simulated and found to be slightly exoergic with a relatively high free energy barrier. To further elucidate the mechanisms of all four PCET reactions, we analyzed the vibronic and electron-proton nonadiabaticities. This analysis suggests that interfacial PCET between Y356 and Y731 is vibronically and electronically nonadiabatic, whereas PCET between Y731 and Y730 and between C439 and GDP is fully adiabatic and PCET between Y730 and C439 is in the intermediate regime. These insights provide guidance for selecting suitable rate constant expressions for these PCET reactions.

The carbonyl-lock mechanism underlying non-aromatic fluorescence in biological matter

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics.

Exploring Proton-Coupled Electron Transfer at Multiple Scales

The coupling of electron and proton transfer is critical for chemical and biological processes spanning a wide range of length and time scales and often occurring in complex environments. Thus, diverse modeling strategies, including analytical theories, quantum chemistry, molecular dynamics, and kinetic modeling, are essential for a comprehensive understanding of such proton-coupled electron transfer reactions. Each of these computational methods provides one piece of the puzzle, and all these pieces must be viewed together to produce the full picture.

Dynamically Reversible Interconversion of Molecular Catalysts for Efficient Electrooxidation of Propylene into Propylene Glycol

For the electrooxidation of propylene into 1,2-propylene glycol (PG), the process involves two key steps of the generation of *OH and the transfer of *OH to the C═C bond in propylene. The strong *OH binding energy (E_B(*OH)) favors the dissociation of H₂O into *OH, whereas the transfer of *OH to propylene will be impeded. The scaling relationship of the E_B(*OH) plays a key role in affecting the catalytic performance toward propylene electrooxidation. Herein, we adopt an immobilized Ag pyrazole molecular catalyst (denoted as AgPz) as the electrocatalyst. The pyrrolic N-H in AgPz could undergo deprotonation to form pyrrolic N (denoted as AgPz-H_vac), which can be protonated reversibly. During propylene electrooxidation, the strong E_B(*OH) on AgPz favors the dissociation of H₂O into *OH. Subsequently, the AgPz transforms into AgPz-H_vac that possesses weak E_B(*OH), benefiting to the further combination of *OH and propylene. The dynamically reversible interconversion between AgPz and AgPz-H_vac accompanied by changeable E_B(*OH) breaks the scaling relationship, thus greatly lowering the reaction barrier. At 2.0 V versus Ag/AgCl electrode, AgPz achieves a remarkable yield rate of 288.9 mmol_PG g_cat^-1 h^-1, which is more than one order of magnitude higher than the highest value ever reported.

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.

Theoretical study on hydrogen transfer in the dissociation of dimethyl disulfide radical cations

Hydrogen transfer (HT) is of crucial importance in biochemistry and atmospheric chemistry. Here, HT processes involved in the dissociation reaction of dimethyl disulfide radical cations (DMDS˙⁺, CH₃SSCH₃˙⁺) are investigated using quantum chemical calculations. Four HTs from the C to S atom and one HT from the S to S atom are observed and the most probable paths are proposed in the dissociation channel from DMDS˙⁺ to CH_nS⁺ (n = 2-4). The mechanisms of all these five HTs are described as hydrogen atom transfer (HAT) and four of them are accompanied by electron transfer (ET). Considering the catalytic effect of water molecules existing in organisms and the atmosphere, five HT processes in the dissociation of the [DMDS + H₂O]˙⁺ complex are further explored, which show lower free energy barriers. With the participation of water molecules acting as a base, two HTs from the C to the S atom, which have the largest decrease in energy barriers, are characterized as concerted proton-coupled electron transfer (cPCET). These results can be extended to understand the mechanism of the HT process during the dissociation of disulfide and help provide a strategy to design a rare cPCET mechanism for the activation of the C-H bond.

17 O Hyperfine Spectroscopy Reveals Hydration Structure of Nitroxide Radicals in Aqueous Solutions

The hydration structure of nitroxide radicals in aqueous solutions is elucidated by advanced ¹⁷ O hyperfine (hf) spectroscopy with support of quantum chemical calculations and MD simulations. A piperidine and a pyrrolidine-based nitroxide radical are compared and show clear differences in the preferred directionality of H-bond formation. We demonstrate that these scenarios are best represented in ¹⁷ O hf spectra, where in-plane coordination over σ ${\sigma }$ -type H-bonding leads to little spin density transfer on the water oxygen and small hf couplings, whereas π ${{\rm \pi }}$ -type perpendicular coordination generates much larger hf couplings. Quantitative analysis of the spectra based on MD simulations and DFT predicted hf parameters is consistent with a distribution of close solvating water molecules, in which directionality is influenced by subtle steric effects of the ring and the methyl group substituents.

Electrochemical and Structural Study of the Buried Tryptophan in Azurin: Effects of Hydration and Polarity on the Redox Potential of W48

Tryptophan serves as an important redox-active amino acid in mediating electron transfer and mitigating oxidative damage in proteins. We previously showed a difference in electrochemical potentials for two tryptophan residues in azurin with distinct hydrogen-bonding environments. Here, we test whether reducing the side chain bulk at position Phe110 to Leu, Ser, or Ala impacts the electrochemical potentials (E°) for tryptophan at position 48. X-ray diffraction confirmed the influx of crystallographically resolved water molecules for both the F110A and F110L tyrosine free azurin mutants. The local environments of W48 in all azurin mutants were further evaluated by UV resonance Raman (UVRR) spectroscopy to probe the impact of mutations on hydrogen bonding and polarity. A correlation between the frequency of the ω17 mode─considered a vibrational marker for hydrogen bonding─and E° is proposed. However, the trend is opposite to the expectation from a previous study on small molecules. Density functional theory calculations suggest that the ω17 mode reflects hydrogen bonding as well as local polarity. Further, the UVRR data reveal different intensity/frequency shifts of the ω9/ω10 vibrational modes that characterize the local H-bonding environments of tryptophan. The cumulative data support that the presence of water increases E° and reveal properties of the protein microenvironment surrounding tryptophan.

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.