Measuring electrostatic fields in both hydrogen-bonding and non-hydrogen-bonding environments using carbonyl vibrational probes

Electrostatic Field-Induced Oscillator Strength Focusing in Molecules

The development of light-harvesting devices based on molecular materials depends critically on the ability to focus the electronic oscillator strength of molecules into the UV-vis spectral window. Typical molecular chromophores have only about 1% of their total electronic oscillator strength in this spectral region and thus perform at a small fraction of their possible effectiveness. This theoretical study finds that the electronic oscillator strength of polyenes in the UV-vis region may be enhanced by 1 order of magnitude using electrostatic fields, motivating specific experimental studies of oscillator strength focusing. We find scaling relationships between the polyene length, the intensity of the applied field, and the field-induced increase in oscillator strength that are useful for the implementation of light-harvesting strategies based on polyenes.

A Preorganized Electric Field Leads to Minimal Geometrical Reorientation in the Catalytic Reaction of Ketosteroid Isomerase

Electrostatic interactions play a pivotal role in enzymatic catalysis and are increasingly modeled explicitly in computational enzyme design; nevertheless, they are challenging to measure experimentally. Using vibrational Stark effect (VSE) spectroscopy, we have measured electric fields inside the active site of the enzyme ketosteroid isomerase (KSI). These studies have shown that these fields can be unusually large, but it has been unclear to what extent they specifically stabilize the transition state (TS) relative to a ground state (GS). In the following, we use crystallography and computational modeling to show that KSI's intrinsic electric field is nearly perfectly oriented to stabilize the geometry of its reaction's TS. Moreover, we find that this electric field adjusts the orientation of its substrate in the ground state so that the substrate needs to only undergo minimal structural changes upon activation to its TS. This work provides evidence that the active site electric field in KSI is preorganized to facilitate catalysis and provides a template for how electrostatic preorganization can be measured in enzymatic systems.

The Volumetric Source Function: Looking Inside van der Waals Interactions

The study of van der Waals interactions plays a central role in the understanding of bonding across a range of biological, chemical and physical phenomena. The presence of van der Waals interactions can be identified through analysis of the reduced density gradient, a fundamental parameter at the core of Density Functional Theory. An extension of Bader’s Quantum Theory of Atoms in Molecules is developed here through combination with the analysis of the reduced density gradient. Through this development, a new quantum chemical topological tool is presented: the volumetric source function. This technique allows insight into the atomic composition of van der Waals interactions, offering the first route towards applying the highly successful source function to these disperse interactions. A new algorithm has been implemented in the open-source code, CRITIC2, and tested on acetone, adipic and maleic acids molecular crystals, each stabilized by van der Waals interactions. This novel technique for studying van der Waals interactions at an atomic level offers unprecedented opportunities in the fundamental study of intermolecular interactions and molecular design for crystal engineering, drug design and bio-macromolecular processes.

Effect of Oriented External Electric Fields on the Photo and Thermal Isomerization of Azobenzene

Azobenzene is a prototype molecule with potential applications in molecular switches, solar thermal batteries, sensors, photoresponsive membranes, molecular electronics, data storage, and nonlinear optics. Photo and thermal isomerization pathways exhibit different charge-transfer character and dipole moments, implying that the use of electric fields can be used to modulate the reactivity of azobenzene. This article examines the differential effect of orientated electric fields on the rotation and inversion thermal and photoisomerization pathways of azobenzene to explore the feasibility of using electric fields in the design of azobenzene-based molecular devices. Our findings demonstrate that the application of orientated electric fields modifies the accessibility of the S₀/S₁ seam of electronic degeneracy, as well as changes the energetically favored relaxation pathway in the branching space to yield different photoproducts. In addition, we observed strong-field dipole-inversion effects that cause a topographical change in the response of the potential energy surface to the applied field and can result in geometric minima that do not exist under field-free conditions. On the S₀ surface, transition barriers can be modified on the order of ±10 kcal mol^-1, enabling control of thermal isomerization rates.

Calcium-Lipid Interactions Observed with Isotope-Edited Infrared Spectroscopy

