Ribonuclease S dynamics measured using a nitrile label with 2D IR vibrational echo spectroscopy

Minimalist IR and fluorescence probes of protein function

Spectroscopic studies of small proteins and peptides, especially those requiring fine spatial and/or temporal resolution, demand synthetic probes that confer the minimal possible steric and functional change on the native properties. Here we review the recent progress in development of minimally disruptive probes for fluorescence and infrared spectroscopies, as well as the methods to efficiently incorporate them into proteins. Advances in spectroscopy on the one hand result in high specialization of synthetic probes for a particular purpose, but on the other hand allow for the same probes be used for different techniques to gather complementary biochemical information.

Two-Dimensional Infrared Study of Vibrational Coupling between Azide and Nitrile Reporters in a RNA Nucleoside

The vibrations in the azide, N3, asymmetric stretching region and nitrile, CN, symmetric stretching region of 2'-azido-5-cyano-2'-deoxyuridine (N3CNdU) are examined by two-dimensional infrared (2D IR) spectroscopy. At earlier waiting times, the 2D IR spectrum shows the presence of both vibrational transitions along the diagonal and off-diagonal cross peaks indicating vibrational coupling. The coupling strength is determined from the off-diagonal anharmonicity to be 66 cm(-1) for the intramolecular distance of ∼7.9 Å, based on a structural map generated for this model system. In addition, the frequency-frequency correlation decay is detected, monitoring the solvent dynamics around each individual probe position. Overall, these vibrational reporters can be utilized in tandem to simultaneously track global structural information and fast structural fluctuations.

Unraveling the dynamics and structure of functionalized self-assembled monolayers on gold using 2D IR spectroscopy and MD simulations

Functionalized self-assembled monolayers (SAMs) are the focus of ongoing investigations because they can be chemically tuned to control their structure and dynamics for a wide variety of applications, including electrochemistry, catalysis, and as models of biological interfaces. Here we combine reflection 2D infrared vibrational echo spectroscopy (R-2D IR) and molecular dynamics simulations to determine the relationship between the structures of functionalized alkanethiol SAMs on gold surfaces and their underlying molecular motions on timescales of tens to hundreds of picoseconds. We find that at higher head group density, the monolayers have more disorder in the alkyl chain packing and faster dynamics. The dynamics of alkanethiol SAMs on gold are much slower than the dynamics of alkylsiloxane SAMs on silica. Using the simulations, we assess how the different molecular motions of the alkyl chain monolayers give rise to the dynamics observed in the experiments.

Probing the effectiveness of spectroscopic reporter unnatural amino acids: a structural study

The X-ray crystal structures of superfolder green fluorescent protein (sfGFP) containing the spectroscopic reporter unnatural amino acids (UAAs) 4-cyano-L-phenylalanine (pCNF) or 4-ethynyl-L-phenylalanine (pCCF) at two unique sites in the protein have been determined. These UAAs were genetically incorporated into sfGFP in a solvent-exposed loop region and/or a partially buried site on the β-barrel of the protein. The crystal structures containing the UAAs at these two sites permit the structural implications of UAA incorporation for the native protein structure to be assessed with high resolution and permit a direct correlation between the structure and spectroscopic data to be made. The structural implications were quantified by comparing the root-mean-square deviation (r.m.s.d.) between the crystal structure of wild-type sfGFP and the protein constructs containing either pCNF or pCCF in the local environment around the UAAs and in the overall protein structure. The results suggest that the selective placement of these spectroscopic reporter UAAs permits local protein environments to be studied in a relatively nonperturbative fashion with site-specificity.

Applications of two-dimensional infrared spectroscopy

