Vibrational Stark Effects on Carbonyl, Nitrile, and Nitrosyl Compounds Including Heme Ligands, CO, CN, and NO, Studied with Density Functional Theory

Hydrogen Bond Blueshifts in Nitrile Vibrational Spectra Are Dictated by Hydrogen Bond Geometry and Dynamics

Vibrational Stark effect (VSE) spectroscopy has become one of the most important experimental approaches to determine the strength of noncovalent, electrostatic interactions in chemistry and biology and to quantify their influence on structure and reactivity. Nitriles (C≡N) have been widely used as VSE probes, but their application has been complicated by an anomalous hydrogen bond (HB) blueshift which is not encompassed within the VSE framework. We present an empirical model describing the anomalous HB blueshift in terms of H-bonding geometry, i.e., as a function of HB distance and angle with respect to the C≡N group. This model is obtained by comparing vibrational observables from density functional theory and electrostatics from the polarizable AMOEBA force field, and it provides a physical explanation for the HB blueshift in terms of underlying multipolar and Pauli repulsion contributions. Additionally, we compare predicted blueshifts with experimental results and find our model provides a useful, direct framework to analyze HB geometry for rigid HBs, such as within proteins or chemical frameworks. In contrast, nitriles in highly dynamic H-bonding environments like protic solvents are no longer a function solely of geometry; this is a consequence of motional narrowing, which we demonstrate by simulating IR spectra. Overall, when HB geometry and dynamics are accounted for, an excellent correlation is found between observed and predicted HB blueshifts. This correlation includes different types of nitriles and HB donors, suggesting that our model is general and can aid in understanding HB blueshifts wherever nitriles can be implemented.

Mechanisms of Copper Selectivity and Release by the Metallochaperone CusF: Insights from CO-Binding, Rapid-Freeze-Quench EXAFS, and Unnatural Amino Acid Substitution

Metallochaperones are small proteins that shuttle essential metal ions such as Cu selectively to their cellular targets. CusF has unusual Cu(I) coordination, bound by two methionines, one histidine and a capping tryptophan residue, W44. Here we compare the CO binding reactivity of the wild type (WT) protein and its W44A, F, and M variants. Fourier-transform infrared (FTIR) indicates that W44A provides unhindered access for CO, while W44M is unreactive. WT is also largely unreactive to CO suggesting that the tryptophan cap is effective in shielding the Cu(I) center from exogenous adduct formation, while the Phe variant shows partial reactivity suggestive of an equilibrium between cap-on and cap-off conformers. Rates of metal transfer to the partner CusB are consistent with the π-cation cap providing both selectivity and redox protection. Unnatural amino acid substitutions of the W44 ligand with cyano-Phe and Br-Phe underpin the conclusion that the Phe ligand is a less effective capping residue. Finally, density functional theory (DFT) calculations validate the CO-binding strategy. Overall, the study suggests that CusF uses the tryptophan cap to protect against exogenous ligand (O2) attack while the mechanism of protein-protein complex formation allows the cap to swing out of the way, and thus have minimal effect on the rates of metal transfer.

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.

Nitrile Vibrational Lifetimes as Probes of Local Electric Fields

Optical measurements of electric fields have wide-ranging applications in the fields of chemistry and biology. Previously, such measurements focused on shifts in intensity or frequency. Here, we show that nitrile vibrational lifetimes can report local electric fields through ultrasensitive picosecond mid-infrared-near-infrared double-resonance fluorescence spectro-microscopy on Rhodamine 800. Using a robust convolution fitting approach, we observe that the nitrile vibrational lifetimes are strongly linearly correlated (R2 = 0.841) with solvent reaction fields. Supported by density functional theory, we rationalize this trend through a doorway model of intramolecular vibrational energy redistribution. This work provides new fundamental insights into the nature of vibrational energy flow in large polyatomic molecular systems and establishes a theoretical basis for electric field sensing with vibrational lifetimes, offering a new experimental dimension for probing intracellular electrostatics.

Modulation of Metal Carbonyl Stretching Frequencies in the Second Coordination Sphere through the Internal Stark Effect

Spectroscopic and computational examination of a homologous series of rhodium(I) pybox carbonyl complexes has revealed a correlation between the conformation of the flanking aryl-substituted oxazoline donors and the carbonyl stretching frequency. This relationship is also observed experimentally for octahedral rhodium(III) and ruthenium(II) variants and cannot be explained through the classical, Dewar-Chatt-Duncanson, interpretation of metal-carbonyl bonding. Instead, these findings are reconciled by local changes in the magnitude of the electric field that is projected along the metal-carbonyl vector: the internal Stark effect.