Calcium ions bind to lipid membranes containing anionic lipids; however, characterizing the specific ion-lipid interactions in multicomponent membranes has remained challenging because it requires nonperturbative lipid-specific probes. Here, using a combination of isotope-edited infrared spectroscopy and molecular dynamics simulations, we characterize the effects of a physiologically relevant (2 mM) Ca²⁺ concentration on zwitterionic phosphatidylcholine and anionic phosphatidylserine lipids in mixed lipid membranes. We show that Ca²⁺ alters hydrogen bonding between water and lipid headgroups by forming a coordination complex involving the lipid headgroups and water. These interactions distort interfacial water orientations and prevent hydrogen bonding with lipid ester carbonyls. We demonstrate, experimentally, that these effects are more pronounced for the anionic phosphatidylserine lipids than for zwitterionic phosphatidylcholine lipids in the same membrane.

A Critical Evaluation of Vibrational Stark Effect (VSE) Probes with the Local Vibrational Mode Theory

Over the past two decades, the vibrational Stark effect has become an important tool to measure and analyze the in situ electric field strength in various chemical environments with infrared spectroscopy. The underlying assumption of this effect is that the normal stretching mode of a target bond such as CO or CN of a reporter molecule (termed vibrational Stark effect probe) is localized and free from mass-coupling from other internal coordinates, so that its frequency shift directly reflects the influence of the vicinal electric field. However, the validity of this essential assumption has never been assessed. Given the fact that normal modes are generally delocalized because of mass-coupling, this analysis was overdue. Therefore, we carried out a comprehensive evaluation of 68 vibrational Stark effect probes and candidates to quantify the degree to which their target normal vibration of probe bond stretching is decoupled from local vibrations driven by other internal coordinates. The unique tool we used is the local mode analysis originally introduced by Konkoli and Cremer, in particular the decomposition of normal modes into local mode contributions. Based on our results, we recommend 31 polyatomic molecules with localized target bonds as ideal vibrational Stark effect probe candidates.

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.

Unified Model for Photophysical and Electro-Optical Properties of Green Fluorescent Proteins

Green fluorescent proteins (GFPs) have become indispensable imaging and optogenetic tools. Their absorption and emission properties can be optimized for specific applications. Currently, no unified framework exists to comprehensively describe these photophysical properties, namely the absorption maxima, emission maxima, Stokes shifts, vibronic progressions, extinction coefficients, Stark tuning rates, and spontaneous emission rates, especially one that includes the effects of the protein environment. In this work, we study the correlations among these properties from systematically tuned GFP environmental mutants and chromophore variants. Correlation plots reveal monotonic trends, suggesting that all these properties are governed by one underlying factor dependent on the chromophore's environment. By treating the anionic GFP chromophore as a mixed-valence compound existing as a superposition of two resonance forms, we argue that this underlying factor is defined as the difference in energy between the two forms, or the driving force, which is tuned by the environment. We then introduce a Marcus-Hush model with the bond length alternation vibrational mode, treating the GFP absorption band as an intervalence charge transfer band. This model explains all of the observed strong correlations among photophysical properties; related subtopics are extensively discussed in the Supporting Information. Finally, we demonstrate the model's predictive power by utilizing the additivity of the driving force. The model described here elucidates the role of the protein environment in modulating the photophysical properties of the chromophore, providing insights and limitations for designing new GFPs with desired phenotypes. We argue that this model should also be generally applicable to both biological and nonbiological polymethine dyes.

Interfacial H-Bond Dynamics in Reverse Micelles: The Role of Surfactant Heterogeneity

Characterizing the hydrogen bond structure and dynamics at surfactant interfaces is essential for understanding how microscopic interactions translate to bulk microemulsion properties. Heterogeneous blends containing tens or hundreds of surfactants are common in the industry, but the most fundamental studies have been carried out on micelles composed of a single surfactant species. Therefore, the effect of surfactant heterogeneity on the interfacial structure and dynamics remains poorly understood. Here, we use ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations to characterize sub-picosecond solvation dynamics as a function of the surfactant composition in ∼120 nm water-in-oil reverse micelles. We probe the ester carbonyl vibrations of nonionic sorbitan surfactants, which are located precisely at the interface between the polar and nonpolar regions, and as such, report on the interfacial water dynamics. We show a 7% increase in hydrogen bond populations together with a 37% slowdown of interfacial hydrogen bond dynamics in heterogeneous mixtures containing hundreds of species, compared to more uniform compositions. Simulations, which are in semiquantitative agreement with experiments, indicate that structural diversity leads to decreased packing efficiency, which in turn drives water further into the otherwise hydrophobic region. Interestingly, this increase in hydration is accompanied by a slowdown of dynamics, indicating that water molecules solvating surfactants are conformationally constrained. These studies demonstrate that the composition and heterogeneity are key factors in determining interfacial properties.

