Simulation of vibrational energy transfer in two-dimensional infrared spectroscopy of amide I and amide II modes in solution

Refining protein amide I spectrum simulations with simple yet effective electrostatic models for local wavenumbers and dipole derivative magnitudes

Analysis of the amide I band of proteins is probably the most wide-spread application of bioanalytical infrared spectroscopy. Although highly desirable for a more detailed structural interpretation, a quantitative description of this absorption band is still difficult. This work optimized several electrostatic models with the aim to reproduce the effect of the protein environment on the intrinsic wavenumber of a local amide I oscillator. We considered the main secondary structures - α-helices, parallel and antiparallel β-sheets - with a maximum of 21 amide groups. The models were based on the electric potential and/or the electric field component along the CO bond at up to four atoms in an amide group. They were bench-marked by comparison to Hessian matrices reconstructed from density functional theory calculations at the BPW91, 6-31G** level. The performance of the electrostatic models depended on the charge set used to calculate the electric field and potential. Gromos and DSSP charge sets, used in common force fields, were not optimal for the better performing models. A good compromise between performance and the stability of model parameters was achieved by a model that considered the electric field at the positions of the oxygen, nitrogen, and hydrogen atoms of the considered amide group. The model describes also some aspects of the local conformation effect and performs similar on its own as in combination with an explicit implementation of the local conformation effect. It is better than a combination of a local hydrogen bonding model with the local conformation effect. Even though the short-range hydrogen bonding model performs worse, it captures important aspects of the local wavenumber sensitivity to the molecular surroundings. We improved also the description of the coupling between local amide I oscillators by developing an electrostatic model for the dependency of the dipole derivative magnitude on the protein environment.

Two-Dimensional Infrared Spectroscopy of Isolated Molecular Ions

Two-dimensional infrared (2D IR) spectroscopy of mass-selected, cryogenically cooled molecular ions is presented. Nonlinear response pathways, encoded in the time-domain photodissociation action response of weakly bound N2 messenger tags, were isolated using pulse shaping techniques following excitation with four collinear ultrafast IR pulses. 2D IR spectra of Re(CO)3(CH3CN)3+ ions capture off-diagonal cross-peak bleach signals between the asymmetric and symmetric carbonyl stretching transitions. These cross peaks display intensity variations as a function of pump-probe delay time due to coherent coupling between the vibrational modes. Well-resolved 2D IR features in the congested fingerprint region of protonated caffeine (C8H10N4O2H+) are also reported. Importantly, intense cross-peak signals were observed at 3 ps waiting time, indicating that tag-loss dynamics are not competing with the measured nonlinear signals. These demonstrations pave the way for more precise studies of molecular interactions and dynamics that are not easily obtainable with current condensed-phase methodologies.

Structural Dependence of Extended Amide III Vibrations in Two-Dimensional Infrared Spectra

Two-dimensional infrared (2D-IR) spectroscopy is a powerful experimental method for probing the structure and dynamics of proteins in aqueous solution. So far, most experimental studies have focused on the amide I vibrations, for which empirical vibrational exciton models provide a means of interpreting such experiments. However, such models are largely lacking for other regions of the vibrational spectrum. To close this gap, we employ an efficient quantum-chemical methodology for the calculation of 2D-IR spectra, which is based on anharmonic theoretical vibrational spectroscopy with localized modes. We apply this approach to explore the potential of 2D-IR spectroscopy in the extended amide III region. Using calculations for a dipeptide model as well as alanine polypeptides, we show that distinct 2D-IR cross-peaks in the extended amide III region can potentially be used to distinguish α-helix and β-strand structures. We propose that the extended amide III region could be a promising target for future 2D-IR experiments.

Simulation of two-dimensional infrared Raman spectroscopy with application to proteins

Two-dimensional infrared Raman spectroscopy is a powerful technique for studying the structure and interaction in molecular and biological systems. Here, we present a new implementation of the simulation of the two-dimensional infrared Raman signals. The implementation builds on the numerical integration of the Schrödinger equation approach. It combines the prediction of dynamics from molecular dynamics with a map-based approach for obtaining Hamiltonian trajectories and response function calculations. The new implementation is tested on the amide-I region for two proteins, where one is dominated by α-helices and the other by β-sheets. We find that the predicted spectra agree well with experimental observations. We further find that the two-dimensional infrared Raman spectra at least of the studied proteins are much less sensitive to the laser polarization used compared to conventional two-dimensional infrared experiments. The present implementation and findings pave the way for future applications for the interpretation of two-dimensional infrared Raman spectra.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

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.

