Vibrational echo correlation spectroscopy probes of hydrogen bond dynamics in water and methanol

Fast Collective Hydrogen-Bond Dynamics in Hexafluoroisopropanol Related to its Chemical Activity

Using fluorinated mono-alcohols, in particular hexafluoro-isopropanol (HFIP), as a solvent can enhance chemical reaction rates in a spectacular manner. Previous work has shown evidence that this enhancement is related to the hydrogen-bond structure of these liquids. Here, we investigate the hydrogen-bond dynamics of HFIP and compare it to that of its non-fluorinated analog, isopropanol. Ultrafast infrared spectroscopy experiments show that the dynamics of individual hydrogen-bonds is about twice as slow in HFIP as in isopropanol. Surprisingly, from dielectric spectroscopy we find the opposite behavior for the dynamics of hydrogen-bonded clusters: collective rearrangements are 3 times faster in HFIP than in isopropanol. This difference indicates that the hydrogen-bonded clusters in HFIP are smaller than in isopropanol. The differences in cluster size can be traced to changes in the hydrogen-bond donor and acceptor strengths upon fluorination. The smaller cluster size can boost reaction rates in HFIP by increasing the concentration of reactive, terminal OH-groups of the clusters, whereas the fast collective dynamics can increase the rate of formation of hydrogen-bonds with the reactants. The longer lifetime of the individual hydrogen-bonds in HFIP can enhance the stability of the hydrogen-bonded clusters, and so increase the probability of reactant-solvent hydrogen-bonding.

Enhanced Chemical Sensing with Multiorder Coherent Raman Scattering Spectroscopic Dephasing

Molecular vibrational spectroscopy is widely used in various sensing and imaging applications, providing intrinsic information at the molecular level. Nonlinear optical interactions using ultrashort laser pulses facilitate the selective coherent excitation of molecular vibrational modes by focusing energy into specific molecular bonds, boosting the signal level for multiple orders of magnitude. The dephasing of such coherence, which is susceptible to the local molecular environment, however, is often neglected. The unique capability of vibrational dephasing dynamics to serve as a unique probe for complex molecular interactions and the effect of local nano- and microenvironments are beyond the reach of conventional, intensity-based spectroscopy. Here, we developed a novel multiorder coherent Raman spectroscopy platform with a special focus on the temporal evolution of molecular vibrational dephasing, termed as time-resolved coherent Raman scattering (T-CRS) spectroscopy. By utilizing a high dynamic range detection, molecular vibrational dynamics and the environmental effects are demonstrated with multidimensional spectroscopic sensing, which promises a new range of applications in biology, materials, and chemical sciences.

Ultrafast vibrational dynamics of a solute correlates with dynamics of the solvent

Two-dimensional infrared (2D-IR) spectroscopy is used to measure the spectral dynamics of the metal carbonyl complex cyclopentadienyl manganese tricarbonyl (CMT) in a series of linear alkyl nitriles. 2D-IR spectroscopy provides direct readout of solvation dynamics through spectral diffusion, probing the decay of frequency correlation induced by fluctuations of the solvent environment. 2D-IR simultaneously monitors intramolecular vibrational energy redistribution (IVR) among excited vibrations, which can also be influenced by the solvent through the spectral density rather than the dynamical friction underlying solvation. Here, we report that the CMT vibrational probe reveals solvent dependences in both the spectral diffusion and the IVR time scales, where each slows with increased alkyl chain length. In order to assess the degree to which solute-solvent interactions can be correlated with bulk solvent properties, we compared our results with low-frequency dynamics obtained from optical Kerr effect (OKE) spectroscopy-performed by others-on the same nitrile solvent series. We find excellent correlation between our spectral diffusion results and the orientational dynamics time scales from OKE. We also find a correlation between our IVR time scales and the amplitudes of the low-frequency spectral densities evaluated at the 90-cm^-1 energy difference, corresponding to the gap between the two strong vibrational modes of the carbonyl probe. 2D-IR and OKE provide complementary perspectives on condensed phase dynamics, and these findings provide experimental evidence that at least at the level of dynamical correlations, some aspects of a solute vibrational dynamics can be inferred from properties of the solvent.