The key to the yellow-to-cyan tuning in the green fluorescent protein family is polarisation

Computational approaches have to date failed to fully capture the large (about 0.4 eV) excitation energy tuning displayed by the nearly identical anionic chromophore in different green fluorescent protein (GFP) variants. Here, we present a thorough comparative study of a set of proteins in this sub-family, including the most red- (phiYFP) and blue-shifted (mTFP0.7) ones. We employ a classical polarisable embedding through induced dipoles and combine it with time-dependent density functional theory and multireference perturbation theory in order to capture both state-specific induction contributions and the coupling of the polarisation of the protein to the chromophore transition density. The obtained results show that only upon inclusion of both these two effects generated by the mutual polarisation between the chromophore and the protein can the full spectral tuning be replicated. We finally discuss how this mutual polarisation affects the correlation between excitation energies, dipole moment variation, and molecular electrostatic field.

4-Oxoproline as a Site-Specific Infrared Probe: Application To Assess Proline Isomerization and Dimer Formation

Due to its unique structure, proline plays important structural and functional roles in proteins. However, this special amino acid lacks an adequate vibrational mode that can be exploited to probe its local electrostatic and hydration status via infrared spectroscopy. Herein, we show that the C═O stretching vibration of a proline derivative, 4-oxoproline, is sensitive to local environment and hence can be used as a site-specific infrared probe. We further validate this notion by applying this unnatural amino acid to assess the thermodynamics of proline cis-trans isomerization in a peptide environment and examine the amino acid dimer formation in concentrated proline and glycine solutions.

Site-Specific 1D and 2D IR Spectroscopy to Characterize the Conformations and Dynamics of Protein Molecular Recognition

Proteins exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires knowledge of the populated states and thus the experimental tools to characterize them. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution, and 2D IR methods that provide richer information are becoming more routine. Although application of IR spectroscopy for investigation of proteins is challenged by spectral congestion, the issue can be overcome by site-specific introduction of amino acid side chains that have IR probe groups with frequency-resolved absorptions, which furthermore enables selective characterization of different locations in proteins. Here, we briefly introduce the biophysical methods and summarize the current progress toward the study of proteins. We then describe our efforts to apply site-specific 1D and 2D IR spectroscopy toward elucidation of protein conformations and dynamics to investigate their involvement in protein molecular recognition, in particular mediated by dynamic complexes: plastocyanin and its binding partner cytochrome f, cytochrome P450s and substrates or redox partners, and Src homology 3 domains and proline-rich peptide motifs. We highlight the advantages of frequency-resolved probes to characterize specific, local sites in proteins and uncover variation among different locations, as well as the advantage of the fast time scale of IR spectroscopy to detect rapidly interconverting states. In addition, we illustrate the greater insight provided by 2D methods and discuss potential routes for further advancement of the field of biomolecular 2D IR spectroscopy.

Measurement of Electric Fields Experienced by Urea Guest Molecules in the 18-Crown-6/Urea (1:5) Host-Guest Complex: An Experimental Reference Point for Electric-Field-Assisted Catalysis

High-resolution synchrotron and neutron single-crystal diffraction data of 18-crown-6/(pentakis)urea measured at 30 K are combined, with the aim of better appreciating the electrostatics associated with intermolecular interactions in condensed matter. With two 18-crown-6 molecules and five different urea molecules in the crystal, this represents the most ambitious combined X-ray/synchrotron and neutron experimental charge density analysis to date on a cocrystal or host-guest system incorporating such a large number of unique molecules. The dipole moments of the five urea guest molecules in the crystal are enhanced considerably compared to values determined for isolated molecules, and 2D maps of the electrostatic potential and electric field show clearly how the urea molecules are oriented with dipole moments aligned along the electric field exerted by their molecular neighbors. Experimental electric fields in the range of 10-19 GV m^-1, obtained for the five different urea environments, corroborate independent measurements of electric fields in the active sites of enzymes and provide an important experimental reference point for recent discussions focused on electric-field-assisted catalysis.

Interspecies Bombolitins Exhibit Structural Diversity upon Membrane Binding, Leading to Cell Specificity

