Water dynamics and sum-frequency generation spectra at electrode/aqueous electrolyte interfaces.

The dynamics of water at interfaces between an electrode and an electrolyte is essential for the transport of redox species and for the kinetics of charge transfer reactions next to the electrode. However, while the effects of electrode potential and ion concentration on the electric double layer structure have been extensively studied, a comparable understanding of dynamical aspects is missing. Interfacial water dynamics presents challenges since it is expected to result from the complex combination of water-water, water-electrode and water-ion interactions. Here we perform molecular dynamics simulations of aqueous NaCl solutions at the interface with graphene electrodes, and examine the impact of both ion concentration and electrode potential on interfacial water reorientational dynamics. We show that for all salt concentrations water dynamics exhibits strongly asymmetric behavior: it slows down at increasingly positively charged electrodes but it accelerates at increasingly negatively charged electrodes. At negative potentials water dynamics is determined mostly by the electrode potential value, but in contrast at positive potentials it is governed both by ion-water and electrode-water interactions. We show how these strikingly different behaviors are determined by the interfacial hydrogen-bond network structure and by the ions' surface affinity. Finally, we indicate how the structural rearrangements impacting water dynamics can be probed via vibrational sum-frequency generation spectroscopy.

Intact carbonic acid is a viable protonating agent for biological bases.

Carbonic acid H2CO3 (CA) is a key constituent of the universal CA/bicarbonate/CO2 buffer maintaining the pH of both blood and the oceans. Here we demonstrate the ability of intact CA to quantitatively protonate bases with biologically-relevant pKas and argue that CA has a previously unappreciated function as a major source of protons in blood plasma. We determine with high precision the temperature dependence of pKa(CA), pKa(T) = -373.604 + 16,500/T + 56.478 ln T. At physiological-like conditions pKa(CA) = 3.45 (I = 0.15 M, 37 °C), making CA stronger than lactic acid. We further demonstrate experimentally that CA decomposition to H2O and CO2 does not impair its ability to act as an ordinary carboxylic acid and to efficiently protonate physiological-like bases. The consequences of this conclusion are far reaching for human physiology and marine biology. While CA is somewhat less reactive than (H+)aq, it is more than 1 order of magnitude more abundant than (H+)aq in the blood plasma and in the oceans. In particular, CA is about 70× more abundant than (H+)aq in the blood plasma, where we argue that its overall protonation efficiency is 10 to 20× greater than that of (H+)aq, often considered to be the major protonating agent there. CA should thus function as a major source for fast in vivo acid-base reactivity in the blood plasma, possibly penetrating intact into membranes and significantly helping to compensate for (H+)aq's kinetic deficiency in sustaining the large proton fluxes that are vital for metabolic processes and rapid enzymatic reactions.

A Model Electron Transfer Reaction in Confined Aqueous Solution.

Liquid water confined within nanometer-sized channels exhibits a strongly reduced local dielectric constant perpendicular to the wall, especially at the interface, and this has been suggested to induce faster electron transfer kinetics at the interface than in the bulk. We study a model electron transfer reaction in aqueous solution confined between graphene sheets with classical molecular dynamics. We show that the solvent reorganization energy is reduced at the interface compared to the bulk, which explains the larger rate constant. However, this facilitated solvent reorganization is due to the partial desolvation by the graphene sheet of the ions involved in the electron transfer and not to a local dielectric constant reduction effect.

Water dynamics at electrified graphene interfaces: a jump model perspective.

The reorientation dynamics of water at electrified graphene interfaces was recently shown [J. Phys. Chem. Lett., 2020, 11, 624-631] to exhibit a surprising and strongly asymmetric behavior: positive electrode potentials slow down interfacial water reorientation, while for increasingly negative potentials water dynamics first accelerates before reaching an extremum and then being retarded for larger potentials. Here we use classical molecular dynamics simulations to determine the molecular mechanisms governing water dynamics at electrified interfaces. We show that changes in water reorientation dynamics with electrode potential arise from the electrified interfaces' impacts on water hydrogen-bond jump exchanges, and can be quantitatively described by the extended jump model. Finally, our simulations indicate that no significant dynamical heterogeneity occurs within the water interfacial layer next to the weakly interacting graphene electrode.

