Orientation of non-spherical protonated water clusters revealed by infrared absorption dichroism

Proton Release Reactions in the Inward H+ Pump NsXeR

Directional ion transport across biological membranes plays a central role in many cellular processes. Elucidating the molecular determinants for vectorial ion transport is key to understanding the functional mechanism of membrane-bound ion pumps. The extensive investigation of the light-driven proton pump bacteriorhodopsin from Halobacterium salinarum(HsBR) enabled a detailed description of outward proton transport. Although the structure of inward-directed proton pumping rhodopsins is very similar to HsBR, little is known about their protonation pathway, and hence, the molecular reasons for the vectoriality of proton translocation remain unclear. Here, we employ a combined experimental and theoretical approach to tracking protonation steps in the light-driven inward proton pump xenorhodopsin from Nanosalina sp. (NsXeR). Time-resolved infrared spectroscopy reveals the transient deprotonation of D220 concomitantly with deprotonation of the retinal Schiff base. Our molecular dynamics simulations support a proton release pathway from the retinal Schiff base via a hydrogen-bonded water wire leading to D220 that could provide a putative gating point for the proton release and with allosteric interactions to the retinal Schiff base. Our findings support the key role of D220 in mediating proton release to the cytoplasmic side and provide evidence that this residue is not the primary proton acceptor of the proton transiently released by the retinal Schiff base.

Proton-transfer spectroscopy beyond the normal-mode scenario

A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200–300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm⁻¹, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm⁻¹ and agrees in location and width with previous experimental mid-infrared spectra.

The spectroscopic signatures of excess protons in HCl solutions are studied by ab initio simulations and THz experiments. Two contributions beyond the normal-mode scenario are identified that reflect proton-waiting and proton-transfer processes.

Ultrafast proton-coupled isomerization in the phototransformation of phytochrome

The biological function of phytochromes is triggered by an ultrafast photoisomerization of the tetrapyrrole chromophore biliverdin between two rings denoted C and D. The mechanism by which this process induces extended structural changes of the protein is unclear. Here we report ultrafast proton-coupled photoisomerization upon excitation of the parent state (Pfr) of bacteriophytochrome Agp2. Transient deprotonation of the chromophore's pyrrole ring D or ring C into a hydrogen-bonded water cluster, revealed by a broad continuum infrared band, is triggered by electronic excitation, coherent oscillations and the sudden electric-field change in the excited state. Subsequently, a dominant fraction of the excited population relaxes back to the Pfr state, while ~35% follows the forward reaction to the photoproduct. A combination of quantum mechanics/molecular mechanics calculations and ultrafast visible and infrared spectroscopies demonstrates how proton-coupled dynamics in the excited state of Pfr leads to a restructured hydrogen-bond environment of early Lumi-F, which is interpreted as a trigger for downstream protein structural changes.

Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin

An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

pH-Dependent Water Clusters in Photoacid Solution: Real-Time Observation by ToF-SIMS at a Submicropore Confined Liquid-Vacuum Interface

Water clusters are ubiquitously formed in aqueous solutions by hydrogen bonding, which is quite sensitive to various environment factors such as temperature, pressure, electrolytes, and pH. Investigation of how the environment has impact on water structure is important for further understanding of the nature of water and the interactions between water and solutes. In this work, pH-dependent water structure changes were studied by monitoring the changes for the size distribution of protonated water clusters by in-situ liquid ToF-SIMS. In combination with a light illumination system, in-situ liquid ToF-SIMS was used to real-time measure the changes of a light-activated organic photoacid under different light illumination conditions. Thus, the proton transfer and pH-mediated water cluster changes were analyzed in real-time. It was found that higher concentration of free protons could lead to a strengthened local hydrogen bonding network as well as relatively larger protonated water clusters in both organic acid and inorganic acid. Besides, the accumulation of protons at the liquid-vacuum interface under light illumination was observed owing to the affinity of organic molecules to the low-pressure gas phase. The application of in-situ liquid ToF-SIMS analysis in combination with in-situ light illumination system opened up an avenue to real-time investigate light-activated reactions. Besides, the results regarding water structure changes in acidic solutions showed important insights in related atmospheric and physiochemical processes.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

S. B. The influence of molecular vicinity (expressed in terms of dielectric constant) on the infrared spectra of embedded species in ices and solid matrices

In this theoretical work we evaluate how the chemical environment influences some features presented in the infrared spectrum, such as band intensities and band location of embedded species in icy matrices. The calculations were performed employing the Polarized Continuum Model (PCM) approach with the second-order Møller–Plesset perturbation theory (MP2) level using the Gaussian 09 package. Here, we simulate the effects of molecular vicinity around embedded species in terms of the effects of the dielectric constant (ε) of the icy and solid samples. Gas phase calculation was also performed for comparison purpose. The investigated embedded single molecules were CO, CO₂, CH₄, NH₃, SO₂ HCOOH, CH₃OH and also H₂O. The results suggest that for most vibrational modes, the strengths of IR bands show an increase with ε, which implies they also decrease with respect to porosity. The frequency shifts showed opposite behavior in relation to the band strengths, with few exceptions. A correlation between calculated band intensities with the band strengths A (taken from literature) was determined and described by a linear function I ∼ 6 × 10¹⁸A [km mol⁻¹], with A in unity of cm per molecule. In addition, an associative exponential function was adjusted to the studied dataset to characterize the evolution of frequency-shift and intensity-shift and band strength ratio as function of the dielectric constant. Since astrophysical ice mantles over cold dust grains can vastly vary in composition in space (having different dielectric constants) they are a challenge to be well characterized. Therefore, this work can help the astrochemistry community to better understand astrophysical ices and its observations in the infrared.