Structural and dynamical changes observed when transitioning from an ionic liquid to a deep eutectic solvent

The microscopic molecular structure and dynamics of a new deep eutectic solvent (DES) composed of an ionic liquid (1-hexyl-3-methylimidazolium chloride) and an amide (trifluoroacetamide) at various molar ratios were investigated using linear and non-linear infrared spectroscopy with a vibrational probe. The use of the ionic liquid allows us to investigate the changes that the system undergoes with the addition of the amide or, equivalently, the changes from an ionic liquid to a DES. Our studies revealed that the vibrational probe in the DES senses a very similar local environment irrespective of the cation chemical structure. In addition, the amide also appears to perceive the same molecular environment. The concentration dependence studies also showed that the amide changes from being isolated from other amides in the ionic liquid environment to an environment where the amide-amide interactions are favored. In the case of the vibrational probe, the addition of the amide produced significant changes in the slow dynamics associated with the making and breaking of the ionic cages but did not affect the rattling-in-cage motions perceived by it. Furthermore, the concentration dependence of slow dynamics showed two regimes which are linked to the changes in the overall structure of the solution. These observations are interpreted in the context of a nanoscopic heterogeneous environment in the DES which, according to the observed dynamical regimes, appears at very large concentrations of the amide (molar ratio of greater than 1:1) since for lower amide molar ratios, the amide appears to be not segregated from the ionic liquid. This proposed molecular picture is supported by small angle x-ray scattering experiments.

Study on hydrogen bonding network in aqueous methanol solution by Raman spectroscopy

The Raman spectra of aqueous methanol solution with various concentrations were measured at room temperature and atmospheric pressure. We found that the CO stretching vibration mode of methanol showed a significant blue shift at V_m (V_m represents the volume fraction of methanol) >0.4, while the CH symmetric and asymmetric stretching vibration modes exhibited red shift under the same conditions. These results illustrate that the variation of hydrogen bond between methanol and water molecules lead to a phase transition of the methanol-water complex at V_m = 0.4. Furthermore, the red shift of the CH vibration mode indicates that there is no hydrogen bond formed between the CH₃ group of methanol and water molecules. In addition, we found that the frequency shift of C-H is affected by the hydrogen bond C-O…H-O formed between methanol and water molecules, and the corresponding theoretical discussion is given. Finally, the phase transition process of methanol-water complex in methanol-water binary solution was given by theoretical analysis.

Histidine Orientation Modulates the Structure and Dynamics of a de Novo Metalloenzyme Active Site

The ultrafast dynamics of a de novo metalloenzyme active site is monitored using two-dimensional infrared spectroscopy. The homotrimer of parallel, coiled coil α-helices contains a His3-Cu(I) metal site where CO is bound and serves as a vibrational probe of the hydrophobic interior of the self-assembled complex. The ultrafast spectral dynamics of Cu-CO reveals unprecedented ultrafast (2 ps) nonequilibrium structural rearrangements launched by vibrational excitation of CO. This initial rapid phase is followed by much slower ∼40 ps vibrational relaxation typical of metal-CO vibrations in natural proteins. To identify the hidden coupled coordinate, small molecule analogues and the full peptide were studied by QM and QM/MM calculations, respectively. The calculations show that variation of the histidines' dihedral angles in coordinating Cu controls the coupling between the CO stretch and the Cu-C-O bending coordinates. Analysis of different optimized structures with significantly different electrostatic field magnitudes at the CO ligand site indicates that the origin of the stretch-bend coupling is not directly due to through-space electrostatics. Instead, the large, ∼3.6 D dipole moments of the histidine side chains effectively transduce the electrostatic environment to the local metal coordination orientation. The sensitivity of the first coordination sphere to the protein electrostatics and its role in altering the potential energy surface of the bound ligands suggests that long-range electrostatics can be leveraged to fine-tune function through enzyme design.

Background-Free Fourth-Order Sum Frequency Generation Spectroscopy