Benzimidazoles as Metal-Free and Recyclable Hydrides for CO2 Reduction to Formate.

We report a novel metal-free chemical reduction of CO2 by a recyclable benzimidazole-based organo-hydride, whose choice was guided by quantum chemical calculations. Notably, benzimidazole-based hydride donors rival the hydride-donating abilities of noble-metal-based hydrides such as [Ru(tpy)(bpy)H]+ and [Pt(depe)2H]+. Chemical CO2 reduction to the formate anion (HCOO-) was carried out in the absence of biological enzymes, a sacrificial Lewis acid, or a base to activate the substrate or reductant. 13CO2 experiments confirmed the formation of H13COO- by CO2 reduction with the formate product characterized by 1H NMR and 13C NMR spectroscopy and ESI-MS. The highest formate yield of 66% was obtained in the presence of potassium tetrafluoroborate under mild conditions. The likely role of exogenous salt additives in this reaction is to stabilize and shift the equilibrium toward the ionic products. After CO2 reduction, the benzimidazole-based hydride donor was quantitatively oxidized to its aromatic benzimidazolium cation, establishing its recyclability. In addition, we electrochemically reduced the benzimidazolium cation to its organo-hydride form in quantitative yield, demonstrating its potential for electrocatalytic CO2 reduction. These results serve as a proof of concept for the electrocatalytic reduction of CO2 by sustainable, recyclable, and metal-free organo-hydrides.

Water Dynamics in the Hydration Shells of Biomolecules.

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Dynamical Disorder in the DNA Hydration Shell.

The reorientation and hydrogen-bond dynamics of water molecules within the hydration shell of a B-DNA dodecamer, which are of interest for many of its biochemical functions, are investigated via molecular dynamics simulations and an analytic jump model, which provide valuable new molecular level insights into these dynamics. Different sources of heterogeneity in the hydration shell dynamics are determined. First, a pronounced spatial heterogeneity is found at the DNA interface and explained via the jump model by the diversity in local DNA interfacial topographies and DNA-water H-bond interactions. While most of the hydration shell is moderately retarded with respect to the bulk, some water molecules confined in the narrow minor groove exhibit very slow dynamics. An additional source of heterogeneity is found to be caused by the DNA conformational fluctuations, which modulate the water dynamics. The groove widening aids the approach of, and the jump to, a new water H-bond partner. This temporal heterogeneity is especially strong in the minor groove, where groove width fluctuations occur on the same time scale as the water H-bond rearrangements, leading to a strong dynamical disorder. The usual simplifying assumption that hydration shell dynamics is much faster than DNA dynamics is thus not valid; our results show that biomolecular conformational fluctuations are essential to facilitate the water motions and accelerate the hydration dynamics in confined groove sites.

Molecules in motion: chemical reaction and allied dynamics in solution and elsewhere.

After my acceptance of the kind invitation from Todd Martínez and Mark Johnson, Co-Editors of this journal, to write this article, I had to decide just how to actually do this, given the existence of a fairly personal and extended autobiographical account of recent vintage detailing my youth, education, and assorted experiences and activities at the University of Colorado, Boulder, and later also at Ecole Normale Supérieure in Paris ( 1 ). In the end, I settled on a differently styled recounting of the adventures with my students, postdocs, collaborators, and colleagues in trying to unravel, comprehend, describe, and occasionally even predict the manifestations and consequences of the myriad assortment of molecular dances that contribute to and govern the rates and mechanisms of chemical reactions in solution (and elsewhere). The result follows.

Non-adiabatic transition probability dependence on conical intersection topography.

We derive a closed form analytical expression for the non-adiabatic transition probability for a distribution of trajectories passing through a generic conical intersection (CI), based on the Landau-Zener equation for the non-adiabatic transition probability for a single straight-line trajectory in the CI's vicinity. We investigate the non-adiabatic transition probability's variation with topographical features and find, for the same crossing velocity, no intrinsic difference in efficiency at promoting non-adiabatic decay between peaked and sloped CIs, a result in contrast to the commonly held view. Any increased efficiency of peaked over sloped CIs is thus due to dynamical effects rather than to any increased transition probability of topographical origin. It is also shown that the transition probability depends in general on the direction of approach to the CI, and that the coordinates' reduced mass can affect the transition probability via its influence on the CI topography in mass-scaled coordinates. The resulting predictions compare well with surface hopping simulation results.

