Vibrational spectroscopic determination of local solvent electric field, solute-solvent electrostatic interaction energy, and their fluctuation amplitudes

Electric fields and potentials in condensed phases

The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.

Empirical S=O stretch vibrational frequency map

Dimethyl sulfoxide (DMSO) water mixtures have been widely studied due to their unique concentration-dependent bulk properties. Here, we present an empirical spectroscopic map for the sulfinyl (S=O) stretching mode. The model can be used to interpret infrared (IR) absorption and ultrafast two-dimensional infrared (2D IR) spectra and quantify hydrogen bond populations and lifetimes by directly connecting spectroscopic measurements with structures and dynamics from molecular dynamics simulations. The electrostatic map is directly parameterized against experimental absorption spectra in the S=O stretching region (980-1100 cm^-1) of dilute DMSO in water. A comparison of center peak frequencies shows that the map performs well across the entire DMSO concentration range, accurately reproducing the ∼10 cm^-1 red-shift per hydrogen bond observed in the experiments. We further benchmark the map by comparing experimental and simulated 2D IR spectra generated by direct numerical integration of the Schrödinger equation. We expect that this empirical frequency map will provide a quantitative platform for investigating intermolecular interactions, microscopic heterogeneity, and ultrafast dynamics in complex liquid mixtures containing DMSO.

Electrostatic frequency maps for amide-I mode of β-peptide: Comparison of molecular mechanics force field and DFT calculations

The spectroscopy of amide-I vibrations has been widely utilized for the understanding of dynamical structure of polypeptides. For the modeling of amide-I spectra, two frequency maps were built for β-peptide analogue (N-ethylpropionamide, NEPA) in a number of solvents within different schemes (molecular mechanics force field based, GM map; DFT calculation based, GD map), respectively. The electrostatic potentials on the amide unit that originated from solvents and peptide backbone were correlated to the amide-I frequency shift from gas phase to solution phase during map parameterization. GM map is easier to construct with negligible computational cost since the frequency calculations for the samples are purely based on force field, while GD map utilizes sophisticated DFT calculations on the representative solute-solvent clusters and brings insight into the electronic structures of solvated NEPA and its chemical environments. The results show that the maps' predicted amide-I frequencies present solvation environmental sensitivities and exhibit their specific characters with respect to the map protocols, and the obtained vibrational parameters are in satisfactory agreement with experimental amide-I spectra of NEPA in solution phase. Although different theoretical schemes based maps have their advantages and disadvantages, the present maps show their potentials in interpreting the amide-I spectra for β-peptides, respectively.

Docking Simulation of the Binding Interactions of Saxitoxin Analogs Produced by the Marine Dinoflagellate Gymnodinium catenatum to the Voltage-Gated Sodium Channel Nav1.4

Saxitoxin (STX) and its analogs are paralytic alkaloid neurotoxins that block the voltage-gated sodium channel pore (Na_v), impeding passage of Na⁺ ions into the intracellular space, and thereby preventing the action potential in the peripheral nervous system and skeletal muscle. The marine dinoflagellate Gymnodinium catenatum produces an array of such toxins, including the recently discovered benzoyl analogs, for which the mammalian toxicities are essentially unknown. We subjected STX and its analogs to a theoretical docking simulation based upon two alternative tri-dimensional models of the Na_v1.4 to find a relationship between the binding properties and the known mammalian toxicity of selected STX analogs. We inferred hypothetical toxicities for the benzoyl analogs from the modeled values. We demonstrate that these toxins exhibit different binding modes with similar free binding energies and that these alternative binding modes are equally probable. We propose that the principal binding that governs ligand recognition is mediated by electrostatic interactions. Our simulation constitutes the first in silico modeling study on benzoyl-type paralytic toxins and provides an approach towards a better understanding of the mode of action of STX and its analogs.

General Applicable Frequency Map for the Amide-I Mode in β-Peptides

