Nitrile Groups as Vibrational Probes: Calculations of the C≡N Infrared Absorption Line Shape of Acetonitrile in Water and Tetrahydrofuran

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.

Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch

The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.

Infrared Spectroscopy of Liquid Solutions as a Benchmarking Tool of Semiempirical QM Methods: The Case of GFN2-xTB

The accurate description of large molecular systems has triggered the development of new computational methods. Due to the computational cost of modeling large systems, the methods usually require a trade-off between accuracy and speed. Therefore, benchmarking to test the accuracy and precision of the method is an important step in their development. The typical gold standard for evaluating these methods is isolated molecules, because of the low computational cost. However, the advent of high-performance computing has made it possible to benchmark computational methods using observables from more complex systems such as liquid solutions. To this end, infrared spectroscopy provides a suitable set of observables (i.e., vibrational transitions) for liquid systems. Here, IR spectroscopy observables are used to benchmark the predictions of the newly developed GFN2-xTB semiempirical method. Three different IR probes (i.e., N-methylacetamide, benzonitrile, and semiheavy water) in solution are selected for this purpose. The work presented here shows that GFN2-xTB predicts central frequencies with errors of less than 10% in all probes. In addition, the method captures detailed properties of the molecular environment such as weak interactions. Finally, the GFN2-xTB correctly assesses the vibrational solvatochromism for N-methylacetamide and semiheavy water but does not have the accuracy needed to properly describe benzonitrile. Overall, the results indicate not only that GFN2-xTB can be used to predict the central frequencies and their dependence on the molecular environment with reasonable accuracy but also that IR spectroscopy data of liquid solutions provide a suitable set of observables for the benchmarking of computational methods.

M←NCCH3, M-η2-(NC)-CH3, and CN-M-CH3 Prepared by Reactions of Ce, Sm, Eu, and Lu Atoms with Acetonitrile: Matrix Infrared Spectra and Theoretical Calculations

The reactions of laser-ablated Ce, Sm, Eu, and Lu atoms with acetonitrile were studied by matrix infrared spectra in a neon matrix, and M←NCCH₃, M-η²-(NC)-CH₃, and CN-M-CH₃ were identified with isotopic substitution and quantum chemical calculations. The major product is the insertion complex (CN-M-CH₃), while the end-on and side-on complexes (M←NCCH₃ and M-η²-(NC)-CH₃) are also trapped in the matrix. The CCN antisymmetric stretching mode for Ce-η²-(NC)-CH₃ was observed at 1536.9 cm^-1, which is much lower than the same modes observed for other lanthanides. NBO analysis reveals that Ce exhibits a remarkable 4f-orbital contribution in bonding to N and to C, reconfirming an active 4f-orbital contribution of cerium in bonding in the side-on complex, while the 4f contributions of Sm and Eu to the M-N and M-C bonds are much lower and the 4f orbital of Lu is not involved in bonding.

Competing Reaction Pathways in Heterogeneously Catalyzed Hydrogenation of Allyl Cyanide: The Chemical Nature of Surface Species

We present a mechanistic study on the formation of an active ligand layer over Pd(111), turning the catalytic surface highly active and selective in partial hydrogenation of an α,β-unsaturated aldehyde acrolein. Specifically, we investigate the chemical composition of a ligand layer consisting of allyl cyanide deposited on Pd(111) and its dynamic changes under the hydrogenation conditions. On pristine surface, allyl cyanide largely retains its chemical structure and forms a layer of molecular species with the CN bond oriented nearly parallel to the underlying metal. In the presence of hydrogen, the chemical composition of allyl cyanide strongly changes. At 100 K, allyl cyanide transforms to unsaturated imine species, containing the C=C and C=N double bonds. At increasing temperatures, these species undergo two competing reaction pathways. First, the C=C bond become hydrogenated and the stable N-butylimine species are produced. In the competing pathway, the unsaturated imine reacts with hydrogen to fully hydrogenate the imine group and produce butylamine. The latter species are unstable under the hydrogenation reaction conditions and desorb from the surface, while the N-butylimine adsorbates formed in the first reaction pathway remain adsorbed and act as an active ligand layer in selective hydrogenation of acrolein.