Vibrational Stark fields in carboxylic acid dimers

Carboxylic acids form exceptionally stable dimers and have been used to model proton and double proton transfer processes. The stabilization energies of the carboxylic acid dimers are very weakly dependent on the nature of substitution. However, the electric field experienced by the OH group of a particular carboxylic acid is dependent more on the nature of the substitution on the dimer partner. In general, the electric field was higher when the partner was substituted with an electron-donating group and lower with an electron-withdrawing substituent on the partner. The Stark tuning rate (Δ) of the O-H stretching vibrations calculated at the MP2/aug-cc-pVDZ level was found to be weakly dependent on the nature of substitution on the carboxylic acid. The average Stark tuning rate of O-H stretching vibrations of a particular carboxylic acid when paired with other acids was 5.7 cm^-1 (MV cm^-1)^-1, while the corresponding average Stark tuning rate of the partner acids due to a particular carboxylic acid was 21.9 cm^-1 (MV cm^-1)^-1. The difference in the Stark tuning rate is attributed to the primary and secondary effects of substitution on the carboxylic acid. The average Stark tuning rate for the anharmonic O-D frequency shifts is about 40-50% higher than the corresponding harmonic O-D frequency shifts calculated at the B3LYP/aug-cc-pVDZ level, much greater than the typical scaling factors used, indicating the strong anharmonicity of O-H/O-D oscillators in carboxylic acid dimers. Finally, the linear correlation observed between pK_a and the electric field was used to estimate the pK_a of fluoroformic acid to be around 0.9.

Quantification of the electric field inside protein active sites and fullerenes

While electrostatic interactions are exceedingly accountable for biological functions, no simple method exists to directly estimate or measure the electrostatic field in protein active sites. The electrostatic field inside the protein is generally inferred from the shift in the vibrational stretching frequencies of nitrile and thionitrile probes at the active sites through several painstaking and time-consuming experiments like vibrational Stark effect spectroscopy (VSS). Here we present a simple, fast, and reliable methodology, which can efficiently predict the vibrational Stark tuning rates (VSRs) of a large variety of probes within 10% error of the reported experimental data. Our methodology is based on geometry optimization and frequency calculations in the presence of an external electric field to predict the accurate VSR of newly designed nitrile/thionitrile probes. A priori information of VSRs is useful for difficult experiments such as catalytic/enzymatic study and in structural biology. We also applied our methodology successfully to estimate the electric field inside fullerenes and nano-onions, which is encouraging for researchers to adopt it for further applications in materials science and supramolecular chemistry.

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.

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

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

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.

Sensitivity of core-level spectroscopy to electrostatic environments of nitrile groups: An ab initio study

Ab initio quantum chemistry calculations have been performed to probe the influence of hydrogen bonding on the electronic structure of hydrogen cyanide (HCN). Our calculations determine the origin of nitrogen-specific Raman spectral features from resonant inelastic X-ray scattering occurring in the presence of a water molecule and an electric dipole field. The similarity of the two interactions in altering the electronic structure of the nitrogen atom differs only in the covalent contributions from the water molecule. The CN stretching mode as a structural probe was also investigated to study the electronic origin of the anomalous frequency shift of the nitrile group when subjected to hydrogen bonding and an electrostatic dipole field. The major changes in the electronic structure of HCN are electrostatic in nature and originate from dipole-dipole interactions. The relative shifts of the CN stretching frequency are in good agreement with those experimentally observed.

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.

Structural dynamics of nitrosylruthenium isomeric complexes studied with steady-state and transient pump-probe infrared spectroscopies

The characteristic nitrosyl stretching (NO) in the region of 1800-1900cm(-1) was used to study the geometric and ligand effect on two nitrosylruthenium complexes, namely [Ru(OAc)(2QN)2NO] (QN=2-chloro-8-quinolinol (H2cqn) or QN=2-methyl-8-quinolinol (H2mqn)). The NO stretching frequency (νNO) was found in the following order: νcis-1 (2cqn)>νcis-2 (2cqn)>νcis-1 (2mqn)>νtrans (2mqn). The results exhibited a spectral sensitivity of the NO mode to both charge distribution and ligand arrangement, which was supported by ab initio computations and natural bond orbital (NBO) analyses. Further, the vibrational population of the vibrationally excited NO stretching mode was found to relax on the order of 7-10ps, showing less than 30% variation from one isomer to another, which were explained on the basis of NO local structures and solute-solvent interactions in these isomeric nitrosylruthenium complexes.

