Vibrational solvatochromism and electrochromism of infrared probe molecules containing C≡O, C≡N, C=O, or C−F vibrational chromophore

Solvatochromic charge model of isonitrile probes for investigating hydrogen-bond dynamics with 2DIR spectroscopy

Isonitrile-derivatized amino acids are emerging as highly effective infrared (IR) probes for investigating the structures and dynamics of hydrogen (H)-bonds. These probes enable the quantification of chemical exchange processes in solute-solvent complexes via two-dimensional IR spectroscopy and hold significant promise for site-specific dynamic studies within proteins. Despite their potential, theoretical models that elucidate the solvatochromism of isonitriles remain underdeveloped. Here, we present the development and validation of a solvatochromic charge model for isonitrile (N≡C) probes. Using density functional theory calculations, we parameterized solvatochromic charges for isonitrile and integrated them into classical molecular dynamics (MD) simulations of β-isocyanoalanine in various solvents, including water and fluorinated alcohols. The model incorporates solvent-induced frequency shifts and accurately reproduces complex experimental line shapes, including asymmetric features from non-Gaussian dynamics. The model successfully reproduced the bimodal distribution of frequency shifts corresponding to free and H-bonded species in alcohols, as well as cross-peaks due to chemical exchange. Achieving reproducibility required long MD trajectories, which were computationally demanding. To manage this, we implemented graphics processing unit acceleration, drastically reducing the computational time and enabling the efficient processing of extensive MD data. While some discrepancies in population ratios suggest the need for refined solvent force field parameters and modeling transition dipole moment variations, the developed solvatochromic model is a reliable tool for studying the solvation dynamics. The model enables more detailed investigations of ultrafast dynamics in solute-solvent complexes and represents important steps toward modeling site-specific dynamics of biomolecules with isonitrile probes.

Evaluating aliphatic CF, CF2, and CF3 groups as vibrational Stark effect reporters

Given the extensive use of fluorination in molecular design, it is imperative to understand the solvation properties of fluorinated compounds and the impact of the C-F bond on electrostatic interactions. Vibrational spectroscopy can provide direct insights into these interactions by using the C-F bond stretching [v(C-F)] as an electric field probe through the vibrational Stark effect (VSE). In this work, we explore the VSE of the three basic patterns of aliphatic fluorination, i.e., mono-, di-, and trifluorination in CF, CF2, and CF3 groups, respectively, and compare their response to the well-studied aromatic v(C-F). Magnitudes (i.e., Stark tuning rates) and orientations of the difference dipole vectors of the v(C-F)-containing normal modes were determined using density functional theory and a molecular dynamics (MD)-assisted solvatochromic analysis of model compounds in solvents of varying polarity. We obtain Stark tuning rates of 0.2-0.8 cm-1/(MV/cm), with smallest and largest electric field sensitivities for CFaliphatic and CF3,aliphatic, respectively. While average electric fields of solvation were oriented along the main symmetry axis of the CFn, and thus along its static dipole, the Stark tuning rate vectors were tilted by up to 87° potentially enabling to map electrostatics in multiple dimensions. We discuss the influence of conformational heterogeneity on spectral shifts and point out the importance of multipolar and/or polarizable MD force fields to describe the electrostatics of fluorinated molecules. The implications of this work are of direct relevance for studies of fluorinated molecules as found in pharmaceuticals, fluorinated peptides, and proteins.

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Vibrational Probe at the Electrochemical Interface: Dependence on Plasmon Coupling and Potential of the Lineshape in Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy of vibrational probes at an electrode surface shows promise for studying the structural dynamics at an active electrochemical interface. This interface is a complex environment where the solution structures in response to the applied potential. A strategy for achieving the necessary monolayer sensitivity is to use a plasmonically active electrode, which enhances the electromagnetic fields that produce the spectroscopic response. Here, we show how the coupling between the plasmon and the vibrations of the molecular monolayer impacts the FTIR and 2D IR spectroscopy, with an emphasis on the electrochemical potential difference spectra. We show how mixing between the vibrational and plasmonic states gives rise to the distortions that are observed in these measurements. This provides an important step toward 2D IR measurements of vibrational probes at the electrochemical interface as a tool for probing the structural dynamics in the double layer.

Effects of the Intramolecular Group and Solvent on Vibrational Coupling Modes and Strengths of Fermi Resonances in Aryl Azides: A DFT Study of 4-Azidotoluene and 4-Azido-N-phenylmaleimide