Probing the Differential Dynamics of the Monomeric and Dimeric Insulin from Amide-I IR Spectroscopy

The monomer-dimer equilibrium for insulin is one of the essential steps in forming the receptor-binding competent monomeric form of the hormone. Despite this importance, the thermodynamic stability, in particular for modified insulins, is quite poorly understood, in part, due to experimental difficulties. This work explores one- and two-dimensional infrared spectroscopy in the range of the amide-I band for the hydrated monomeric and dimeric wild-type hormone. It is found that for the monomer the frequency fluctuation correlation function (FFCF) and the one-dimensional infrared spectra are position sensitive. The spectra of the -CO probes at the dimerization interface (residues Phe24, Phe25, and Tyr26) split and indicate an asymmetry despite the overall (formal) point symmetry of the dimer structure. Also, the decay times of the FFCF for the same -CO probe in the monomer and the dimer can differ by up to 1 order of magnitude, for example, for residue PheB24, which is solvent exposed for the monomer but at the interface for the dimer. The spectroscopic shifts correlate approximately with the average number of hydration waters and the magnitude of the FFCF at time zero. However, this correlation is only qualitative due to the heterogeneous and highly dynamical environment. Based on density functional theory calculations, the dominant contribution for solvent-exposed -CO is found to arise from the surrounding water (∼75%), whereas the protein environment contributes considerably less. The results suggest that infrared spectroscopy is a positionally sensitive probe of insulin dimerization, in particular in conjunction with isotopic labeling of the probe.

Development of Vibrational Frequency Maps for Nucleobases

Vibrational spectroscopy provides a powerful tool to probe the structure and dynamics of nucleic acids because specific normal modes, particularly the base carbonyl stretch modes, are highly sensitive to the hydrogen bonding patterns and stacking configurations in these biomolecules. In this work, we develop vibrational frequency maps for the C═O and C═C stretches in nucleobases that allow the calculations of their site frequencies directly from molecular dynamics simulations. We assess the frequency maps by applying them to nucleobase derivatives in aqueous solutions and nucleosides in organic solvents and demonstrate that the predicted infrared spectra are in good agreement with experimental measurements. The frequency maps can be readily used to model the linear and nonlinear vibrational spectroscopy of nucleic acids and elucidate the molecular origin of the experimentally observed spectral features.

Hydrophobic Collapse in N-Methylacetamide-Water Mixtures

Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump-probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with "frozen" amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.

Simple Quantum Dynamics with Thermalization

In this paper, we introduce two simple quantum dynamics methods. One is based on the popular surface-hopping method, and the other is based on rescaling of the propagation on the bath ground-state potential surface. The first method is special, as it avoids specific feedback from the simulated quantum system to the bath and can be applied for precalculated classical trajectories. It is based on the equipartition theorem to determine if hops between different potential energy surfaces are allowed. By comparing with the formally exact Hierarchical Equations Of Motion approach for four model systems we find that the method generally approximates the quantum dynamics toward thermal equilibrium very well. The second method is based on rescaling of the nonadiabatic coupling and also neglect the effect of the state of the quantum system on the bath. By the nature of the approximations, they cannot reproduce the effect of bath relaxation following excitation. However, the methods are both computationally more tractable than the conventional fewest switches surface hopping, and we foresee that the methods will be powerful for simulations of quantum dynamics in systems with complex bath dynamics, where the system-bath coupling is not too strong compared to the thermal energy.

The dynamics of peptide-water interactions in dialanine: An ultrafast amide I 2D IR and computational spectroscopy study

We present a joint experimental and computational study of the dynamic interactions of dialanine (Ala-Ala) with water, comparing the results of ultrafast 2D IR and infrared transient absorption spectroscopy of its amide I vibration with spectra modeled from molecular dynamics (MD) simulations. The experimental data are analyzed to describe vibrational frequency fluctuations, vibrational energy relaxation, and chemical exchange processes. The origin of these processes in the same underlying fluctuating forces allows a common description in terms of the fluctuations and conformational dynamics of the peptide and associated solvent. By comparing computational spectroscopy from MD simulations with multiple force fields and water models, we describe how the dynamics of water hydrogen bond fluctuations and switching processes act as a source of friction that governs the dephasing and vibrational relaxation, and provide a description of coupled water and peptide motions that give rise to spectroscopic exchange processes.