Electric Fields Influence Intramolecular Vibrational Energy Relaxation and Line Widths

Intramolecular vibrational energy relaxation (IVR) is fundamentally important to chemical dynamics. We show that externally applied electric fields affect IVR and vibrational line widths by changing the anharmonic couplings and frequency detunings between modes. We demonstrate this effect in benzonitrile for which prior experimental results show a decrease in vibrational line width as a function of applied electric field. We identify three major channels for IVR that depend on electric field. In the dominant channel, the electric field affects the frequency detuning, while in the other two channels, variation of anharmonic couplings as a function of field is the underlying mechanism. Consistent with experimental results, we show that the combination of all channels gives rise to reduced line widths with increasing electric field in benzonitrile. Our results are relevant for controlling IVR with external or internal fields and for gaining a more complete interpretation of line widths of vibrational Stark probes.

An ultrafast vibrational study of dynamical heterogeneity in the protic ionic liquid ethyl-ammonium nitrate. I. Room temperature dynamics

Using ultrafast two-dimensional infrared spectroscopy (2D-IR), a vibrational probe (thiocyanate, SCN^-) was used to investigate the hydrogen bonding network of the protic ionic liquid ethyl-ammonium nitrate (EAN) in comparison to H₂O. The 2D-IR experiments were performed in both parallel (⟨ZZZZ⟩) and perpendicular (⟨ZZXX⟩) polarizations at room temperature. In EAN, the non-Gaussian lineshape in the FTIR spectrum of SCN^- suggests two sub-ensembles. Vibrational relaxation rates extracted from the 2D-IR spectra provide evidence of the dynamical differences between the two sub-ensembles. We support the interpretation of two sub-ensembles with response function simulations of two overlapping bands with different vibrational relaxation rates and, otherwise, similar dynamics. The measured rates for spectral diffusion depend on polarization, indicating reorientation-induced spectral diffusion (RISD). A model of restricted molecular rotation (wobbling in a cone) fully describes the observed spectral diffusion in EAN. In H₂O, both RISD and structural spectral diffusion contribute with similar timescales. This complete characterization of the dynamics at room temperature provides the basis for the temperature-dependent measurements in Paper II of this series.

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.

Ultrafast vibrational dynamics of a trigonal planar anionic probe in ionic liquids (ILs): A two-dimensional infrared (2DIR) spectroscopic investigation

A major impediment limiting the widespread application of ionic liquids (ILs) is their high shear viscosity. Incorporation of a tricyanomethanide (TCM^-) anion in ILs leads to low shear viscosity and improvement of several characteristics suitable for large scale applications. However, properties including interactions of TCM^- with the local environment and dynamics of TCM^- have not been thoroughly investigated. Herein, we have studied the ultrafast dynamics of TCM^- in several imidazolium ILs using linear IR and two-dimensional infrared spectroscopy techniques. The spectral diffusion dynamics of the CN stretching modes of TCM^- in all ILs exhibit a nonexponential behavior with a short time component of ∼2 ps and a long time component spanning ∼9 ps to 14 ps. The TCM^- vibrational probe reports a significantly faster relaxation of ILs compared to those observed previously using linear vibrational probes, such as thiocyanate and selenocyanate. Our results indicate a rapid relaxation of the local ion-cage structure embedding the vibrational probe in the ILs. The faster relaxation suggests that the lifetime of the local ion-cage structure decreases in the presence of TCM^- in the ILs. Linear IR spectroscopic results show that the hydrogen-bonding interaction between TCM^- and imidazolium cations in ILs is much weaker. Shorter ion-cage lifetimes together with weaker hydrogen-bonding interactions account for the low shear viscosity of TCM^- based ILs compared to commonly used ILs. In addition, this study demonstrates that TCM^- can be used as a potential vibrational reporter to study the structure and dynamics of ILs and other molecular systems.

Inhomogeneity of Interfacial Electric Fields at Vibrational Probes on Electrode Surfaces