The high transition dipole strength of the azide asymmetric stretch makes aryl azides good candidates as vibrational probes (VPs). However, aryl azides have complex absorption profiles due to Fermi resonances (FRs). Understanding the origin and the vibrational modes involved in FRs of aryl azides is critically important toward developing them as VPs for studies of protein structures and structural changes in response to their surroundings. As such, we studied vibrational couplings in 4-azidotoluene and 4-azido-N-phenylmaleimide in two solvents, N,N-dimethylacetamide and tetrahydrofuran, to explore the origin and the effects of intramolecular group and solvent on the FRs of aryl azides using density functional theory (DFT) calculations with the B3LYP functional and seven basis sets, 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), and 6-311++G(df,pd). Two combination bands consisting of the azide symmetric stretch and another mode form strong FRs with the azide asymmetric stretch for both molecules. The FR profile was altered by replacing the methyl group with maleimide. Solvents change the relative peak position and intensity more significantly for 4-azido-N-phenylmaleimide, which makes it a more sensitive VP. Furthermore, the DFT results indicate that a comparison among the results from different basis sets can be used as a means to predict more reliable vibrational spectra.

Vibrational Solvatochromism Study of the C-H···O Improper Hydrogen Bond

The improper C-H···O hydrogen bond is an important weak interaction, with broad implications for protein and nucleic acid structure, molecular recognition, enzyme catalysis, and drug interaction. Despite its wide identification in crystal structures, the general existence of C-H···O hydrogen bonds remains elusive especially for natural C-H groups in bulk aqueous solutions at room temperature. Vibrational spectroscopy is a promising methodology to tackle this challenge, as formation of C-H···O hydrogen bonds usually causes shifts of the C-H stretch frequency. Yet, prior observations are inconclusive, as they are all based on a simple blue-shift in aqueous solution and cannot distinguish if it is an effect caused by solvent reorganization or a specific hydrogen-bonding interaction. In this work, we used vibrational solvatochromism as a calibration of the solvent reorganization effect and identified a specific H-bonding interaction. We performed vibrational solvatochromism study of C-H(D) of multiple alcohol molecules including the CH mode of CD₃CH(OH)CD₃ and the CD₃ modes of CD₃OH, CD₃CH₂OH, and CD₃CH(OH)CD₃ in a series of solvents. We found an abnormal blue-shift of the Raman frequency of the C-H and C-D bonds at both the C_α and C_β positions of alcohols in water, which lies in an opposite direction to the expected trend due to vibrational solvatochromism. This experimental evidence supports that the improper C-H···O hydrogen bonds might generally exist between nonpolarized C-H and water in liquid solutions at room temperature.

Theoretical Sum Frequency Generation Spectra of Protein Amide with Surface-Specific Velocity-Velocity Correlation Functions

Vibrational sum frequency generation (vSFG) spectroscopy is widely used to probe the protein structure at interfaces. Because protein vSFG spectra are complex, they can only provide detailed structural information if combined with computer simulations of protein molecular dynamics and spectra calculations. We show how vSFG spectra can be accurately modeled using a surface-specific velocity-velocity scheme based on ab initio normal modes. Our calculated vSFG spectra show excellent agreement with the experimental sum frequency spectrum of LTα14 peptide and provide insight into the origin of the characteristic α-helical amide I peak. Analysis indicates that the peak shape can be explained largely by two effects: (1) the uncoupled response of amide groups located on opposite sides of the α-helix will have different orientations with respect to the interface and therefore different local environments affecting the local mode vibrations and (2) vibrational splitting from nearest neighbor coupling evaluated as inter-residue vibrational correlation. The conclusion is consistent with frequency mapping techniques with an empirically based ensemble of peptide structures, thus showing how time correlation approaches and frequency mapping techniques can give independent yet complementary molecular descriptions of protein vSFG. These models reveal the sensitive relationship between protein structure and their amide I response, allowing exploitation of the complicated molecular vibrations and their interference to derive the structures of proteins under native conditions at interfaces.

Solvent Organization and Electrostatics Tuned by Solute Electronic Structure: Amide versus Non-Amide Carbonyls

The ability to exploit carbonyl groups to measure electric fields in enzymes and other complex reactive environments by using the vibrational Stark effect has inspired growing interest in how these fields can be measured, tuned, and ultimately designed. Previous studies have concentrated on the role of the solvent in tuning the fields exerted on the solute. Here, we explore instead the role of the solute electronic structure in modifying the local solvent organization and electric field exerted on the solute. By measuring the infrared absorption spectra of amide-containing molecules, as prototypical peptides, and contrasting them with non-amide carbonyls in a wide range of solvents, we show that these solutes experience notable differences in their frequency shifts in polar solvents. Using vibrational Stark spectroscopy and molecular dynamics simulations, we demonstrate that while some of these differences can be rationalized by using the distinct intrinsic Stark tuning rates of the solutes, the larger frequency shifts for amides and dimethylurea primarily result from the larger solvent electric fields experienced by their carbonyl groups. These larger fields arise due to their stronger p-π conjugation, which results in larger C═O bond dipole moments that further induce substantial solvent organization. Using electronic structure calculations, we decompose the electric fields into contributions from solvent molecules that are in the first solvation shell and those from the bulk and show that both of these contributions are significant and become larger with enhanced conjugation in solutes. These results show that structural modifications of a solute can be used to tune both the solvent organization and electrostatic environment, indicating the importance of a solute-centric paradigm in modulating and designing the electrostatic environment in condensed-phase chemical processes.