The recently developed 2D sum frequency generation spectroscopy offers new possibilities to analyze the structure and structural dynamics of interfaces in a surface-specific manner. Its implementation, however, has so far remained limited to the pump-probe geometry, with its inherent restrictions. Here we present 2D SFG experiments utilizing a novel noncollinear geometry of four incident laser pulses generating a 2D SFG response, analogous to the triangle geometry applied in bulk-sensitive 2D infrared spectroscopy. This approach allows for background-free measurements of fourth-order nonlinear signals, which is demonstrated by measuring the fourth-order material response from a GaAs (110) surface. The implementation of phase-sensitive detection and broadband excitation pulses allows for both highest possible time resolution and high spectral resolution of the pump axis of a measured 2D SFG spectrum. To reduce the noise in our spectra, we employ a referencing procedure, for which we use noncollinear pathways and individual focusing for the signal and local oscillator beams. The 2D spectra recorded from the GaAs (110) surface show nonzero responses for the real and imaginary component, pointing to contributions from resonant electronic pathways to the χ((4)) response.

A novel setup for femtosecond pump-repump-probe IR spectroscopy with few cycle CEP stable pulses

We present a three-color mid-IR setup for vibrational pump-repump-probe experiments with a temporal resolution well below 100 fs and a freely selectable spectral resolution of 20 to 360 cm(-1) for the pump and repump. The usable probe range without optical realignment is 900 cm(-1). The experimental design employed is greatly simplified compared to the widely used setups, highly robust and includes a novel means for generation of tunable few-cycle pulses with stable carrier-envelope phase. A Ti:sapphire pump system operating with 1 kHz and a modest 150 fs pulse duration supplies the total pump energy of just 0.6 mJ. The good signal-to-noise ratio of the setup allows the determination of spectrally resolved transient probe changes smaller than 6·10(-5) OD at 130 time delays in just 45 minutes. The performance of the spectrometer is demonstrated with transient IR spectra and decay curves of HDO molecules in lithium nitrate trihydrate and ice and a first all MIR pump-repump-probe measurement.

Simplified and economical 2D IR spectrometer design using a dual acousto-optic modulator

Over the last decade two-dimensional infrared (2D IR) spectroscopy has proven to be a very useful extension of infrared spectroscopy, yet the technique remains restricted to a small group of specialized researchers because of its experimental complexity and high equipment cost. We report on a spectrometer that is compact, mechanically robust, and is much less expensive than previous designs because it uses a single pixel MCT detector rather than an array detector. Moreover, each axis of the spectrum can be collected in either the time or frequency domain via computer programming. We discuss pulse sequences for scanning the probe axis, which were not previously possible. We present spectra on metal carbonyl compounds at 5 µm and a model peptide at 6 µm. Data collection with a single pixel MCT takes longer than using an array detector, but publishable quality data are still achieved with only a few minutes of averaging.

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

A nitrile-labeled amino acid, p-cyanophenylalanine, is introduced near the active site of the semisynthetic enzyme ribonuclease S to serve as a probe of protein dynamics and fluctuations. Ribonuclease S is the limited proteolysis product of subtilisin acting on ribonuclease A, and consists of a small fragment including amino acids 1-20, the S-peptide, and a larger fragment including residues 21-124, the S-protein. A series of two-dimensional vibrational echo experiments performed on the nitrile-labeled S-peptide and the RNase S are described. The time-dependent changes in the two-dimensional infrared vibrational echo line shapes are analyzed using the center line slope method to obtain the frequency-frequency correlation function (FFCF). The observations show that the nitrile probe in the S-peptide has dynamics that are similar to, but faster than, those of the single amino acid p-cyanophenylalanine in water. In contrast, the dynamics of the nitrile label when the peptide is bound to form ribonuclease S are dominated by homogeneous dephasing (motionally narrowed) contributions with only a small contribution from very fast inhomogeneous structural dynamics. The results provide insights into the nature of the structural dynamics of the ribonuclease S complex. The equilibrium dynamics of the nitrile labeled S-peptide and the ribonuclease S complex are also investigated by molecular dynamics simulations. The experimentally determined FFCFs are compared to the FFCFs obtained from the molecular dynamics simulations, thereby testing the capacity of simulations to determine the amplitudes and time scales of protein structural fluctuations on fast time scales under thermal equilibrium conditions.

Mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon-tetrachloride mixtures II: excited state hydrogen bonding structure and dynamics, infrared emission spectrum, and excited state lifetime

We present a mixed quantum-classical molecular dynamics study of the hydrogen-bonding structure and dynamics of a vibrationally excited hydroxyl stretch in methanol/carbon-tetrachloride mixtures. The adiabatic Hamiltonian of the quantum-mechanical hydroxyl is diagonalized on-the-fly to obtain the ground and first-excited adiabatic energy levels and wave functions which depend parametrically on the instantaneous configuration of the classical degrees of freedom. The dynamics of the classical degrees of freedom are determined by Hellmann-Feynman forces obtained by taking the expectation value of the force with respect to the ground or excited vibrational wave functions. Polarizable force fields are used which were previously shown to reproduce the experimental infrared absorption spectrum rather well, for different isotopomers and over a wide composition range [Kwac, K.; Geva, E. J. Phys. Chem. B 2011, 115, 9184]. We show that the agreement of the absorption spectra with experiment can be further improved by accounting for the dependence of the dipole moment derivatives on the configuration of the classical degrees of freedom. We find that the propensity of a methanol molecule to form hydrogen bonds increases upon photoexcitation of its hydroxyl stretch, thereby leading to a sizable red-shift of the corresponding emission spectrum relative to the absorption spectrum. Treating the relaxation from the first excited to the ground state as a nonadiabatic process, and calculating its rate within the framework of Fermi's golden rule and the harmonic-Schofield quantum correction factor, we were able to predict a lifetime which is of the same order of magnitude as the experimental value. The experimental dependence of the lifetime on the transition frequency is also reproduced. Nonlinear mapping relations between the hydroxyl transition frequency and bond length in the excited state and the electric field along the hydroxyl bond axis are established. These mapping relations make it possible to reduce the computational cost of the mixed quantum-classical treatment to that of a fully classical treatment.

Structure dependent energy transport: relaxation-assisted 2DIR measurements and theoretical studies

Vibrational energy relaxation and transport in a molecule that is far from thermal equilibrium can affect its chemical reactivity. Understanding the energy transport dynamics in such molecules is also important for measuring molecular structural constraints via relaxation-assisted two-dimensional infrared (RA 2DIR) spectroscopy. In this paper we investigated vibrational relaxation and energy transport in the ortho, meta, and para isomers of acetylbenzonitrile (AcPhCN) originated from excitation of the CN stretching mode. The amplitude of the cross-peak among the CN and CO stretching modes served as an indicator for the energy transport from the CN group toward the CO group. A surprisingly large difference is observed in both the lifetimes of the CN mode and in the energy transport rates for the three isomers. The anharmonic DFT calculations and energy transport modeling performed to understand the origin of the differences and to identify the main cross-peak contributors in these isomers described well the majority of the experimental results including mode excited-state lifetimes and the energy transport dynamics. The strong dependence of the energy transport on molecular structure found in this work could be useful for recognizing different isomers of various compounds via RA 2DIR spectroscopy.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

Conformational switching between protein substates studied with 2D IR vibrational echo spectroscopy and molecular dynamics simulations