Bombolitins, a class of peptides produced by bees of the genus Bombus, target and disrupt cellular membranes, leading to lysis. Antimicrobial peptides exhibit various mechanisms of action resulting from the interplay between peptide structure, lipid composition, and cellular target membrane selectivity. Herein, two bombolitins displaying significant amino-acid-sequence similarity, BII and BL6, were assessed for antimicrobial activity as well as correlated dodecylphosphocholine (DPC) micelle binding and membrane-induced peptide conformational changes. Infrared and circular dichroism spectroscopies were used to assess the structure-function relationship of each bombolitin, and the results indicate that BII forms a rigid and helically ordered secondary structure upon binding to DPC micelles, whereas BL6 largely lacks secondary structural order. Moreover, the binding affinity of each peptide to DPC micelles was determined, revealing that BL6 displayed a difference in binding affinity by over two orders of magnitude. Further investigations into the growth-inhibitory activity of the two bombolitins were performed against Escherichia coli and Saccharomyces cerevisiae. Interestingly, BII specifically targeted S. cerevisiae, whereas BL6 more effectively inhibited E. coli growth. Overall, the antimicrobial selectivity and specificity of BII and BL6 are largely dependent on the primary as well as secondary structural content of the peptides and the membrane composition.

Local and Global Electric Field Asymmetry in Photosynthetic Reaction Centers

The origin of unidirectional electron transfer in photosynthetic reaction centers (RCs) has been widely discussed. Despite the high level of structural similarity between the two branches of pigments that participate in the initial electron transfer steps of photosynthesis, electron transfer only occurs along one branch. One possible explanation for this functional asymmetry is the differences in the electrostatic environment between the active and the inactive branches arising from the charges and dipoles of the organized protein structure. We present an analysis of electric fields in the RC of the purple bacterium Rhodobacter sphaeroides using the intrinsic carbonyl groups of the pigments as vibrational reporters whose vibrational frequency shifts can be converted into electric fields based on the vibrational Stark effect and also provide Stark effect data for plant pigments that can be used in future studies. The carbonyl stretches of the isolated pigments show pronounced Stark effects. We use these data, solvatochromism, molecular dynamics simulations, and data in the literature from IR and Raman spectra to evaluate differences in fields at symmetry-related positions, in particular at the 9-keto and 2-acetyl positions of the pigments involved in primary charge separation.

Surface-Driven Keto-Enol Tautomerization: Atomistic Insights into Enol Formation and Stabilization Mechanisms

Tautomerisation of simple carbonyl compounds to their enol counterparts on metal surfaces is envisaged to enable an easier route for hydrogenation of the C=O bond in heterogeneously catalyzed reactions. To understand the mechanisms of enol formation and stabilization over catalytically active metal surfaces, we performed a mechanistic study on keto-enol tautomerization of a monocarbonyl compound acetophenon over Pt(111) surface. By employing infrared reflection adsorption spectroscopy in combination with scanning tunneling microscopy, we found that enol can be formed by building a ketone-enol dimer, in which one molecule in the enol form is stabilized through hydrogen bonding to the carbonyl group of the second ketone molecule. Based on the investigations of the co-adsorption behavior of acetophenone and hydrogen, we conclude that keto-enol tautomerization occurs in the intramolecular process and does not involve hydrogen transfer through the surface hypothesized previously.

Quantifying conformations of ester vibrational probes with hydrogen-bond-induced Fermi resonances

Solvatochromic shifts of local vibrational probes report on the strength of the surrounding electric fields and the probe's hydrogen bonding status. Stretching vibrational mode of the ester carbonyl group is a popular solvatochromic reporter used in the studies of peptides and proteins. Small molecules, used to calibrate the response of the vibrational probes, sometimes involve Fermi resonances (FRs) induced by inter-molecular interactions. In the present work, we focus on the scenario where FR does not appear in the infrared spectrum of the ester carbonyl stretching mode in aprotic solvents; however, it is intensified when a hydrogen bond with the reporter is established. When two molecules form hydrogen bonds to the same carbonyl oxygen atom, FR leads to strong hybridization of the involved modes and splitting of the absorption peak. Spectral overlap between the Fermi doublets associated with singly and doubly hydrogen-bonded carbonyl groups significantly complicates quantifying different hydrogen-bonded conformations. We employed a combination of linear and third-order (2DIR) infrared spectroscopy with chemometrics analysis to reveal the individual line shapes and to estimate the occupations of the hydrogen-bonded conformations in methyl acetate, a model small molecule. We identified a hydrogen-bond-induced FR in complexes of methyl acetate with alcohols and water and found that FR is lifted in larger molecules used for control experiments-cholesteryl stearate and methyl cyanoacetate. Applying this methodology to analyze acetonitrile-water solutions revealed that when dissolved in neat water, methyl acetate occupies a single hydrogen-bonding conformation, which is in contrast to the conclusions of previous studies. Our approach can be generally used when FRs prevent direct quantification of the hydrogen bonding status of the vibrational probe.

