Peptide conformational heterogeneity revealed from nonlinear vibrational spectroscopy and molecular-dynamics simulations

Intrinsic α-helical and β-sheet conformational preferences: a computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.

The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development

The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.

A coupled two-dimensional main chain torsional potential for protein dynamics: generation and implementation

Main chain torsions of alanine dipeptide are parameterized into coupled 2-dimensional Fourier expansions based on quantum mechanical (QM) calculations at M06 2X/aug-cc-pvtz//HF/6-31G** level. Solvation effect is considered by employing polarizable continuum model. Utilization of the M06 2X functional leads to precise potential energy surface that is comparable to or even better than MP2 level, but with much less computational demand. Parameterization of the 2D expansions is against the full main chain torsion space instead of just a few low energy conformations. This procedure is similar to that for the development of AMBER03 force field, except unique weighting factor was assigned to all the grid points. To avoid inconsistency between quantum mechanical calculations and molecular modeling, the model peptide is further optimized at molecular mechanics level with main chain dihedral angles fixed before the calculation of the conformational energy on molecular mechanical level at each grid point, during which generalized Born model is employed. Difference in solvation models at quantum mechanics and molecular mechanics levels makes this parameterization procedure less straightforward. All force field parameters other than main chain torsions are taken from existing AMBER force field. With this new main chain torsion terms, we have studied the main chain dihedral distributions of ALA dipeptide and pentapeptide in aqueous solution. The results demonstrate that 2D main chain torsion is effective in delineating the energy variation associated with rotations along main chain dihedrals. This work is an implication for the necessity of more accurate description of main chain torsions in the future development of ab initio force field and it also raises a challenge to the development of quantum mechanical methods, especially the quantum mechanical solvation models.

pH-Independence of trialanine and the effects of termini blocking in short peptides: a combined vibrational, NMR, UVCD, and molecular dynamics study

Several lines of evidence now well establish that unfolded peptides in general, and alanine in specific, have an intrinsic preference for the polyproline II (pPII) conformation. Investigation of local order in the unfolded state is, however, complicated by experimental limitations and the inherent dynamics of the system, which has in some cases yielded inconsistent results from different types of experiments. One method of studying these systems is the use of short model peptides, and specifically short alanine peptides, known for predominantly sampling pPII structure in aqueous solution. Recently, He et al. ( J. Am. Chem. Soc. 2012 , 134 , 1571 - 1576 ) proposed that unblocked tripeptides may not be suitable models for studying conformational propensities in unfolded peptides due to the presence of end effect, that is, electrostatic interactions between investigated amino acid residues and terminal charges. To determine whether changing the protonation states of the N- and C-termini influence the conformational manifold of the central amino acid residue in tripeptides, we have examined the pH-dependence of unblocked trialanine and the conformational preferences of alanine in the alanine dipeptide. To this end, we measured and globally analyzed amide I' band profiles and NMR J-coupling constants. We described conformational distributions as the superposition of two-dimensional Gaussian distributions assignable to specific subspaces of the Ramachandran plot. Results show that the conformational ensemble of trialanine as a whole, and the pPII content (χpPII = 0.84) in particular, remains practically unaffected by changing the protonation state. We found that compared to trialanine, the alanine dipeptide has slightly lower pPII content (χpPII = 0.74) and an ensemble more reminiscent of the unblocked Gly-Ala-Gly model peptide. In addition, a two-state thermodynamic analysis of the conformational sensitive Δε(T) and (3)J(H(N)H(α))(T) data obtained from electronic circular dichroism and H NMR spectra indicate that the free energy landscape of trialanine is similar in all protonation states. MD simulations for the investigated peptides corroborate this notion and show further that the hydration shell around unblocked trialanine is unaffected by the protonation/deprotonation of the C-terminal group. In contrast, the alanine dipeptide shows a reduced water density around the central residue as well as a less ordered hydration shell, which decreases the pPII propensity and reduces the lifetime of sampled conformations.

Disorder and order in unfolded and disordered peptides and proteins: a view derived from tripeptide conformational analysis. I. Tripeptides with long and predominantly hydrophobic side chains