Bis(cyano) Iron(III) Porphyrinates: What Is the Ground State?

The synthesis of six new bis(cyano) iron(III) porphyrinate derivatives is reported. The anionic porphyrin complexes utilized tetraphenylporphyrin, tetramesitylporphyrin, and tetratolylporphyrin as the porphyrin ligand. The potassium salts of Kryptofix-222 and 18-C-6 were used as the cations. These complexes have been characterized by X-ray structure analysis, solid-state Mössbauer spectroscopy, and EPR spectroscopy, both in frozen CH2Cl2 solution and in the microcrystalline state. These data show that these anionic complexes can exist in either the (dxz,dyz)(4)(dxy)(1) or the (dxy)(2)(dxz,dyz)(3) electronic configuration and some can clearly readily interconvert. This is a reflection that these two states can be very close in energy. In addition to the effects of varying the porphyrin ligand, subtle effects of the cyanide ligand environment in the solid state and in solution are sufficient to shift the balance between the two electronic states.

Calculations of the electric fields in liquid solutions

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

Computational infrared and two-dimensional infrared photon echo spectroscopy of both wild-type and double mutant myoglobin-CO proteins

The CO stretching mode of both wild-type and double mutant ( T67R / S92D ) MbCO (carbonmonoxymyoglobin) proteins is an ideal infrared (IR) probe for studying the local electrostatic environment inside the myoglobin heme pocket. Recently, to elucidate the conformational switching dynamics between two distinguishable states, extensive IR absorption, IR pump-probe, and two-dimensional (2D) IR spectroscopic studies for various mutant MbCO's have been performed by the Fayer group. They showed that the 2D IR spectroscopy of the double mutant, which has a peroxidase enzyme activity, reveals a rapid chemical exchange between two distinct states, whereas that of the wild-type does not. Despite the fact that a few simulation studies on these systems were already performed and reported, such complicated experimental results have not been fully reproduced nor described in terms of conformational state-to-state transition processes. Here, we first develop a distributed vibrational solvatochromic charge model for describing the CO stretch frequency shift reflecting local electric potential changes. Then, by carrying out molecular dynamic simulations of the two MbCO's and examining their CO frequency trajectories, it becomes possible to identify a proper reaction coordinate consisting of His64 imidazole ring rotation and its distance to the CO ligand. From the 2D surfaces of the resulting potential of mean forces, the spectroscopically distinguished A1 and A3 states of the wild-type as well as two more substates of the double mutant are identified and their vibrational frequencies and distributions are separately examined. Our simulated IR absorption and 2D IR spectra of the two MbCO's are directly compared with the previous experimental results reported by the Fayer group. The chemical exchange rate constants extracted from the two-state kinetic analyses of the simulated 2D IR spectra are in excellent agreement with the experimental values. On the basis of the quantitative agreement between the simulated spectra and experimental ones, we further examine the conformational differences in the heme pockets of the two proteins and show that the double mutation, T67R / S92D , suppresses the A1 population, restricts the imidazole ring rotation, and increases hydrogen-bond strength between the imidazole Nε-H and the oxygen atom of the CO ligand. It is believed that such delicate change of distal His64 imidazole ring dynamics induced by the double mutation may be responsible for its enhanced peroxidase catalytic activity as compared to the wild-type myoglobin.

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

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

Dynamics and the mesomorphic properties of a novel antiferroelectric liquid crystalline thiobenzoate MHPSBO10: thermal, optical and dielectric spectroscopy study