New Insights into Quinine-DNA Binding Using Raman Spectroscopy and Molecular Dynamics Simulations

Quinine's ability to bind DNA and potentially inhibit transcription and translation has been examined as a mode of action for its antimalarial activity. UV absorption and fluorescence-based studies have lacked the chemical specificity to develop an unambiguous molecular-level picture of the binding interaction. To address this, we use Raman spectroscopy and molecular dynamics (MD) to investigate quinine-DNA interactions. We demonstrate that quinine's strongest Raman band in the fingerprint region, which derives from a symmetric stretching mode of the quinoline ring, is highly sensitive to the local chemical environment and pH. The frequency shifts observed for this mode in solvents of varying polarity can be explained in terms of the Stark effect using a simple Onsager solvation model, indicating that the vibration reports on the local electrostatic environment. However, specific chemical interactions between the quinoline ring and its environment, such as hydrogen bonding and π-stacking, perturb the frequency of this mode in a more complicated but predictable manner. We use this vibration as a spectroscopic probe to investigate the binding interaction between quinine and DNA. We find that, when the quinoline ring is protonated, quinine weakly intercalates into DNA by forming π-stacking interactions with the base pairs. The Raman spectra indicate that quinine can intercalate into DNA with a ratio reaching up to roughly one molecule per 25 base pairs. Our results are confirmed by MD simulations, which also show that the quinoline ring adopts a t-shaped π-stacking geometry with the DNA base pairs, whereas the quinuclidine head group weakly interacts with the phosphate backbone in the minor groove. We expect that the spectral correlations determined here will enable future studies to probe quinine's antimalarial activities, such as disrupting hemozoin biocrystallization, which is hypothesized to be, among other things, one of its primary modes of action against Plasmodium parasites.

An Ab Initio QM/MM Study of the Electrostatic Contribution to Catalysis in the Active Site of Ketosteroid Isomerase

The electric field in the hydrogen-bond network of the active site of ketosteroid isomerase (KSI) has been experimentally measured using vibrational Stark effect (VSE) spectroscopy, and utilized to study the electrostatic contribution to catalysis. A large gap was found in the electric field between the computational simulation based on the Amber force field and the experimental measurement. In this work, quantum mechanical (QM) calculations of the electric field were performed using an ab initio QM/MM molecular dynamics (MD) simulation and electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method. Our results demonstrate that the QM-derived electric field based on the snapshots from QM/MM MD simulation could give quantitative agreement with the experiment. The accurate calculation of the electric field inside the protein requires both the rigorous sampling of configurations, and a QM description of the electrostatic field. Based on the direct QM calculation of the electric field, we theoretically confirmed that there is a linear correlation relationship between the activation free energy and the electric field in the active site of wild-type KSI and its mutants (namely, D103N, Y16S, and D103L). Our study presents a computational protocol for the accurate simulation of the electric field in the active site of the protein, and provides a theoretical foundation that supports the link between electric fields and enzyme catalysis.

Physiological Calcium Concentrations Slow Dynamics at the Lipid-Water Interface

Phospholipids can interact strongly with ions at physiological concentrations, and these interactions can alter membrane properties. Here, we describe the effects of calcium ions on the dynamics in phospholipid membranes. We used a combination of time-resolved ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations. We found that millimolar Ca²⁺ concentrations lead to slower fluctuations in the local environment at the lipid-water interface of membranes with phosphatidylserine. The effect was only observed in bilayers containing anionic phosphatidylserine; membranes composed of only zwitterionic phosphatidylcholine did not experience a slowdown. Local water dynamics were measured using the ester groups as label-free probes and were found to be up to 50% slower with 2.5 mM Ca²⁺. Molecular dynamics simulations show that Ca²⁺ primarily binds to the carboxylate group of phosphatidylserines. These findings have implications for apoptotic and diseased cells in which phosphatidylserine is exposed to extracellular calcium and for the biophysical effects of divalent cations on lipid bilayers.