In this work, a general applicable amide-I vibrational frequency map (GA map) for β-peptides in a number of common solvents was constructed, based on a peptide derivative, N-ethylpropionamide (NEPA). The map utilizes force fields at the ab initio computational level to accurately describe molecular structure and solute-solvent interactions, and also force fields at the molecular mechanics level to take into account long-range solute-solvent interactions. The results indicate that the GA map works reasonably for mapping the vibrational frequencies of the amide-I local-modes for β-peptides, holding promises for understanding the complicated infrared spectra of the amide-I mode in β-polypeptides.

Solvent effect on the vibrational spectrum of Michler's ketone. Experimental and theoretical investigations

We examined solvent effect on the IR and Raman spectra of MK in several solvents of different polarity and proticity, for understanding of intermolecular interactions, focusing on solvent effect in detail. It has been found that change of solvent polarity has an ambiguous influence on solvatochromism of MK. We have observed that not only vibrations of carbonyl group are affected by the solvent polarity, but also mode ν(CN) and ν(CC) in IR and Raman spectra of MK. Experimental investigations have been supported by the quantum-mechanical computations to gain more insight into the solvatochromic behavior of Michler's ketone. Calculations have been carried using Kohn-Sham formulation of Density Functional Theory (DFT) and the Polarizable Continuum Model (PCM) was employed to account for solute solvent interactions.

General strategy for the bioorthogonal incorporation of strongly absorbing, solvation-sensitive infrared probes into proteins

A high-sensitivity metal-carbonyl-based IR probe is described that can be incorporated into proteins or other biomolecules in very high yield via Click chemistry. A two-step strategy is demonstrated. First, a methionine auxotroph is used to incorporate the unnatural amino acid azidohomoalanine at high levels. Second, a tricarbonyl (η(5)-cyclopentadienyl) rhenium(I) probe modified with an alkynyl linkage is coupled via the Click reaction. We demonstrate these steps using the C-terminal domain of the ribosomal protein L9 as a model system. An overall incorporation level of 92% was obtained at residue 109, which is a surface-exposed residue. Incorporation of the probe into a surface site is shown not to perturb the stability or structure of the target protein. Metal carbonyls are known to be sensitive to solvation and protein electrostatics through vibrational lifetimes and frequency shifts. We report that the frequencies and lifetimes of this probe also depend on the isotopic composition of the solvent. Comparison of the lifetimes measured in H2O versus D2O provides a probe of solvent accessibility. The metal carbonyl probe reported here provides an easy and robust method to label very large proteins with an amino-acid-specific tag that is both environmentally sensitive and a very strong absorber.

2D IR spectra of cyanide in water investigated by molecular dynamics simulations

Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN(-) solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN(-) molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN(-) and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T1 times are sensitive to the van der Waals ranges on the CN(-) is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm(-1) vs. 14.9 cm(-1)) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.

Amide I IR probing of core and shell hydrogen-bond structures in reverse micelles

The properties of N-methylacetamide (NMA) molecules encapsulated in the reverse micelles (RMs) formed by anionic surfactant aerosol OT (AOT), are studied with vibrational spectroscopy and computation. Vibrational spectra of the amide I′ mode of the fully deuterated NMA-d7 show gradual increase of peak frequencies and line broadening as the size of RMs decreases. Analyses of the spectral features reveal the presence of three states of NMA-d7 that correspond to NMA located in the core of water phase (absorption frequency of 1606 cm–1) and two types of interfacial NMA near the surfactant layer (1620 and 1644 cm–1). In larger RMs with water content w0 = [D2O]/[AOT] ≥ 10, only the first two states are observed, whereas in smaller RMs, the population of the third state grows up to 25 % at w0 = 2. These results indicate the general validity of the two-state core/shell model for the confined aqueous solution of NMA, with small modifications due to the system-dependent solute-interface interaction. However, simulations of small RM systems with w0 ≤ 15 show continuous variations of the population, frequency shifts, and the solute-solvent interaction strengths at solute-interface distance less than 4 Å. Thus, the distinction of solute core/shell states tends to be blurred in small RMs but is still effective in interpreting the average spectroscopic observables.