We performed a conformational analysis of the central residues of three tripeptides glycyl-L-isoleucyl-glycine (GIG), glycyl-L-tyrosyl-glycine (GYG) and glycyl-L-arginyl-glycine (GRG) in aqueous solution, based on a global analysis of amide I' band profiles and NMR J-coupling constants. The results are compared with recently reported distributions of GVG, GFG and GEG. For GIG and GYG, we found that even though the polyproline II (pPII) fraction is below 0.5, it is still the most populated conformation, whereas GVG and GFG show both a larger β-strand fraction. For GRG, we observed a clear dominance of pPII over β-strand, reminiscent of observations for GEG and GKG. This finding indicates that terminal charges on otherwise hydrophobic residue side chains stabilize pPII over β-strand conformations. For all peptides investigated we found that a variety of compact and turn-like conformations constitute nearly 20 percent of their conformational distributions. Attempts to analyze our data with a simple two-state pPII-->/<--β model therefore do not yield any satisfactory reproduction of experimental results. A comparison of the obtained GxG ensembles with conformational distributions of GxG segments in truncated coil libraries (helices and sheets omitted) revealed a much larger fraction of type II β(i+2) and type III β like conformations for the latter. Thus, a comparison of conformational distributions of unfolded peptide segments in solution and in coil libraries reveal interesting information on how the interplay between intrinsic propensities of amino acid residues and non-local interactions in polypeptide chains determine the conformations of loop segments in proteins.

Hybrid quantum/classical simulations of the vibrational relaxation of the amide I mode of N-methylacetamide in D2O solution

Hybrid quantum/classical molecular dynamics (MD) is applied to simulate the vibrational relaxation (VR) of the amide I mode of deuterated N-methylacetamide (NMAD) in aqueous (D(2)O) solution. A novel version of the vibrational molecular dynamics with quantum transitions (MDQT) treatment is developed in which the amide I mode is treated quantum mechanically while the remaining degrees of freedom are treated classically. The instantaneous normal modes of the initially excited NMAD molecule (INM(0)) are used as internal coordinates since they provide a proper initial partition of the system in quantum and classical subsystems. The evolution in time of the energy stored in each individual normal mode is subsequently quantified using the hybrid quantum-classical instantaneous normal modes (INM(t)). The identities of both the INM(0)s and the INM(t)s are tracked using the equilibrium normal modes (ENMs) as templates. The results extracted from the hybrid MDQT simulations show that the quantum treatment of the amide I mode accelerates the whole VR process versus pure classical simulations and gives better agreement with experiments. The relaxation of the amide I mode is found to be essentially an intramolecular vibrational redistribution (IVR) process with little contribution from the solvent, in agreement with previous theoretical and experimental studies. Two well-defined relaxation mechanisms are identified. The faster one accounts for ≈40% of the total vibrational energy that flows through the NMAD molecule and involves the participation of the lowest frequency vibrations as short-life intermediate modes. The second and slower mechanism accounts for the remaining ≈60% of the energy released and is associated to the energy flow through specific mid-range and high-frequency modes.

A peptide's perspective of water dynamics

This Perspective is focused on amide groups of peptides interacting with water. The 2D IR spectroscopy has already enabled structural aspects of the peptide backbone to be determined through its ability to measure the coupling between different amide-I modes. Here we describe why nonlinear IR is emerging as the method of choice to examine the fast components of the water dynamics near peptides and how isotopically edited peptide links can be used to probe the local water at a residue level in proteins. This type of research necessarily involves an intimate mix of theory and experiment. The description of the results is underpinned by relatively well established quantum-statistical theories that describe the important manifestations of peptide vibrational frequency fluctuations.

Specific interactions of sodium salts with alanine dipeptide and tetrapeptide in water: insights from molecular dynamics

We examine computationally the dipeptide and tetrapeptide of alanine in pure water and solutions of sodium chloride (NaCl) and iodide (NaI), with salt concentrations up to 3 M. Enhanced sampling of the configuration space is achieved by the replica exchange method. In agreement with other works, we observe preferential sodium interactions with the peptide carbonyl groups, which are enhanced in the NaI solutions due to the increased affinity of the less hydrophilic iodide anion for the peptide methyl side-chains and terminal blocking groups. These interactions have been associated with a decrease in the helicities of more complex peptides. In our simulations, both salts have a small effect on the dipeptide, but consistently stabilize the intramolecular hydrogen-bonding interactions and "α-helical" conformations of the tetrapeptide. This behavior, and an analysis of the intermolecular interaction energies show that ion-peptide interactions, or changes in the peptide hydration due to salts, are not sufficient determining factors of the peptide conformational preferences. Additional simulations suggest that the observed stabilizing effect is not due to the employed force-field, and that it is maintained in short peptides but is reversed in longer peptides. Thus, the peptide conformational preferences are determined by an interplay of energetic and entropic factors, arising from the peptide sequence and length and the composition of the solution.