Long-Range Modulations of Electric Fields in Proteins

Electrostatic interactions are essential for controlling the protein structure and function. Whereas so far experimental and theoretical efforts focused on the effect of local electrostatics, this work aims at elucidating the long-range modulation of electric fields in proteins upon binding to charged surfaces. The study is based on cytochrome c (Cytc) variants carrying nitrile reporters for the vibrational Stark effect that are incorporated into the protein via genetic engineering and chemical modification. The Cytc variants were thoroughly characterized with respect to possible structural perturbations due to labeling. For the proteins in solution, the relative hydrogen bond occupancy and the calculated electric fields, both obtained from molecular dynamics (MD) simulations, and the experimental nitrile stretching frequencies were used to develop a relationship for separating hydrogen-bonding and non-hydrogen-bonding electric field effects. This relationship provides an excellent description for the stable Cytc variants in solution. For the proteins bound to Au electrodes coated with charged self-assembled monolayers (SAMs), the underlying MD simulations can only account for the electric field changes Δ E_ads due to the formation of the electrostatic SAM-Cytc complexes but not for the additional contribution, Δ E_int, representing the consequences of the potential drops over the electrode/SAM/protein interfaces. Both Δ E_ads and Δ E_int, determined at distances between 20 and 30 Å with respect to the SAM surface, are comparable in magnitude to the non-hydrogen-bonding electric field in the unbound protein. This long-range modulation of the internal electric field may be of functional relevance for proteins in complexes with partner proteins (Δ E_ads) and attached to membranes (Δ E_ads + Δ E_int).

Exploring local solvation environments of a heme protein using the spectroscopic reporter 4-cyano-l-phenylalanine

The vibrational reporter unnatural amino acid (UAA) 4-cyano-l-phenylalanine (pCNF) was genetically incorporated individually at three sites (5, 36, and 78) in the heme protein Caldanaerobacter subterraneus H-NOX to probe local hydration environments. The UAA pCNF was incorporated site-specifically using an engineered, orthogonal tRNA synthetase in E. coli. The ability of all of the pCNF-containing H-NOX proteins to form the ferrous CO, NO, or O₂ ligated and unligated states was confirmed with UV-Vis spectroscopy. The solvation state at each site of the three sites of pCNF incorporation was assessed using temperature-dependent infrared spectroscopy. Specifically, the frequency-temperature line slope (FTLS) method was utilized to show that the nitrile group at site 36 was fully solvated and the nitrile group at site 78 was de-solvated (buried) in the heme pocket. The nitrile group at site 5 was found to be partially solvated suggesting that the nitrile group was involved in moderate strength hydrogen bonds. These results were confirmed by the determination of the X-ray crystal structure of the H-NOX protein construct containing pCNF at site 5.

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

An efficient computational approach for modeling protein electrostatic is developed according to static point-charge model distributions based on the linear-scaling EE-GMFCC (electrostatically embedded generalized molecular fractionation with conjugate caps) quantum mechanical (QM) method. In this approach, the Electrostatic-Potential atomic charges are obtained from ab initio calculation of protein, both polarization and charge transfer effect are taken into consideration. This approach shows a significant improvement in the description of electrostatic potential and solvation energy of proteins comparing with current popular molecular mechanics (MM) force fields. Therefore, it has gorgeous prospect in many applications, including accurate calculations of electric field or vibrational Stark spectroscopy in proteins and predicting protein-ligand binding affinity. It can also be applied in QM/MM calculations or electronic embedding method of ONIOM to provide a better electrostatic environment.

Electrochemical and Electrostatic Cleavage of Alkoxyamines

Alkoxyamines are heat-labile molecules, widely used as an in situ source of nitroxides in polymer and materials sciences. Here we show that the one-electron oxidation of an alkoxyamine leads to a cation radical intermediate that even at room temperature rapidly fragments, releasing a nitroxide and carbocation. Digital simulations of experimental voltammetry and current-time transients suggest that the unimolecular decomposition which yields the "unmasked" nitroxide (TEMPO) is exceedingly rapid and irreversible. High-level quantum computations indicate that the collapse of the alkoxyamine cation radical is likely to yield a neutral nitroxide radical and a secondary phenylethyl cation. However, this fragmentation is predicted to be slow and energetically very unfavorable. To attain qualitative agreement between the experimental kinetics and computational modeling for this fragmentation step, the explicit electrostatic environment within the double layer must be accounted for. Single-molecule break-junction experiments in a scanning tunneling microscope using solvent of low dielectric (STM-BJ technique) corroborate the role played by electrostatic forces on the lysis of the alkoxyamine C-ON bond. This work highlights the electrostatic aspects played by charged species in a chemical step that follows an electrochemical reaction, defines the magnitude of this catalytic effect by looking at an independent electrical technique in non-electrolyte systems (STM-BJ), and suggests a redox on/off switch to guide the cleavage of alkoxyamines at an electrified interface.