The complementary studies of the mesomorphic properties of a novel antiferroelectric liquid crystal (AFLC) (S)-2-octile 4-S-(4'decyloxybiphenyl-4-tiocarboxy)benzoate, known under MHPSBO10 acronym have been undertaken. The polymorphism has been complementary studied in details by Differential Scanning Calorimetry (DSC), Transmitted Light Intensity (TLI) and Polarization Microscopy (POM). The switching characteristics along with multiple macroscopic parameters describing the mesomorphic properties were determined by using electro-optic measurements, both upon cooling and heating. Frequency domain dielectric spectroscopy (DS), covering a wide frequency range, has been applied to characterize the molecular motions. Several collective modes, including the low frequency processes in the condensed hexatic phase were detected, analyzed in details and followed with the temperature. The presented studies deliver a wide report of the phase transitions, molecular dynamics and the macroscopic properties of the novel antiferroelectric thiobenzoate.

Experimental quantification of electrostatics in X-H···π hydrogen bonds

Hydrogen bonds are ubiquitous in chemistry and biology. The physical forces that govern hydrogen-bonding interactions have been heavily debated, with much of the discussion focused on the relative contributions of electrostatic vs quantum mechanical effects. In principle, the vibrational Stark effect, the response of a vibrational mode to electric field, can provide an experimental method for parsing such interactions into their electrostatic and nonelectrostatic components. In a previous study we showed that, in the case of relatively weak O-H···π hydrogen bonds, the O-H bond displays a linear response to an electric field, and we exploited this response to demonstrate that the interactions are dominated by electrostatics (Saggu, M.; Levinson, N. M.; Boxer, S. G. J. Am. Chem. Soc.2011, 133, 17414-17419). Here we extend this work to other X-H···π interactions. We find that the response of the X-H vibrational probe to electric field appears to become increasingly nonlinear in the order O-H < N-H < S-H. The observed effects are consistent with differences in atomic polarizabilities of the X-H groups. Nonetheless, we find that the X-H stretching vibrations of the model compounds indole and thiophenol report quantitatively on the electric fields they experience when complexed with aromatic hydrogen-bond acceptors. These measurements can be used to estimate the electrostatic binding energies of the interactions, which are found to agree closely with the results of energy calculations. Taken together, these results highlight that with careful calibration vibrational probes can provide direct measurements of the electrostatic components of hydrogen bonds.

A solvatochromic model calibrates nitriles' vibrational frequencies to electrostatic fields

Electrostatic interactions provide a primary connection between a protein's three-dimensional structure and its function. Infrared probes are useful because vibrational frequencies of certain chemical groups, such as nitriles, are linearly sensitive to local electrostatic field and can serve as a molecular electric field meter. IR spectroscopy has been used to study electrostatic changes or fluctuations in proteins, but measured peak frequencies have not been previously mapped to total electric fields, because of the absence of a field-frequency calibration and the complication of local chemical effects such as H-bonds. We report a solvatochromic model that provides a means to assess the H-bonding status of aromatic nitrile vibrational probes and calibrates their vibrational frequencies to electrostatic field. The analysis involves correlations between the nitrile's IR frequency and its (13)C chemical shift, whose observation is facilitated by a robust method for introducing isotopes into aromatic nitriles. The method is tested on the model protein ribonuclease S (RNase S) containing a labeled p-CN-Phe near the active site. Comparison of the measurements in RNase S against solvatochromic data gives an estimate of the average total electrostatic field at this location. The value determined agrees quantitatively with molecular dynamics simulations, suggesting broader potential for the use of IR probes in the study of protein electrostatics.

Quantitative, directional measurement of electric field heterogeneity in the active site of ketosteroid isomerase

Understanding the electrostatic forces and features within highly heterogeneous, anisotropic, and chemically complex enzyme active sites and their connection to biological catalysis remains a longstanding challenge, in part due to the paucity of incisive experimental probes of electrostatic properties within proteins. To quantitatively assess the landscape of electrostatic fields at discrete locations and orientations within an enzyme active site, we have incorporated site-specific thiocyanate vibrational probes into multiple positions within bacterial ketosteroid isomerase. A battery of X-ray crystallographic, vibrational Stark spectroscopy, and NMR studies revealed electrostatic field heterogeneity of 8 MV/cm between active site probe locations and widely differing sensitivities of discrete probes to common electrostatic perturbations from mutation, ligand binding, and pH changes. Electrostatic calculations based on active site ionization states assigned by literature precedent and computational pK(a) prediction were unable to quantitatively account for the observed vibrational band shifts. However, electrostatic models of the D40N mutant gave qualitative agreement with the observed vibrational effects when an unusual ionization of an active site tyrosine with a pK(a) near 7 was included. UV-absorbance and (13)C NMR experiments confirmed the presence of a tyrosinate in the active site, in agreement with electrostatic models. This work provides the most direct measure of the heterogeneous and anisotropic nature of the electrostatic environment within an enzyme active site, and these measurements provide incisive benchmarks for further developing accurate computational models and a foundation for future tests of electrostatics in enzymatic catalysis.

Direct measurements of electric fields in weak OH···π hydrogen bonds

Hydrogen bonds and aromatic interactions are of widespread importance in chemistry, biology, and materials science. Electrostatics play a fundamental role in these interactions, but the magnitude of the electric fields that support them has not been quantified experimentally. Phenol forms a weak hydrogen bond complex with the π-cloud of benzene, and we used this as a model system to study the role of electric fields in weak OH···π hydrogen bonds. The effects of complex formation on the vibrational frequency of the phenol OH or OD stretches were measured in a series of benzene-based aromatic solvents. Large shifts are observed and these can be converted into electric fields via the measured vibrational Stark effect. A comparison of the measured fields with quantum chemical calculations demonstrates that calculations performed in the gas phase are surprisingly effective at capturing the electrostatics observed in solution. The results provide quantitative measurements of the magnitude of electric fields and electrostatic binding energies in these interactions and suggest that electrostatics dominate them. The combination of vibrational Stark effect (VSE) measurements of electric fields and high-level quantum chemistry calculations is a general strategy for quantifying and characterizing the origins of intermolecular interactions.

Enzyme reactivation by hydrogen peroxide in heme-based tryptophan dioxygenase

An intriguing mystery about tryptophan 2,3-dioxygenase is its hydrogen peroxide-triggered enzyme reactivation from the resting ferric oxidation state to the catalytically active ferrous form. In this study, we found that such an odd Fe(III) reduction by an oxidant depends on the presence of L-Trp, which ultimately serves as the reductant for the enzyme. In the peroxide reaction with tryptophan 2,3-dioxygenase, a previously unknown catalase-like activity was detected. A ferryl species (δ = 0.055 mm/s and ΔE(Q) = 1.755 mm/s) and a protein-based free radical (g = 2.0028 and 1.72 millitesla linewidth) were characterized by Mössbauer and EPR spectroscopy, respectively. This is the first compound ES-type of ferryl intermediate from a heme-based dioxygenase characterized by EPR and Mössbauer spectroscopy. Density functional theory calculations revealed the contribution of secondary ligand sphere to the spectroscopic properties of the ferryl species. In the presence of L-Trp, the reactivation was demonstrated by enzyme assays and by various spectroscopic techniques. A Trp-Trp dimer and a monooxygenated L-Trp were both observed as the enzyme reactivation by-products by mass spectrometry. Together, these results lead to the unraveling of an over 60-year old mystery of peroxide reactivation mechanism. These results may shed light on how a metalloenzyme maintains its catalytic activity in an oxidizing environment.

Molecular-level mechanisms of vibrational frequency shifts in a polar liquid

A molecular-level analysis of the origins of the vibrational frequency shifts of the CN stretching mode in neat liquid acetonitrile is presented. The frequency shifts and infrared spectrum are calculated using a perturbation theory approach within a molecular dynamics simulation and are in good agreement with measured values reported in the literature. The resulting instantaneous frequency of each nitrile group is decomposed into the contributions from each molecule in the liquid and by interaction type. This provides a detailed picture of the mechanisms of frequency shifts, including the number of surrounding molecules that contribute to the shift, the relationship between their position and relative contribution, and the roles of electrostatic and van der Waals interactions. These results provide insight into what information is contained in infrared (IR) and Raman spectra about the environment of the probed vibrational mode.

Ligand migration and binding in nonsymbiotic hemoglobins of Arabidopsis thaliana

We have studied carbon monoxide (CO) migration and binding in the nonsymbiotic hemoglobins AHb1 and AHb2 of Arabidopsis thaliana using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures. Both proteins have similar amino acid sequences but display pronounced differences in ligand binding properties, at both physiological and cryogenic temperatures. Near neutral pH, the distal HisE7 side chain is close to the heme-bound ligand in the majority of AHb1-CO molecules, as indicated by a low CO stretching frequency at 1921 cm(-1). In this fraction, two CO docking sites can be populated, the primary site B and the secondary site C. When the pH is lowered, a high-frequency stretching band at approximately 1964 cm(-1) grows at the expense of the low-frequency band, indicating that HisE7 protonates and, concomitantly, moves away from the bound ligand. Geminate rebinding barriers are markedly different for the two conformations, and docking site C is not accessible in the low-pH conformation. Rebinding of NO ligands was observed only from site B of AHb1, regardless of conformation. In AHb2, the HisE7 side chain is removed from the bound ligand; rebinding barriers are low, and CO molecules can populate only primary docking site B. These results are interpreted in terms of differences in the active site structures and physiological functions.

Residue-specific, real-time characterization of lag-phase species and fibril growth during amyloid formation: a combined fluorescence and IR study of p-cyanophenylalanine analogs of islet amyloid polypeptide

Amyloid formation normally exhibits a lag phase followed by a growth phase, which leads to amyloid fibrils. Characterization of the species populated during the lag phase is experimentally challenging, but is critical since the most toxic entities may be pre-fibrillar species. p-Cyanophenylalanine (F(C[triple bond]N)) fluorescence is used to probe the nature of lag-phase species populated during the formation of amyloid by human islet amyloid polypeptide. The polypeptide contains two phenylalanines at positions 15 and 23 and a single tyrosine located at the C-terminus. Each aromatic residue was separately replaced by F(C[triple bond]N). The substitutions do not perturb amyloid formation relative to wild-type islet amyloid polypeptide as detected using thioflavin T fluorescence and electron microscopy. F(C[triple bond]N) fluorescence is high when the cyano group is hydrogen bonded and low when it is not. It can also be quenched via Förster resonance energy transfer to tyrosine. Fluorescence intensity was monitored in real time and revealed that all three positions remained exposed to solvent during the lag phase but less exposed than unstructured model peptides. The time course of amyloid formation as monitored by thioflavin T fluorescence and F(C[triple bond]N) fluorescence is virtually identical. Fluorescence quenching experiments confirmed that each residue remains exposed during the lag phase. These results place significant constraints on the nature of intermediates that are populated during the lag phase and indicate that significant sequestering of the aromatic side chains does not occur until beta-structure sufficient to bind thioflavin T has developed. Seeding studies and analysis of maximum rates confirm that sequestering of the cyano groups occurs concomitantly with the development of thioflavin T binding capability. Overall, the process of amyloid formation and growth appears to be remarkably homogenous in terms of side-chain ordering. F(C[triple bond]N) also provides information about fibril structure. Fluorescence emission measurements, infrared measurements, and quenching studies indicate that the aromatic residues are differentially exposed in the fibril state with Phe15 being the most exposed. F(C[triple bond]N) is readily accommodated into proteins; thus, the approach should be broadly applicable.

5-Cyanotryptophan as an Infrared Probe of Local Hydration Status of Proteins

The nitrile (C≡N) stretching vibration is sensitive to environment, making nitrile-derivatized amino acids an increasingly utilized tool to study various biological processes. Herein, we show that the bandwidth of the C≡N stretching vibration of 5-cyanotryptophan is particularly sensitive to water, rendering it an attractive infrared probe of local hydration status. We confirm the utility of this probe in biological applications by using it to examine how the hydration status of individual tryptophan sidechains of an antimicrobial peptide, indolicidin, changes upon peptide binding to model membranes. Furthermore, we show that p-cyanophenylalanine and 5-cyanotryptophan constitute a useful fluorescence energy transfer pair.

Probing the role of hydration in the unfolding transitions of carbonmonoxy myoglobin and apomyoglobin

We show that the equilibrium unfolding transition of horse carbonmonoxy myoglobin monitored by the stretching vibration of the CO ligand, a local environmental probe, is very sharp and, thus, quite different from those measured by global conformational reporters. In addition, the denatured protein exhibits an A(0)-like CO band. We hypothesize that this sharp transition reports penetration of water into the heme pocket of the protein. Parallel experiments on horse apomyoglobin, wherein an environment-sensitive fluorescent probe, nile red, was used, also reveals a similar putative hydration event. Given the importance of dehydration in protein folding and also the recent debate over the interpretation of probe-dependent unfolding transitions, these results have strong implications on the mechanism of protein folding.

Vibrational stark effect probes for nucleic acids

The vibrational Stark effect (VSE) has proven to be an effective method for the study of electric fields in proteins via the use of infrared probes. To explore the use of VSE in nucleic acids, we investigated the Stark spectroscopy of nine structurally diverse nucleosides. These nucleosides contained nitrile or azide probes in positions that correspond to both the major and minor grooves of DNA. The nitrile probes showed better characteristics and exhibited absorption frequencies over a broad range; that is, from 2253 cm-1 for 2'-O-cyanoethyl ribonucleosides 8 and 9 to 2102 cm(-1) for a 13C-labeled 5-thiocyanatomethyl-2'-deoxyuridine 3c. The largest Stark tuning rate observed was |Deltamu| = 1.1 cm(-1)/(MV/cm) for both 5-cyano-2'-deoxyuridine 1 and N2-nitrile-2'-deoxyguanosine 7. The latter is a particularly attractive probe because of its high extinction coefficient (epsilon = 412 M-1cm-1) and ease of incorporation into oligomers.

Electric field effects on aromatic and aliphatic hydrocarbons: a density-functional study

The influence of a uniform static external electric field on some aliphatic and aromatic molecular species is studied within the density functional theory (DFT) employing the 6-311++G(2d,2p) basis set with B3LYP exchange-correlation prescription. The electric field perturbs the molecular geometry but drastically alters the dipole moments and engenders, to a varying degree, the molecular vibrational Stark effect, i.e., shifts in the infrared (IR) vibrational frequencies accompanied by spectral intensity redistribution. For polar molecules, significant negative ("red") and positive ("blue") frequency shifts are observed for field orientations both parallel and antiparallel to their permanent dipole moments. Further, a selective reordering of frontier orbitals is observed to be brought about by moderately intense fields. In particular, molecules having a lowest unoccupied molecular orbital (LUMO) with predominant pi character possess a threshold field beyond which energy gap between the highest occupied molecular orbital (HOMO) and LUMO diminishes rapidly. A time-dependent (TD) DFT analysis reveals that an increase in the applied field strength by and large increases the excitation energies corresponding to significant electronic transitions among frontier MOs with a concomitant decrease in their oscillator strengths.

Molecular Dynamics Simulations of CN− Dynamics and Spectroscopy in Myoglobin

The vibrational dynamics of the cyanide anion and the heme group in MbCN (CN complexed to Myoglobin) are investigated using molecular dynamics simulations. A previously calculated quantum-chemical heme-ligand potential-energy surface together with a three-center charge model for the iron-ligand center that captures both polarization and ligand-to-metal charge transfer allows for a detailed description of the interactions around the active site. It is found that the CN binding orientation (Fe--CN or Fe--NC) to the heme affects the stretching frequency of the ligand, with a 25 cm-1 difference in the fundamental wavenumber between the two orientations as well as a change in bond length. The charge model also captures such crucial interactions as the possible hydrogen bond between the ligand and the His64 residue. This interaction is weakened when the ligand binds in the Fe--NC conformation but is also sensitive to the protonation state of His64. The structural changes around the active site, the observation of water penetration for the Fe--NC conformation, the computed IR spectrum, and the energetics suggest that the Fe--CN conformation with Hisepsilon64 is the most likely one. The water accessibility of the active site is also found to be related to the protonation state of His64. The presence of water in the active site could also affect the IR band of the C--N stretch mode. Thus, IR spectroscopy of the C--N stretch is a potentially useful reporter of ligand isomers and active-site structure.

New insights on the electronic and molecular structure of cyanide-ligated iron(III) porphyrinates

The preparation and characterization of several new cyano-ligated six-coordinate low-spin iron(III) porphyrinates are reported. The synthesis and structure of the new bis(cyanide) derivative K(222)][Fe(TMP)(CN)2] (TMP = tetramesitylporphyrinate) is described. Three mixed-ligand species of the general form [Fe(Porph)(CN)(L)], where L = 1-methylimidazole or pyridine, have also been prepared and structurally characterized. All complexes have been studied with EPR spectroscopy in frozen solution and in microcrystalline form. In some cases, especially those of the bis(cyanide) derivative above and the previously reported [Fe(TPP)(CN)2](-), there are significant differences in the EPR spectra as a result of the state change. These spectral differences can be correlated with changes in the electron configuration that are the likely result of a differing environment of the coordinated cyanide ligands; the core conformation and electronic structure of the porphyrin ligand are unlikely to play a role. All four new complexes and [Fe(TPP)(CN)2](-) have been studied by Mössbauer spectroscopy with variable-temperature and applied magnetic-field measurements. The sign of the quadrupole splitting value has been established as negative. These measurements have allowed us to give estimates of the energy difference between the two close-lying dpi (dxz and dyz) orbitals. These splitting values range from approximately 267 cm-1 for [Fe(TPP)(CN)2](-) to approximately 614 cm(-1) for [Fe(TPP)(CN)(Py)].

