Search

A bioinspired naphthoquinone model of the quinones in photosynthetic reaction centers but bearing an intramolecular hydrogen-bonded carboxylic acid has been synthesized and characterized electrochemically, spectroscopically, and computationally to provide mechanistic insight into the role of proton-coupled electron transfer (PCET) of quinone reduction in photosynthesis. The reduction potential of this construct is 370 mV more positive than the unsubstituted naphthoquinone. In addition to the reversible cyclic voltammetry, infrared spectroelectrochemistry confirms that the naphthoquinone/naphthoquinone radical anion couple is fully reversible. Calculated redox potentials agree with the experimental trends arising from the intramolecular hydrogen bond. Molecular electrostatic potentials illustrate the reversible proton transfer driving forces, and analysis of the computed vibrational spectra supports the possibility of a combination of electron transfer and PCET processes. The significance of PCET, reversibility, and redox potential management relevant to the design of artificial photosynthetic assemblies involving PCET processes is discussed., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

A Co II -porphyrin complex (1) with an appended aza-crown ether for Lewis acid (LA) binding was synthesized and characterized. NMR spectroscopy and electrochemistry show that cationic group I and II LAs ( i.e. , Li + , Na + , K + , Ca 2+ , Sr 2+ , and Ba 2+ ) bind to the aza-crown ether group of 1. The binding constant for Li + is comparable to that observed for a free aza-crown ether. LA binding causes an anodic shift in the Co II /Co I couple of between 10 and 40 mV and also impacts the Co III /Co II couple. The magnitude of the anodic shift of the Co II /Co I couple varies linearly with the strength of the LA as determined by the p K a of the corresponding metal-aqua complex, with dications giving larger shifts than monocations. The extent of the anodic shift of the Co II /Co I couple also increases as the ionic strength of the solution decreases. This is consistent with electric field effects being responsible for the changes in the redox properties of 1 upon LA binding and provides a novel method to tune the reduction potential. Density functional theory calculations indicate that the bound LA is 5.6 to 6.8 Å away from the Co II ion, demonstrating that long-range electrostatic effects, which do not involve changes to the primary coordination sphere, are responsible for the variations in redox chemistry. Compound 1 was investigated as a CO 2 reduction electrocatalyst and shows high activity but rapid decomposition., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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 α 3 C protein (forming 2MP-α 3 C and 4MP-α 3 C) 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 C 32 inside the protein. With 2MP-α 3 C and 4MP-α 3 C 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-C 32 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 p K a values associated with the MP-α 3 C proteins compared to α 3 Y provided a sufficient change in PT driving force to alter the PCET mechanism. The PCET mechanism for 2MP-α 3 C and 4MP-α 3 C 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 α 3 Y. 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., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear-electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. These calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.

Surface immobilized catalysts are highly promising candidates for a range of energy conversion reactions, and atomistic mechanistic understanding is essential for their rational design. Cobalt tetraphenylporphyrin (CoTPP) nonspecifically adsorbed on a graphitic surface has been shown to undergo concerted proton-coupled electron transfer (PCET) in aqueous solution. Herein, density functional theory calculations on both cluster and periodic models representing π-stacked interactions or axial ligation to a surface oxygenate are performed. As the electrode surface is charged due to applied potential, the adsorbed molecule experiences the electrical polarization of the interface and nearly the same electrostatic potential as the electrode, regardless of the adsorption mode. PCET occurs by electron abstraction from the surface to the CoTPP concerted with protonation to form a cobalt hydride, thereby circumventing Co(II/I) redox. Specifically, the Co(II) d-state localized orbital interacts with a proton from solution and an electron from the delocalized graphitic band states to produce a Co(III)-H bonding orbital below the Fermi level, corresponding to redistribution of electrons from the band states to the bonding states. These insights have broad implications for electrocatalysis by chemically modified electrodes and surface immobilized catalysts., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

We develop an electrostatic map for the vibrational NH stretch (amide A) of the protein backbone with a focus on vibrational chiral sum frequency generation spectroscopy (chiral SFG). Chiral SFG has been used to characterize protein secondary structure at interfaces using the NH stretch and to investigate chiral water superstructures around proteins using the OH stretch. Interpretation of spectra has been complicated because the NH stretch and OH stretch overlap spectrally. Although an electrostatic map for water OH developed by Skinner and co-workers was used previously to calculate the chiral SFG response of water structures around proteins, a map for protein NH that is directly responsive to biological complexity has yet to be developed. Here, we develop such a map, linking the local electric field to vibrational frequencies and transition dipoles. We apply the map to two protein systems and achieve much better agreement with experiment than was possible in our previous studies. We show that couplings between NH and OH vibrations are crucial to the line shape, which informs the interpretation of chiral SFG spectra, and that the chiral NH stretch response is sensitive to small differences in structure. This work increases the utility of the NH stretch in biomolecular spectroscopy.

