Origins of the Sensitivity of Molecular Vibrations to Electric Fields: Carbonyl and Nitrosyl Stretches in Model Compounds and Proteins

Perturbation of Fermi Resonance on Hydrogen-Bonded > C═O: 2D IR Studies of Small Ester Probes

We utilized linear and 2D infrared spectroscopy to analyze the carbonyl stretching modes of small esters in different solvents. Particularly noteworthy were the distinct carbonyl spectral line shapes in aqueous solutions, prompting our investigation of the underlying factors responsible for these differences. Through our experimental and theoretical calculations, we identified the presence of the hydrogen-bond-induced Fermi resonance as the primary contributor to the varied line shapes of small esters in aqueous solutions. Furthermore, our findings revealed that the skeletal deformation mode plays a crucial role in the Fermi resonance for all small esters. Specifically, the first overtone band of the skeletal deformation mode intensifies when hydrogen bonds form with the carbonyl group of esters, whereas such coupling is rare in aprotic organic solvents. These spectral insights carry significant implications for the utilization of esters as infrared probes in both biological and chemical systems.

Nitrile Infrared Intensities Characterize Electric Fields and Hydrogen Bonding in Protic, Aprotic, and Protein Environments

Nitriles are widely used vibrational probes; however, the interpretation of their IR frequencies is complicated by hydrogen bonding (H-bonding) in protic environments. We report a new vibrational Stark effect (VSE) that correlates the electric field projected on the -C≡N bond to the transition dipole moment and, by extension, the nitrile peak area or integrated intensity. This linear VSE applies to both H-bonding and non-H-bonding interactions. It can therefore be generally applied to determine electric fields in all environments. Additionally, it allows for semiempirical extraction of the H-bonding contribution to the blueshift of the nitrile frequency. Nitriles were incorporated at H-bonding and non-H-bonding protein sites using amber suppression, and each nitrile variant was structurally characterized at high resolution. We exploited the combined information available from variations in frequency and integrated intensity and demonstrate that nitriles are a generally useful probe for electric fields.

TfNN15N: A γ-15N-Labeled Diazo-Transfer Reagent for the Synthesis of β-15N-Labeled Azides

Azides are infrared (IR) probes that are important for structure and dynamics studies of proteins. However, they often display complex IR spectra owing to Fermi resonances and multiple conformers. Isotopic substitution of azides weakens the Fermi resonance, allowing more accurate IR spectral analysis. Site-specifically ¹⁵N-labeled aromatic azides, but not aliphatic azides, are synthesized through nitrosation. Both ¹⁵N-labeled aromatic and aliphatic azides are synthesized through nucleophilic substitution or diazo-transfer reaction but as an isotopomeric mixture. We present the synthesis of TfNN¹⁵N, a γ-¹⁵N-labeled diazo-transfer reagent, and its use to prepare β-¹⁵N-labeled aliphatic as well as aromatic azides.

Geometrical, electrical, and energetic parameters of hetero-disubstituted cumulenes and polyynes in the presence and absence of the external electric field

Cumulenes and polyynes have the potential to be applied as linear, sp-hybridized, one-dimensional all-carbon nanowires in molecular electronics and optoelectronics. The delocalization and conductivity descriptors of the two π-conjugated systems, heterodisubstituted with the NO2, CN, NH2, and OH groups, were studied using the B3LYP, B3LYP/D3, CAM-B3LYP, and ωB97XD DFT functionals, combined with the aug-cc-pVTZ basis set. Three independent types of molecular descriptors, based on geometry (the HOMA index), electrical properties (trace of the polarizability tensor), and energetic (the HOMO-LUMO energy gap) were shown to be mutually correlated and provided concordant indication that communication through the cumulene chain was considerably better than through the polyyne one. The communication can be tuned by using substituents of significantly different π-electron donor-acceptor properties as well as by the external electric field directed along the carbon chain.

Recent Progress in the Manipulation of Molecules with DC Electric Fields