Application of two-dimensional infrared spectroscopy to benchmark models for the amide I band of proteins

In this paper, we present a novel benchmarking method for validating the modelling of vibrational spectra for the amide I region of proteins. We use the linear absorption spectra and two-dimensional infrared spectra of four experimentally well-studied proteins as a reference and test nine combinations of molecular dynamics force fields, vibrational frequency mappings, and coupling models. We find that two-dimensional infrared spectra provide a much stronger test of the models than linear absorption does. The best modelling approach in the present study still leaves significant room for future improvement. The presented benchmarking scheme, thus, provides a way of validating future protocols for modelling the amide I band in proteins.

Isotope-enriched protein standards for computational amide I spectroscopy

We present a systematic isotope labeling study of the protein G mutant NuG2b as a step toward the production of reliable, structurally stable, experimental standards for amide I infrared spectroscopic simulations. By introducing isotope enriched amino acids into a minimal growth medium during bacterial expression, we induce uniform labeling of the amide bonds following specific amino acids, avoiding the need for chemical peptide synthesis. We use experimental data to test several common amide I frequency maps and explore the influence of various factors on map performance. Comparison of the predicted absorption frequencies for the four maps tested with empirical assignments to our experimental spectra yields a root-mean-square error of 6-12 cm(-1), with outliers of at least 12 cm(-1) in all models. This means that the predictions may be useful for predicting general trends such as changes in hydrogen bonding configuration; however, for finer structural constraints or absolute frequency assignments, the models are unreliable. The results indicate the need for careful testing of existing literature maps and shed light on possible next steps for the development of quantitative spectral maps.

Distinguishing gramicidin D conformers through two-dimensional infrared spectroscopy of vibrational excitons

Gramicidin D is a short peptide which dimerizes to form helical pores, adopting one of two conformations in the process. These conformations differ primarily in number of residues per turn and the hydrogen-bond registry between rungs of the helix. Using amide I 2D infrared (IR) and FTIR, we have demonstrated that it is possible to distinguish between the different conformers of gramicidin D in solution. We show that the spectra observed for this helical peptide bear no resemblance to the spectra of α- or 310-helices and that while the FTIR spectra appear similar to spectra of β-sheets, 2D IR reveals that the observed resonances arise from vibrational modes unlike those observed in β-sheets. We also present an idealized model which reproduces the experimental data with high fidelity. This model is able to explain the polarization-dependence of the experimental 2D IR data. Using this model, we show the coupling between the rungs of the helix dominates the spectra, and as a consequence of this, the number of residues per turn can greatly influence the amide I spectra of gramicidin D.

Visualizing KcsA conformational changes upon ion binding by infrared spectroscopy and atomistic modeling

The effect of ion binding in the selectivity filter of the potassium channel KcsA is investigated by combining amide I Fourier-transform infrared spectroscopy with structure-based spectral modeling. Experimental difference IR spectra between K(+)-bound KcsA and Na(+)-bound KcsA are in good qualitative agreement with spectra modeled from structural ensembles generated from molecular dynamics simulations. The molecular origins of the vibrational modes contributing to differences in these spectra are determined not only from structural differences in the selectivity filter but also from the pore helices surrounding this region. Furthermore, the coordination of K(+) or Na(+) to carbonyls in the selectivity filter effectively decouples the vibrations of those carbonyls from the rest of the protein, creating local probes of the electrostatic environment. The results suggest that it is necessary to include the influence of the surrounding helices in discussing selectivity and transport in KcsA and, on a more general level, that IR spectroscopy offers a nonperturbative route to studying the structure and dynamics of ion channels.

Vibrational ladder climbing in carboxy-hemoglobin: effects of the protein environment

We present simulations on vibrational ladder climbing in carboxy-hemoglobin. Motivated by recent experiments, we study the influence of different realistic pump probe parameters. To allow for a direct comparison with experimental results, transient absorption spectra obtained by a weak probe pulse following the strong, shaped pump pulse are calculated. The influence of the protein fluctuations is taken into account using a recently developed microscopic model. This model consists of a quantum Hamiltonian describing the CO vibration in carboxy-hemoglobin, together with a fluctuating potential, which is obtained by electronic structure calculation based on a large number of protein configurations. Using realistic pulse parameters, vibrational excitations to very high-lying states are possible, in qualitative agreement with experimental observations.