Myoglobin is an important protein for the study of structure and dynamics. Three conformational substates have been identified for the carbonmonoxy form of myoglobin (MbCO). These are manifested as distinct peaks in the IR absorption spectrum of the CO stretching mode. Ultrafast 2D IR vibrational echo chemical exchange experiments are used to observed switching between two of these substates, A(1) and A(3), on a time scale of <100 ps for two mutants of wild-type Mb. The two mutants are a single mutation of Mb, L29I, and a double mutation, T67R/S92D. Molecular dynamics (MD) simulations are used to model the structural differences between the substates of the two MbCO mutants. The MD simulations are also employed to examine the substate switching in the two mutants as a test of the ability of MD simulations to predict protein dynamics correctly for a system in which there is a well-defined transition over a significant potential barrier between two substates. For one mutant, L29I, the simulations show that translation of the His64 backbone may differentiate the two substates. The simulations accurately reproduce the experimentally observed interconversion time for the L29I mutant. However, MD simulations exploring the same His64 backbone coordinate fail to display substate interconversion for the other mutant, T67R/S92D, thus pointing to the likely complexity of the underlying protein interactions. We anticipate that understanding conformational dynamics in MbCO via ultrafast 2D IR vibrational echo chemical exchange experiments can help to elucidate fast conformational switching processes in other proteins.

The role of intercalated water in multilayered graphene oxide

A detailed in situ infrared spectroscopy analysis of single layer and multilayered graphene oxide (GO) thin films reveals that the normalized infrared absorption in the carbonyl region is substantially higher in multilayered GO upon mild annealing. These results highlight the fact that the reduction chemistry of multilayered GO is dramatically different from the single layer GO due to the presence of water molecules confined in the ∼1 nm spacing between sheets. IR spectroscopy, XPS analysis, and DFT calculations all confirm that the water molecules play a significant role interacting with basal plane etch holes through passivation, via evolution of CO(2) leading to the formation of ketone and ester carbonyl groups. Displacement of water from intersheet spacing with alcohol significantly changes the chemistry of carbonyl formation with temperature.

Relaxation-assisted two-dimensional infrared (RA 2DIR) method: accessing distances over 10 A and measuring bond connectivity patterns

Development of new approaches for measuring three-dimensional structures and dynamics of structural changes is important for a number of natural sciences, including structural biology, where it can lead to understanding the physical bases of molecular recognition and catalysis. A two-dimensional infrared (2DIR) spectroscopy method permits measuring pairwise interactions among vibrational modes in molecules providing a molecular scale ruler for delivering structural constraints, such as the distances between the vibrational modes, angles between their transition dipoles, and the energy-transfer rates between them. While there is a large variety of systems that have recently been interrogated using 2DIR, questions remain of how to measure structural features of larger molecules. The challenges of working with larger molecules, such as proteins, include very congested vibrational spectra, a small range of distances accessible by the 2DIR method, and sensitivity issues. This Account describes the efforts of our laboratory to overcome some of these challenges. First, we discuss the dual-frequency 2DIR approach, which provides the highest selectivity to a particular pair of vibrational reporters and highest sensitivity. Second, we describe our steps in developing vibrational labels, novel for 2DIR, such as C identical withN and C-D stretching modes that have frequencies in the water transparency region, as well as the modes in the fingerprint region. The schemes suitable for labeling amino acids are discussed. Next, we describe the novel relaxation-assisted 2DIR (RA 2DIR) method, developed in our laboratory. The method uses vibrational relaxation and vibrational energy transport in molecules and the thermalization process on a molecular scale, to generate stronger cross-peaks. An 18-fold cross-peak amplification was observed for the modes separated by about 11 A using the RA 2DIR method, and larger amplifications are expected for larger distances between the modes. Large amplification provided by the RA 2DIR method enhances the sensitivity of 2DIR spectroscopy and permits longer range structural measurements. In addition to generating stronger cross-peaks, a correlation of the energy transport time with the intermode distance is demonstrated. This correlation permits measurements of mode-connectivity patterns in molecules much similar to those available in total correlation spectroscopy (TOCSY) and heteronuclear multiple-bond correlation (HMBC) methods of 2D nuclear magnetic resonance (NMR) spectroscopy. It is our hope that, with a proper calibration, the RA 2DIR method will permit speedy assessments of distances and the bond connectivity patterns in molecules and reach the level of an analytical method.

How to turn your pump-probe instrument into a multidimensional spectrometer: 2D IR and Vis spectroscopies via pulse shaping