The structure and reactivity of a molecule in the condensed phase are governed by its intermolecular interactions with the surrounding environment. The multipole expansion of each molecule in the condensed phase indicates that the intermolecular interactions are essentially electrostatic (e.g., ion-dipole, dipole-dipole, dipole-quadrupole, dipole-induced dipole). The electrostatic field is a fundamental language of intermolecular communications. Therefore, understanding the influence of the electrostatic field on a molecule, that is, the mechanisms by which an electrostatic field manipulates a molecule, from the perspective of molecular structure, energy states, and dynamics is indispensable for illustrating and, by extension, controlling the chemistry in molecular systems.In this Account, we describe the recent progress made in manipulation of molecular processes using an external DC electrostatic field. An electrostatic field with unprecedentedly high strength (≤4 × 10⁸ V/m) was applied in a controlled manner across a molecular film sample using the ice film nanocapacitor method. This field strength is comparable in magnitude to that of weak intermolecular interactions such as van der Waals interactions in the condensed phases. The samples were prepared using a thin film growing technique in vacuum to obtain the desired chemically tailored molecular systems. The examples of prepared systems included small molecules and molecular clusters isolated in cryogenic Ar matrices, frozen molecular films in amorphous or crystalline phase, and interfaces of multilayered molecular films. The response of the molecules to the external field was monitored by reflection-absorption infrared spectroscopy. This approach allowed us to investigate a variety of molecular systems with various intermolecular strength and environments under the influence of strong electrostatic fields. The range of observed molecular behaviors includes the manipulation of molecular orientation, intramolecular dynamics, and proton transfer reactions as an example of stereodynamic control of chemical reactivity. These observations improve our understanding of molecular behaviors in strong electric fields and broaden our perspective on electrostatic manipulation of molecules. This information is also relevant to a variety of research topics in physical and biological sciences where electric fields play a role in molecular and biological functions.

Two-dimensional IR spectroscopy reveals a hidden Fermi resonance band in the azido stretch spectrum of β-azidoalanine

Azido stretch modes in a variety of azido-derivatized nonnatural amino acids and nucleotides have been used as a site-specific infrared (IR) probe for monitoring changes in their conformations and local electrostatic environments. The vibrational bands of azide probes are often accompanied by complex line shapes with shoulder peaks, which may arise either from incomplete background subtraction, Fermi resonance, or multiple conformers. The isotope substitution in the infrared probe has thus been introduced to remove Fermi resonances without causing a significant perturbation to the structure. Here, we synthesized and labeled the mid-N atoms of aliphatic azide derivatives with 15N to study the effects of isotope labelling on their vibrational properties. The FT-IR spectra of the aliphatic azide with asymmetric lineshape became a single symmetric band upon isotope substitution, which might be an indication of the removal of the hidden Fermi resonance from the system. We also noticed that the 2D-IR spectrum of unlabeled aliphatic azide has cross-peaks, even though it is not apparently identifiable. The 1D slice spectra obtained from the 2D-IR spectra reveal the existence of a hidden Fermi resonance peak. Furthermore, we show that this weak Fermi resonance does not produce discernible oscillatory beating patterns in the IR pump-probe spectrum, which has been used as evidence of the Fermi resonance. Therefore, we confirm that isotope labelling combined with 2D-IR spectroscopy is the most efficient and incisive way to identify the origin of small shoulder peaks in the linear and nonlinear vibrational spectra of various IR probe molecules.

Caught in the Loop: Binding of the [Ru(phen)2 (dppz)]2+ Light-Switch Compound to Quadruplex DNA in Solution Informed by Time-Resolved Infrared Spectroscopy