Two-dimensional infrared (2D IR) spectroscopy has recently emerged as a powerful tool with applications in many areas of scientific research. The inherent high time resolution coupled with bond-specific spatial resolution of IR spectroscopy enable direct characterization of rapidly interconverting species and fast processes, even in complex systems found in chemistry and biology. In this minireview, we briefly outline the fundamental principles and experimental procedures of 2D IR spectroscopy. Using illustrative example studies, we explain the important features of 2D IR spectra and their capability to elucidate molecular structure and dynamics. Primarily, this minireview aims to convey the scope and potential of 2D IR spectroscopy by highlighting select examples of recent applications including the use of innate or introduced vibrational probes for the study of nucleic acids, peptides/proteins, and materials.

Line shape analysis of two-dimensional infrared spectra

Ultrafast two-dimensional infrared (2D IR) spectroscopy probes femtosecond to picosecond time scale dynamics ranging from solvation to protein motions. The frequency-frequency correlation function (FFCF) is the quantitative measure of the spectral diffusion that reports those dynamics and, within certain approximations, can be extracted directly from 2D IR line shapes. A variety of methods have been developed to extract the FFCF from 2D IR spectra, which, in principle, should give the same FFCF parameters, but the complexity of real experimental systems will affect the results of these analyses differently. Here, we compare five common analysis methods using both simulated and experimental 2D IR spectra to understand the effects of apodization, anharmonicity, phasing errors, and finite signal-to-noise ratios on the results of each of these analyses. Our results show that although all of the methods can, in principle, yield the FFCF under idealized circumstances, under more realistic experimental conditions they behave quite differently, and we find that the centerline slope analysis yields the best compromise between the effects we test and is most robust to the distortions that they cause.

Fast infrared spectroscopy of protein dynamics: advancing sensitivity and selectivity

2D-IR spectroscopy has matured to a powerful technique to study the structure and dynamics of peptides, but its extension to larger proteins is still in its infancy, the major limitations being sensitivity and selectivity. Site-selective information requires measuring single vibrational probes at sub-millimolar concentrations where most proteins are still stable, which is a severe challenge for conventional (FT)IR spectroscopy. Besides its ultrafast time-resolution, a so far largely underappreciated potential of 2D-IR spectroscopy lies in its sensitivity gain. The present paper sets the goals and outlines strategies how to use that sensitivity gain together with properly designed vibrational labels to make IR spectroscopy a versatile tool to study a wide class of proteins.

Site-specific infrared probes of proteins

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

Complexation dynamics of CH3SCN and Li+ in acetonitrile studied by two-dimensional infrared spectroscopy

A chemical exchange 2DIR study of ion–molecule complexation dynamics in electrolyte solutions.

β-Isocyanoalanine as an IR probe: comparison of vibrational dynamics between isonitrile and nitrile-derivatized IR probes

An infrared (IR) probe based on isonitrile (NC)-derivatized alanine 1 was synthesized and the vibrational properties of its NC stretching mode were investigated using FTIR and femtosecond IR pump–probe spectroscopy.

Site-selective Characterization of Src Homology 3 Domain Molecular Recognition with Cyanophenylalanine Infrared Probes

Local heterogeneity of microenvironments in proteins is important in biological function, but difficult to characterize experimentally. One approach is the combination of infrared (IR) spectroscopy and site-selective incorporation of probe moieties with spectrally resolved IR absorptions that enable characterization within inherently congested protein IR spectra. We employed this method to study molecular recognition of a Src homology 3 (SH3) domain from the yeast protein Sho1 for a peptide containing the proline-rich recognition sequence of its physiological binding partner Pbs2. Nitrile IR probes were introduced at four distinct sites in the protein by selective incorporation of p-cyanophenylalanine via the amber codon suppressor method and characterized by IR spectroscopy. Variation among the IR absorption bands reports on heterogeneity in local residue environments dictated by the protein structure, as well as on residue-dependent changes upon peptide binding. The study informs on the molecular recognition of SH3 ^Sho1 and illustrates the speed and simplicity of this approach for characterization of select microenvironments within proteins.

Synthesis and Protein Incorporation of Azido-Modified Unnatural Amino Acids