Anharmonic vibrational modes of nucleic acid bases revealed by 2D IR spectroscopy

Polarization-dependent two-dimensional infrared (2D IR) spectra of the purine and pyrimadine base vibrations of five nucleotide monophosphates (NMPs) were acquired in D(2)O at neutral pH in the frequency range 1500-1700 cm(-1). The distinctive cross-peaks between the ring deformations and carbonyl stretches of NMPs indicate that these vibrational modes are highly coupled, in contrast with the traditional peak assignment, which is based on a simple local mode picture such as C═O, C═N, and C═C double bond stretches. A model of multiple anharmonically coupled oscillators was employed to characterize the transition energies, vibrational anharmonicities and couplings, and transition dipole strengths and orientations. No simple or intuitive structural correlations are found to readily assign the spectral features, except in the case of guanine and cytosine, which contain a single local CO stretching mode. To help interpret the nature of these vibrational modes, we performed density functional theory (DFT) calculations and found that multiple ring vibrations are coupled and delocalized over the purine and pyrimidine rings. Generally, there is close correspondence between the experimental and computational results, provided that the DFT calculations include explicit waters solvating hydrogen-bonding sites. These results provide direct experimental evidence of the delocalized nature of the nucleotide base vibrations via a nonperturbative fashion and will serve as building blocks for constructing a structure-based model of DNA and RNA vibrational spectroscopy.

Temperature dependent equilibrium native to unfolded protein dynamics and properties observed with IR absorption and 2D IR vibrational echo experiments

Dynamic and structural properties of carbonmonoxy (CO)-coordinated cytochrome c(552) from Hydrogenobacter thermophilus (Ht-M61A) at different temperatures under thermal equilibrium conditions were studied with infrared absorption spectroscopy and ultrafast two-dimensional infrared (2D IR) vibrational echo experiments using the heme-bound CO as the vibrational probe. Depending on the temperature, the stretching mode of CO shows two distinct bands corresponding to the native and unfolded proteins. As the temperature is increased from low temperature, a new absorption band for the unfolded protein grows in and the native band decreases in amplitude. Both the temperature-dependent circular dichroism and the IR absorption area ratio R(A)(T), defined as the ratio of the area under the unfolded band to the sum of the areas of the native and unfolded bands, suggest a two-state transition from the native to the unfolded protein. However, it is found that the absorption spectrum of the unfolded protein increases its inhomogeneous line width and the center frequency shifts as the temperature is increased. The changes in line width and center frequency demonstrate that the unfolding does not follow simple two-state behavior. The temperature-dependent 2D IR vibrational echo experiments show that the fast dynamics of the native protein are virtually temperature independent. In contrast, the fast dynamics of the unfolded protein are slower than those of the native protein, and the unfolded protein fast dynamics and at least a portion of the slower dynamics of the unfolded protein change significantly, becoming faster as the temperature is raised. The temperature dependence of the absorption spectrum and the changes in dynamics measured with the 2D IR experiments confirm that the unfolded ensemble of conformers continuously changes its nature as unfolding proceeds, in contrast to the native state, which displays a temperature-independent distribution of structures.

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.

UV resonance Raman investigation of the conformations and lowest energy allowed electronic excited states of tri- and tetraalanine: charge transfer transitions

UV resonance Raman excitation profiles and Raman depolarization ratios were measured for trialanine and tetraalanine between 198 and 210 nm. Excitation within the pi --> pi* electronic transitions of the peptide bond results in UVRR spectra dominated by amide peptide bond vibrations. In addition to the resonance enhancement of the normal amide vibrations, we find enhancement of the symmetric terminal COO(-) vibration. The Ala(3) UVRR AmIII(3) band frequencies indicate that poly-proline II and 2.5(1) helix conformations and type II turns are present in solution. We also find that the conformation of the interior peptide bond of Ala(4) is predominantly poly-proline-II-like. The Raman excitation profiles of both Ala(3) and Ala(4) reveal a charge transfer electronic transition at 202 nm, where electron transfer occurs from the terminal nonbonding carboxylate orbital to the adjacent peptide bond pi* orbital. Raman depolarization ratio measurements support this assignment. An additional electronic transition is found in Ala(4) at 206 nm.