Electric fields control chemical reactivity in a wide range of systems, including enzymes and electrochemical interfaces. Characterizing the electric fields at electrode-solution interfaces is critical for understanding heterogeneous catalysis and associated energy conversion processes. To address this challenge, recent experiments have probed the response of the nitrile stretching frequency of 4-mercaptobenzonitrile (4-MBN) attached to a gold electrode to changes in the solvent and applied electrode potential. Herein, this system is modeled with periodic density functional theory using a multilayer dielectric continuum treatment of the solvent and at constant applied potentials. The impact of the solvent dielectric constant and the applied electrode potential on the nitrile stretching frequency computed with a grid-based method is in qualitative agreement with the experimental data. In addition, the interfacial electrostatic potentials and electric fields as a function of applied potential were calculated directly with density functional theory. Substantial spatial inhomogeneity of the interfacial electric fields was observed, including oscillations in the region of the molecular probe attached to the electrode. These simulations highlight the microscopic inhomogeneity of the electric fields and the role of molecular polarizability at electrode-solution interfaces, thereby demonstrating the limitations of mean-field models and providing insights relevant to the interpretation of vibrational Stark effect experiments.

Solvation properties of protic ionic liquid–molecular solvent mixtures

In this study, we have investigated the solvation properties of binary mixtures of PILs with molecular solvents. The selected binary solvent systems are the PILs ethylammonium nitrate (EAN) and propylammonium nitrate (PAN) combined with either water, methanol, acetonitrile or DMSO.

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.

Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: III. Dynamics and Spectroscopy

In recent years, interest in carbon capture and sequestration has led to numerous investigations of the ability of ionic liquids to act as recyclable CO₂-sorbent materials. Herein, we investigate the structure and dynamics of a model physisorbing ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate ([C₄C₁Im][PF₆]), from the perspective of CO₂ using two-dimensional (2D) IR spectroscopy and molecular dynamics simulations. A direct comparison of experimentally measured and calculated 2D IR line shapes confirms the validity of the simulations and spectroscopic calculations. Taken together, the simulations and experiments reveal new insights into the interactions of a CO₂ solute with the surrounding ionic liquid and how these interactions manifest in the 2D IR spectra. In particular, higher CO₂ asymmetric stretch vibrational frequencies are associated with softer, less populated solvent cages and lower frequencies are associated with tighter, more highly populated solvent cages. The CO₂ interacts most strongly with the anions, and these interactions persist for more than 1 ns. The second strongest interactions are with the imidazolium cation ring that last 100 ps, and the weakest interactions are with the cation butyl tail that persist for 10 ps. The principal contributors to spectral diffusion of the CO₂ asymmetric stretch vibrational frequency due to the dynamical evolution of the solvent are through Lennard-Jones interactions at short times and electrostatics at long times.

Synthesis of 5-Cyano-Tryptophan as a Two-Dimensional Infrared Spectroscopic Reporter of Structure

A concise synthesis of protected 5-cyano-l-tryptophan (Trp5CN ) has been developed for 2D IR spectroscopic investigations within either peptides or proteins. To assess the potential of differently substituted cyano-tryptophans, several model cyano-indole systems were characterized using IR spectroscopy. Upon assessment of their spectroscopic properties, Trp5CN was integrated into a model peptide sequence, Trp5CN -Gly-Phe4CN , to elucidate its structure. This peptide demonstrates the capability of this probe to capture structural information by 2D IR spectroscopy. The 2D IR spectrum of the peptide in water was simulated to reveal a unique spectral signature resulting from the presence of dipolar coupling. The coupling strength between cyano labels was determined to be 1.4 cm-1 by matching the slopes along the max contour for the simulated and experimental spectrum. Using transition dipole coupling, a distance between the two probes of 13 Å was calculated.

Residue-Specific Dynamics and Local Environmental Changes in Aβ40 Oligomer and Fibril Formation