Ultrafast time-resolved infrared (TRIR) is used to report on the binding site of the [Ru(phen)₂ (dppz)]²⁺ "light-switch" complex with both bimolecular (Oxytricha nova telomere) and intramolecular (human telomere) guanine-quadruplex structures in both K⁺ and Na⁺ containing solutions. TRIR permits the simultaneous monitoring both of the "dark" and "bright" states of the complex and of the quadruplex nucleobase bases, the latter via a Stark effect induced by the excited state of the complex. These data are used to establish the contribution of guanine base stacking and loop interactions to the binding site of this biologically relevant DNA structure in solution. A particularly striking observation is the strong thymine signal observed for the Na⁺ form of the human telomere sequence, which is expected to be in the anti-parallel conformation.

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.

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

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

Local and Global Electric Field Asymmetry in Photosynthetic Reaction Centers

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

Current and future prospects for the use of pulsed electric field in the meat industry

Pulsed electric field (PEF) is a novel non-thermal technology that has recently attracted the attention of meat scientists and technologists due to its ability to modify membrane structure and enhance mass transfer. Several studies have confirmed the potential of pulsed electric field for improving meat tenderness in both pre-rigor and post-rigor muscles during aging. However, there is a high degree of variability between studies and the underlying mechanisms are not clearly understood. While some studies have suggested physical disruption as the main cause of PEF induced tenderness, enzymatic nature of the tenderization seems to be the most plausible mechanism. Several studies have suggested the potential of PEF to mediate the tenderization process due to its membrane altering properties causing early release of calcium ions and early activation of the calpain proteases. However, experimental research is yet to confirm this postulation. Recent studies have also reported increased post-mortem proteolysis in PEF treated muscles during aging. PEF has also been reported to accelerate curing, enhance drying and reduce the numbers of both pathogens and spoilage organisms in meat, although that demands intense processing conditions. While tenderization, meat safety and accelerated curing appears to be the areas where PEF could provide attractive options in meat processing, further research is required before the application of PEF becomes a commercial reality in the meat industry. It needs to deal with carcasses which vary biochemically and in composition (muscle, fat, and bones). This review critically evaluates the published reports on the topic with the aim of reaching a clear understanding of the possible applications of PEF in the meat sector in addition to providing some insight on critical issues that need to be addressed for the technology to be a practical option for the meat industry.

Electric Fields and Enzyme Catalysis

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

Solvent-Independent Anharmonicity for Carbonyl Oscillators

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

An expanded genetic code for probing the role of electrostatics in enzyme catalysis by vibrational Stark spectroscopy

BACKGROUND

To find experimental validation for electrostatic interactions essential for catalytic reactions represents a challenge due to practical limitations in assessing electric fields within protein structures.

SCOPE OF REVIEW

This review examines the applications of non-canonical amino acids (ncAAs) as genetically encoded probes for studying the role of electrostatic interactions in enzyme catalysis.

MAJOR CONCLUSIONS

ncAAs constitute sensitive spectroscopic probes to detect local electric fields by exploiting the vibrational Stark effect (VSE) and thus have the potential to map the protein electrostatics.

GENERAL SIGNIFICANCE

Mapping the electrostatics in proteins will improve our understanding of natural catalytic processes and, in beyond, will be helpful for biocatalyst engineering. This article is part of a Special Issue entitled "Biochemistry of Synthetic Biology - Recent Developments" Guest Editor: Dr. Ilka Heinemann and Dr. Patrick O'Donoghue.

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

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

Direct observation by time-resolved infrared spectroscopy of the bright and the dark excited states of the [Ru(phen)2(dppz)]2+ light-switch compound in solution and when bound to DNA