Vibrational solvatochromism. II. A first-principle theory of solvation-induced vibrational frequency shift based on effective fragment potential method

Vibrational solvatochromism is a solvation-induced effect on fundamental vibrational frequencies of molecules in solutions. Here we present a detailed first-principle coarse-grained theory of vibrational solvatochromism, which is an extension of our previous work [B. Błasiak, H. Lee, and M. Cho, J. Chem. Phys. 139(4), 044111 (2013)] by taking into account electrostatic, exchange-repulsion, polarization, and charge-transfer interactions. By applying our theory to the model N-methylacetamide-water clusters, solute-solvent interaction-induced effects on amide I vibrational frequency are fully elucidated at Hartree-Fock level. Although the electrostatic interaction between distributed multipole moments of solute and solvent molecules plays the dominant role, the contributions from exchange repulsion and induced dipole-electric field interactions are found to be of comparable importance in short distance range, whereas the charge-transfer effect is negligible. The overall frequency shifts calculated by taking into account the contributions of electrostatics, exchange-repulsion, and polarization terms are in quantitative agreement with ab initio results obtained at the Hartree-Fock level of theory.

Vibrational energy transport in molecules studied by relaxation-assisted two-dimensional infrared spectroscopy

This review presents an overview of the relaxation-assisted two-dimensional infrared (RA 2DIR) spectroscopy method for measuring structures and energy transport dynamics in molecules. The method strongly enhances the range of accessible distances compared to traditional 2DIR and offers new structural reporters, such as the energy transport time, cross-peak amplification factors, and connectivity patterns. The use of the method for assigning vibrational modes with various levels of delocalization is illustrated. RA 2DIR relies on vibrational energy transport in molecules; as such, the transport mechanism can be conveniently studied by the method. Applications to identify diffusive and ballistic energy transport are demonstrated.

Parallel β-sheet vibrational couplings revealed by 2D IR spectroscopy of an isotopically labeled macrocycle: quantitative benchmark for the interpretation of amyloid and protein infrared spectra

Infrared spectroscopy is playing an important role in the elucidation of amyloid fiber formation, but the coupling models that link spectra to structure are not well tested for parallel β-sheets. Using a synthetic macrocycle that enforces a two stranded parallel β-sheet conformation, we measured the lifetimes and frequency for six combinations of doubly (13)C═(18)O labeled amide I modes using 2D IR spectroscopy. The average vibrational lifetime of the isotope labeled residues was 550 fs. The frequencies of the labels ranged from 1585 to 1595 cm(-1), with the largest frequency shift occurring for in-register amino acids. The 2D IR spectra of the coupled isotope labels were calculated from molecular dynamics simulations of a series of macrocycle structures generated from replica exchange dynamics to fully sample the conformational distribution. The models used to simulate the spectra include through-space coupling, through-bond coupling, and local frequency shifts caused by environment electrostatics and hydrogen bonding. The calculated spectra predict the line widths and frequencies nearly quantitatively. Historically, the characteristic features of β-sheet infrared spectra have been attributed to through-space couplings such as transition dipole coupling. We find that frequency shifts of the local carbonyl groups due to nearest neighbor couplings and environmental factors are more important, while the through-space couplings dictate the spectral intensities. As a result, the characteristic absorption spectra empirically used for decades to assign parallel β-sheet secondary structure arises because of a redistribution of oscillator strength, but the through-space couplings do not themselves dramatically alter the frequency distribution of eigenstates much more than already exists in random coil structures. Moreover, solvent exposed residues have amide I bands with >20 cm(-1) line width. Narrower line widths indicate that the amide I backbone is solvent protected inside the macrocycle. This work provides calculated and experimentally verified couplings for parallel β-sheets that can be used in structure-based models to simulate and interpret the infrared spectra of β-sheet containing proteins and protein assemblies, such as amyloid fibers.

Solvent and conformation dependence of amide I vibrations in peptides and proteins containing proline