Elucidating local dynamics of protein aggregation is crucial for understanding the mechanistic details of protein amyloidogenesis. Herein, we studied the residue-specific dynamics and local environmental changes of Aβ40 along the course of aggregation by using para-cyanophenylalanine (Phe_CN ) as a fluorescent and vibrational probe. Our results show that the Phe_CN residues introduced at various positions all exhibited an immediate decay of fluorescence intensity, indicating a relatively synergistic process in early oligomer formation. The fast decreases in the fluorescence intensities of residues 19 and 20 in the central hydrophobic core region and residue 10 in the N-terminal region suggest that they play crucial roles in the formation of the oligomeric core. The Phe_CN 4 residue exhibits a remarkably slower decrease in fluorescence intensity, implicating its dynamic conformational characteristics in oligomer and fibril formation. Our results also suggest that the N-terminal residues in fibrils are surrounded by a relatively hydrophobic local environment, as opposed to being solvated.

Cyanylated Cysteine Reports Site-Specific Changes at Protein-Protein-Binding Interfaces Without Perturbation

To investigate the cyanylated cysteine vibrational probe group’s ability to report on binding-induced changes along a protein–protein interface, the probe group was incorporated at several sites in a peptide of the calmodulin (CaM)-binding domain of skeletal muscle myosin light chain kinase. Isothermal titration calorimetry was used to determine the binding thermodynamics between calmodulin and each peptide. For all probe positions, the binding affinity was nearly identical to that of the unlabeled peptide. The CN stretching infrared band was collected for each peptide free in solution and bound to calmodulin. Binding-induced shifts in the IR spectral frequencies were correlated with estimated solvent accessibility based on molecular dynamics simulations. This work generally suggests (1) that site-specific incorporation of this vibrational probe group does not cause major perturbations to its local structural environment and (2) that this small probe group might be used quite broadly to map dynamic protein-binding interfaces. However, site-specific perturbations due to artificial labeling groups can be somewhat unpredictable and should be evaluated on a site-by-site basis through complementary measurements. A fully quantitative, simulation-based interpretation of the rich probe IR spectra is still needed but appears to be possible given recent advances in simulation techniques.

A Direct, Quantitative Connection between Molecular Dynamics Simulations and Vibrational Probe Line Shapes

A quantitative connection between molecular dynamics simulations and vibrational spectroscopy of probe-labeled systems would enable direct translation of experimental data into structural and dynamical information. To constitute this connection, all-atom molecular dynamics (MD) simulations were performed for two SCN probe sites (solvent-exposed and buried) in a calmodulin-target peptide complex. Two frequency calculation approaches with substantial nonelectrostatic components, a quantum mechanics/molecular mechanics (QM/MM)-based technique and a solvatochromic fragment potential (SolEFP) approach, were used to simulate the infrared probe line shapes. While QM/MM results disagreed with experiment, SolEFP results matched experimental frequencies and line shapes and revealed the physical and dynamic bases for the observed spectroscopic behavior. The main determinant of the CN probe frequency is the exchange repulsion between the probe and its local structural neighbors, and there is a clear dynamic explanation for the relatively broad probe line shape observed at the "buried" probe site. This methodology should be widely applicable to vibrational probes in many environments.

Vibrational Stark spectroscopy for assessing ligand-binding strengths in a protein

Nitrile groups are potentially useful spectroscopic probes in the infrared to characterize the binding and dynamics of ligands in proteins. This opens the possibility of locating and determining the binding mode of suitably labelled ligands in proteins based on optical spectroscopy, without the need for determining an X-ray structure. However, relating structure and spectroscopy requires means to accurately compute infrared spectra. This is investigated for benzonitrile (PhCN) in water, wild type (WT) and two lysozyme mutants in solution. The force field is validated by comparing with experimental data for benzonitrile in water which is the basis for computing the Stark shift and time scale for spectral diffusion of PhCN in WT and the L99A and L99G mutants of T4 lysozyme. The 1-d spectra for PhCN in WT and the two mutant proteins differ in their maximum absorption by up to 4 cm^-1, which reflects the modified electrostatic environments in the three proteins. It is also tested whether extending from 1-d to 2-d infrared spectroscopy provides further discrimination in the ligand-binding modes. First, for PhCN in solution the frequency fluctuation correlation function (FFCF) decays to zero at short times whereas in the protein a pronounced static inhomogeneous component is found. Secondly, the decay time of the FFCF for the mutant to which PhCN binds most strongly has the longest decay time. It is demonstrated explicitly that the ligand-binding free energy with respect to the three protein variants correlates with the Stark shift. This makes 1-d infrared spectroscopy together with computations a valuable tool for characterizing binding modes and potentially binding locations in proteins.