We have recently developed a new and simple way of collecting 2D infrared and visible spectra that utilizes a pulse shaper and a partly collinear beam geometry. 2D IR and Vis spectroscopies are powerful tools for studying molecular structures and their dynamics. They can be used to correlate vibrational or electronic eigenstates, measure energy transfer rates, and quantify the dynamics of lineshapes, for instance, all with femtosecond time-resolution. As a result, they are finding use in systems that exhibit fast dynamics, such as sub-millisecond chemical and biological dynamics, and in hard-to-study environments, such as in membranes. While powerful, these techniques have been difficult to implement because they require a series of femtosecond pulses to be spatially and temporally overlapped with precise time-resolution and interferometric phase stability. However, many of the difficulties associated with implementing 2D spectroscopies are eliminated by using a pulse shaper and a simple beam geometry, which substantially lowers the technical barriers required for researchers to enter this exciting field while simultaneously providing many new capabilities. The aim of this paper is to provide an overview of the methods for collecting 2D spectra so that an outsider considering using 2D spectroscopy in their own research can judge which approach would be most suitable for their research aims. This paper focuses primarily on 2D IR spectroscopy, but also includes our recent work on adapting this technology to collecting 2D Vis spectra. We review work that has already been published as well as cover several topics that we have not reported previously, including phase cycling methods to remove background signals, eliminate unwanted scatter, and shift data collection into the rotating frame.

Phasing problem of heterodyne-detected two-dimensional infrared spectroscopy

A rigorous method is presented to measure and adjust the phase difference between the two pulse pairs used in heterodyne-detected 2D-IR spectroscopy with an accuracy better than |Deltaphi|=0.1 rad. The method, which can easily be automated, avoids the otherwise tedious measurement of the much weaker pump-probe spectrum as a reference, which is the commonly used approach to phase 2D-IR spectra in a postprocessing step.

Native and unfolded cytochrome c--comparison of dynamics using 2D-IR vibrational echo spectroscopy

Unfolded vs native CO-coordinated horse heart cytochrome c (h-cyt c) and a heme axial methionine mutant cyt c552 from Hydrogenobacter thermophilus ( Ht-M61A) are studied by IR absorption spectroscopy and ultrafast 2D-IR vibrational echo spectroscopy of the CO stretching mode. The unfolding is induced by guanidinium hydrochloride (GuHCl). The CO IR absorption spectra for both h-cyt c and Ht-M61A shift to the red as the GuHCl concentration is increased through the concentration region over which unfolding occurs. The spectra for the unfolded state are substantially broader than the spectra for the native proteins. A plot of the CO peak position vs GuHCl concentration produces a sigmoidal curve that overlays the concentration-dependent circular dichroism (CD) data of the CO-coordinated forms of both Ht-M61A and h-cyt c within experimental error. The coincidence of the CO peak shift curve with the CD curves demonstrates that the CO vibrational frequency is sensitive to the structural changes induced by the denaturant. 2D-IR vibrational echo experiments are performed on native Ht-M61A and on the protein in low- and high-concentration GuHCl solutions. The 2D-IR vibrational echo is sensitive to the global protein structural dynamics on time scales from subpicosecond to greater than 100 ps through the change in the shape of the 2D spectrum with time (spectral diffusion). At the high GuHCl concentration (5.1 M), at which Ht-M61A is essentially fully denatured as judged by CD, a very large reduction in dynamics is observed compared to the native protein within the approximately 100 ps time window of the experiment. The results suggest the denatured protein may be in a glassy-like state involving hydrophobic collapse around the heme.

Solute-solvent complex kinetics and thermodynamics probed by 2D-IR vibrational echo chemical exchange spectroscopy