Tutorials in vibrational sum frequency generation spectroscopy. III. Collecting, processing, and analyzing vibrational sum frequency generation spectra

In this Tutorial series, we aim to provide an accessible introduction to vibrational sum frequency generation (VSFG) spectroscopy, targeted toward people entering the VSFG world without a rigorous formal background in optical physics or nonlinear spectroscopy. In this article, we discuss in detail the processes of collecting and processing VSFG data, and user-friendly processing software (sfgtools) is provided for use by people new to the field. Some discussion of analyzing VSFG spectra is also given, specifically with a discussion of fitting homodyne VSFG spectra, and a discussion of what can be learned (both qualitatively and quantitatively) from VSFG spectra.

From Networked to Isolated: Observing Water Hydrogen Bonds in Concentrated Electrolytes with Two-Dimensional Infrared Spectroscopy

Superconcentrated electrolytes have emerged as a promising class of materials for energy storage devices, with evidence that high voltage performance is possible even with water as the solvent. Here, we study the changes in the water hydrogen bonding network induced by the dissolution of lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) in concentrations ranging from the dilute to the superconcentrated regimes. Using time-resolved two-dimensional infrared spectroscopy, we observe the progressive disruption of the water-water hydrogen bond network and the appearance of isolated water molecules interacting only with ions, which can be identified and spectroscopically isolated through the intermolecular cross-peaks between the water and the TFSI^- ions. Analyzing the vibrational relaxation of excitations of the H₂O stretching mode, we observe a transition in the dominant relaxation path as the bulk-like water vanishes and is replaced by ion-solvation water with the rapid single-step relaxation of delocalized stretching vibrations into the low frequency modes being replaced by multistep relaxation through the intramolecular H₂O bend and into the TFSI^- high frequency modes prior to relaxing to the low frequency structural degrees of freedom. These results definitively demonstrate the absence of vibrationally bulk-like water in the presence of high concentrations of LiTFSI and especially in the superconcentrated regime, while additionally revealing aspects of the water hydrogen bond network that have been difficult to discern from the vibrational spectroscopy of the neat liquid.

The Solvation-Induced Onsager Reaction Field Rather than the Double-Layer Field Controls CO2 Reduction on Gold

The selectivity and activity of the carbon dioxide reduction (CO₂R) reaction are sensitive functions of the electrolyte cation. By measuring the vibrational Stark shift of in situ-generated CO on Au in the presence of alkali cations, we quantify the total electric field present at catalytic active sites and deconvolute this field into contributions from (1) the electrochemical Stern layer and (2) the Onsager (or solvation-induced) reaction field. Contrary to recent theoretical reports, the CO₂R kinetics does not depend on the Stern field but instead is closely correlated with the strength of the Onsager reaction field. These results show that in the presence of adsorbed (bent) CO₂, the Onsager field greatly exceeds the Stern field and is primarily responsible for CO₂ activation. Additional measurements of the cation-dependent water spectra using vibrational sum frequency generation spectroscopy show that interfacial solvation strongly influences the CO₂R activity. These combined results confirm that the cation-dependent interfacial water structure and its associated electric field must be explicitly considered for accurate understanding of CO₂R reaction kinetics.

Interfacial dynamics in inverted-headgroup lipid membranes

Inverted-headgroup (choline-phosphate) lipids are synthetic lipids that are not found in nature and are used as model systems to understand the role of headgroup dipole orientation. Recently, studies revealed that the net orientation of interfacial water strongly depends on the headgroup electrostatics, i.e., the charges and dipole generated by the phosphate and the choline groups. In order to characterize interfacial H-bond dynamics, we measured two-dimensional infrared spectra of the ester carbonyl band and performed molecular dynamics simulations in fully hydrated 1,2-dioleoyl-sn-glycero-3-phosphocholine and 2-((2,3-bis(oleoyloxy)propyl)-dimethyl-ammonio)ethyl ethyl phosphate (DOCPe) lipid bilayers. The experiments and simulations suggest that the reverse dipole generated by the inverted-headgroup in DOCPe does not affect the carbonyl H-bond populations or the interfacial water H-bond dynamics. However, while phosphate-associated waters in both lipids appear to show a similar H-bond structure, carbonyl-associated waters are characterized by a slightly disrupted H-bond structure in the DOCPe bilayer, especially within the second hydration shell. Our findings show that changes in net water orientation perturb the water H-bonds at the linker region between the headgroup and the lipid tail, although this perturbation is weak.