The [Ru(phen)2(dppz)]2+ complex (1) is non-emissive in water but is highly luminescent in organic solvents or when bound to DNA, making it a useful probe for DNA binding. To date, a complete mechanistic explanation for this "light-switch" effect is still lacking. With this in mind we have undertaken an ultrafast time resolved infrared (TRIR) study of 1 and directly observe marker bands between 1280-1450 cm-1, which characterise both the emissive "bright" and the non-emissive "dark" excited states of the complex, in CD3CN and D2O respectively. These characteristic spectral features are present in the [Ru(dppz)3]2+ solvent light-switch complex but absent in [Ru(phen)3]2+, which is luminescent in both solvents. DFT calculations show that the vibrational modes responsible for these characteristic bands are predominantly localised on the dppz ligand. Moreover, they reveal that certain vibrational modes of the "dark" excited state couple with vibrational modes of two coordinating water molecules, and through these to the bulk solvent, thus providing a new insight into the mechanism of the light-switch effect. We also demonstrate that the marker bands for the "bright" state are observed for both Λ- and Δ-enantiomers of 1 when bound to DNA and that photo-excitation of the complex induces perturbation of the guanine and cytosine carbonyl bands. This perturbation is shown to be stronger for the Λ-enantiomer, demonstrating the different binding site properties of the two enantiomers and the ability of this technique to determine the identity and nature of the binding site of such intercalators.

Monitoring one-electron photo-oxidation of guanine in DNA crystals using ultrafast infrared spectroscopy

To understand the molecular origins of diseases caused by ultraviolet and visible light, and also to develop photodynamic therapy, it is important to resolve the mechanism of photoinduced DNA damage. Damage to DNA bound to a photosensitizer molecule frequently proceeds by one-electron photo-oxidation of guanine, but the precise dynamics of this process are sensitive to the location and the orientation of the photosensitizer, which are very difficult to define in solution. To overcome this, ultrafast time-resolved infrared (TRIR) spectroscopy was performed on photoexcited ruthenium polypyridyl-DNA crystals, the atomic structure of which was determined by X-ray crystallography. By combining the X-ray and TRIR data we are able to define both the geometry of the reaction site and the rates of individual steps in a reversible photoinduced electron-transfer process. This allows us to propose an individual guanine as the reaction site and, intriguingly, reveals that the dynamics in the crystal state are quite similar to those observed in the solvent medium.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

Measuring electric fields and noncovalent interactions using the vibrational stark effect

Over the past decade, we have developed a spectroscopic approach to measure electric fields inside matter with high spatial (<1 Å) and field (<1 MV/cm) resolution. The approach hinges on exploiting a physical phenomenon known as the vibrational Stark effect (VSE), which ultimately provides a direct mapping between observed vibrational frequencies and electric fields. Therefore, the frequency of a vibrational probe encodes information about the local electric field in the vicinity around the probe. The VSE method has enabled us to understand in great detail the underlying physical nature of several important biomolecular phenomena, such as drug-receptor selectivity in tyrosine kinases, catalysis by the enzyme ketosteroid isomerase, and unidirectional electron transfer in the photosynthetic reaction center. Beyond these specific examples, the VSE has provided a conceptual foundation for how to model intermolecular (noncovalent) interactions in a quantitative, consistent, and general manner. The starting point for research in this area is to choose (or design) a vibrational probe to interrogate the particular system of interest. Vibrational probes are sometimes intrinsic to the system in question, but we have also devised ways to build them into the system (extrinsic probes), often with minimal perturbation. With modern instruments, vibrational frequencies can increasingly be recorded with very high spatial, temporal, and frequency resolution, affording electric field maps correspondingly resolved in space, time, and field magnitude. In this Account, we set out to explain the VSE in broad strokes to make its relevance accessible to chemists of all specialties. Our intention is not to provide an encyclopedic review of published work but rather to motivate the underlying framework of the methodology and to describe how we make and interpret the measurements. Using certain vibrational probes, benchmarked against computer models, it is possible to use the VSE to measure absolute electric fields in arbitrary environments. The VSE approach provides an organizing framework for thinking generally about intermolecular interactions in a quantitative way and may serve as a useful conceptual tool for molecular design.

Spontaneous electric fields in solid carbon monoxide

Reflection–absorption infrared spectroscopy (RAIRS) is shown to provide a means of observing the spontelectric phase of matter, the defining characteristic of which is the occurrence of a spontaneous and powerful static electric field within a film of material.

Spontaneously electrical solids in a new light