The formation and dissociation kinetics of a series of triethylsilanol/solvent weakly hydrogen bonding complexes with enthalpies of formation ranging from -1.4 to -3.3 kcal/mol are measured with ultrafast two-dimensional infrared (2D IR) chemical exchange spectroscopy in liquid solutions at room temperature. The correlation between the complex enthalpies of formation and dissociation rate constants can be expressed with an equation similar to the Arrhenius equation. The experimental results are in accord with previous observations on eight phenol/solvent complexes with enthalpies of formation from -0.6 to -2.5 kcal/mol. It was found that the inverse of the solute-solvent complex dissociation rate constant is linearly related to exp(-DeltaH0/RT) where DeltaH0 is the complex enthalpy of formation. It is shown here, that the triethylsilanol-solvent complexes obey the same relationship with the identical proportionality constant, that is, all 13 points, five silanol complexes and eight phenol complexes, fall on the same line. In addition, features of 2D IR chemical exchange spectra at long reaction times (spectral diffusion complete) are explicated using the triethylsilanol systems. It is shown that the off-diagonal chemical exchange peaks have shapes that are a combination (outer product) of the absorption line shapes of the species that give rise to the diagonal peaks.

Disulfide bond influence on protein structural dynamics probed with 2D-IR vibrational echo spectroscopy

Intramolecular disulfide bonds are understood to play a role in regulating protein stability and activity. Because disulfide bonds covalently link different components of a protein, they influence protein structure. However, the effects of disulfide bonds on fast (subpicosecond to approximately 100 ps) protein equilibrium structural fluctuations have not been characterized experimentally. Here, ultrafast 2D-IR vibrational echo spectroscopy is used to examine the constraints an intramolecular disulfide bond places on the structural fluctuations of the protein neuroglobin (Ngb). Ngb is a globin family protein found in vertebrate brains that binds oxygen reversibly. Like myoglobin (Mb), Ngb has the classical globin fold and key residues around the heme are conserved. Furthermore, the heme-ligated CO vibrational spectra of Mb (Mb-CO) and Ngb (Ngb-CO) are virtually identical. However, in contrast to Mb, human Ngb has an intramolecular disulfide bond that affects its oxygen affinity and protein stability. By using 2D-IR vibrational echo spectroscopy, we investigated the equilibrium protein dynamics of Ngb-CO by observing the CO spectral diffusion (time dependence of the 2D-IR line shapes) with and without the disulfide bond. Despite the similarity of the linear FTIR spectra of Ngb-CO with and without the disulfide bond, 2D-IR measurements reveal that the equilibrium sampling of different protein configurations is accelerated by disruption of the disulfide bond. The observations indicate that the intramolecular disulfide bond in Ngb acts as an inhibitor of fast protein dynamics even though eliminating it does not produce significant conformational change in the protein's structure.

Frequency-frequency correlation functions and apodization in two-dimensional infrared vibrational echo spectroscopy: a new approach

Ultrafast two-dimensional infrared (2D-IR) vibrational echo spectroscopy can probe structural dynamics under thermal equilibrium conditions on time scales ranging from femtoseconds to approximately 100 ps and longer. One of the important uses of 2D-IR spectroscopy is to monitor the dynamical evolution of a molecular system by reporting the time dependent frequency fluctuations of an ensemble of vibrational probes. The vibrational frequency-frequency correlation function (FFCF) is the connection between the experimental observables and the microscopic molecular dynamics and is thus the central object of interest in studying dynamics with 2D-IR vibrational echo spectroscopy. A new observable is presented that greatly simplifies the extraction of the FFCF from experimental data. The observable is the inverse of the center line slope (CLS) of the 2D spectrum. The CLS is the inverse of the slope of the line that connects the maxima of the peaks of a series of cuts through the 2D spectrum that are parallel to the frequency axis associated with the first electric field-matter interaction. The CLS varies from a maximum of 1 to 0 as spectral diffusion proceeds. It is shown analytically to second order in time that the CLS is the T(w) (time between pulses 2 and 3) dependent part of the FFCF. The procedure to extract the FFCF from the CLS is described, and it is shown that the T(w) independent homogeneous contribution to the FFCF can also be recovered to yield the full FFCF. The method is demonstrated by extracting FFCFs from families of calculated 2D-IR spectra and the linear absorption spectra produced from known FFCFs. Sources and magnitudes of errors in the procedure are quantified, and it is shown that in most circumstances, they are negligible. It is also demonstrated that the CLS is essentially unaffected by Fourier filtering methods (apodization), which can significantly increase the efficiency of data acquisition and spectral resolution, when the apodization is applied along the axis used for obtaining the CLS and is symmetrical about tau=0. The CLS is also unchanged by finite pulse durations that broaden 2D spectra.