Are there dynamical effects in enzyme catalysis? Some thoughts concerning the enzymatic chemical step.

We offer some thoughts on the much debated issue of dynamical effects in enzyme catalysis, and more specifically on their potential role in the acceleration of the chemical step. Since the term 'dynamics' has been used with different meanings, we find it useful to first return to the Transition State Theory rate constant, its assumptions and the choices it involves, and detail the various sources of deviations from it due to dynamics (or not). We suggest that much can be learned about the key current questions for enzyme catalysis from prior extensive studies of dynamical and other effects in the case of reactions in solution. We analyze dynamical effects both in the neighborhood of the transition state and far from it, together with the situation when quantum nuclear motion is central to the reaction, and we illustrate our discussion with various examples of enzymatic reactions.

Reduction of CO2 to methanol catalyzed by a biomimetic organo-hydride produced from pyridine.

We use quantum chemical calculations to elucidate a viable mechanism for pyridine-catalyzed reduction of CO2 to methanol involving homogeneous catalytic steps. The first phase of the catalytic cycle involves generation of the key catalytic agent, 1,2-dihydropyridine (PyH2). First, pyridine (Py) undergoes a H(+) transfer (PT) to form pyridinium (PyH(+)), followed by an e(-) transfer (ET) to produce pyridinium radical (PyH(0)). Examples of systems to effect this ET to populate PyH(+)'s LUMO (E(0)(calc) ∼ -1.3 V vs SCE) to form the solution phase PyH(0) via highly reducing electrons include the photoelectrochemical p-GaP system (E(CBM) ∼ -1.5 V vs SCE at pH 5) and the photochemical [Ru(phen)3](2+)/ascorbate system. We predict that PyH(0) undergoes further PT-ET steps to form the key closed-shell, dearomatized (PyH2) species (with the PT capable of being assisted by a negatively biased cathode). Our proposed sequential PT-ET-PT-ET mechanism for transforming Py into PyH2 is analogous to that described in the formation of related dihydropyridines. Because it is driven by its proclivity to regain aromaticity, PyH2 is a potent recyclable organo-hydride donor that mimics important aspects of the role of NADPH in the formation of C-H bonds in the photosynthetic CO2 reduction process. In particular, in the second phase of the catalytic cycle, which involves three separate reduction steps, we predict that the PyH2/Py redox couple is kinetically and thermodynamically competent in catalytically effecting hydride and proton transfers (the latter often mediated by a proton relay chain) to CO2 and its two succeeding intermediates, namely, formic acid and formaldehyde, to ultimately form CH3OH. The hydride and proton transfers for the first of these reduction steps, the homogeneous reduction of CO2, are sequential in nature (in which the formate to formic acid protonation can be assisted by a negatively biased cathode). In contrast, these transfers are coupled in each of the two subsequent homogeneous hydride and proton transfer steps to reduce formic acid and formaldehyde.

Biomolecular hydration dynamics: a jump model perspective.

The dynamics of water molecules within the hydration shell surrounding a biomolecule can have a crucial influence on its biochemical function. Characterizing their properties and the extent to which they differ from those of bulk water have thus been long-standing questions. Following a tutorial approach, we review the recent advances in this field and the different approaches which have probed the dynamical perturbation experienced by water in the vicinity of proteins or DNA. We discuss the molecular factors causing this perturbation, and describe how they change with temperature. We finally present more biologically relevant cases beyond the dilute aqueous situation. A special focus is on the jump model for water reorientation and hydrogen bond rearrangement.

Water jump reorientation: from theoretical prediction to experimental observation.