We present a mixed quantum-classical model for studying the amide I vibrational dynamics (predominantly CO stretching) in peptides and proteins containing proline. There are existing models developed for determining frequencies of and couplings between the secondary amide units. However, these are not applicable to proline because this amino acid has a tertiary amide unit. Therefore, a new parametrization is required for infrared-spectroscopic studies of proteins that contain proline, such as collagen, the most abundant protein in humans and animals. Here, we construct the electrostatic and dihedral maps accounting for solvent and conformation effects on frequency and coupling for the proline unit. We examine the quality and the applicability of these maps by carrying out spectral simulations of a number of peptides with proline in D(2)O and compare with experimental observations.

Frequency distribution of the amide-I vibration sorted by residues in amyloid fibrils revealed by 2D-IR measurements and simulations

The infrared optical response of amyloid fibrils Aβ(1-40) is investigated. Simulations of two models corresponding to different protonation states are compared with experiment. The simulations reveal that vibrational frequency distributions inside the fibrils are dominated by side chain fluctuations. We further confirm earlier suggestions based on 2D-IR measurements that water molecules can be trapped inside the fibrils.

Redistribution of carbonyl stretch mode energy in isolated and solvated N-methylacetamide: kinetic energy spectral density analyses

The vibrational energy transfer from the excited carbonyl stretch mode in N-deuterated N-methylacetamide (NMA-d), both in isolation and in a heavy water cluster, is studied with nonequilibrium molecular dynamics (NEMD) simulations, employing a quantum mechanical/molecular mechanical (QM∕MM) force field at the semiempirical PM3 level. The nonequilibrium ensemble of vibrationally excited NMA-d is prepared by perturbing the positions and velocities of the carbonyl C and O atoms and its NEMD trajectories are obtained with a leap-frog algorithm properly modified for the initial perturbation. In addition to the time-domain analysis of the kinetic and potential energies, a novel method for the spectral analysis of the atomic kinetic energies is developed, in terms of the spectral density of kinetic energy, which provides the time-dependent changes of the frequency-resolved kinetic energies without the complications of normal mode analysis at every MD time step. Due to the QM description of the solute electronic structure, the couplings among the normal modes are captured more realistically than with classical force fields. The energy transfer in the isolated NMA-d is found to proceed first from the carbonyl bond to other modes with time scales of 3 ps or less, and then among the other modes over 3-21 ps. In the solvated NMA-d, most of the excess energy is first transferred to other intramolecular modes within 5 ps, which is subsequently dissipated to solvent with 7-19 ps time scales. The contribution of the direct energy transfer from the carbonyl bond to solvent was only 5% with ~7 ps time scale. Solvent reorganization that leads to destabilization of the electrostatic interactions is found to be crucial in the long time relaxation of the excess energy, while the water intramolecular modes do not contribute significantly. Detailed mode-specific energy transfer pathways are deduced for the isolated and solvated NMA-d and they show that the energy transfer in NMA-d is a highly cooperative process among the intramolecular modes and there is no single dominant pathway with more than 30% of transient contribution.

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.

Quantum decoherence in finite size exciton-phonon systems

Based on the operatorial formulation of the perturbation theory, the properties of a confined exciton coupled with phonons in thermal equilibrium is revisited. Within this method, the dynamics is governed by an effective Hamiltonian which accounts for exciton-phonon entanglement. The exciton is dressed by a virtual phonon cloud whereas the phonons are clothed by virtual excitonic transitions. Special attention is thus paid for describing the time evolution of the excitonic coherences at finite temperature. As in an infinite lattice, temperature-enhanced quantum decoherence takes place. However, it is shown that the confinement softens the decoherence. The coherences are very sensitive to the excitonic states so that the closer to the band center the state is located, the slower the coherence decays. In particular, for odd lattice sizes, the coherence between the vacuum state and the one-exciton state exactly located at the band center survives over an extremely long time scale. A superimposition involving the vacuum and this specific one-exciton state behaves as an ideal qubit insensitive to its environment.

Development and validation of transferable amide I vibrational frequency maps for peptides

Infrared (IR) spectroscopy of the amide I band has been widely utilized for the analysis of peptides and proteins. Theoretical modeling of IR spectra of proteins requires an accurate and efficient description of the amide I frequencies. In this paper, amide I frequency maps for protein backbone and side chain groups are developed from experimental spectra and vibrational lifetimes of N-methylacetamide and acetamide in different solvents. The frequency maps, along with established nearest-neighbor frequency shift and coupling schemes, are then applied to a variety of peptides in aqueous solution and reproduce experimental spectra well. The frequency maps are designed to be transferable to different environments; therefore, they can be used for heterogeneous systems, such as membrane proteins.