Ultrafast vibrational spectroscopy (2D-IR) of CO2 in ionic liquids: Carbon capture from carbon dioxide's point of view

The CO2ν3 asymmetric stretching mode is established as a vibrational chromophore for ultrafast two-dimensional infrared (2D-IR) spectroscopic studies of local structure and dynamics in ionic liquids, which are of interest for carbon capture applications. CO2 is dissolved in a series of 1-butyl-3-methylimidazolium-based ionic liquids ([C4C1im][X], where [X](-) is the anion from the series hexafluorophosphate (PF6 (-)), tetrafluoroborate (BF4 (-)), bis-(trifluoromethyl)sulfonylimide (Tf2N(-)), triflate (TfO(-)), trifluoroacetate (TFA(-)), dicyanamide (DCA(-)), and thiocyanate (SCN(-))). In the ionic liquids studied, the ν3 center frequency is sensitive to the local solvation environment and reports on the timescales for local structural relaxation. Density functional theory calculations predict charge transfer from the anion to the CO2 and from CO2 to the cation. The charge transfer drives geometrical distortion of CO2, which in turn changes the ν3 frequency. The observed structural relaxation timescales vary by up to an order of magnitude between ionic liquids. Shoulders in the 2D-IR spectra arise from anharmonic coupling of the ν2 and ν3 normal modes of CO2. Thermal fluctuations in the ν2 population stochastically modulate the ν3 frequency and generate dynamic cross-peaks. These timescales are attributed to the breakup of ion cages that create a well-defined local environment for CO2. The results suggest that the picosecond dynamics of CO2 are gated by local diffusion of anions and cations.

Site-specific infrared probes of proteins

Infrared spectroscopy has played an instrumental role in the study of a wide variety of biological questions. However, in many cases, it is impossible or difficult to rely on the intrinsic vibrational modes of biological molecules of interest, such as proteins, to reveal structural and environmental information in a site-specific manner. To overcome this limitation, investigators have dedicated many recent efforts to the development and application of various extrinsic vibrational probes that can be incorporated into biological molecules and used to site-specifically interrogate their structural or environmental properties. In this review, we highlight recent advancements in this rapidly growing research area.

Site-specific dynamics of amyloid formation and fibrillar configuration of Aβ1–23 using an unnatural amino acid

Distinct local dynamics of Aβ_1–23 amyloid formation are characterized using an unnatural amino acid p-cyanophenylalanine as a spectroscopic probe.

Probing the electrostatics of active site microenvironments along the catalytic cycle for Escherichia coli dihydrofolate reductase

Electrostatic interactions play an important role in enzyme catalysis by guiding ligand binding and facilitating chemical reactions. These electrostatic interactions are modulated by conformational changes occurring over the catalytic cycle. Herein, the changes in active site electrostatic microenvironments are examined for all enzyme complexes along the catalytic cycle of Escherichia coli dihydrofolate reductase (ecDHFR) by incorporation of thiocyanate probes at two site-specific locations in the active site. The electrostatics and degree of hydration of the microenvironments surrounding the probes are investigated with spectroscopic techniques and mixed quantum mechanical/molecular mechanical (QM/MM) calculations. Changes in the electrostatic microenvironments along the catalytic environment lead to different nitrile (CN) vibrational stretching frequencies and ¹³C NMR chemical shifts. These environmental changes arise from protein conformational rearrangements during catalysis. The QM/MM calculations reproduce the experimentally measured vibrational frequency shifts of the thiocyanate probes across the catalyzed hydride transfer step, which spans the closed and occluded conformations of the enzyme. Analysis of the molecular dynamics trajectories provides insight into the conformational changes occurring between these two states and the resulting changes in classical electrostatics and specific hydrogen-bonding interactions. The electric fields along the CN axes of the probes are decomposed into contributions from specific residues, ligands, and solvent molecules that make up the microenvironments around the probes. Moreover, calculation of the electric field along the hydride donor–acceptor axis, along with decomposition of this field into specific contributions, indicates that the cofactor and substrate, as well as the enzyme, impose a substantial electric field that facilitates hydride transfer. Overall, experimental and theoretical data provide evidence for significant electrostatic changes in the active site microenvironments due to conformational motion occurring over the catalytic cycle of ecDHFR.