Liquid water is remarkably labile in reorganizing its hydrogen-bond (HB) network through the breaking and forming of HBs. This rapid restructuring, which occurs on the picosecond time scale, is critical not only for many of the pure liquid's special features but also for a range of aqueous media phenomena, including chemical reactions and protein activity. An essential part of the HB network reorganization is water molecule reorientation, which has long been described as Debye rotational diffusion characterized by very small angular displacements. Recent theoretical work, however, has presented a starkly contrasting picture: a sudden, large-amplitude jump mechanism, in which the reorienting water molecule rapidly exchanges HB partners in what amounts to an activated chemical reaction. In this Account, we first briefly review the jump mechanism and then discuss how it is supported by a series of experiments. These studies range from indirect indications to direct characterization of the jumps through pioneering two-dimensional infrared spectroscopy (2D-IR), the power of which accords it a special focus here. The scenarios in which experimental signatures of the jump mechanism are sought increase in complexity throughout the Account, beginning with pure water. Here 2D-IR in combination with theory can give a glimpse of the jumps, but the tell-tale markers are not pronounced. A more fruitful arena is provided by aqueous ionic solutions. The difference between water-water and water-anion HB strengths provides the experimental handle of differing OH stretch frequencies; in favorable cases, the kinetic exchange of a water between these two sites can be monitored. Sole observation of this exchange, however, is insufficient to establish the jump mechanism. Fortunately, 2D-IR with polarized pulses has demonstrated that HB exchange is accompanied by significant angular displacement, with an estimated jump angle similar to theoretical estimates. The Janus-like character of amphiphilic solutes, with their hydrophobic and hydrophilic faces, presents a special challenge for theory and experiment. Here a consensus on the 2D-IR interpretation has not yet been achieved; this lack of accord impedes the understanding of, for example, biochemical solutes and interfaces. We argue that the influence of hydrophobic groups on water jumps is only modest and well accounted for by an excluded volume effect in the HB exchange process. Conversely, hydrophilic groups have an important influence when their HB strength with water differs significantly from that of the water-water HB. The power of 2D-IR is argued to be accompanied by subtleties that can lead to just the opposite and, in our view, erroneous conclusion. We close with a prediction that a hydrophobic surface offers an arena in which the dynamics of "dangling" water OHs, bereft of a HB, could provide a 2D-IR confirmation of water jumps.

Roles of the Lewis acid and base in the chemical reduction of CO2 catalyzed by frustrated Lewis pairs.

We employ quantum chemical calculations to discover how frustrated Lewis pairs (FLP) catalyze the reduction of CO2 by ammonia borane (AB); specifically, we examine how the Lewis acid (LA) and Lewis base (LB) of an FLP activate CO2 for reduction. We find that the LA (trichloroaluminum, AlCl3) alone catalyzes hydride transfer (HT) to CO2 while the LB (trimesitylenephosphine, PMes3) actually hinders HT; inclusion of the LB increases the HT barrier by ∼8 kcal/mol relative to the reaction catalyzed by LAs only. The LB hinders HT by donating its lone pair to the LUMO of CO2, increasing the electron density on the C atom and thus lowering its hydride affinity. Although the LB hinders HT, it nonetheless plays a crucial role by stabilizing the active FLP·CO2 complex relative to the LA dimer, free CO2, and free LB. This greatly increases the concentration of the reactive complex in the form FLP·CO2 and thus increases the rate of reaction. We expect that the principles we describe will aid in understanding other catalytic CO2 reductions.

GreenLight Model 960.

The importance of food safety has resulted in a demand for a more rapid, high-throughput method for total viable count (TVC). The industry standard for TVC determination (ISO 4833:2003) is widely used but presents users with some drawbacks. The method is materials- and labor-intensive, requiring multiple agar plates per sample. More importantly, the method is slow, with 72 h typically required for a definitive result. Luxcel Biosciences has developed the GreenLight Model 960, a microtiter plate-based assay providing a rapid high-throughput method of aerobic bacterial load assessment through analysis of microbial oxygen consumption. Results are generated in 1-12 h, depending on microbial load. The mix and measure procedure allows rapid detection of microbial oxygen consumption and equates oxygen consumption to microbial load (CFU/g), providing a simple, sensitive means of assessing the microbial contamination levels in foods (1). As bacteria in the test sample grow and respire, they deplete O2, which is detected as an increase in the GreenLight probe signal above the baseline level (2). The time required to reach this increase in signal can be used to calculate the CFU/g of the original sample, based on a predetermined calibration. The higher the initial microbial load, the earlier this threshold is reached (1).