Sulfoxide stretching mode as a structural reporter via dual-frequency two-dimensional infrared spectroscopy

The S=O stretching mode in sulfoxides, having a frequency in the 950-1150 cm(-1) range, is tested as a structural label via dual-frequency two-dimensional infrared (2DIR) spectroscopy. The properties of this structural reporter are studied in several compounds, including (4,4(')-dimethyl-2,2(')-bipyridyl)(o-methylsulfinylbenzoate) ruthenium II, [Ru(dmb)(2)(BzSO)](+), (RuBzSO), octylsulfinylpropionic acid (OSPA), and o- and p-methylsulfinyl-benzoic acid (oMSBA and pMSBA). The mode assignment in the fingerprint region for these compounds is made using a combination of density functional theory calculations and 2DIR and relaxation-assisted 2DIR (RA 2DIR) spectroscopies. The SO stretching mode frequency and IR intensity demonstrate substantial sensitivity to the molecular structure. Multiple cross peaks of the C=O and S=O stretching modes with modes in the fingerprint region (930-1450 cm(-1)) were recorded. The 2DIR and RA 2DIR spectra focusing at interactions of a high-frequency mode of a ligand with the modes in the fingerprint region provide a spectral fingerprint of a compound and help mode assignment in the often congested fingerprint region. The cross-peak amplitudes in oMSBA, pMSBA, and OSPA were compared with the theoretical predictions based on the computed values for the off-diagonal anharmonicities and a reasonable match is found. The SO stretching mode provides means for assigning other modes in the fingerprint region and constitutes a promising structural reporter for the 2DIR and RA 2DIR spectroscopy measurements.

Modeling the vibrational dynamics and nonlinear infrared spectra of coupled amide I and II modes in peptides

The amide vibrational modes play an important role in energy transport and relaxation in polypeptides and proteins and provide us with spectral markers for structure and structural dynamics of these macromolecules. Here, we present a detailed model to describe the dynamic properties of the amide I and amide II modes and the resulting linear and nonlinear spectra. These two modes have large oscillator strengths, and their mutual coupling plays an important role in their relaxation. Using first-principles calculations of NMA-d(7) and a dipeptide in a fluctuating bath described by molecular dynamics simulations, we model the frequencies of the local vibrations as well as the coupling between them. Both the coherent couplings and the fluctuations induced by contact with their environment are taken into account. We apply the resulting model of interacting fluctuating oscillators to study the collective vibrations and the partially coherent transport of vibrational energy through a model α-helix. We find that the instantaneous vibrations are delocalized over a few (up to four) amide units, while the coherences in the helix survive for 0.5-1 ps, leading to coherent transport on a similar time scale.

Amide I'-II' 2D IR spectroscopy provides enhanced protein secondary structural sensitivity

We demonstrate how multimode 2D IR spectroscopy of the protein amide I' and II' vibrations can be used to distinguish protein secondary structure. Polarization-dependent amide I'-II' 2D IR experiments on poly-l-lysine in the beta-sheet, alpha-helix, and random coil conformations show that a combination of amide I' and II' diagonal and cross peaks can effectively distinguish between secondary structural content, where amide I' infrared spectroscopy alone cannot. The enhanced sensitivity arises from frequency and amplitude correlations between amide II' and amide I' spectra that reflect the symmetry of secondary structures. 2D IR surfaces are used to parametrize an excitonic model for the amide I'-II' manifold suitable to predict protein amide I'-II' spectra. This model reveals that the dominant vibrational interaction contributing to this sensitivity is a combination of negative amide II'-II' through-bond coupling and amide I'-II' coupling within the peptide unit. The empirically determined amide II'-II' couplings do not significantly vary with secondary structure: -8.5 cm(-1) for the beta sheet, -8.7 cm(-1) for the alpha helix, and -5 cm(-1) for the coil.

Toward detecting the formation of a single helical turn by 2D IR cross peaks between the amide-I and -II modes