Classical Molecular Electrostatics: Recognition of Ligands in Proteins and the Vibrational Stark Effect

It is shown that classical electrostatics quantitatively describes both the binding of the diatomic ligands XO (X = C, N, O) to the heme group in myoglobin and the dependence of their vibrational frequencies upon an external field, the vibrational Stark effect. The key is a proper treatment of induced dipoles. The results suggest that ligand binding occurs via an "electrostatic bond", a generalization of the standard ionic bond to include induction, and, more generally, that classical electrostatics can replace quantum mechanics for a considerable simplification of some complex problems.

Infrared absorption study of the heme pocket dynamics of carbonmonoxyheme proteins

The temperature dependencies of the infrared absorption CO bands of carboxy complexes of horseradish peroxidase (HRP(CO)) in glycerol/water mixture at pH 6.0 and 9.3 are interpreted using the theory of optical absorption bandshape. The bands' anharmonic behavior is explained assuming that there is a higher-energy set of conformational substates (CSS(h)), which are populated upon heating and correspond to the protein substates with disordered water molecules in the heme pocket. Analysis of the second moments of the CO bands of the carboxy complexes of myoglobin (Mb(CO)) and hemoglobin (Hb(CO)), and of HRP(CO) with benzohydroxamic acid (HRP(CO)+BHA), shows that the low energy CSS(h) exists also in the open conformation of Mb(CO), where the heme pocket is spacious enough to accommodate a water molecule. In the HRP(CO)+BHA and closed conformations of Mb(CO) and Hb(CO), the heme pocket is packed with BHA and different amino acids, the CSS(h) has much higher energy and is hardly populated even at the highest temperatures. Therefore only motions of these amino acids contribute to the band broadening. These motions are linked to the protein surface and frozen in the glassy matrix, whereas in the liquid solvent they are harmonic. Thus the second moment of the CO band is temperature-independent in glass and is proportional to the temperature in liquid. The temperature dependence of the second moment of the CO peak of HRP(CO) in the trehalose glass exhibits linear coupling to an oscillator. This oscillator can be a moving water molecule locked in the heme pocket in the whole interval of temperatures or a trehalose molecule located in the heme pocket.