Characterization of Acetonitrile Isotopologues as Vibrational Probes of Electrolytes

Acetonitrile has emerged as a solvent candidate for novel electrolyte formulations in metal-ion batteries and supercapacitors. It features a bright local C≡N stretch vibrational mode whose infrared (IR) signature is sensitive to battery-relevant cations (Li⁺, Mg²⁺, Zn²⁺, Ca²⁺) both in pure form and in the presence of water admixture across a full possible range of concentrations from the dilute to the superconcentrated regime. Stationary and time-resolved IR spectroscopy thus emerges as a natural tool to study site-specific intermolecular interactions from the solvent perspective without introducing an extrinsic probe that perturbs solution morphology and may not represent the intrinsic dynamics in these electrolytes. The metal-coordinated acetonitrile, water-separated metal-acetonitrile pair, and free solvent each have a distinct vibrational signature that allows their unambiguous differentiation. The IR band frequency of the metal-coordinated acetonitrile depends on the ion charge density. To study the ion transport dynamics, it is necessary to differentiate energy-transfer processes from structural interconversions in these electrolytes. Isotope labeling the solvent is a necessary prerequisite to separate these processes. We discuss the design principles and choice of the CD₃¹³CN label and characterize its vibrational spectroscopy in these electrolytes. The Fermi resonance between ¹³C≡N and C-D stretches complicates the spectral response but does not prevent its effective utilization. Time-resolved two-dimensional (2D) IR spectroscopy can be performed on a mixture of acetonitrile isotopologues and much can be learned about the structural dynamics of various species in these formulations.

A fixed multi-site interaction charge model for an accurate prediction of the QM/MM interactions

A fixed multi-site interaction charge (FMIC) model was proposed for the accurate prediction of intermolecular electrostatic interactions based on the quantum mechanical linear response of a molecule to an external electric field. In such a model, some additional off-center interaction sites were added for capturing multipole interactions for a given molecule. By multivariate least-square fitting analysis of the calculated QM/MM interactions of a given molecule with the electrostatic environment and the electrostatic potentials of the environment at the pre-defined distributed interaction sites, the FMIC of the molecule was obtained. The model system of CO in myoglobin (Mb) was utilized to demonstrate the derivation of the FMIC. The accuracy of FMIC in predicting the electrostatic interactions between CO and the Mb environment was investigated using 10 000 different Mb-CO configurations generated from the 400 ps QM/MM MD simulation. In comparison to the QM/MM calculations at the B3LYP/aug-cc-pVTZ/ff99SB level, the mean unsigned error (MUE) of the results based on the FMIC model was merely 0.10 kcal mol^-1, and the root mean square error (RMSE) was only 0.13 kcal mol^-1, which are significantly lower than the results predicted by the ESP charge model (MUE = 1.45 kcal mol^-1, and RMSE = 1.7 kcal mol^-1, respectively). The transferability of FMIC was tested by applying the obtained FMIC in the wild type Mb-CO system to the mutants of V68F and H64L Mb-CO systems. The MUEs of the obtained results for 10 000 different configurations are both smaller than 0.2 kcal mol^-1 for the V68F and H64L Mb-CO systems in comparison to the B3LYP/aug-cc-pVTZ/ff99SB calculations, and the RMSEs are also lower than 0.2 kcal mol^-1 for both mutants. The applications of FMIC were extended to model the electrostatic interactions between a hydrogen fluoride molecule and 492 waters in a truncated octahedron box; our study showed that the FMIC could give satisfactory results with a MUE of 0.12 kcal mol^-1 and a RMSE of 0.16 kcal mol^-1 in comparison to the B3LYP/aug-cc-pVDZ/TIP3P calculations for 10 000 different configurations generated using the 10 ns classical MD simulation. Therefore, the FMIC method provides an accurate and efficient tool for predicting intermolecular electrostatic interactions, which can be utilized in the future development of molecular force fields.

Dynamic Duo: Vibrational Sum Frequency Scattering Investigation of pH-Switchable Carboxylic Acid/Carboxylate Surfactants on Nanodroplet Surfaces