Multistep drug intercalation: molecular dynamics and free energy studies of the binding of daunomycin to DNA.

Atomic-scale molecular dynamics and free energy calculations in explicit aqueous solvent are used to study the complex mechanism by which a molecule can intercalate between successive base pairs of the DNA double helix. We have analyzed the intercalation pathway for the anticancer drug daunomycin using two different methods: metadynamics and umbrella sampling. The resulting free energy pathways are found to be consistent with one another and point, within an equilibrium free energy context, to a three-step process. Daunomycin initially binds in the minor groove of DNA. An activated step then leads to rotation of the drug, coupled with DNA deformation that opens a wedge between the base pairs, bends DNA toward the major groove, and forms a metastable intermediate that resembles structures seen within the interfaces between DNA and minor-groove-binding proteins. Finally, crossing a small free energy barrier leads to further rotation of daunomycin and full intercalation of the drug, reestablishing stacking with the flanking base pairs and straightening the double helix.

Reorientation and allied dynamics in water and aqueous solutions.

The reorientation of a water molecule is important for a host of phenomena, ranging over--in an only partial listing--the key dynamic hydrogen-bond network restructuring of water itself, aqueous solution chemical reaction mechanisms and rates, ion transport in aqueous solution and membranes, protein folding, and enzymatic activity. This review focuses on water reorientation and related dynamics in pure water, and for aqueous solutes with hydrophobic, hydrophilic, and amphiphilic character, ranging from tetra-methylurea to halide ions and amino acids. Attention is given to the application of theory, simulation, and experiment in the probing of these dynamics, in usefully describing them, and in assessing the description. Special emphasis is placed on a novel sudden, large-amplitude jump mechanism for water reorientation, which contrasts with the commonly assumed Debye rotational diffusion mechanism, characterized by small-amplitude angular motion. Some open questions and directions for further research are also discussed.

Tracking energy transfer from excited to accepting modes: application to water bend vibrational relaxation.

We extend, via a reformulation in terms of Poisson brackets, the method developed previously (Rey et al., J. Phys. Chem. A, 2009, 113, 8949) allowing analysis of the pathways of an excited molecule's ultrafast vibrational relaxation in terms of intramolecular and intermolecular contributions. In particular we show how to ascertain, through the computation of power and work, which portion of an initial excess molecular energy (e.g. vibrational) is transferred to various degrees of freedom (e.g. rotational, translational) of the excited molecule itself and its neighbors. The particular case of bend excess energy relaxation in pure water is treated in detail, completing the picture reported in the work cited above. It is shown explicitly, within a classical description, that almost all of the initial water bend excitation energy is transferred-either indirectly, via Fermi resonance centrifugal coupling to the bend-excited water's rotation, or directly via intermolecular coupling- to local water librations, only involving molecules in the first two hydration shells of the vibrationally excited water molecule. Finally, it is pointed out that the Poisson bracket formulation can also be applied to elucidate the microscopic character of solvation and rotational dynamics, and should prove useful in developing a quantum treatment for energy flow in condensed phases.

Photoisomerization for a model protonated Schiff base in solution: sloped/peaked conical intersection perspective.

The topographical character of conical intersections (CIs)--either sloped or peaked--has played a fundamental and important role in the discussion of the efficiency of CIs as photochemical "funnels." Here this perspective is employed in connection with a recent study of a model protonated Schiff base (PSB) cis to trans photoisomerization in solution [Malhado et al., J. Phys. Chem. A 115, 3720 (2011)]. In that study, the calculated reduced photochemical quantum yield for the successful production of trans product versus cis reactant in acetonitrile solvent compared to water was interpreted in terms of a dynamical solvent effect related to the dominance, for the acetonitrile case, of S(1) to S(0) nonadiabatic transitions prior to the reaching the seam of CIs. The solvent influence on the quantum yield is here re-examined in the sloped/peaked CI topographical perspective via conversion of the model's two PSB internal coordinates and a nonequilibrium solvent coordinate into an effective branching space description, which is then used to re-analyze the generalized Langevin equation/surface hopping results. The present study supports the original interpretation and enriches it in terms of topographical detail.