Vibrational solvatochromism: towards systematic approach to modeling solvation phenomena

Vibrational solvatochromic frequency shift of IR probe is an effect of interaction between local electric field and IR probe in condensed phases. Despite prolonged efforts to develop empirical maps for vibrational frequency shifts and transition dipoles of IR probes, a systematic approach to ab initio calculation of vibrational solvatochromic charges and multipoles has not been developed. Here, we report on density functional theory (DFT) calculations of N-methylacetamide (NMA) frequency shifts using implicit and coarse-grained models. The solvatochromic infrared spectral shifts are estimated based on the distributed multipole analysis of electronic densities calculated for gas-phase equilibrium structure of NMA. Thus obtained distributed solvatochromic multipole parameters are used to calculate the amide I vibrational frequency shifts of NMA in water clusters that mimic the instantaneous configurations of the liquid water. Our results indicate that the spectral shifts are primarily electrostatic in nature and can be quantitatively reproduced using the proposed model with semi-quantitative accuracy when compared to the corresponding DFT results.

Predicting mutation-induced Stark shifts in the active site of a protein with a polarized force field

The electric field inside a protein has a significant effect on the protein structure, function, and dynamics. Recent experimental developments have offered a direct approach to measure the electric field by utilizing a nitrile-containing inhibitor as a probe that can deliver a unique vibration to the specific site of interest in the protein. The observed frequency shift of the nitrile stretching vibration exhibits a linear dependence on the electric field at the nitrile site, thus providing a direct measurement of the relative electric field. In the present work, molecular dynamics simulations were carried out to compute the electric field shift in human aldose reductase (hALR2) using a polarized protein-specific charge (PPC) model derived from fragment-based quantum-chemistry calculations in implicit solvent. Calculated changes of electric field in the active site of hALR2 between the wild type and mutants were directly compared with measured vibrational frequency shifts (Stark shifts). Our study demonstrates that the Stark shifts calculated using the PPC model are in much better agreement with the experimental data than widely used nonpolarizable force fields, indicating that the electronic polarization effect is important for the accurate prediction of changes in the electric field inside proteins.

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.

Calculation of vibrational shifts of nitrile probes in the active site of ketosteroid isomerase upon ligand binding

The vibrational Stark effect provides insight into the roles of hydrogen bonding, electrostatics, and conformational motions in enzyme catalysis. In a recent application of this approach to the enzyme ketosteroid isomerase (KSI), thiocyanate probes were introduced in site-specific positions throughout the active site. This paper implements a quantum mechanical/molecular mechanical (QM/MM) approach for calculating the vibrational shifts of nitrile (CN) probes in proteins. This methodology is shown to reproduce the experimentally measured vibrational shifts upon binding of the intermediate analogue equilinen to KSI for two different nitrile probe positions. Analysis of the molecular dynamics simulations provides atomistic insight into the roles that key residues play in determining the electrostatic environment and hydrogen-bonding interactions experienced by the nitrile probe. For the M116C-CN probe, equilinen binding reorients an active-site water molecule that is directly hydrogen-bonded to the nitrile probe, resulting in a more linear C≡N--H angle and increasing the CN frequency upon binding. For the F86C-CN probe, equilinen binding orients the Asp103 residue, decreasing the hydrogen-bonding distance between the Asp103 backbone and the nitrile probe and slightly increasing the CN frequency. This QM/MM methodology is applicable to a wide range of biological systems and has the potential to assist in the elucidation of the fundamental principles underlying enzyme catalysis.

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.