Reflection–absorption infrared spectroscopy (RAIRS) of nitrous oxide (N₂O) thin films is shown to provide an independent means of observing the spontelectric state; the first new structural phase of matter, with unique electrical properties, to have emerged in decades.

Fourier transform infrared spectroscopy study of ligand photodissociation and migration in inducible nitric oxide synthase

Inducible nitric oxide synthase (iNOS) is a homodimeric heme enzyme that catalyzes the formation of nitric oxide (NO) from dioxygen and L-arginine (L-Arg) in a two-step process. The produced NO can either diffuse out of the heme pocket into the surroundings or it can rebind to the heme iron and inhibit enzyme action. Here we have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy at cryogenic temperatures, using the carbon monoxide (CO) and NO stretching bands as local probes of the active site of iNOS. Characteristic changes were observed in the spectra of the heme-bound ligands upon binding of the cofactors. Unlike photolyzed CO, which becomes trapped in well-defined orientations, as indicated by sharp photoproduct bands, photoproduct bands of NO photodissociated from the ferric heme iron were not visible, indicating that NO does not reside in the protein interior in a well-defined location or orientation. This may be favorable for NO release from the enzyme during catalysis because it reduces self-inhibition. Moreover, we used temperature derivative spectroscopy (TDS) with FTIR monitoring to explore the dynamics of NO and carbon monoxide (CO) inside iNOS after photodissociation at cryogenic temperatures. Only a single kinetic photoproduct state was revealed, but no secondary docking sites as in hemoglobins. Interestingly, we observed that intense illumination of six-coordinate ferrous iNOS oxy-NO ruptures the bond between the heme iron and the proximal thiolate to yield five-coordinate ferric iNOS oxy-NO, demonstrating the strong trans effect of the heme-bound NO.

Fourier transform infrared spectroscopy study of ligand photodissociation and migration in inducible nitric oxide synthase

Inducible nitric oxide synthase (iNOS) is a homodimeric heme enzyme that catalyzes the formation of nitric oxide (NO) from dioxygen and L-arginine (L-Arg) in a two-step process. The produced NO can either diffuse out of the heme pocket into the surroundings or it can rebind to the heme iron and inhibit enzyme action. Here we have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy at cryogenic temperatures, using the carbon monoxide (CO) and NO stretching bands as local probes of the active site of iNOS. Characteristic changes were observed in the spectra of the heme-bound ligands upon binding of the cofactors. Unlike photolyzed CO, which becomes trapped in well-defined orientations, as indicated by sharp photoproduct bands, photoproduct bands of NO photodissociated from the ferric heme iron were not visible, indicating that NO does not reside in the protein interior in a well-defined location or orientation. This may be favorable for NO release from the enzyme during catalysis because it reduces self-inhibition. Moreover, we used temperature derivative spectroscopy (TDS) with FTIR monitoring to explore the dynamics of NO and carbon monoxide (CO) inside iNOS after photodissociation at cryogenic temperatures. Only a single kinetic photoproduct state was revealed, but no secondary docking sites as in hemoglobins. Interestingly, we observed that intense illumination of six-coordinate ferrous iNOS oxy-NO ruptures the bond between the heme iron and the proximal thiolate to yield five-coordinate ferric iNOS oxy-NO, demonstrating the strong trans effect of the heme-bound NO.

Raman spectroscopic signatures of echovirus 1 uncoating