Surfactants containing pH-switchable, carboxylic acid moieties are utilized in a variety of environmental, industrial, and biological applications that require controlled stability of hydrophobic droplets in water. For nanoemulsions, kinetically stable oil droplets in water, surface adsorption of the anionic form of the carboxylic acid surfactant stabilizes the droplet, whereas a dominant surface presence of the neutral form leads to destabilization. Through the use of dynamic light scattering, ζ-potential, and vibrational sum frequency scattering spectroscopy (VSFSS), we investigate this mechanism and the relative surface population of the neutral and charged species as pH is adjusted. We find that the relative population of the two surfactant species at the droplet surface is distinctly different than their bulk equilibrium concentrations. The ζ-potential measurements show that the surface concentration of the charged surfactant stays nearly constant throughout the stabilizing pH range. In contrast, VSFSS shows that the neutral carboxylic acid form increasingly adsorbs to the surface with increased acidity. The spectral features of the headgroup vibrational modes confirm this behavior and go further to reveal additional molecular details of their adsorption. A significant hydrogen-bonding interaction occurs between the headgroups that, along with hydrophobic chain-chain interactions, assists in drawing more carboxylic acid surfactant to the interface. The charged surfactant provides the stabilizing force for these droplets, while the neutral surfactant introduces complexity to the interfacial structure as the pH is lowered. The results are significantly different than what has been found for the planar oil/water studies where stabilization of the interface is not a factor.

Testing the Limitations of MD-Based Local Electric Fields Using the Vibrational Stark Effect in Solution: Penicillin G as a Test Case

Noncovalent interactions underlie nearly all molecular processes in the condensed phase from solvation to catalysis. Their quantification within a physically consistent framework remains challenging. Experimental vibrational Stark effect (VSE)-based solvatochromism can be combined with molecular dynamics (MD) simulations to quantify the electrostatic forces in solute-solvent interactions for small rigid molecules and, by extension, when these solutes bind in enzyme active sites. While generalizing this approach toward more complex (bio)molecules, such as the conformationally flexible and charged penicillin G (PenG), we were surprised to observe inconsistencies in MD-based electric fields. Combining synthesis, VSE spectroscopy, and computational methods, we provide an intimate view on the origins of these discrepancies. We observe that the electric fields are correlated to conformation-dependent effects of the flexible PenG side chain, including both the local solvation structure and solute conformational sampling in MD. Additionally, we identified that MD-based electric fields are consistently overestimated in three-point water models in the vicinity of charged groups; this cannot be entirely ameliorated using polarizable force fields (AMOEBA) or advanced water models. This work demonstrates the value of the VSE as a direct method for experiment-guided refinements of MD force fields and establishes a general reductionist approach to calibrating vibrational probes for complex (bio)molecules.

An electrostatic energy-based charge model for molecular dynamics simulation

The interactions of the polar chemical bonds such as C=O and N-H with an external electric field were investigated, and a linear relationship between the QM/MM interaction energies and the electric field along the chemical bond is established in the range of weak to intermediate electrical fields. The linear relationship indicates that the electrostatic interactions of a polar group with its surroundings can be described by a simple model of a dipole with constant moment under the action of an electric field. This relationship is employed to develop a general approach to generating an electrostatic energy-based charge (EEC) model for molecules containing single or multiple polar chemical bonds. Benchmark test studies of this model were carried out for (CH₃)₂-CO and N-methyl acetamide in explicit water, and the result shows that the EEC model gives more accurate electrostatic energies than those given by the widely used charge model based on fitting to the electrostatic potential (ESP) in direct comparison to the energies computed by the QM/MM method. The MD simulations of the electric field at the active site of ketosteroid isomerase based on EEC demonstrated that EEC gave a better representation of the electrostatic interaction in the hydrogen-bonding environment than the Amber14SB force field by comparison with experiment. The current study suggests that EEC should be better suited for molecular dynamics study of molecular systems with polar chemical bonds such as biomolecules than the widely used ESP or RESP (restrained ESP) charge models.

Quantification of the Electrostatic Effect on Redox Potential by Positive Charges in a Catalyst Microenvironment

Charged functional groups in the secondary coordination sphere (SCS) of a heterogeneous nanoparticle or homogeneous electrocatalyst are of growing interest due to enhancements in reactivity that derive from specific interactions that stabilize substrate binding or charged intermediates. At the same time, accurate benchmarking of electrocatalyst systems most often depends on the development of linear free-energy scaling relationships. However, the thermodynamic axis in those kinetic-thermodynamic correlations is most often obtained by a direct electrochemical measurement of the catalyst redox potential and might be influenced by electrostatic effects of a charged SCS. In this report, we systematically probe positive charges in a SCS and their electrostatic contributions to the electrocatalyst redox potential. A series of 11 iron carbonyl clusters modified with charged and uncharged ligands was probed, and a linear correlation between the ν_CO absorption band energy and electrochemical redox potentials is observed except where the SCS is positively charged.