A novel fluorescent probe for protein binding and folding studies:p-cyano-phenylalanine

Recently, it is has been shown that the C=N stretching vibration of a non-natural amino acid, p-cyano-phenylalanine (PheCN), could be used as an infrared reporter of local environment. Here, we further showed that the fluorescence emission of PheCN is also sensitive to solvent and, therefore, could be used as a novel optical probe for protein binding and folding studies. Moreover, we found that the fluorescence quantum yield of PheCN is nearly five times larger than that of phenylalanine and, more importantly, can be selectively excited even when other aromatic amino acids are present, thus making it a more versatile fluorophore. To test the feasibility of using PheCN as a practical fluorescent probe, we studied the binding of calmodulin (CaM) to a peptide derived from the CaM-binding domain of skeletal muscle myosin light chain kinase (MLCK). The peptide (MLCK3CN) contains a single PheCN residue and has been shown to bind to CaM with high affinity. As expected, addition of CaM into a MLCK3CN solution resulted in quenching of the PheCN fluorescence. A series of stochiometric titrations further allowed us to determine the binding affinity (Kd) of this peptide to CaM. Taken together, these results indicated that the PheCN fluorescence is sensitive to environment and could be applicable to a wide variety of biological problems.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

Internal dynamics and protein–matrix coupling in trehalose-coated proteins