In recent decades, Raman spectroscopy has entered the biological and medical fields. It enables nondestructive analysis of structural details at the molecular level and has been used to study viruses and their constituents. Here, we used Raman spectroscopy to study echovirus 1 (EV1), a small, nonenveloped human pathogen, in two different uncoating states induced by heat treatments. Raman signals of capsid proteins and RNA genome were observed from the intact virus, the uncoating intermediate, and disrupted virions. Transmission electron microscopy data revealed general structural changes between the studied particles. Compared to spectral characteristics of proteins in the intact virion, those of the proteins of the heat-treated particles indicated reduced α-helix content with respect to β-sheets and coil structures. Changes observed in tryptophan and tyrosine signals suggest an increasingly hydrophilic environment around these residues. RNA signals revealed a change in the environment of the genome and in its conformation. The ionized-carbonyl vibrations showed small changes between the intact virion and the uncoating intermediate, which points to cleavage of salt bridges in the protein structure during the uncoating process. In conclusion, our data reveal distinguishable Raman signatures of the intact, intermediate, and disrupted EV1 particles. These changes indicate structural, chemical, and solute-solvent alterations in the genome and in the capsid proteins and lay the essential groundwork for investigating the uncoating of EV1 and related viruses in real time.

IMPORTANCE

In order to combat virus infection, we need to know the details of virus uncoating. We present here the novel Raman signatures for opened and intact echovirus 1. This gives hope that the signatures may be used in the near future to evaluate the ambient conditions in endosomes leading to virus uncoating using, e.g., coherent anti-Stokes Raman spectroscopy (CARS) imaging. These studies will complement structural studies on virus uncoating. In addition, Raman/CARS imaging offers the possibility of making dynamic live measurements in vitro and in cells which are impossible to measure by, for example, cryo-electron tomography. Furthermore, as viral Raman spectra can be overwhelmed with various contaminants, our study is highly relevant in demonstrating the importance of sample preparation for Raman spectroscopy in the field of virology.

Solvatochromism and the solvation structure of benzophenone

Many complex molecular phenomena, including macromolecular association, protein folding, and chemical reactivity, are determined by the nuances of their electrostatic landscapes. The measurement of such electrostatic effects is nonetheless difficult, and is typically accomplished by exploiting a spectroscopic probe within the system of interest, such as through the vibrational Stark effect. Raman spectroscopy and solvatochromism afford an alternative to this method, circumventing the limitations of infrared spectroscopy, providing a lower detection limit, and permitting measurement in a native chemical environment. To explore this possibility, the solvatochromism of the C=O and aromatic C-H stretching modes of benzophenone are investigated using Raman spectroscopy. In conjunction with density functional theory calculations, these observations are sufficient to determine the probe electrostatic environment as well as contributions from halogen and hydrogen bonding. Further analysis using a detailed Kubo-Anderson lineshape model permits the detailed assignment of distinct hydrogen bonding configurations for water in the benzophenone solvation shell. These observations reinforce the use of benzophenone as an effective electrostatic probe for complex chemical systems.

New approach to Tolman's electronic parameter based on local vibrational modes

Tolman's electronic parameter (TEP) derived from the A1-symmetrical CO stretching frequency of nickel-phosphine-tricarbonyl complexes, R3PNi(CO)3, is brought to a new, improved level by replacing normal with local vibrational frequencies. CO normal vibrational frequencies are always flawed by mode-mode coupling especially with metal-carbon stretching modes, which leads to coupling frequencies as large as 100 cm(-1) and can become even larger when the transition metal and the number of ligands is changed. Local TEP (LTEP) values, being based on local CO stretching force constants rather than normal mode frequencies, no longer suffer from mode coupling and mass effects. For 42 nickel complexes of the type LNi(CO)3, it is shown that LTEP values provide a different ordering of ligand electronic effects as previously suggested by TEP and CEP values. The general applicability of the LTEP concept is demonstrated.

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.

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

Vibrational probes can provide a direct readout of the local electrostatic field in complex molecular environments, such as protein binding sites and enzyme active sites. This information provides an experimental method to explore the underlying physical causes of important biomolecular processes such as binding and catalysis. However, specific chemical interactions such as hydrogen bonds can have complicated effects on vibrational probes and confound simple electrostatic interpretations of their frequency shifts. We employ vibrational Stark spectroscopy along with infrared spectroscopy of carbonyl probes in different solvent environments and in ribonuclease S to understand the sensitivity of carbonyl frequencies to electrostatic fields, including those due to hydrogen bonds. Additionally, we carried out molecular dynamics simulations to calculate ensemble-averaged electric fields in solvents and in ribonuclease S and found excellent correlation between calculated fields and vibrational frequencies. These data enabled the construction of a robust field-frequency calibration curve for the C═O vibration. The present results suggest that carbonyl probes are capable of quantitatively assessing the electrostatics of hydrogen bonding, making them promising for future study of protein function.

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.