Spectroscopic Study of Sol-Gel Entrapped Triruthenium Dodecacarbonyl Catalyst Reveals Hydride Formation

Triruthenium dodecacarbonyl exhibits increased catalytic activity toward hydrogenation reactions when encapsulated in alumina sol-gels. In this study, we demonstrate structural and electronic changes induced by the encapsulation process. Fourier transform infrared (FTIR) spectroscopy reveals that the carbonyl vibrational modes dramatically red shift during aging in the sol-gel glass. These shifts are attributed to the formation of the metal hydride: [HRu₃(CO)₁₁]^-. A comparison to the FTIR spectrum of synthesized [NEt₄][HRu₃(CO)₁₁] confirms this assignment. XPS studies show that the Ru 3d_5/2 peak of [HRu₃(CO)₁₁]^- also shifts to lower binding energy, consistent with an increased electron density on the Ru nuclei compared to Ru₃(CO)₁₂ and confirmed by density functional calculations. This study should open the door to further investigations into the hydride's role in the previously observed catalytic activity. To the best of our knowledge, this is the first study to identify the presence of [HRu₃(CO)₁₁]^- in the alumina sol-gel.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Photo-initiated rupture of azobenzene micelles to enable the spectroscopic analysis of antimicrobial peptide dynamics

Antimicrobial peptides (AMPs) show promise for the treatment of bacterial infections, but many have undesired hemolytic activities. The AMP MP1 not only has broad spectrum bactericidal activity, but has been shown to have antitumor activity. The interaction between AMPs and cellular membranes gives rise to a peptide's cell-specificity and activity. However, direct analysis of the biophysical interactions between peptides and membrane is complex, in part due to the nature of membrane environments as well as structural changes in the peptide that occurs upon binding to the membrane. In order to investigate the interplay between cell selectivity, activity, and secondary structural changes involved in antimicrobial peptide activity, we sought to implement photolizable membrane mimics to assess the types of information available from infrared spectroscopic measurements that follow from photoinitiated peptide dynamics. Azo-surfactants (APEG) form micelles containing a photolizable azobenzene core, which upon irradiation can induce membrane deformation resulting in breakdown of micelles. Spectroscopic analysis of membrane deformation may provide insights into the physical behavior associated with unfolding and dissociation of antimicrobial peptides within a membrane environment. Herein, we synthesized and characterized two new azo-surfactants, APEGTMG and APEGNEt2MeI. Furthermore, we demonstrate the viability of azosurfactants as membrane mimics by examining both the membrane binding and dissociation induced secondary structural changes of the antimicrobial peptide, MP1.

Optical mapping of biological water in single live cells by stimulated Raman excited fluorescence microscopy

Water is arguably the most common and yet least understood material on Earth. Indeed, the biophysical behavior of water in crowded intracellular milieu is a long-debated issue. Understanding of the spatial and compositional heterogeneity of water inside cells remains elusive, largely due to a lack of proper water-sensing tools with high sensitivity and spatial resolution. Recently, stimulated Raman excited fluorescence (SREF) microscopy was reported as the most sensitive vibrational imaging in the optical far field. Herein we develop SREF into a water-sensing tool by coupling it with vibrational solvatochromism. This technique allows us to directly visualize spatially-resolved distribution of water states inside single mammalian cells. Qualitatively, our result supports the concept of biological water and reveals intracellular water heterogeneity between nucleus and cytoplasm. Quantitatively, we unveil a compositional map of the water pool inside living cells. Hence we hope SREF will be a promising tool to study intracellular water and its relationship with cellular activities.

Mapping the amide-I vibrations of model dipeptides with secondary structure sensitivity and amino acid residue specificity, and its application to amyloid β-peptide in aqueous solution

Vibrational spectroscopy has been known as particularly well-suited for deciphering the polypeptide's structure. To decode structural information encoded in IR spectra, we developed amide-I frequency maps on the basis of model dipeptides to correlate the amide-I frequency of interest to the combination of the calculated secondary structure dependent amide-I frequency by using DFT method and the electrostatic potentials that projected onto the amide unit from the micro-environment within molecular mechanics force field. The constructed maps were applied to model dipeptides and amyloid β-peptide fragment (Aβ_25-35). The dipeptide specified map (DS map) and the hybrid map (HYB map) predicted amide-I bands of Aβ_25-35 in solution satisfactorily reproduce experimental observation, and indicate the preference of forming β-sheet and random coil structure for Aβ_25-35 in D₂O just as the results of cluster analysis suggested. These maps with secondary structural sensitivity and amino acid residue specificity open up a way for the interpretation of amide-I vibrations and show their potentials in the understanding of molecular structure of polypeptides in solution.

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.

Vibrational Approach to the Dynamics and Structure of Protein Amyloids