We review recent studies on the role played by non-liquid, water-containing matrices on the dynamics and structure of embedded proteins. Two proteins were studied, in water-trehalose matrices: a water-soluble protein (carboxy derivative of horse heart myoglobin) and a membrane protein (reaction centre from Rhodobacter sphaeroides). Several experimental techniques were used: Mossbauer spectroscopy, elastic neutron scattering, FTIR spectroscopy, CO recombination after flash photolysis in carboxy-myoglobin, kinetic optical absorption spectroscopy following pulsed and continuous photoexcitation in Q(B) containing or Q(B) deprived reaction centre from R. sphaeroides. Experimental results, together with the outcome of molecular dynamics simulations, concurred to give a picture of how water-containing matrices control the internal dynamics of the embedded proteins. This occurs, in particular, via the formation of hydrogen bond networks that anchor the protein surface to the surrounding matrix, whose stiffness increases by lowering the sample water content. In the conclusion section, we also briefly speculate on how the protein-matrix interactions observed in our samples may shed light on the protein-solvent coupling also in liquid aqueous solutions.

Temperature excursion infrared (TEIR) spectroscopy used to study hydrogen bonding between water and biomolecules

Water is a highly polar molecule that is capable of making four H-bonding linkages. Stability and specificity of folding of water-soluble protein macromolecules are determined by the interplay between water and functional groups of the protein. Yet, under some conditions, water can be replaced with sugar or other polar protic molecules with retention of protein structure. Infrared (IR) spectroscopy allows one to probe groups on the protein that interact with solvent, whether the solvent is water, sugar or glycerol. The basis of the measurement is that IR spectral lines of functional groups involved in H-bonding show characteristic spectral shifts with temperature excursion, reflecting the dipolar nature of the group and its ability to H-bond. For groups involved in H-bonding to water, the stretching mode absorption bands shift to lower frequency, whereas bending mode absorption bands shift to higher frequency as temperature decreases. The results indicate increasing H-bonding and decreasing entropy occurring as a function of temperature, even at cryogenic temperatures. The frequencies of the amide group modes are temperature dependent, showing that as temperature decreases, the amide group H-bonds to water strengthen. These results are relevant to protein stability as a function of temperature. The influence of solvent relaxation is demonstrated for tryptophan fluorescence over the same temperature range where the solvent was examined by infrared spectroscopy.