Two new azidophenylalanine residues (3 and 4) have been synthesized and, in combination with 4-azido-L-phenylalanine (1) and 4-azidomethyl-L-phenylalanine (2), form a series of unnatural amino acids (UAAs) containing the azide vibrational reporter at varying distances from the aromatic ring of phenylalanine. These UAAs were designed to probe protein hydration with high spatial resolution by utilizing the large extinction coefficient and environmental sensitivity of the azide asymmetric stretch vibration. The sensitivity of the azide reporters was investigated in solvents that mimic distinct local protein environments. Three of the four azido-modified phenylalanine residues were successfully genetically incorporated into a surface site in superfolder green fluorescent protein (sfGFP) utilizing an engineered, orthogonal aminoacyl-tRNA synthetase in response to an amber codon with high efficiency and fidelity. SDS-PAGE and ESI-Q-TOF mass analysis verified the site-specific incorporation of these UAAs. The observed azide asymmetric stretch in the linear IR spectra of these UAAs incorporated into sfGFP indicated that the azide groups were hydrated in the protein.

The non-uniform early structural response of globular proteins to cold denaturing conditions: a case study with Yfh1

The mechanism of cold denaturation in proteins is often incompletely understood due to limitations in accessing the denatured states at extremely low temperatures. Using atomistic molecular dynamics simulations, we have compared early (nanosecond timescale) structural and solvation properties of yeast frataxin (Yfh1) at its temperature of maximum stability, 292 K (Ts), and the experimentally observed temperature of complete unfolding, 268 K (Tc). Within the simulated timescales, discernible "global" level structural loss at Tc is correlated with a distinct increase in surface hydration. However, the hydration and the unfolding events do not occur uniformly over the entire protein surface, but are sensitive to local structural propensity and hydrophobicity. Calculated infrared absorption spectra in the amide-I region of the whole protein show a distinct red shift at Tc in comparison to Ts. Domain specific calculations of IR spectra indicate that the red shift primarily arises from the beta strands. This is commensurate with a marked increase in solvent accessible surface area per residue for the beta-sheets at Tc. Detailed analyses of structure and dynamics of hydration water around the hydrophobic residues of the beta-sheets show a more bulk water like behavior at Tc due to preferential disruption of the hydrophobic effects around these domains. Our results indicate that in this protein, the surface exposed beta-sheet domains are more susceptible to cold denaturing conditions, in qualitative agreement with solution NMR experimental results.

Probing the phosphopantetheine arm conformations of acyl carrier proteins using vibrational spectroscopy

Acyl carrier proteins (ACPs) are universal and highly conserved domains central to both fatty acid and polyketide biosynthesis. These proteins tether reactive acyl intermediates with a swinging 4'-phosphopantetheine (Ppant) arm and interact with a suite of catalytic partners during chain transport and elongation while stabilizing the growing chain throughout the biosynthetic pathway. The flexible nature of the Ppant arm and the transient nature of ACP-enzyme interactions impose a major obstacle to obtaining structural information relevant to understanding polyketide and fatty acid biosynthesis. To overcome this challenge, we installed a thiocyanate vibrational spectroscopic probe on the terminal thiol of the ACP Ppant arm. This site-specific probe successfully reported on the local environment of the Ppant arm of two ACPs previously characterized by solution NMR, and was used to determine the solution exposure of the Ppant arm of an ACP from 6-deoxyerythronolide B synthase (DEBS). Given the sensitivity of the probe's CN stretching band to conformational distributions resolved on the picosecond time scale, this work lays a foundation for observing the dynamic action-related structural changes of ACPs using vibrational spectroscopy.

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

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

2D IR spectroscopy reveals the role of water in the binding of channel-blocking drugs to the influenza M2 channel