Cryogenic ion vibrational spectra of D 2 -tagged cyanobenzoate (CBA) derivatives are obtained and analyzed to characterize the intrinsic spectroscopic responses of the -CO 2 - headgroup to its location on the ring in both the isolated anions and the cationic complexes with divalent metal ions, M 2+ (M = Mg, Ca, Sr). The benzonitrile functionality establishes the different ring isomers ( para , meta , ortho ) according to the location of the carboxylate and provides an additional reporter on the molecular response to the proximal charge center. The aromatic carboxylates display shifts slightly smaller than those observed for a related aliphatic system upon metal ion complexation. Although the CBA anions display very similar band patterns for all three ring positions, upon complexation with metal ions, the ortho isomer yields dramatically different spectral responses in both the -CO 2 - moiety and the CN group. This behavior is traced to the emergence of a tridentate binding motif unique to the ortho isomer in which the metal ions bind to both the oxygen atoms of the carboxylate group and the N atom of the cyano group. In that configuration, the -CO 2 - moiety is oriented perpendicular to the phenyl ring, and the CN stretching fundamental is both strong and red-shifted relative to its behavior in the isolated neutral. The behaviors of the metal-bound ortho complexes occur in contrast to the usual blue shifts associated with "Lewis" type binding of metal ions end-on to -CN. The origins of these spectroscopic features are analyzed with the aid of electronic structure calculations, which also explore differences expected for complexation of monovalent cations to the ortho carboxylate. The resulting insights have implications for understanding the balance between electrostatic and steric interactions at metal binding sites in chemical and biological systems.

Understanding the role of water in biological processes remains a central challenge in the life sciences. Water structures in hydration shells of biomolecules are difficult to study in situ due to overwhelming background from aqueous environments. Biological interfaces introduce additional complexity because biomolecular hydration differs at interfaces compared to bulk solution. Here, we perform experimental and computational studies of chiral sum frequency generation (chiral SFG) spectroscopy to probe chirality transfer from a protein to the surrounding water molecules. This work reveals that chiral SFG probes the first hydration shell around the protein almost exclusively. We explain the selectivity to the first hydration shell in terms of the asymmetry induced by the protein structure and specific protein-water hydrogen-bonding interactions. This work establishes chiral SFG as a powerful technique for studying hydration shell structures around biomolecules at interfaces, presenting new possibilities to address grand research challenges in biology, including the molecular origins of life., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

Accurate modeling of important nuclear quantum effects, such as nuclear delocalization, zero-point energy, and tunneling, as well as non-Born-Oppenheimer effects, requires treatment of both nuclei and electrons quantum mechanically. The nuclear-electronic orbital (NEO) method provides an elegant framework to treat specified nuclei, typically protons, on the same level as the electrons. In conventional electronic structure theory, finding a converged ground state can be a computationally demanding task; converging NEO wavefunctions, due to their coupled electronic and nuclear nature, is even more demanding. Herein, we present an efficient simultaneous optimization method that uses the direct inversion in the iterative subspace method to simultaneously converge wavefunctions for both the electrons and quantum nuclei. Benchmark studies show that the simultaneous optimization method can significantly reduce the computational cost compared to the conventional stepwise method for optimizing NEO wavefunctions for multicomponent systems.

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.

Our recent experimental and theoretical investigations have shown that fluorene C-H bonds can be activated through a mechanism in which the proton and electron are transferred from the C-H bond to a separate base and oxidant in a concerted, elementary step. This multisite proton-coupled electron transfer (MS-PCET) mechanism for C-H bond activation was shown to be facilitated by shorter proton donor-acceptor distances. With the goal of intentionally modulating this donor-acceptor distance, we have now studied C-H MS-PCET in the 3-methyl-substituted fluorenyl benzoate (2-Flr-3-Me-BzO - ). This derivative was readily oxidized by ferrocenium oxidants by initial C-H MS-PCET, with rate constants that were 6-21 times larger than those for 2-Flr-BzO - with the same oxidants. Structural comparisons by X-ray crystallography and by computations showed that addition of the 3-methyl group caused the expected steric compression; however, the relevant C···O - proton donor-acceptor distance was longer, due to a twist of the carboxylate group. The structural changes induced by the 3-Me group increased the basicity of the carboxylate, weakened the C-H bond, and reduced the reorganization energy for C-H bond cleavage. Thus, the rate enhancement for 2-Flr-3-Me-BzO - was due to effects on the thermochemistry and kinetic barrier, rather than from compression of the C···O - proton donor-acceptor distance. These results highlight both the challenges of controlling molecules on the 0.1 Å length scale and the variety of parameters that affect PCET rate constants.