In this theoretical work we evaluate how the chemical environment influences some features presented in the infrared spectrum, such as band intensities and band location of embedded species in icy matrices.

Global and local minima of protonated acetonitrile clusters

Potential energy surfaces of protonated acetonitrile clusters have been explored to locate global and local minima energy structures. The structures are stabilized by strong hydrogen bonds, anti-parallel dimers, dipole–dipole and CH⋯N interactions.

Characterization of Intact Eukaryotic Cells with Subcellular Spatial Resolution by Photothermal-Induced Resonance Infrared Spectroscopy and Imaging

Photothermal-induced resonance (PTIR) spectroscopy and imaging with infrared light has seen increasing application in the molecular spectroscopy of biological samples. The appeal of the technique lies in its capability to provide information about IR light absorption at a spatial resolution better than that allowed by light diffraction, typically below 100 nm. In the present work, we tested the capability of the technique to perform measurements with subcellular resolution on intact eukaryotic cells, without drying or fixing. We demonstrate the possibility of obtaining PTIR images and spectra from the nucleus and multiple organelles with high resolution, better than that allowed by diffraction with infrared light. We obtain particularly strong signal from bands typically assigned to acyl lipids and proteins. We also show that while a stronger signal is obtained from some subcellular structures, other large subcellular components provide a weaker or undetectable PTIR response. The mechanism that underlies such variability in response is presently unclear. We propose and discuss different possibilities, addressing thermomechanical, geometrical, and electrical properties of the sample and the presence of cellular water, from which the difference in response may arise.

Non-Markovian barrier crossing with two-time-scale memory is dominated by the faster memory component

We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times, [Formula: see text] and [Formula: see text] . Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier-crossing time becomes independent of the longer memory time if at least one of the two memory times is larger than the intrinsic diffusion time. When we associate memory effects with coupled degrees of freedom that are orthogonal to a one-dimensional reaction coordinate, this counterintuitive result shows that the faster orthogonal degrees of freedom dominate barrier-crossing kinetics in the non-Markovian limit and that the slower orthogonal degrees become negligible, quite contrary to the standard time-scale separation assumption and with important consequences for the proper setup of coarse-graining procedures in the non-Markovian case. By asymptotic matching and symmetry arguments, we construct a crossover formula for the barrier crossing time that is valid for general multi-exponential memory kernels. This formula can be used to estimate barrier-crossing times for general memory functions for high friction, i.e. in the overdamped regime, as well as for low friction, i.e. in the inertial regime. Typical examples where our results are important include protein folding in the high-friction limit and chemical reactions such as proton-transfer reactions in the low-friction limit.

Acid solvation versus dissociation at "stardust conditions": Reaction sequence matters

Chemical reactions at ultralow temperatures are of fundamental importance to primordial molecular evolution as it occurs on icy mantles of dust nanoparticles or on ultracold water clusters in dense interstellar clouds. As we show, studying reactions in a stepwise manner in ultracold helium nanodroplets by mass-selective infrared (IR) spectroscopy provides an avenue to mimic these "stardust conditions" in the laboratory. In our joint experimental/theoretical study, in which we successively add H₂O molecules to HCl, we disclose a unique IR fingerprint at 1337 cm^-1 that heralds hydronium (H₃O⁺) formation and, thus, acid dissociation generating solvated protons. In stark contrast, no reaction is observed when reversing the sequence by allowing HCl to interact with preformed small embryonic ice-like clusters. Our ab initio simulations demonstrate that not only reaction stoichiometry but also the reaction sequence needs to be explicitly considered to rationalize ultracold chemistry.

Self-Assembled Monolayer Enables Slurry-Coating of Li Anode

Li metal has long been considered as the ultimate anodic material for high-energy-density batteries. Protection of Li metal in electrochemical cycling and in the manufacturing environment is critical for practical applications. Here, we present the passivation of the Li metal-carbon nanotube (CNT) composite with molecular self-assembly of a long-chain aliphatic phosphonic acid. The dynamics of the self-assembly process is investigated with sum-frequency generation spectroscopy (SFG). The aliphatic phosphonic acid molecules self-assemble on the Li metal surface via the lithium phosphate bonding, while the well-aligned long chains of the molecules help to prevent corrosion of lithium by oxygen and water in the air. As a result, the self-assembled monolayer (SAM) passivated Li-CNT composite displays excellent stability in dry or even humid air, and could be slurry-coated with organic solvents. The resulting slurry-coated Li anode exhibits a high Coulombic efficiency of 98.8% under a 33% depth of discharge (DOD) at a 1C rate in full battery cycling. The concept of molecular self-assembly on Li metal and the stability of the resulting SAM layer open vast possibilities of designed reagents for surface passivation of Li, which may eventually pave the way for practical application of Li metal in secondary batteries.

Salt Bridge in Aqueous Solution: Strong Structural Motifs but Weak Enthalpic Effect

Salt bridges are elementary motifs of protein secondary and tertiary structure and are commonly associated with structural driving force that increases stability. Often found on the interface to the solvent, they are highly susceptible to solvent–solute interactions, primarily with water but also with other cosolvents (especially ions). We have investigated the interplay of an Arginine–Aspartic acid salt bridge with simple salt ions in aqueous solution by means of molecular dynamics simulations. Besides structural and dynamical features at equilibrium, we have computed the mean force along the dissociation pathway of the salt bridge. We demonstrate that solvated ions influence the behavior of the salt bridge in a very specific and local way, namely the formation of tight ionic pairs Li⁺/Na⁺–Asp⁻. Moreover, our findings show that the enthalpic relevance of the salt bridge is minor, regardless of the presence of solvated ions.