Water is an integral part of the homotetrameric M2 proton channel of the influenza A virus, which not only assists proton conduction but could also play an important role in stabilizing channel-blocking drugs. Herein, we employ two dimensional infrared (2D IR) spectroscopy and site-specific IR probes, i.e., the amide I bands arising from isotopically labeled Ala30 and Gly34 residues, to probe how binding of either rimantadine or 7,7-spiran amine affects the water dynamics inside the M2 channel. Our results show, at neutral pH where the channel is non-conducting, that drug binding leads to a significant increase in the mobility of the channel water. A similar trend is also observed at pH 5.0 although the difference becomes smaller. Taken together, these results indicate that the channel water facilitates drug binding by increasing its entropy. Furthermore, the 2D IR spectral signatures obtained for both probes under different conditions collectively support a binding mechanism whereby amantadine-like drugs dock in the channel with their ammonium moiety pointing toward the histidine residues and interacting with a nearby water cluster, as predicted by molecular dynamics simulations. We believe these findings have important implications for designing new anti-influenza drugs.

Microscopic insights into the protein-stabilizing effect of trimethylamine N-oxide (TMAO)

Although it is widely known that trimethylamine N-oxide (TMAO), an osmolyte used by nature, stabilizes the folded state of proteins, the underlying mechanism of action is not entirely understood. To gain further insight into this important biological phenomenon, we use the C≡N stretching vibration of an unnatural amino acid, p-cyano-phenylalanine, to directly probe how TMAO affects the hydration and conformational dynamics of a model peptide and a small protein. By assessing how the lineshape and spectral diffusion properties of this vibration change with cosolvent conditions, we are able to show that TMAO achieves its protein-stabilizing ability through the combination of (at least) two mechanisms: (i) It decreases the hydrogen bonding ability of water and hence the stability of the unfolded state, and (ii) it acts as a molecular crowder, as suggested by a recent computational study, that can increase the stability of the folded state via the excluded volume effect.

Vibrational dynamics and solvatochromism of the label SCN in various solvents and hemoglobin by time dependent IR and 2D-IR spectroscopy

The vibrational label SCN is used to report on local structural dynamics in a protein revealing spectral diffusion on a picosecond scale. The SCN spectra are compared to the response of methylthiocyanate in solvents with different polarity and hydrogen-bonding capabilities.

Two-dimensional spectroscopy of coupled vibrations with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical representation of the nonlinear vibrational response function used to compute multidimensional infrared spectra. In this method, response functions are calculated from a sequence of classical trajectories linked by discontinuities representing the effects of radiation-matter interactions, thus providing an approximation to quantum dynamics using classical inputs. This approach was previously formulated and assessed numerically for a single anharmonic degree of freedom. Our previous work is generalized here in two respects. First, the derivation of the OMT is extended to any number of coupled anharmonic vibrations by determining semiclassical approximations for pairs of double-sided Feynman diagrams. Second, an efficient numerical procedure is developed for calculating two-dimensional infrared spectra of coupled anharmonic vibrations in the OMT approximation. The OMT approximation is shown to reproduce the fundamental features of the quantum response function including both coherence and population dynamics.

Dynamics in the isotropic phase of nematogens using 2D IR vibrational echo measurements on natural-abundance 13CN and extended lifetime probes