Perspective: THz-driven nuclear dynamics from solids to molecules

Recent years have seen dramatic developments in the technology of intense pulsed light sources in the THz frequency range. Since many dipole-active excitations in solids and molecules also lie in this range, there is now a tremendous potential to use these light sources to study linear and nonlinear dynamics in such systems. While several experimental investigations of THz-driven dynamics in solid-state systems have demonstrated a variety of interesting linear and nonlinear phenomena, comparatively few efforts have been made to drive analogous dynamics in molecular systems. In the present Perspective article, we discuss the similarities and differences between THz-driven dynamics in solid-state and molecular systems on both conceptual and practical levels. We also discuss the experimental parameters needed for these types of experiments and thereby provide design criteria for a further development of this new research branch. Finally, we present a few recent examples to illustrate the rich physics that may be learned from nonlinear THz excitations of phonons in solids as well as inter-molecular vibrations in liquid and gas-phase systems.

Electric Fields and Enzyme Catalysis

What happens inside an enzyme's active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists' attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme's active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

Two-Dimensional Infrared Spectroscopy Reveals Cosolvent-Composition-Dependent Crossover in Intermolecular Hydrogen-Bond Dynamics

Cosolvents have versatile composition-dependent applications in chemistry and biology. The simultaneous presence of hydrophobic and hydrophilic groups in dimethyl sulfoxide (DMSO), an industrially important amphiphilic cosolvent, when combined with the unique properties of water, plays key roles in the diverse fields of pharmacology, cryoprotection, and cell biology. Moreover, molecules dissolved in aqueous DMSO exhibit an anomalous concentration-dependent nonmonotonic behavior in stability and activity near a critical DMSO mole fraction of 0.15. An experimental identification of the origin of this anomaly can lead to newer chemical and biological applications. We report a direct spectroscopic observation of the anomalous behavior using ultrafast two-dimensional infrared spectroscopy experiments. Our results demonstrate the cosolvent-concentration-dependent nonmonotonicity arises from nonidentical mechanisms in ultrafast hydrogen-bond-exchange dynamics of water above and below the critical cosolvent concentration. Comparison of experimental and theoretical results provides a molecular-level mechanistic understanding: a distinct difference in the stabilization of the solute through dynamic solute-solvent interactions is the key to the anomalous behavior.

Solvent-Independent Anharmonicity for Carbonyl Oscillators

The physical origins of vibrational frequency shifts have been extensively studied in order to understand noncovalent intermolecular interactions in the condensed phase. In the case of carbonyls, vibrational solvatochromism, MD simulations, and vibrational Stark spectroscopy suggest that the frequency shifts observed in simple solvents arise predominately from the environment's electric field due to the vibrational Stark effect. This is contrary to many previously invoked descriptions of vibrational frequency shifts, such as bond polarization, whereby the bond's force constant and/or partial nuclear charges are altered due to the environment, often illustrated in terms of favored resonance structures. Here we test these hypotheses using vibrational solvatochromism as measured using 2D IR to assess the solvent dependence of the bond anharmonicity. These results indicate that the carbonyl bond's anharmonicity is independent of solvent as tested using hexanes, DMSO, and D₂O and is supported by simulated 2D spectra. In support of the linear vibrational Stark effect, these 2D IR measurements are consistent with the assertion that the Stark tuning rate is unperturbed by the electric field generated by both hydrogen and non-hydrogen bonding environments and further extends the general applicability of carbonyl probes for studying intermolecular interactions.

Correction: Vibrational dynamics and solvatochromism of the label SCN in various solvents and hemoglobin by time dependent IR and 2D-IR spectroscopy

Correction for 'Vibrational dynamics and solvatochromism of the label SCN in various solvents and hemoglobin by time dependent IR and 2D-IR spectroscopy' by Luuk J. G. W. van Wilderen et al., Phys. Chem. Chem. Phys., 2014, 16, 19643-19653.