Amyloid diseases, including neurodegenerative diseases such as Alzheimer’s and Parkinson’s, are linked to a poorly understood progression of protein misfolding and aggregation events that culminate in tissue-selective deposition and human pathology. Elucidation of the mechanistic details of protein aggregation and the structural features of the aggregates is critical for a comprehensive understanding of the mechanisms of protein oligomerization and fibrillization. Vibrational spectroscopies, such as Fourier transform infrared (FTIR) and Raman, are powerful tools that are sensitive to the secondary structure of proteins and have been widely used to investigate protein misfolding and aggregation. We address the application of the vibrational approaches in recent studies of conformational dynamics and structural characteristics of protein oligomers and amyloid fibrils. In particular, introduction of isotope labelled carbonyl into a peptide backbone, and incorporation of the extrinsic unnatural amino acids with vibrational moieties on the side chain, have greatly expanded the ability of vibrational spectroscopy to obtain site-specific structural and dynamic information. The applications of these methods in recent studies of protein aggregation are also reviewed.

Site-specific 2D IR spectroscopy: a general approach for the characterization of protein dynamics with high spatial and temporal resolution

The conformational heterogeneity and dynamics of protein side chains contribute to function, but investigating exactly how is hindered by experimental challenges arising from the fast timescales involved and the spatial heterogeneity of protein structures. The potential of two-dimensional infrared (2D IR) spectroscopy for measuring conformational heterogeneity and dynamics with unprecedented spatial and temporal resolution has motivated extensive effort to develop amino acids with functional groups that have frequency-resolved absorptions to serve as probes of their protein microenvironments. We demonstrate the full advantage of the approach by selective incorporation of the probe p-cyanophenylalanine at six distinct sites in a Src homology 3 domain and the application of 2D IR spectroscopy to site-specifically characterize heterogeneity and dynamics and their contribution to cognate ligand binding. The approach revealed a wide range of microenvironments and distinct responses to ligand binding, including at the three adjacent, conserved aromatic residues that form the recognition surface of the protein. Molecular dynamics simulations performed for all the labeled proteins provide insight into the underlying heterogeneity and dynamics. Similar application of 2D IR spectroscopy and site-selective probe incorporation will allow for the characterization of heterogeneity and dynamics of other proteins, how heterogeneity and dynamics are affected by solvation and local structure, and how they might contribute to biological function.

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.

Strategy for Modeling the Electrostatic Responses of the Spectroscopic Properties of Proteins

For better understanding and more efficient use of the spectroscopic probes (vibrational and NMR) of the local electrostatic situations inside proteins, appropriate modeling of the properties of those probes is essential. The present study is devoted to examining the strategy for constructing such models. A more well-founded derivation than the ones in previous studies is given in constructing the models. Theoretical analyses are conducted on two representative example cases related to proteins, i.e., the peptide group of the main chains and the CO and NO ligands to the Fe²⁺ ion of heme, with careful treatment of the behavior of electrons in the electrostatic responses and with verification of consistency with observable quantities. It is shown that, for the stretching frequencies and NMR chemical shifts, it is possible to construct reasonable electrostatic interaction models that encompass the situations of hydration and uniform electric field environment and thus are applicable also to the cases of nonuniform electrostatic situations, which are highly expected for inside of proteins.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

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.

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.

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.

Simple method to introduce an ester infrared probe into proteins

The ester carbonyl stretching vibration has recently been shown to be a sensitive and convenient infrared (IR) probe of protein electrostatics due to the linear dependence of its frequency on local electric field. While an ester moiety can be easily incorporated into peptides via solid-phase synthesis, currently there is no method available to site-specifically incorporate it into a large protein. Herein, we show that it is possible to use a cysteine alkylation reaction to achieve this goal and demonstrate the feasibility of this simple method by successfully incorporating a methyl ester group (CH₂ COOCH₃ ) into a model peptide (YGGCGG), two amyloid-forming peptides derived from the insulin B chain and Aβ, and bovine serum albumin (BSA). IR results obtained with those peptide and protein systems further confirm the utility of this vibrational probe in monitoring, for example, the structural integrity of amyloid fibrils and ligand binding-induced changes in protein local hydration status.

Observation of Electronic Excitation Transfer Through Light Harvesting Complex II Using Two-Dimensional Electronic-Vibrational Spectroscopy

Light-harvesting complex II (LHCII) serves a central role in light harvesting for oxygenic photosynthesis and is arguably the most important photosynthetic antenna complex. In this work, we present two-dimensional electronic-vibrational (2DEV) spectra of LHCII isolated from spinach, demonstrating the possibility of using this technique to track the transfer of electronic excitation energy between specific pigments within the complex. We assign the spectral bands via comparison with the 2DEV spectra of the isolated chromophores, chlorophyll a and b, and present evidence that excitation energy between the pigments of the complex are observed in these spectra. Finally, we analyze the essential components of the 2DEV spectra using singular value decomposition, which makes it possible to reveal the relaxation pathways within this complex.