The long time scale orientational relaxation of nematogens in the isotropic phase is associated with the randomization of pseudonematic domains, which have a correlation length that grows as the isotropic-to-nematic phase transition temperature is approached from above. Here we begin to address the fast dynamics of the nematogen molecules within the domains using two-dimensional infrared (2D IR) vibrational echo experiments. The problems of performing ultrafast IR experiments in pure liquids are discussed, and solutions are presented. In addition, the issue of short vibrational lifetimes, which limit the ability of 2D IR experiments to examine dynamics over a wide range of times, is addressed. The experiments were performed on the nematogen 4-cyano-4'-pentylbiphenyl (5CB), with the CN stretch initially used as the vibrational probe. Although the CN stretch has a small transition dipole, because the sample is a pure liquid it is necessary to use an exceedingly thin sample to perform the experiments. The small sample volume leads to massive heating effects that distort the results. In addition, the high concentration in the pure liquid can result in vibrational excitation transfer that interferes with the measurements of structural dynamics, and the CN vibrational lifetime is very short (3.6 ps). These problems were overcome by performing the experiments on the natural-abundance (13)CN stretch (5(13)CB), which greatly reduced the absorbance, eliminating the heating problems; also, this stretch has a longer lifetime (7.9 ps). Experiments were also performed on benzonitrile, which showed that the heating problems associated with pure liquids are not unique to 5CB. Again, the problems were eliminated by conducting measurements on the (13)CN stretch, which has an even longer lifetime (20.2 ps) compared with the (12)CN stretch (5.6 ps). Finally, to extend the range of the dynamical measurements, 4-pentyl-4'-thiocyanobiphenyl (5SCB) was synthesized and studied as a dilute solute in 5CB. The CN stretch of 5SCB has a vibrational lifetime of 103 ps, which permits dynamical measurements to 200 ps, revealing the full range of fast structural dynamics in the isotropic phase of 5CB. It is shown that the 5SCB probe reports essentially the same dynamics as 5(13)CB on the short time scale that is observable with the 5(13)CB vibrational probe.

Sensitive, site-specific, and stable vibrational probe of local protein environments: 4-azidomethyl-L-phenylalanine

We have synthesized the unnatural amino acid (UAA), 4-azidomethyl-L-phenylalanine (pN₃CH₂Phe), to serve as an effective vibrational reporter of local protein environments. The position, extinction coefficient, and sensitivity to local environment of the azide asymmetric stretch vibration of pN₃CH₂Phe are compared to the vibrational reporters: 4-cyano-L-phenylalanine (pCNPhe) and 4-azido-L-phenylalanine (pN₃Phe). This UAA was genetically incorporated in a site-specific manner utilizing an engineered, orthogonal aminoacyl-tRNA synthetase in response to an amber codon with high efficiency and fidelity into two distinct sites in superfolder green fluorescent protein (sfGFP). This allowed for the dependence of the azide asymmetric stretch vibration of pN₃CH₂Phe to different protein environments to be measured. The photostability of pN₃CH₂Phe was also measured relative to the photoreactive UAA, pN₃Phe.

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

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

An optimized semiclassical approximation for vibrational response functions

The observables of multidimensional infrared spectroscopy may be calculated from nonlinear vibrational response functions. Fully quantum dynamical calculations of vibrational response functions are generally impractical, while completely classical calculations are qualitatively incorrect at long times. These challenges motivate the development of semiclassical approximations to quantum mechanics, which use classical mechanical information to reconstruct quantum effects. The mean-trajectory (MT) approximation is a semiclassical approach to quantum vibrational response functions employing classical trajectories linked by deterministic transitions representing the effects of the radiation-matter interaction. Previous application of the MT approximation to the third-order response function R(3)(t3, t2, t1) demonstrated that the method quantitatively describes the coherence dynamics of the t3 and t1 evolution times, but is qualitatively incorrect for the waiting-time t2 period. Here we develop an optimized version of the MT approximation by elucidating the connection between this semiclassical approach and the double-sided Feynman diagrams (2FD) that represent the quantum response. Establishing the direct connection between 2FD and semiclassical paths motivates a systematic derivation of an optimized MT approximation (OMT). The OMT uses classical mechanical inputs to accurately reproduce quantum dynamics associated with all three propagation times of the third-order vibrational response 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.

Dynamics in the interior of AOT lamellae investigated with two-dimensional infrared spectroscopy