We have combined two-dimensional infrared (2D IR) spectroscopy and isotope substitutions to reveal the vibrational couplings between a pair of amide-I and -II modes that are several residues away but directly connected through a hydrogen bond in a helical peptide. This strategy is demonstrated on a 3(10)-helical hexapeptide, Z-Aib-L-Leu-(Aib)2-Gly-Aib-OtBu, and its 13C=18O-Leu monolabeled and 13C=18O-Leu/15N-Gly bis-labeled isotopomers in CDCl3. The isotope-dependent amide-I/II cross peaks clearly show that the second and fourth peptide linkages are vibrationally coupled as they are in proximity, forming a 3(10)-helical turn. The experimental spectra are compared to simulations based on a vibrational exciton Hamiltonian model that fully takes into account the amide-I and -II modes. The amide-II local mode frequency is evaluated by a new model based on the effects of hydrogen-bond geometry and sites. Ab initio nearest-neighbor coupling maps of the amide-I/I, -I/II, -II/I and -II/II modes are generated by isotopically isolating the local modes of N-acetyl-glycine N'-methylamide (AcGlyNHMe). Longer range couplings are modeled by transition charge interactions. The effects of the capping groups are incorporated and isotope effects are analyzed based on ab initio calculations of six model compounds. The main features of the 2D IR spectra are reproduced by this modeling. The conformational sensitivity of the isotope-dependent amide-I/II cross peaks is discussed in comparison with the calculated spectra for a semiextended structure. Our experimental and theoretical study demonstrates that the combination of 2D IR and 13C=18O/15N labeling is a useful structural method for detecting helical turn formation with residue-level specificity.

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.

Waiting time dynamics in two-dimensional infrared spectroscopy

We review recent work on the waiting time dynamics of coherent two-dimensional infrared (2DIR) spectroscopy. This dynamics can reveal chemical and physical processes that take place on the femto- and picosecond time scale, which is faster than the time scale that may be probed by, for example, nuclear magnetic resonance spectroscopy. A large number of chemically relevant processes take place on this time scale. Such processes range from forming and breaking hydrogen bonds and proton transfer to solvent exchange and vibrational population transfer. In typical 2DIR spectra, multiple processes contribute to the waiting time dynamics and the spectra are often congested. This makes the spectra challenging to interpret, and the aid of theoretical models and simulations is often needed. To be useful, such models need to account for all dynamical processes in the sample simultaneously. The numerical integration of the Schrodinger equation (NISE) method has proven to allow for a very general treatment of the dynamical processes. It accounts for both the motional narrowing resulting from solvent-induced frequency fluctuations and population transfer between coupled vibrations. At the same time, frequency shifts arising from chemical-exchange reactions and changes of the transition dipoles because of either non-Condon effects or molecular reorientation are included in the treatment. This method therefore allows for the disentanglement of all of these processes. The NISE method has thus far been successfully applied to study chemical-exchange processes. It was demonstrated that 2DIR is not only sensitive to reaction kinetics but also to the more detailed reaction dynamics. NISE has also been applied to the study of population transfer within the amide I band (CO stretch) and between the amide I and amide II bands (CN stretch and NH bend) in polypeptides. From the amide I studies, it was found that the population transfer can be used to enhance cross-peaks that act as structural markers for beta-sheet structure in proteins. From the amide I/II investigation, it was found that the amide II band and the hydrogen-bond stretch vibration are important parts of the relaxation pathway for the amide I vibration. With the development of simple approximations, it becomes possible to apply the NISE method even to very big systems, such as the OH stretch of bulk water, which can only be described well when large numbers of coupled vibrations are taken into account.

2D-IR experiments and simulations of the coupling between amide-I and ionizable side chains in proteins: application to the Villin headpiece

The carboxylate side chains of Asp and Glu have significant coupling with the amide states of the backbone of the Villin headpiece. In two-dimensional spectroscopy, cross peaks are observed between these side chains and the main amide-I band. To model the absorption of the side chains, the electric field variations of vibrational frequencies of a carboxylic acid group (neutral form, CH(3)-COOH) and a carboxylate group (ionized form, CH(3)-COO(-)) are parametrized by means of density functional theory calculations. Simulations indicate that the side chains significantly couple to only one or two amide-I modes out of all of the amino acid residues which makes them useful as spectroscopic markers, providing information about the local structural behavior of the protein. Both experiment and simulations show that the cross peaks between the carboxylate and the amide-I bands are significantly diminished above the melting temperature.