CO, NO and O2 as Vibrational Probes of Heme Protein Interactions

The gaseous XO molecules (X = C, N or O) bind to the heme prosthetic group of heme proteins, and thereby activate or inhibit key biological processes. These events depend on interactions of the surrounding protein with the FeXO adduct, interactions that can be monitored via the frequencies of the Fe-X and X-O bond stretching modes, νFeX and νXO. The frequencies can be determined by vibrational spectroscopy, especially resonance Raman spectroscopy. Backbonding, the donation of Fe dπ electrons to the XO π* orbitals, is a major bonding feature in all the FeXO adducts. Variations in backbonding produce negative νFeX/νXO correlations, which can be used to gauge electrostatic and H-bonding effects in the protein binding pocket. Backbonding correlations have been established for all the FeXO adducts, using porphyrins with electron donating and withdrawing substituents. However the adducts differ in their response to variations in the nature of the axial ligand, and to specific distal interactions. These variations provide differing vantages for evaluating the nature of protein-heme interactions. We review experimental studies that explore these variations, and DFT computational studies that illuminate the underlying physical mechanisms.

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.

Ligand binding to heme proteins: a comparison of cytochrome c variants with globins

We have studied the binding of carbon monoxide (CO) in mutants of Cyt c having its methionine at position 80 replaced by alanine, aspartate, and arginine, so that the sixth coordination is available for ligand binding. We have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy to examine interactions of the heme-bound and photolyzed CO (and also nitric oxide, NO) in the small heme pocket created by the mutations. By using FTIR temperature derivative spectroscopy (TDS) and nanosecond flash photolysis, the enthalpy barrier distributions for CO rebinding were determined. In flash photolysis experiments, the majority of ligands rebind to the heme iron on picosecond time scales so that only the high-barrier tail of the distributions is visible on the nanosecond scale. By continuous wave excitation prior to TDS characterization of the barriers, however, each Cyt c molecule is photoexcited multiple times and complete photodissociation can be achieved, which likely arises from a rotation of the CO within the heme pocket so that the oxygen faces the heme iron. Apparently, reorientation prior to rebinding constitutes an additional and significant contribution to the rebinding barrier. Our experiments reveal that the compact, rigid structure of Cyt c offers no alternative binding sites for photodissociated ligands in the protein matrix. A comparison of ligand binding in these Cyt c mutants and hemoglobins underscores the importance of internal ligand docking sites and ligand migration routes for conveying a ligand binding function to heme proteins.

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

Changes in the matrix electric field in a protein, due for example to mutations or structural fluctuations, can be correlated with changes in the vibrational transition frequencies of suitable chromophores measured by IR spectroscopy through the Stark tuning rate. To make this correlation, the Stark tuning rate must be known from experiment or theory. In this paper, density functional theory at the B3LYP/TZV level of theory is used to compute the Stark tuning rate of adducts of heme porphyrin, namely, -CO, -CN, and -NO+ compounds. The results are compared with the corresponding vibrational frequencies-field dependencies from vibrational Stark spectroscopy of heme-proteins. The zero-field computed Stark tuning rate of 1.3 cm-1/MV/cm for heme-CO is in agreement with experiment, where typically the rate is 2.4/f cm-1/MV/cm for myoglobin, where f is a local field correction between 1.1 and 1.4. Several small nitrile, carbonyl, and dinitrile molecules were studied to rationalize the findings for the heme-adducted models. Here, the higher B3LYP/6-311++G(2d,2p) level of theory could be used so the agreement with recent experimental results is even better than for heme-adducted groups.