The dynamics inside the organic regions of aerosol-OT (AOT)/water mixtures in the lamellar mesophase, bicontinuous cubic (BC) phase, and in an analogous molecule without the charged sulfonate headgroup are investigated by observing spectral diffusion, orientational relaxation and population relaxation using ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy and IR pump-probe experiments on the asymmetric CO stretch of a vibrational probe, tungsten hexacarbonyl (W(CO)6). The water layer thickness between the bilayer planes in the lamellar phase was varied. For comparison, the dynamics of W(CO)6 in the normal liquid bis(2-ethylhexyl) succinate (EHS), which is analogous to AOT but has no charged sulfonate headgroup, were also studied. The 2D IR experiments measure spectral diffusion, which results from the structural evolution of the system. Spectral diffusion is quantified by the frequency-frequency correlation function (FFCF). In addition to a homogeneous component, the FFCFs are biexponential decays with fast and slow time components of ∼12.5 and ∼150 ps in the lamellar phase. Both components of the FFCF are independent of the number of water molecules per headgroup for the lamellae, but they slow somewhat in the BC phase. The dynamics in the ordered phases are in sharp contrast to the dynamics in EHS, which displays fast and slow components of the FFCF of 5 and 80 ps, respectively. As the hydration level of AOT increases, vibrational lifetime decreases, suggesting some change in the local environment of W(CO)6 with water content.

Vibrational solvatochromism and electrochromism. II. Multipole analysis

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

Ligand binding studied by 2D IR spectroscopy using the azidohomoalanine label

We explore the capability of the azidohomoalanine (Aha) as a vibrational label for 2D IR spectroscopy to study the binding of the target peptide to the PDZ2 domain. The Aha label responds sensitively to its local environment and its peak extinction coefficient of 350-400 M(-1) cm(-1) is high enough to routinely measure it in the low millimolar concentration regime. The central frequency, inhomogeneous width and spectral diffusion times deduced from the 2D IR line shapes of the Aha label at various positions in the peptide sequence is discussed in relationship to the known X-ray structure of the peptide bound to the PDZ2 domain. The results suggest that the Aha label introduces only a small perturbation to the overall structure of the peptide in the binding pocket. Finally, Aha is a methionine analog that can be incorporated also into larger proteins at essentially any position using protein expression. Altogether, Aha thus fulfills the requirements a versatile label should have for studies of protein structure and dynamics by 2D IR spectroscopy.

Folding and self-assembly of a small protein complex

The synthetic homotetrameric ββα (BBAT1) protein possesses a stable quaternary structure with a ββα fold. Because of its small size (a total of 84 residues), the homotetramer is an excellent model system with which to study the self-assembly and protein-protein interactions. We find from replica exchange molecular dynamics simulations with the coarse-grain UNRES force field that the folding and association pathway consists of three well-separated steps, where that association to a tetramer precedes and facilitates folding of the four chains. At room temperature the tetramer exists in an ensemble of diverse structures. The crystal structure becomes energetically favored only when the molecule is put in a dense and crystal-like environment. The observed picture of folding promoted by association may mirror the mechanism according to which intrinsically unfolded proteins assume their functional structure.

Fast dynamics of HP35 for folded and urea-unfolded conditions

The changes in fast dynamics of HP35 with a double CN vibrational dynamics label (HP35-P(2)) as a function of the extent of denaturation by urea were investigated with two-dimensional infrared (2D IR) vibrational echo spectroscopy. Cyanophenylalanine (PheCN) replaces the native phenylalanine at two residues in the hydrophobic core of HP35, providing vibrational probes. NMR data show that HP35-P(2) maintains the native folded structure similar to wild type and that both PheCN residues share essentially the same environment within the peptide. A series of time-dependent 2D IR vibrational echo spectra were obtained for the folded peptide and the increasingly unfolded peptide. Analysis of the time dependence of the 2D spectra yields the system's spectral diffusion, which is caused by the sampling of accessible structures of the peptide under thermal equilibrium conditions. The structural dynamics become faster as the degree of unfolding is increased.

Conformational dynamics and stability of HP35 studied with 2D IR vibrational echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.

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.