Vibrational Stark Effects of Carbonyl Probes Applied to Reinterpret IR and Raman Data for Enzyme Inhibitors in Terms of Electric Fields at the Active Site

IR and Raman frequency shifts have been reported for numerous probes of enzyme transition states, leading to diverse interpretations. In the case of the model enzyme ketosteroid isomerase (KSI), we have argued that IR spectral shifts for a carbonyl probe at the active site can provide a connection between the active site electric field and the activation free energy (Fried et al. Science 2014, 346, 1510-1514). Here we generalize this approach to a much broader set of carbonyl probes (e.g., oxoesters, thioesters, and amides), first establishing the sensitivity of each probe to an electric field using vibrational Stark spectroscopy, vibrational solvatochromism, and MD simulations, and then applying these results to reinterpret data already in the literature for enzymes such as 4-chlorobenzoyl-CoA dehalogenase and serine proteases. These results demonstrate that the vibrational Stark effect provides a general framework for estimating the electrostatic contribution to the catalytic rate and may provide a metric for the design or modification of enzymes. Opportunities and limitations of the approach are also described.

Conformational Heterogeneity in the Michaelis Complex of Lactate Dehydrogenase: An Analysis of Vibrational Spectroscopy Using Markov and Hidden Markov Models

Lactate dehydrogenase (LDH) catalyzes the interconversion of pyruvate and lactate. Recent isotope-edited IR spectroscopy suggests that conformational heterogeneity exists within the Michaelis complex of LDH, and this heterogeneity affects the propensity toward the on-enzyme chemical step for each Michaelis substate. By combining molecular dynamics simulations with Markov and hidden Markov models, we obtained a detailed kinetic network of the substates of the Michaelis complex of LDH. The ensemble-average electric fields exerted onto the vibrational probe were calculated to provide a direct comparison with the vibrational spectroscopy. Structural features of the Michaelis substates were also analyzed on atomistic scales. Our work not only clearly demonstrates the conformational heterogeneity in the Michaelis complex of LDH and its coupling to the reactivities of the substates, but it also suggests a methodology to simultaneously resolve kinetics and structures on atomistic scales, which can be directly compared with the vibrational spectroscopy.

Accurate Calculation of Electric Fields Inside Enzymes

The specific electric field generated by a protease at its active site is considered as an important source of the catalytic power. Accurate calculation of electric field at the active site of an enzyme has both fundamental and practical importance. Measuring site-specific changes of electric field at internal sites of proteins due to, eg, mutation, has been realized by using molecular probes with CO or CN groups in the context of vibrational Stark effect. However, theoretical prediction of change in electric field inside a protein based on a conventional force field, such as AMBER or OPLS, is often inadequate. For such calculation, quantum chemical approach or quantum-based polarizable or polarized force field is highly preferable. Compared with the result from conventional force field, significant improvement is found in predicting experimentally measured mutation-induced electric field change using quantum-based methods, indicating that quantum effect such as polarization plays an important role in accurate description of electric field inside proteins. In comparison, the best theoretical prediction comes from fully quantum mechanical calculation in which both polarization and inter-residue charge transfer effects are included for accurate prediction of electrostatics in proteins.

Site-Specific Characterization of Cytochrome P450cam Conformations by Infrared Spectroscopy

Conformational changes are central to protein function but challenging to characterize with both high spatial and temporal precision. The inherently fast time scale and small chromophores of infrared (IR) spectroscopy are well-suited for characterization of potentially rapidly fluctuating environments, and when frequency-resolved probes are incorporated to overcome spectral congestion, enable characterization of specific sites in proteins. We selectively incorporated p-cyanophenylalanine (CNF) as a vibrational probe at five distinct locations in the enzyme cytochrome P450cam and used IR spectroscopy to characterize the environments in substrate and/or ligand complexes reflecting those in the catalytic cycle. Molecular dynamics (MD) simulations were performed to provide a structural basis for spectral interpretation. Together the experimental and simulation data suggest that the CN frequencies are sensitive to both long-range influences, resulting from the particular location of a residue within the enzyme, as well as short-range influences from hydrogen bonding and packing interactions. The IR spectra demonstrate that the environments and effects of substrate and/or ligand binding are different at each position probed and also provide evidence that a single site can experience multiple environments. This study illustrates how IR spectroscopy, when combined with the spectral decongestion and spatial selectivity afforded by CNF incorporation, provides detailed information about protein structural changes that underlie function.