Chiral vibrational sum frequency generation (SFG) spectroscopy probes the structure of the solvation shell around chiral macromolecules. The dominant theoretical framework for understanding the origin of chiral SFG signals is based on the analysis of molecular symmetry, which assumes no interaction between molecules. However, water contains strong intermolecular interactions that significantly affect its properties. Here, the role of intermolecular vibrational coupling in the chiral SFG response of the O-H stretch of water surrounding an antiparallel β-sheet at the vacuum-water interface is investigated. Both intramolecular and intermolecular couplings between O-H groups are required to simulate the full lineshape of the chiral SFG signal. This dependence is also observed for a chiral water dimer, illustrating that this phenomenon is not specific to larger systems. We also find that a dimer of C 3v molecules predicted to be chirally SFG-inactive by the symmetry-based theory can generate a chiral SFG signal when intermolecular couplings are considered, suggesting that even highly symmetric solvent molecules may produce chiral SFG signals when interacting with a chiral solute. The consideration of intermolecular couplings extends the prevailing theory of the chiral SFG response to structures larger than individual molecules and provides guidelines for future modeling.

Small ribozymes cleave their RNA phosphodiester backbone by catalyzing a transphosphorylation reaction wherein a specific O2' functions as the nucleophile. While deprotonation of this alcohol through its acidification would increase its nucleophilicity, little is known about the p K a 's are not readily accessible experimentally. Herein, we turn to molecular dynamics to calculate the p K a 's are not readily accessible experimentally. Herein, we turn to molecular dynamics to calculate the p K a of the nucleophilic O2' in the hairpin ribozyme and to study interactions within the active site that may impact its value. We estimate the p K a of the nucleophilic O2' in the wild-type hairpin ribozyme to be 18.5 ± 0.8, which is higher than the reference compound, and identify a correlation between proper positioning of the O2' for nucleophilic attack and elevation of its p K a . We find that monovalent ions may play a role in depression of the O2' p K a , while the exocyclic amine appears to be important for organizing the ribozyme active site. Overall, this study suggests that the p K a of the O2' is raised in the ground state and lowers during the course of the reaction owing to positioning and metal ion interactions.