Neuroglobin dynamics observed with ultrafast 2D-IR vibrational echo spectroscopy

Neuroglobin (Ngb), a protein in the globin family, is found in vertebrate brains. It binds oxygen reversibly. Compared with myoglobin (Mb), the amino acid sequence has limited similarity, but key residues around the heme and the classical globin fold are conserved in Ngb. The CO adduct of Ngb displays two CO absorption bands in the IR spectrum, referred to as N(3) (distal histidine in the pocket) and N(0) (distal histidine swung out of the pocket), which have absorption spectra that are almost identical with the Mb mutants L29F and H64V, respectively. The Mb mutants mimic the heme pocket structures of the corresponding Ngb conformers. The equilibrium protein dynamics for the CO adduct of Ngb are investigated by using ultrafast 2D-IR vibrational echo spectroscopy by observing the CO vibration's spectral diffusion (2D-IR spectra time dependence) and comparing the results with those for the Mb mutants. Although the heme pocket structure and the CO FTIR peak positions of Ngb are similar to those of the mutant Mb proteins, the 2D-IR results demonstrate that the fast structural fluctuations of Ngb are significantly slower than those of the mutant Mbs. The results may also provide some insights into the nature of the energy landscape in the vicinity of the folded protein free energy minimum.

Dynamics around solutes and solute-solvent complexes in mixed solvents

Ultrafast 2D-IR vibrational echo experiments, IR pump-probe experiments, and FT-IR spectroscopy of the hydroxyl stretch of phenol-OD in three solvents, CCl4, mesitylene (1, 3, 5 trimethylbenzene), and the mixed solvent of mesitylene and CCl4 (0.83 mole fraction CCl4), are used to study solute-solvent dynamics via observation of spectral diffusion. Phenol forms a complex with Mesitylene. In the mesitylene solution, there is only complexed phenol; in the CCl4 solution, there is only uncomplexed phenol; and in the mixed solvent, both phenol species are present. Dynamics of the free phenol in CCl4 or the mixed solvent are very similar, and dynamics of the complex in mesitylene and in the mixed solvent are very similar. However, there are differences in the slowest time scale dynamics between the pure solvents and the mixed solvents. The mixed solvent produces slower dynamics that are attributed to first solvent shell solvent composition variations. The composition variations require a longer time to randomize than is required in the pure solvents, where only density variations occur. The experimental results and recent MD simulations indicate that the solvent structure around the solute may be different from the mixed solvent's mole fraction.

Automated 2D IR spectroscopy using a mid-IR pulse shaper and application of this technology to the human islet amyloid polypeptide

The capability of 2D IR spectroscopy to elucidate time-evolving structures is enhanced by a programmable mid-IR pulse shaper that greatly improves the ease, speed, and accuracy of data collection. Traditional ways of collecting 2D IR spectra are difficult to implement, cause distorted peak shapes, and result in poor time resolution and/or phase problems. We report on several methods for collecting 2D IR spectra by using a computer-controlled germanium acoustooptic modulator that overcomes the above problems. The accuracy and resolution of each method is evaluated by using model metal carbonyl compounds that have well defined lineshapes. Furthermore, phase cycling can now be employed to largely alleviate background scatter from heterogeneous samples. With these methods in hand, we apply 2D IR spectroscopy to study the structural diversity in amyloid fibers of aggregated human islet amyloid polypeptide (hIAPP), which is involved with type 2 diabetes. The 2D IR spectra reveal that the beta-sheet fibers have a large structural distribution, as evidenced by an inhomogeneously broadened beta-sheet peak and strong coupling to random coil conformations. Structural diversity is an important characteristic of hIAPP because it may be that partly folded peptides cause the disease. This experiment on hIAPP is an example of how computer generation of 2D IR pulse sequences is a key step toward automating 2D IR spectroscopy, so that new pulse sequences can be implemented quickly and a diverse range of systems can be studied more easily.