Molecular mechanics force field-based general map for the solvation effect on amide I probe of peptide in different micro-environments

A general electrostatic potential map based on molecular mechanics force field for modeling the amide I frequency is presented. This map is applied to N-methylacetamide (NMA) and designed to be transferable in different micro-environments. The electrostatic potentials from solvent and peptide side chain are projected on the amide unit of NMA to induce the frequency shift of amide I mode. It is shown that the predicted amide I frequency reproduces the experimental data satisfactorily, especially when NMA in polar solvents. The amide I frequency shift is largely determined by the solvents in aqueous solution while it is dominated by the local structure of peptide in other solvent environments. The map parameters are further applied on NMA-MeOH system and the obtained IR spectra show doublet peak profile with negligible deviation from the experimental data, suggesting the usefulness of this general map for providing information about vibrational parameters of amide motions of peptide in different environments.

Calculations of the electric fields in liquid solutions

The electric field created by a condensed-phase environment is a powerful and convenient descriptor for intermolecular interactions. Not only does it provide a unifying language to compare many different types of interactions, but it also possesses clear connections to experimental observables, such as vibrational Stark effects. We calculate here the electric fields experienced by a vibrational chromophore (the carbonyl group of acetophenone) in an array of solvents of diverse polarities using molecular dynamics simulations with the AMOEBA polarizable force field. The mean and variance of the calculated electric fields correlate well with solvent-induced frequency shifts and band broadening, suggesting Stark effects as the underlying mechanism of these key solution-phase spectral effects. Compared to fixed-charge and continuum models, AMOEBA was the only model examined that could describe nonpolar, polar, and hydrogen bonding environments in a consistent fashion. Nevertheless, we found that fixed-charge force fields and continuum models were able to replicate some results of the polarizable simulations accurately, allowing us to clearly identify which properties and situations require explicit polarization and/or atomistic representations to be modeled properly, and to identify for which properties and situations simpler models are sufficient. We also discuss the ramifications of these results for modeling electrostatics in complex environments, such as proteins.

A strongly absorbing class of non-natural labels for probing protein electrostatics and solvation with FTIR and 2D IR spectroscopies

A series of non-natural infrared probes is reported that consist of a metal-tricarbonyl modified with a -(CH2)n- linker and cysteine-specific leaving group. They can be site-specifically attached to proteins using mutagenesis and similar protocols for EPR spin labels, which have the same leaving group. We characterize the label's frequencies and lifetimes using 2D IR spectroscopy in solvents of varying dielectric. The frequency range spans 10 cm(-1), and the variation in lifetimes ranges from 6 to 19 ps, indicating that these probes are very sensitive to their environments. Also, we attached probes with -(CH2)-, -(CH2)3-, and -(CH2)4- linkers to ubiquitin at positions 6 and 63 and collected spectra in aqueous buffer. The frequencies and lifetimes were correlated for 3C and 4C linkers, as they were in the solvents, but did not correlate for the 1C linker. We conclude that lifetime measures solvation, whereas frequency reflects the electrostatics of the environment, which in the case of the 1C linker is a measure of the protein electrostatic field. We also labeled V71C α-synuclein in buffer and membrane-bound. Unlike most other infrared labels, this label has extremely strong cross sections and thus can be measured with 2D IR spectroscopy at sub-millimolar concentrations. We expect that these labels will find use in studying the structure and dynamics of membrane-bound, aggregated, and kinetically evolving proteins for which high signal-to-noise at low protein concentrations is imperative.

Effect of solvent polarity on the vibrational dephasing dynamics of the nitrosyl stretch in an Fe(II) complex revealed by 2D IR spectroscopy