Proton-coupled electron transfer (PCET) reactions depend on the hydrogen-bond connectivity between sites of proton donors and acceptors. The 2-(2'-hydroxyphenyl) benzimidazole (BIP) based systems, which mimic the natural Tyr Z -His190 pair of Photosystem II, have been useful for understanding the associated PCET process triggered by one-electron oxidation of the phenol. Substitution of the benzimidazole by an appropriate terminal proton acceptor (TPA) group allows for two-proton translocations. However, the prototropic properties of substituted benzimidazole rings and rotation around the bond linking the phenol and the benzimidazole can lead to isomers that interrupt the intramolecular hydrogen-bonded network and thereby prevent a second proton translocation. Herein, a strategic symmetrization of a benzimidazole based system with two identical TPAs yields an uninterrupted network of intramolecular hydrogen bonds regardless of the isomeric form. NMR data confirms the presence of a single isomeric form in the disubstituted system but not in the monosubstituted system in certain solvents. Infrared spectroelectrochemistry demonstrates a two-proton transfer process associated with the oxidation of the phenol occurring at a lower redox potential in the disubstituted system relative to its monosubstituted analogue. Computational studies support these findings and show that the disubstituted system stabilizes the oxidized two-proton transfer product through the formation of a bifurcated hydrogen bond. Considering the prototropic properties of the benzimidazole heterocycle in the context of multiple PCET will improve the next generation of novel, bioinspired constructs built by concatenated units of benzimidazoles, thus allowing proton translocations at nanoscale length., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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 α 3 Y. 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.

Biomolecular hydration is fundamental to biological functions. Using phase-resolved chiral sum-frequency generation spectroscopy (SFG), we probe molecular architectures and interactions of water molecules around a self-assembling antiparallel β-sheet protein. We find that the phase of the chiroptical response from the O-H stretching vibrational modes of water flips with the absolute chirality of the (l-) or (d-) antiparallel β-sheet. Therefore, we can conclude that the (d-) antiparallel β-sheet organizes water solvent into a chiral supermolecular structure with opposite handedness relative to that of the (l-) antiparallel β-sheet. We use molecular dynamics to characterize the chiral water superstructure at atomic resolution. The results show that the macroscopic chirality of antiparallel β-sheets breaks the symmetry of assemblies of surrounding water molecules. We also calculate the chiral SFG response of water surrounding (l-) and (d-) LK 7 β to confirm the presence of chiral water structures. Our results offer a different perspective as well as introduce experimental and computational methodologies for elucidating hydration of biomacromolecules. The findings imply potentially important but largely unexplored roles of water solvent in chiral selectivity of biomolecular interactions and the molecular origins of homochirality in the biological world., Competing Interests: The authors declare no competing interest.

Blue light using flavin (BLUF) photoreceptor proteins are critical for many light-activated biological processes and are promising candidates for optogenetics because of their modular nature and long-range signaling capabilities. Although the photocycle of the Slr1694 BLUF domain has been characterized experimentally, the identity of the light-adapted state following photoexcitation of the bound flavin remains elusive. Herein hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations of this photocycle provide a nonequilibrium dynamical picture of a possible mechanism for the formation of the light-adapted state. Photoexcitation of the flavin induces a forward proton-coupled electron transfer (PCET) process that leads to the formation of an imidic acid tautomer of Gln50. The calculations herein show that the subsequent rotation of Gln50 allows a reverse PCET process that retains this tautomeric form. In the resulting purported light-adapted state, the glutamine tautomer forms a hydrogen bond with the flavin carbonyl group. Additional ensemble-averaged QM/MM calculations of the dark-adapted and purported light-adapted states demonstrate that the light-adapted state with the imidic acid glutamine tautomer reproduces the experimentally observed spectroscopic signatures. Specifically, the calculations reproduce the red shifts in the flavin electronic absorption and carbonyl stretch infrared spectra in the light-adapted state. Further hydrogen-bonding analyses suggest the formation of hydrogen-bonding interactions between the flavin and Arg65 in the light-adapted state, providing a plausible explanation for the experimental observation of faster photoinduced PCET in this state. These characteristics of the light-adapted state may also be essential for the long-range signaling capabilities of this photoreceptor protein., Competing Interests: The authors declare no competing interest.

The coupling between electrons and protons and the long-range transport of protons play important roles throughout biology. Biomimetic systems derived from benzimidazole-phenol (BIP) constructs have been designed to undergo proton-coupled electron transfer (PCET) upon electrochemical or photochemical oxidation. Moreover, these systems can transport protons along hydrogen-bonded networks or proton wires through multiproton PCET. Herein, the nonequilibrium dynamics of both single and double proton transfer in BIP molecules initiated by oxidation are investigated with first-principles molecular dynamics simulations. Although these processes are concerted in that no thermodynamically stable intermediate is observed, the simulations predict that they are predominantly asynchronous on the ultrafast time scale. For both systems, the first proton transfer typically occurs ∼100 fs after electron transfer. For the double proton transfer system, typically the second proton transfer occurs hundreds of femtoseconds after the initial proton transfer. A machine learning algorithm was used to identify the key molecular vibrational modes essential for proton transfer: a slow, in-plane bending mode that dominates the overall inner-sphere reorganization, the proton donor-acceptor motion that leads to vibrational coherence, and the faster donor-hydrogen stretching mode. The asynchronous double proton transfer mechanism can be understood in terms of a significant mode corresponding to the two anticorrelated proton donor-acceptor motions, typically decreasing only one donor-acceptor distance at a time. Although these PCET processes appear concerted on the time scale of typical electrochemical experiments, attaching these BIP constructs to photosensitizers may enable the detection of the asynchronicity of the electron and multiple proton transfers with ultrafast two-dimensional spectroscopy. Understanding the fundamental PCET mechanisms at this level will guide the design of PCET systems for catalysis and energy conversion processes., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

Designing molecular platforms for controlling proton and electron movement in artificial photosynthetic systems is crucial to efficient catalysis and solar energy conversion. The transfer of both protons and electrons during a reaction is known as proton-coupled electron transfer (PCET) and is used by nature in myriad ways to provide low overpotential pathways for redox reactions and redox leveling, as well as to generate bioenergetic proton currents. Herein, we describe theoretical and electrochemical studies of a series of bioinspired benzimidazole-phenol (BIP) derivatives and a series of dibenzimidazole-phenol (BI 2 P) analogs with each series bearing the same set of terminal proton-accepting (TPA) groups. The set of TPAs spans more than 6 p K units. These compounds have been designed to explore the role of the bridging benzimidazole(s) in a one-electron oxidation process coupled to intramolecular proton translocation across either two (the BIP series) or three (the BI a units. These compounds have been designed to explore the role of the bridging benzimidazole(s) in a one-electron oxidation process coupled to intramolecular proton translocation across either two (the BIP series) or three (the BI 2 P series) acid/base sites. These molecular constructs feature an electrochemically active phenol connected to the TPA group through a benzimidazole-based bridge, which together with the phenol and TPA group form a covalent framework supporting a Grotthuss-type hydrogen-bonded network. Infrared spectroelectrochemistry demonstrates that upon oxidation of the phenol, protons translocate across this well-defined hydrogen-bonded network to a TPA group. The experimental data show the benzimidazole bridges are non-innocent participants in the PCET process in that the addition of each benzimidazole unit lowers the redox potential of the phenoxyl radical/phenol couple by 60 mV, regardless of the nature of the TPA group. Using a series of hypothetical thermodynamic steps, density functional theory calculations correctly predicted the dependence of the redox potential of the phenoxyl radical/phenol couple on the nature of the final protonated species and provided insight into the thermodynamic role of dibenzimidazole units in the PCET process. This information is crucial for developing molecular "dry proton wires" with these moieties, which can transfer protons via a Grotthuss-type mechanism over long distances without the intervention of water molecules., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Electric fields control chemical reactivity in a wide range of systems, including enzymes and electrochemical interfaces. Characterizing the electric fields at electrode-solution interfaces is critical for understanding heterogeneous catalysis and associated energy conversion processes. To address this challenge, recent experiments have probed the response of the nitrile stretching frequency of 4-mercaptobenzonitrile (4-MBN) attached to a gold electrode to changes in the solvent and applied electrode potential. Herein, this system is modeled with periodic density functional theory using a multilayer dielectric continuum treatment of the solvent and at constant applied potentials. The impact of the solvent dielectric constant and the applied electrode potential on the nitrile stretching frequency computed with a grid-based method is in qualitative agreement with the experimental data. In addition, the interfacial electrostatic potentials and electric fields as a function of applied potential were calculated directly with density functional theory. Substantial spatial inhomogeneity of the interfacial electric fields was observed, including oscillations in the region of the molecular probe attached to the electrode. These simulations highlight the microscopic inhomogeneity of the electric fields and the role of molecular polarizability at electrode-solution interfaces, thereby demonstrating the limitations of mean-field models and providing insights relevant to the interpretation of vibrational Stark effect experiments., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

The selective reduction of O 2 , typically with the goal of forming H 2 O, represents a long-standing challenge in the field of catalysis. Macrocyclic transition-metal complexes, and cobalt porphyrins in particular, have been the focus of extensive study as catalysts for this reaction. Here, we show that the mononuclear Co-tetraarylporphyrin complex, Co(por OMe ) (por OMe = meso-tetra(4-methoxyphenyl)porphyrin), catalyzes either 2e - /2H + or 4e - /4H + reduction of O 2 with high selectivity simply by changing the identity of the Brønsted acid in dimethylformamide (DMF). The thermodynamic potentials for O 2 reduction to H 2 O 2 or H 2 O in DMF are determined and exhibit a Nernstian dependence on the acid p K a , while the Co III/II redox potential is independent of the acid p K a . The reaction product, H 2 O or H 2 O 2 , is defined by the relationship between the thermodynamic potential for O 2 reduction to H 2 O 2 and the Co III/II redox potential: selective H 2 O 2 formation is observed when the Co III/II potential is below the O 2 /H 2 O 2 potential, while H 2 O formation is observed when the Co III/II potential is above the O 2 /H 2 O 2 potential. Mechanistic studies reveal that the reactions generating H 2 O 2 and H 2 O exhibit different rate laws and catalyst resting states, and these differences are manifested as different slopes in linear free energy correlations between the log(rate) versus p K a and log(rate) versus effective overpotential for the reactions. This work shows how scaling relationships may be used to control product selectivity, and it provides a mechanistic basis for the pursuit of molecular catalysts that achieve low overpotential reduction of O 2 to H 2 O., Competing Interests: The authors declare no competing financial interest.

Photoreceptor proteins control vital cellular responses to light. The photocycle of the Slr1694 blue light using flavin photoreceptor is initiated by photoexcitation to a locally excited state within the flavin, followed by electron transfer from Tyr8 to the flavin and a proton relay from Tyr8 to the flavin via an intervening glutamine. Herein, the two-dimensional excited state potential energy surfaces associated with this double proton-transfer reaction are computed using the complete active space self-consistent-field method and multiconfigurational perturbation theory, including the protein and solvent environment with electrostatic embedding. The double proton-transfer reaction was found to be energetically unfavorable in the ground state and locally excited state but energetically favorable in the charge-transfer state corresponding to electron transfer from Tyr8 to the flavin. These results indicate that the proton-coupled electron transfer process is sequential, with electron transfer preceding double proton transfer, and that the double proton-transfer reaction is also sequential, with proton transfer from Tyr8 to Gln50 followed by proton transfer from Gln50 to the flavin. The barrier is lower for the first proton-transfer reaction, and both barriers are significantly influenced by geometrical changes within the active site, particularly the proton donor-acceptor distance as well as the protein environment. These calculations provide insight into the impact of protein reorganization and electrostatics on the excited electronic states prior to and during the double proton-transfer reaction. This interplay between excited states and the environment has implications for other photoreceptor proteins.

Nature employs a Tyr Z -His pair as a redox relay that couples proton transfer to the redox process between P680 and the water oxidizing catalyst in photosystem II. Artificial redox relays composed of different benzimidazole-phenol dyads (benzimidazole models His and phenol models Tyr) with substituents designed to simulate the hydrogen bond network surrounding the Tyr Z -His pair have been prepared. When the benzimidazole substituents are strong proton acceptors such as primary or tertiary amines, theory predicts that a concerted two proton transfer process associated with the electrochemical oxidation of the phenol will take place. Also, theory predicts a decrease in the redox potential of the phenol by ∼300 mV and a small kinetic isotope effect (KIE). Indeed, electrochemical, spectroelectrochemical, and KIE experimental data are consistent with these predictions. Notably, these results were obtained by using theory to guide the rational design of artificial systems and have implications for managing proton activity to optimize efficiency at energy conversion sites involving water oxidation and reduction.

Human DNA polymerase η (Pol η) plays a vital role in protection against skin cancer caused by damage from ultraviolet light. This enzyme rescues stalled replication forks at cyclobutane thymine-thymine dimers (TTDs) by inserting nucleotides opposite these DNA lesions. Residue R61 is conserved in the Pol η enzymes across species, but the corresponding residue, as well as its neighbor S62, is different in other Y-family polymerases, Pol ι and Pol κ. Herein, R61 and S62 are mutated to their Pol ι and Pol κ counterparts. Relative binding free energies of dATP to mutant Pol η•DNA complexes with and without a TTD were calculated using thermodynamic integration. The binding free energies of dATP to the Pol η•DNA complex with and without a TTD are more similar for all of these mutants than for wild-type Pol η, suggesting that these mutations decrease the ability of this enzyme to distinguish between a TTD lesion and undamaged DNA. Molecular dynamics simulations of the mutant systems provide insights into the molecular level basis for the changes in relative binding free energies. The simulations identified differences in hydrogen-bonding, cation-π, and π-π interactions of the side chains with the dATP and the TTD or thymine-thymine (TT) motif. The simulations also revealed that R61 and Q38 act as a clamp to position the dATP and the TTD or TT and that the mutations impact the balance among the interactions related to this clamp. Overall, these calculations suggest that R61 and S62 play key roles in the specificity and effectiveness of Pol η for bypassing TTD lesions during DNA replication. Understanding the basis for this specificity is important for designing drugs aimed at cancer treatment.

NiFe oxyhydroxide materials are highly active electrocatalysts for the oxygen evolution reaction (OER), an important process for carbon-neutral energy storage. Recent spectroscopic and computational studies increasingly support iron as the site of catalytic activity but differ with respect to the relevant iron redox state. A combination of hybrid periodic density functional theory calculations and spectroelectrochemical experiments elucidate the electronic structure and redox thermodynamics of Ni-only and mixed NiFe oxyhydroxide thin-film electrocatalysts. The UV/visible light absorbance of the Ni-only catalyst depends on the applied potential as metal ions in the film are oxidized before the onset of OER activity. In contrast, absorbance changes are negligible in a 25% Fe-doped catalyst up to the onset of OER activity. First-principles calculations of proton-coupled redox potentials and magnetizations reveal that the Ni-only system features oxidation of Ni 2+ to Ni 3+ , followed by oxidation to a mixed Ni 3+/4+ state at a potential coincident with the onset of OER activity. Calculations on the 25% Fe-doped system show the catalyst is redox inert before the onset of catalysis, which coincides with the formation of Fe 4+ and mixed Ni oxidation states. The calculations indicate that introduction of Fe dopants changes the character of the conduction band minimum from Ni-oxide in the Ni-only to predominantly Fe-oxide in the NiFe electrocatalyst. These findings provide a unified experimental and theoretical description of the electrochemical and optical properties of Ni and NiFe oxyhydroxide electrocatalysts and serve as an important benchmark for computational characterization of mixed-metal oxidation states in heterogeneous catalysts.

Blue light using flavin adenine dinucleotide (BLUF) proteins are essential for the light regulation of a variety of physiologically important processes and serve as a prototype for photoinduced proton-coupled electron transfer (PCET). Free-energy simulations elucidate the active site conformations in the AppA (activation of photopigment and puc expression) BLUF domain before and following photoexcitation. The free-energy profile for interconversion between conformations with either Trp104 or Met106 closer to the flavin, denoted Trp in /Met out and Trp out /Met in , reveals that both conformations are sampled on the ground state, with the former thermodynamically favorable by ∼3 kcal/mol. These results are consistent with the experimental observation of both conformations. To analyze the proton relay from Tyr21 to the flavin via Gln63, the free-energy profiles for Gln63 rotation were calculated on the ground state, the locally excited state of the flavin, and the charge-transfer state associated with electron transfer from Tyr21 to the flavin. For the Trp in /Met out conformation, the hydrogen-bonding pattern conducive to the proton relay is not thermodynamically favorable on the ground state but becomes more favorable, corresponding to approximately half of the configurations sampled, on the locally excited state. The calculated energy gaps between the locally excited and charge-transfer states suggest that electron transfer from Tyr21 to the flavin is more facile for configurations conducive to proton transfer. When the active site conformation is not conducive to PCET from Tyr21, Trp104 can directly compete with Tyr21 for electron transfer to the flavin through a nonproductive pathway, impeding the signaling efficiency.

Although the binding of a positron to a neutral atom has not been directly observed experimentally, high-level theoretical methods have predicted that a positron will bind to a neutral atom. In the present study, the binding energies of a positron to lithium, sodium, beryllium, and magnesium, as well as the electron-positron annihilation rates for these systems, are calculated using the reduced explicitly correlated Hartree-Fock (RXCHF) method within the nuclear-electronic orbital (NEO) framework. Due to the lack of explicit electron-positron correlation, NEO Hartree-Fock and full configuration interaction calculations with reasonable electronic and positronic basis sets do not predict positron binding to any of these atoms. In contrast, the RXCHF calculations predict positron binding energies and electron-positron annihilation rates in qualitative agreement with previous highly accurate but computationally expensive stochastic variational method calculations. These results illustrate that the RXCHF method can successfully describe the binding of a positron to a neutral species with no dipole moment. Moreover, the RXCHF method will be computationally tractable for calculating positron binding to molecular systems. The RXCHF approach offers a balance of accuracy and computational tractability for studying these types of positronic systems.

Nonadiabatic dynamics simulations of photoinduced proton-coupled electron transfer (PCET) in a phenol-amine complex in solution were performed. The electronic potential energy surfaces were generated on-the-fly with a hybrid quantum mechanical/molecular mechanical approach that described the solute with a multiconfigurational method in a bath of explicit solvent molecules. The transferring hydrogen nucleus was represented as a quantum mechanical wave function calculated with grid-based methods, and surface hopping trajectories were propagated on the adiabatic electron-proton vibronic surfaces. Following photoexcitation to the excited S1 electronic state, the overall decay to the ground vibronic state was found to be comprised of relatively fast decay from a lower proton vibrational state of S1 to a highly excited proton vibrational state of the ground S0 electronic state, followed by vibrational relaxation within the S0 state. Proton transfer can occur either on the highly excited proton vibrational states of S0 due to small environmental fluctuations that shift the delocalized vibrational wave functions or on the low-energy proton vibrational states of S1 due to solvent reorganization that alters the asymmetry of the proton potential and reduces the proton transfer barrier. The isotope effect arising from replacing the transferring hydrogen with deuterium is predicted to be negligible because hydrogen and deuterium behave similarly in both types of proton transfer processes. Although an isotope effect could be observed for other systems, in general the absence of an isotope effect does not imply the absence of proton transfer in photoinduced PCET systems. This computational approach is applicable to a wide range of other photoinduced PCET processes.

The development of more effective energy conversion processes is critical for global energy sustainability. The design of molecular electrocatalysts for the hydrogen evolution reaction is an important component of these efforts. Proton-coupled electron transfer (PCET) reactions, in which electron transfer is coupled to proton transfer, play an important role in these processes and can be enhanced by incorporating proton relays into the molecular electrocatalysts. Herein nickel porphyrin electrocatalysts with and without an internal proton relay are investigated to elucidate the hydrogen evolution mechanisms and thereby enable the design of more effective catalysts. Density functional theory calculations indicate that electrochemical reduction leads to dearomatization of the porphyrin conjugated system, thereby favoring protonation at the meso carbon of the porphyrin ring to produce a phlorin intermediate. A key step in the proposed mechanisms is a thermodynamically favorable PCET reaction composed of intramolecular electron transfer from the nickel to the porphyrin and proton transfer from a carboxylic acid hanging group or an external acid to the meso carbon of the porphyrin. The C-H bond of the active phlorin acts similarly to the more traditional metal-hydride by reacting with acid to produce H2. Support for the theoretically predicted mechanism is provided by the agreement between simulated and experimental cyclic voltammograms in weak and strong acid and by the detection of a phlorin intermediate through spectroelectrochemical measurements. These results suggest that phlorin species have the potential to perform unique chemistry that could prove useful in designing more effective electrocatalysts.

Rate constant expressions for vibronically nonadiabatic proton transfer and proton-coupled electron transfer reactions are presented and analyzed. The regimes covered include electronically adiabatic and nonadiabatic reactions, as well as high-frequency and low-frequency proton donor-acceptor vibrational modes. These rate constants differ from previous rate constants derived with the cumulant expansion approach in that the logarithmic expansion of the vibronic coupling in terms of the proton donor-acceptor distance includes a quadratic as well as a linear term. The analysis illustrates that inclusion of this quadratic term in the framework of the cumulant expansion framework may significantly impact the rate constants at high temperatures for proton transfer interfaces with soft proton donor-acceptor modes that are associated with small force constants and weak hydrogen bonds. The effects of the quadratic term may also become significant in these regimes when using the vibronic coupling expansion in conjunction with a thermal averaging procedure for calculating the rate constant. In this case, however, the expansion of the coupling can be avoided entirely by calculating the couplings explicitly for the range of proton donor-acceptor distances sampled. The effects of the quadratic term for weak hydrogen-bonding systems are less significant for more physically realistic models that prevent the sampling of unphysical short proton donor-acceptor distances. Additionally, the rigorous relation between the cumulant expansion and thermal averaging approaches is clarified. In particular, the cumulant expansion rate constant includes effects from dynamical interference between the proton donor-acceptor and solvent motions and becomes equivalent to the thermally averaged rate constant when these dynamical effects are neglected. This analysis identifies the regimes in which each rate constant expression is valid and thus will be important for future applications to proton transfer and proton-coupled electron transfer in chemical and biological processes.

With recent developments in simulating nonadiabatic systems to high accuracy, it has become possible to determine how much energy is attributed to nuclear quantum effects beyond zero-point energy. In this work, we calculate the non-relativistic ground-state energies of atomic and molecular systems without the Born-Oppenheimer approximation. For this purpose, we utilize the fixed-node diffusion Monte Carlo method, in which the nodes depend on both the electronic and ionic positions. We report ground-state energies for all systems studied, ionization energies for the first-row atoms and atomization energies for the first-row hydrides. We find the ionization energies of the atoms to be nearly independent of the Born-Oppenheimer approximation, within the accuracy of our results. The atomization energies of molecular systems, however, show small effects of the nonadiabatic coupling between electrons and nuclei.

The nuclear electronic orbital (NEO) reduced explicitly correlated Hartree-Fock (RXCHF) approach couples select electronic orbitals to the nuclear orbital via Gaussian-type geminal functions. This approach is extended to enable the use of a restricted basis set for the explicitly correlated electronic orbitals and an open-shell treatment for the other electronic orbitals. The working equations are derived and the implementation is discussed for both extensions. The RXCHF method with a restricted basis set is applied to HCN and FHF(-) and is shown to agree quantitatively with results from RXCHF calculations with a full basis set. The number of many-particle integrals that must be calculated for these two molecules is reduced by over an order of magnitude with essentially no loss in accuracy, and the reduction factor will increase substantially for larger systems. Typically, the computational cost of RXCHF calculations with restricted basis sets will scale in terms of the number of basis functions centered on the quantum nucleus and the covalently bonded neighbor(s). In addition, the RXCHF method with an odd number of electrons that are not explicitly correlated to the nuclear orbital is implemented using a restricted open-shell formalism for these electrons. This method is applied to HCN(+), and the nuclear densities are in qualitative agreement with grid-based calculations. Future work will focus on the significance of nonadiabatic effects in molecular systems and the further enhancement of the NEO-RXCHF approach to accurately describe such effects.

The nuclear-electronic orbital (NEO) approach treats select nuclei quantum mechanically on the same level as the electrons and includes nonadiabatic effects between the electrons and the quantum nuclei. The practical implementation of this approach is challenging due to the significance of electron-nucleus dynamical correlation. Herein, we present a general extension of the previously developed reduced NEO explicitly correlated Hartree-Fock (RXCHF) approach, in which only select electronic orbitals are explicitly correlated to each quantum nuclear orbital via Gaussian-type geminal functions. Approximations of the electronic exchange between the geminal-coupled electronic orbitals and the other electronic orbitals are also explored. This general approach enables computationally tractable yet accurate calculations on molecular systems with quantum protons. The RXCHF method is applied to the hydrogen cyanide (HCN) and FHF(-) systems, where the proton and all electrons are treated quantum mechanically. For the HCN system, only the two electronic orbitals associated with the CH covalent bond are geminal-coupled to the proton orbital. For the FHF(-) system, only the four electronic orbitals associated with the two FH covalent bonds are geminal-coupled to the proton orbital. For both systems, the RXCHF method produces qualitatively accurate nuclear densities, in contrast to mean field-based NEO approaches. The development and implementation of the RXCHF method provide the framework to perform calculations on systems such as proton-coupled electron transfer reactions, where electron-proton nonadiabatic effects are important.