Probing the folding of mini-protein Beta3s by two-dimensional infrared spectroscopy; simulation study

We propose to use infrared coherent two-dimensional correlation spectroscopy (2DCS) to characterize the folding mechanism of the mini-protein Beta3s. In this study Beta3s was folded by molecular dynamics (MD) simulation and intermediate conformational ensembles were identified. The one and two-dimensional correlation spectrum was calculated for the intermediate and native states of the mini-protein. A direct structure-spectra relationship was determined by analysis of conformational properties and specific residue contributions. We identified the structural origin of diagonal and off-diagonal peaks in the 2DCS spectra for the native and intermediate conformational ensembles in the folding mechanism. This work supports the implementation of computational techniques in conjunction with experimental 2DCS to study the folding mechanism of proteins. In addition to exploring the folding mechanism the work presented here can be applied in combination with experiment to refine and validate current molecular dynamics force fields.

PACS Codes: 87.15.Cc, 87.15.hm, 87.15.hp

Interactions of tyrosine in Leu-enkephalin at a membrane-water interface: an ultrafast two-dimensional infrared study combined with density functional calculations and molecular dynamics simulations

The interactions of neuropeptides and membranes play an important role in peptide hormone function. Our current understanding of peptide-membrane interactions remains limited due to the paucity of experimental techniques capable of probing such interactions. In this work, we study the nature of opioid peptide-membrane interactions using ultrafast two-dimensional infrared (2D IR) spectroscopy. The high temporal resolution of 2D IR is particularly suited for studying highly flexible opioid peptides. We investigate the location of the tyrosine (Tyr) side chain of leucine-enkephalin (Lenk) in lipid bilayer membranes by measuring spectral diffusion of the phenolic ring vibrational mode in three different systems: Lenk in lipid bilayer membranes (bicelles), Lenk in deuterated water, and p-cresol in deuterated water. Frequency-frequency correlation functions obtained from waiting-time-dependent 2D IR spectra reveal an ultrafast decaying component with an approximately 1 ps time constant that is common for all three systems. On the basis of density functional theory calculations and molecular dynamics simulations, this spectral diffusion component is attributed to hydrogen-bond dynamics of the phenolic hydroxyl group interacting with bulk water. Unlike p-cresol in water, both Lenk systems exhibit static spectral inhomogeneity, which can be attributed to conformational distributions of Lenk that do not interconvert within 4 ps. Our results suggest that the Tyr side chain of Lenk in bicelles is located at the water-abundant region at the membrane-water interface and not embedded into the hydrophobic core.

Purely absorptive three-dimensional infrared spectroscopy

We demonstrate a method to collect purely absorptive three-dimensional (3D) fifth-order vibrational spectra on the model system CO(2) in H(2)O. The six beam interferometer is described, as well as a method to experimentally determine the phase of the 3D spectrum. The measured spectra agree very well with simulations of the data based on the cumulant expansion. There are five peaks corresponding to different paths up and down the vibrational ladder. The positions, signs, and amplitudes of the peaks agree with theoretical predictions, and the intensities of the peaks scale linearly with concentration. Based on the concentration dependence and agreement between the simulations and measurements, we conclude that cascaded lower order signals contribute negligibly to the observed signal.

Applications of 2D IR spectroscopy to peptides, proteins, and hydrogen-bond dynamics

Following a survey of 2D IR principles, this article describes recent experiments on the hydrogen-bond dynamics of small ions, amide-I modes, nitrile probes, peptides, reverse transcriptase inhibitors, and amyloid fibrils.

Multidimensional Langevin modeling of biomolecular dynamics

A systematic computational approach to describe the conformational dynamics of biomolecules in reduced dimensionality is presented. The method is based on (i) the decomposition of a high-dimensional molecular dynamics trajectory into a few "system" and (many) "bath" degrees of freedom and (ii) a Langevin simulation of the resulting model. Employing principal component analysis, the dimension of the system is chosen such that it contains all slow large-amplitude motions of the molecule, while the bath coordinates only account for its high-frequency fluctuations. It is shown that a sufficiently large dimension of the model is essential to ensure a clear time scale separation of system and bath variables, which warrants the validity of the memory-free Langevin equation. Applying methods from nonlinear time series analysis, a practical Langevin algorithm is presented which performs a local estimation of the multidimensional Langevin vector fields describing deterministic drift and stochastic driving. Adopting a 800 ns molecular dynamics simulation of the folding of heptaalanine in explicit water, it is shown that a five-dimensional Langevin model correctly reproduces the structure and conformational dynamics of the system. The virtues and limits of the approach are discussed in some detail.