The vibrational dephasing dynamics of the nitrosyl stretching vibration (ν(NO)) in sodium nitroprusside (SNP, Na2[Fe(CN)5NO]·2H2O) are investigated using two-dimensional infrared (2D IR) spectroscopy. The ν(NO) in SNP acts as a model system for the nitrosyl ligand found in metalloproteins which play an important role in the transportation and detection of nitric oxide (NO) in biological systems. We perform a 2D IR line shape study of the ν(NO) in the following solvents: water, deuterium oxide, methanol, ethanol, ethylene glycol, formamide, and dimethyl sulfoxide. The frequency of the ν(NO) exhibits a large vibrational solvatochromic shift of 52 cm(-1), ranging from 1884 cm(-1) in dimethyl sulfoxide to 1936 cm(-1) in water. The vibrational anharmonicity of the ν(NO) varies from 21 to 28 cm(-1) in the solvents used in this study. The frequency-frequency correlation functions (FFCFs) of the ν(NO) in SNP in each of the seven solvents are obtained by fitting the experimentally obtained 2D IR spectra using nonlinear response theory. The fits to the 2D IR line shape reveal that the spectral diffusion time scale of the ν(NO) in SNP varies from 0.8 to 4 ps and is negatively correlated with the empirical solvent polarity scales. We compare our results with the experimentally determined FFCFs of other charged vibrational probes in polar solvents and in the active sites of heme proteins. Our results suggest that the vibrational dephasing dynamics of the ν(NO) in SNP reflect the fluctuations of the nonhomogeneous electric field created by the polar solvents around the nitrosyl and cyanide ligands. The solute solvent interactions occurring at the trans-CN ligand are sensed through the π-back-bonding network along the Fe-NO bond in SNP.

Vibrational solvatochromism and electrochromism. II. Multipole analysis

Small infrared probe molecules have been widely used to study local electrostatic environment in solutions and proteins. Using a variety of time- and frequency-resolved vibrational spectroscopic methods, one can accurately measure the solvation-induced vibrational frequency shifts and the timescales and amplitudes of frequency fluctuations of such IR probes. Since the corresponding frequency shifts are directly related to the local electric field and its spatial derivatives of the surrounding solvent molecules or amino acids in proteins, one can extract information on local electric field around an IR probe directly from the vibrational spectroscopic results. Here, we show that, carrying out a multipole analysis of the solvatochromic frequency shift, the solvatochromic dipole contribution to the frequency shift is not always the dominant factor. In the cases of the nitrile-, thiocyanato-, and azido-derivatized molecules, the solvatochromic quadrupole contributions to the corresponding stretch mode frequency shifts are particularly large and often comparable to the solvatochromic dipole contributions. Noting that the higher multipole moment-solvent electric field interactions are short range effects in comparison to the dipole interaction, the H-bonding interaction-induced vibrational frequency shift can be caused by such short-range multipole-field interaction effects. We anticipate that the present multipole analysis method specifically developed to describe the solvatochromic vibrational frequency shifts will be useful to understand the intermolecular interaction-induced vibrational property changes and to find out a relationship between vibrational solvatochromism and electrochromism of IR probes in condensed phases.

Solute's perspective on how trimethylamine oxide, urea, and guanidine hydrochloride affect water's hydrogen bonding ability

While the thermodynamic effects of trimethylamine oxide (TMAO), urea, and guanidine hydrochloride (GdnHCl) on protein stability are well understood, the underlying mechanisms of action are less well characterized and, in some cases, even under debate. Herein, we employ the stretching vibration of two infrared (IR) reporters, i.e., nitrile (C≡N) and carbonyl (C═O), to directly probe how these cosolvents mediate the ability of water to form hydrogen bonds with the solute of interest, e.g., a peptide. Our results show that these three agents, despite having different effects on protein stability, all act to decrease the strength of the hydrogen bonds formed between water and the infrared probe. While the behavior of TMAO appears to be consistent with its protein-protecting ability, those of urea and GdnHCl are inconsistent with their role as protein denaturants. The latter is of particular interest as it provides strong evidence indicating that although urea and GdnHCl can perturb the hydrogen-bonding property of water their protein-denaturing ability does not arise from a simple indirect mechanism.