Coherent multidimensional vibrational spectroscopy of biomolecules: concepts, simulations, and challenges

The response of complex molecules to sequences of femtosecond infrared pulses provides a unique window into their structure, dynamics, and fluctuating environments. Herein we survey the basic principles of modern two-dimensional infrared (2DIR) spectroscopy, which analogous to those of multidimensional NMR spectroscopy. The perturbative approach for computing the nonlinear optical response of coupled localized chromophores is introduced and applied to the amide backbone transitions of proteins, liquid water, membrane lipids, and amyloid fibrils. The signals are analyzed using classical molecular dynamics simulations combined with an effective fluctuating Hamiltonian for coupled localized anharmonic vibrations whose dependence on the local electrostatic environment is parameterized by an ab initio map. Several simulation methods, (cumulant expansion of Gaussian fluctuation, quasiparticle scattering, the stochastic Liouville equations, direct numerical propagation) are surveyed. Chirality-induced techniques which dramatically enhance the resolution are demonstrated. Signatures of conformational and hydrogen-bonding fluctuations, protein folding, and chemical-exchange processes are discussed.

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.

Taking apart the two-dimensional infrared vibrational echo spectra: more information and elimination of distortions

Ultrafast two-dimensional infrared (2D-IR) vibrational echo spectroscopy can probe the fast structural evolution of molecular systems under thermal equilibrium conditions. Structural dynamics are tracked by observing the time evolution of the 2D-IR spectrum, which is caused by frequency fluctuations of vibrational mode(s) excited during the experiment. However, there are a variety of effects that can produce line shape distortions and prevent the correct determination of the frequency-frequency correlation function (FFCF), which describes the frequency fluctuations and connects the experimental observables to a molecular level depiction of dynamics. In addition, it can be useful to analyze different parts of the 2D spectrum to determine if dynamics are different for subensembles of molecules that have different initial absorption frequencies in the inhomogeneously broadened absorption line. Here, an important extension to a theoretical method for extraction of the FFCF from 2D-IR spectra is described. The experimental observable is the center line slope (CLSomega(m)) of the 2D-IR spectrum. The CLSomega(m) is obtained by taking slices through the 2D spectrum parallel to the detection frequency axis (omega(m)). Each slice is a spectrum. The slope of the line connecting the frequencies of the maxima of the sliced spectra is the CLSomega(m). The change in slope of the CLSomega(m) as a function of time is directly related to the FFCF and can be used to obtain the complete FFCF. CLSomega(m) is immune to line shape distortions caused by destructive interference between bands arising from vibrational echo emission, from the 0-1 vibrational transition (positive), and from the 1-2 vibrational transition (negative) in the 2D-IR spectrum. The immunity to the destructive interference enables the CLSomega(m) method to compare different parts of the bands as well as comparing the 0-1 and 1-2 bands. Also, line shape distortions caused by solvent background absorption and finite pulse durations do not affect the determination of the FFCF with the CLSomega(m) method. The CLSomega(m) can also provide information on the cross correlation between frequency fluctuations of the 0-1 and 1-2 vibrational transitions.

Two-dimensional infrared spectroscopy of photoswitchable peptides

We present a detailed discussion of the complimentary fields of the application of two-dimensional infrared (2D-IR) spectroscopy in comparison with two-dimensional nuclear magnetic resonance (2D-NMR) spectroscopy. Transient 2D-IR (T2D-IR) spectroscopy of nonequilibrium ensembles is probably one of the most promising strengths of 2D-IR spectroscopy, as the possibilities of 2D-NMR spectroscopy are limited in this regime. T2D-IR spectroscopy uniquely combines ultrafast time resolution with microscopic structural resolution. In this article we summarize our recent efforts to investigate the ultrafast structural dynamics of small peptides, such as the unfolding of peptide secondary structure motifs. The work requires two ingredients: 2D-IR spectroscopy and the possibility of triggering a structural transition of a peptide on an ultrafast timescale using embedded or intrinsic photoswitches. Several photoswitches have been tested, and we discuss our progress in merging these two pathways of research.

Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed alpha-helix and extended beta and P(II) conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the P(II) conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is alpha(R). Polarizable CHARMM force field gives almost exclusively P(II) conformation and other force fields predict that both alpha-helical and extended (beta and P(II)) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the P(II) conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.