N-Methylacetamide and Its Hydrogen-Bonded Water Molecules Are Vibrationally Coupled

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Intermolecular Protein-Water Coupling Impedes the Coupling Between the Amide A and Amide I Mode in Interfacial Proteins

The coupling between different vibrational modes in proteins is essential for chemical dynamics and biological functions and is linked to the propagation of conformational changes and pathways of allosteric communication. However, little is known about the influence of intermolecular protein-H2O coupling on the vibrational coupling between amide A (NH) and amide I (C═O) bands. Here, we investigate the NH/CO coupling strength in various peptides with different secondary structures at the lipid cell membrane/H2O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy (SFG-VS) in which a femtosecond infrared pump is used to excite the amide A band, and SFG-VS is used to probe transient spectral evolution in the amide A and amide I bands. Our results reveal that the NH/CO coupling strength strongly depends on the bandwidth of the amide I mode and the coupling of proteins with water molecules. A large extent of protein-water coupling significantly reduces the delocalization of the amide I mode along the peptide chain and impedes the NH/CO coupling strength. A large NH/CO coupling strength is found to show a strong correlation with the high energy transfer rate found in the light-harvesting proteins of green sulfur bacteria, which may understand the mechanism of energy transfer through a molecular system and assist in controlling vibrational energy transfer by engineering the molecular structures to achieve high energy transfer efficiency.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

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.

Ultrafast energy relaxation dynamics of amide I vibrations coupled with protein-bound water molecules

The influence of hydration water on the vibrational energy relaxation in a protein holds the key to understand ultrafast protein dynamics, but its detection is a major challenge. Here, we report measurements on the ultrafast vibrational dynamics of amide I vibrations of proteins at the lipid membrane/H₂O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy. We find that the relaxation time of the amide I mode shows a very strong dependence on the H₂O exposure, but not on the D₂O exposure. This observation indicates that the exposure of amide I bond to H₂O opens up a resonant relaxation channel and facilitates direct resonant vibrational energy transfer from the amide I mode to the H₂O bending mode. The protein backbone motions can thus be energetically coupled with protein-bound water molecules. Our findings highlight the influence of H₂O on the ultrafast structure dynamics of proteins.

Crystalline modification of a rare earth nucleating agent for isotactic polypropylene based on its self-assembly

In this paper, the crystalline modification of a rare earth nucleating agent (WBG) for isotactic polypropylene (PP) based on its supramolecular self-assembly was investigated by differential scanning calorimetry, wide-angle X-ray diffraction and polarized optical microscopy. In addition, the relationship between the self-assembly structure of the nucleating agent and the crystalline structure, as well as the possible reason for the self-assembly behaviour, was further studied. The structure evolution of WBG showed that the self-assembly structure changed from a needle-like structure to a dendritic structure with increase in the content of WBG. When the content of WBG exceeded a critical value (0.4 wt%), it self-assembled into a strip structure. This revealed that the structure evolution of WBG contributed to the K_β and the crystallization morphology of PP with different content of WBG. In addition, further studies implied that the behaviour of self-assembly was a liquid-solid transformation of WBG, followed by a liquid-liquid phase separation of molten isotactic PP and WBG. The formation of the self-assembly structure was based on the free molecules by hydrogen bond dissociation while being heated, followed by aggregation into another structure by hydrogen bond association while being cooled. Furthermore, self-assembly behaviour depends largely on the interaction between WBG themselves.

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.

Strategy for Modeling the Electrostatic Responses of the Spectroscopic Properties of Proteins

For better understanding and more efficient use of the spectroscopic probes (vibrational and NMR) of the local electrostatic situations inside proteins, appropriate modeling of the properties of those probes is essential. The present study is devoted to examining the strategy for constructing such models. A more well-founded derivation than the ones in previous studies is given in constructing the models. Theoretical analyses are conducted on two representative example cases related to proteins, i.e., the peptide group of the main chains and the CO and NO ligands to the Fe²⁺ ion of heme, with careful treatment of the behavior of electrons in the electrostatic responses and with verification of consistency with observable quantities. It is shown that, for the stretching frequencies and NMR chemical shifts, it is possible to construct reasonable electrostatic interaction models that encompass the situations of hydration and uniform electric field environment and thus are applicable also to the cases of nonuniform electrostatic situations, which are highly expected for inside of proteins.

The Role of Backbone Hydration of Poly(N-isopropyl acrylamide) Across the Volume Phase Transition Compared to its Monomer

Thermo-responsive polymers undergo a reversible coil-to-globule transition in water after which the chains collapse and aggregate into bigger globules when passing to above its lower critical solution temperature (LCST). The hydrogen bonding with the amide groups in the side chains has to be contrasted with the hydration interaction of the hydrophobic main-chain hydrocarbons. In the present investigation we study molecular changes in the polymer poly(N-isopropyl acrylamide) (PNIPAM) and in its monomer N-isopropyl acrylamide (NIPAM) in solution across the LCST transition. Employing Fourier-transform infrared spectroscopy we probe changes in conformation and hydrogen bonding. We observe a nearly discontinuous shift of the peak frequencies and areas of vibrational bands across the LCST transition for PNIPAM whereas NIPAM exhibits a continuous linear change with temperature. This supports the crucial role of the polymer backbone with respect to hydration changes in the amide group in combination with cooperative interactions of bound water along the backbone chain.

Tracking the amyloidogenic core of IAPP amyloid fibrils: Insights from micro-Raman spectroscopy

Human islet amyloid polypeptide (hIAPP) is the major protein component of extracellular amyloid deposits, located in the islets of Langerhans, a hallmark of type II diabetes. The underlying mechanisms of IAPP aggregation have not yet been clearly defined, although the highly amyloidogenic sequence of the protein has been extensively studied. Several segments have been highlighted as aggregation-prone regions (APRs), with much attention focused on the central 8-17 and 20-29 stretches. In this work, we employ micro-Raman spectroscopy to identify specific regions that are contributing to or are excluded from the amyloidogenic core of IAPP amyloid fibrils. Our results demonstrate that both the N-terminal region containing a conserved disulfide bond between Cys residues at positions 2 and 7, and the C-terminal region containing the only Tyr residue are excluded from the amyloid core. Finally, by performing detailed aggregation assays and molecular dynamics simulations on a number of IAPP variants, we demonstrate that point mutations within the central APRs contribute to the reduction of the overall amyloidogenic potential of the protein but do not completely abolish the formation of IAPP amyloid fibrils.

Local Mode Analysis: Decoding IR Spectra by Visualizing Molecular Details

Integration of experimental and computational approaches to investigate chemical reactions in proteins has proven to be very successful. Experimentally, time-resolved FTIR difference-spectroscopy monitors chemical reactions at atomic detail. To decode detailed structural information encoded in IR spectra, QM/MM calculations are performed. Here, we present a novel method which we call local mode analysis (LMA) for calculating IR spectra and assigning spectral IR-bands on the basis of movements of nuclei and partial charges from just a single QM/MM trajectory. Through LMA the decoding of IR spectra no longer requires several simulations or optimizations. The novel approach correlates the motions of atoms of a single simulation with the corresponding IR bands and provides direct access to the structural information encoded in IR spectra. Either the contributions of a particular atom or atom group to the complete IR spectrum of the molecule are visualized, or an IR-band is selected to visualize the corresponding structural motions. Thus, LMA decodes the detailed information contained in IR spectra and provides an intuitive approach for structural biologists and biochemists. The unique feature of LMA is the bidirectional analysis connecting structural details to spectral features and vice versa spectral features to molecular motions.

Solvent-dependent dynamics of hydrogen bonding structure 5-(methylthio)-1, 3, 4-thiadiazole-2(3H)-thione as determined by Raman spectroscopy and theoretical calculation

The vibration spectra of 5-(methylthio)-1,3,4-thiadiazole-2(3H)-thione (MTTN) in acetonitrile (CH₃CN), methanol (CH₃OH) and water (H₂O) solvents were collected and evaluated via deuterium isotopic substitution Raman spectroscopic experiments. These experiments were combined with the quantum chemical theoretical calculations using the PCM solvent model and normal mode analysis. The results confirmed that the MTTN in CH₃CN, CH₃OH and H₂O have hydrogen bonding (H-bonding) MTTN(solvent)_n clusters that produce significantly different Raman intensity patterns in different solvents. Combined with the normal Raman assignment, most resonance Raman spectra were assigned to the vibration modes of the H-bonding MTTN(CH₃CN), MTTN(CH₃OH)₃ and MTTN(H₂O)₃ clusters in CH₃CN, CH₃OH and H₂O. The theoretically-predicted frequencies and intensities in different surrounding environments enabled reliable assignments of Raman bands. The intermolecular>NH⋯O and >NH⋯N H-bonding interactions are key constituents of stable thione structures in MTTN. This underlines the significant structural differences of MTTN in CH₃CN, CH₃OH and H₂O. H-bonding perturbation of MTTN reveal important insights about the intermolecular excited state proton transfer (ESPT) reaction mechanisms in the Franck-Condon region structural dynamics of the thione→thiol tautomer in CH₃OH and H₂O.

Polarized Raman Spectroscopy for Determining the Orientation of di-D-phenylalanine Molecules in a Nanotube

Self-assembly of short peptides into nanostructures has become an important strategy for the bottom-up fabrication of nanomaterials. Significant interest to such peptide-based building blocks is due to the opportunity to control the structure and properties of well-structured nanotubes, nanofibrils, and hydrogels. X-ray crystallography and solution NMR, two major tools of structural biology, have significant limitations when applied to peptide nanotubes because of their non-crystalline structure and large weight. Polarized Raman spectroscopy was utilized for structural characterization of well-aligned D-Diphenylalanine nanotubes. The orientation of selected chemical groups relative to the main axis of the nanotube was determined. Specifically, the C-N bond of CNH₃⁺groups is oriented parallel to the nanotube axis, the peptides' carbonyl groups are tilted at approximately 54° from the axis and the COO^- groups run perpendicular to the axis. The determined orientation of chemical groups allowed the understanding of the orientation of D-diphenylalanine molecule that is consistent with its equilibrium conformation. The obtained data indicate that there is only one orientation of D-diphenylalanine molecules with respect to the nanotube main axis.

A polarizable QM/MM approach to the molecular dynamics of amide groups solvated in water

The infrared (IR) spectra of polypeptides are dominated by the so-called amide bands. Because they originate from the strongly polar and polarizable amide groups (AGs) making up the backbone, their spectral positions sensitively depend on the local electric fields. Aiming at accurate computations of these IR spectra by molecular dynamics (MD) simulations, which derive atomic forces from a hybrid quantum and molecular mechanics (QM/MM) Hamiltonian, here we consider the effects of solvation in bulk liquid water on the amide bands of the AG model compound N-methyl-acetamide (NMA). As QM approach to NMA we choose grid-based density functional theory (DFT). For the surrounding MM water, we develop, largely based on computations, a polarizable molecular mechanics (PMM) model potential called GP6P, which features six Gaussian electrostatic sources (one induced dipole, five static partial charge distributions) and, therefore, avoids spurious distortions of the DFT electron density in hybrid DFT/PMM simulations. Bulk liquid GP6P is shown to have favorable properties at the thermodynamic conditions of the parameterization and beyond. Lennard-Jones (LJ) parameters of the DFT fragment NMA are optimized by comparing radial distribution functions in the surrounding GP6P liquid with reference data obtained from a "first-principles" DFT-MD simulation. Finally, IR spectra of NMA in GP6P water are calculated from extended DFT/PMM-MD trajectories, in which the NMA is treated by three different DFT functionals (BP, BLYP, B3LYP). Method-specific frequency scaling factors are derived from DFT-MD simulations of isolated NMA. The DFT/PMM-MD simulations with GP6P and with the optimized LJ parameters then excellently predict the effects of aqueous solvation and deuteration observed in the IR spectra of NMA. As a result, the methods required to accurately compute such spectra by DFT/PMM-MD also for larger peptides in aqueous solution are now at hand.

Infrared Spectroscopy of N-Methylacetamide Revisited by ab Initio Molecular Dynamics Simulations

The density functional theory based molecular dynamics simulation method ("Car-Parrinello") was applied in a numerical study of the electronic properties, hydrogen bonding, and infrared spectroscopy of the trans and cis isomer of N-methylacetamide in aqueous solution. A detailed analysis of the electronic structure of the solvated molecules, in terms of localized Wannier functions and Born atomic charges, is presented. Two schemes for the computation of the solute infrared absorption spectrum are investigated:  In the first method the spectrum is determined by Fourier transforming the time correlation function of the solute dipole as determined from the Wannier function analysis. The second method uses instead the molecular current-current correlation function computed from the Born charges and atomic velocities. The resulting spectral properties of trans- and cis-NMA are carefully compared to each other and to experimental results. We find that the two solvated isomers can be clearly distinguished by their infrared spectral profile in the 1000-2000 cm(-)(1) range.

Amide I Vibrational Properties Affected by Hydrogen Bonding Out-of-Plane of the Peptide Group

The amide I vibrational properties of a peptide-water complex in various intermolecular configurations are analyzed theoretically to see whether a water molecule with a weak out-of-plane hydrogen bond really induces a large low-frequency shift. It is shown that the frequency shift strongly depends on the C═O···H angle, with a larger low-frequency shift as the C═O···H becomes more bent, suggesting that the so-called hydrated helix with a rather low amide I frequency has an additional water molecule located out-of-plane of the peptide group as compared with a typical α-helix. The infrared intensity also depends on the angular position of water. A new model parameter set (that can be combined with molecular dynamics) is developed for a more correct representation of the hydration-induced frequency shift. The question regarding the scalar and vectorial nature of the molecular properties related to the frequency shift is also discussed.

Vibronic structure in the far-UV electronic circular dichroism spectra of proteins

The Franck–Condon effect is considered and the vibrational structure of the π_nbπ* transition of the peptide backbone is incorporated into matrix method calculations of the electronic circular dichroism (CD) spectra of proteins in the far-ultraviolet. We employ the state-averaged CASPT2 method to calculate the ground and π_nbπ* excited state geometries and frequencies of N-methylacetamide (NMA), which represents the peptide chromophore. The results of these calculations are used to incorporate vibronic levels of the excited states into the matrix method calculation. The CD spectra of a set of 49 proteins, comprising a range of structural types, are calculated to assess the influence of the vibrational structure. The calculated spectra of α-helical proteins are better resolved using the vibronic parameters and correlation between the experimental and the calculated intensity of less regular β structure proteins improves over most wavelengths in the far-UV. No obvious improvement is observed in the calculated spectra of regular β-sheet proteins. Our high-level ab initio calculations of the vibronic structure of the π_nbπ* transition in NMA have provided some further insight into the physical origins of the nature of protein CD spectra in the far-UV.

Hydration Effect on Amide I Infrared Bands in Water: An Interpretation Based on an Interaction Energy Decomposition Scheme

The sensitivity of some infrared bands to the local environment can be exploited to shed light on the structure and the dynamics of biological systems. In particular, the amide I band, which is specifically related to vibrations within the peptide bonds, can give information on the ternary structure of proteins, and can be used as a probe of energy transfer. In this work, we propose a model to quantitatively interpret the frequency shift on the amide I band of a model peptide induced by the formation of hydrogen bonds in the first solvation shell. This method allows us to analyze to what extent the electrostatic interaction, electronic polarization and charge transfer affect the position of the amide I band. The impact of the anharmoniticy of the pontential energy surface on the hydration induced shift is elucidated as well.

Polarized Raman Spectroscopy of Aligned Insulin Fibrils

Amyloid fibrils are associated with many neurodegenerative diseases. The application of conventional techniques of structural biology, X-ray crystallography and solution NMR, for fibril characterization is limited because of the non-crystalline and insoluble nature of the fibrils. Here, polarized Raman spectroscopy was used to determine the orientation of selected chemical groups in aligned insulin fibrils, specifically of peptide carbonyls. The methodology is solely based on the measurement of the change in Raman scattered intensity as a function of the angle between the incident laser polarization and the aligned fibrils. The order parameters 〈 P2 〉 and 〈 P4 〉 of the orientation distribution function were obtained, and the most probable distribution of C=O group orientation was calculated. The results indicate that the peptides' carbonyl groups are oriented at an angle of 13±5° from the fibril axis, which is in consistent with previously reported qualitative descriptions of an almost parallel orientation of the C=O groups relative to the main fibril axis.

pH-sensitive transformation of the peptidic bolaamphiphile self-assembly: exploitation for the pH-triggered chemical reaction

Control of the macroscopic self-assembled structure of the amphiphilic molecule has been a challenging issue in micro/nanotechnologies. In this study, the microtubular self-assembly of a novel peptidic bolaamphiphile, bis(N-α-amido-glycylglycine)-1,10-decene dicarboxylate, which undergoes reversible structural transformation between microtubes and precipitates as a function of pH, was exploited for pH-triggered chemical release. At neutral and basic conditions above a pH of 6, the peptidic bolaamphiphilic molecule self-assembled to form tubular structures several hundreds of micrometers in length. When the solution became acidic below a pH of 4, the tubular assembly disintegrated to form aggregated precipitates. The reversible transformation of precipitate to microtube was achieved by raising the pH above 6. From the Raman spectroscopy results, it was revealed that the hydrogen bonds of the amide group and carboxylate were enhanced under the acidic conditions. These variation of the hydrogen bonds resulted in precipitation of peptidic bolaamphiphilic molecules while rupturing the microtubular structure. The pH-sensitive microscopic structural transformation was exploited for release of a chemical in which the pH-triggered release of a model chemical, a fluorescence dye of ANS, was demonstrated. The ANS dye was released gradually with decreasing pH, which suggests gradual disintegration of the microtubular self-assembly. Furthermore, this pH-triggered release of a chemical was exploited for the chemical reaction of gold ion reduction to produce solid clusters. This study demonstrated the reversible transformation of peptidic bolaamphiphile and its application as a pH-sensitive host matrix.

Observation of persistent α-helical content and discrete types of backbone disorder during a molten globule to ordered peptide transition via deep-UV resonance Raman spectroscopy

The molten globule state can aide in the folding of a protein to a functional structure and is loosely defined as an increase in structural disorder with conservation of the ensemble secondary structure content. Simultaneous observation of persistent secondary structure content with increased disorder has remained experimentally problematic. As a consequence, modeling how the molten globule state remains stable and how it facilitates proper folding remains difficult due to a lack of amenable spectroscopic techniques to characterize this class of partially unfolded proteins. Previously, deep-UV resonance Raman (dUVRR) spectroscopy has proven useful in the resolution of global and local structural fluctuations in the secondary structure of proteins. In this work, dUVRR was employed to study the molten globule to ordered transition of a model four-helix bundle protein, HP7. Both the average ensemble secondary structure and types of local disorder were monitored, without perturbation of the solvent, pH, or temperature. The molten globule to ordered transition is induced by stepwise coordination of two heme molecules. Persistent dUVRR spectral features in the amide III region at 1295-1301 and 1335-1338 cm-1 confirm previous observations that HP7 remains predominantly helical in the molten globule versus the fully ordered state. Additionally, these spectra represent the first demonstration of conserved helical content in a molten globule protein. With successive heme binding significant losses are observed in the spectral intensity of the amide III3 and S regions (1230-1260 and 1390 cm-1, respectively), which are known to be sensitive to local disorder. These observations indicate that there is a decrease in the structural populations able to explore various extended conformations, with successive heme binding events. DUVRR spectra indicate that the first heme coordination between two helical segments diminishes exploration of more elongated backbone structural conformations in the inter-helical regions. A second heme coordination by the remaining two helices further restricts protein motion.

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.

Structural analysis of unfolded peptides by Raman spectroscopy

Raman spectroscopy has positioned itself as an invaluable tool in the study of complex biological systems, consistently being used to obtain information illustrating a vast array of fundamental properties. Of primary interest, with respect to the focus of this chapter, are conformational changes of peptide backbones. For short peptides to larger biological systems this understanding can be extended to local hydrogen bonding interactions and the probing of other structural or organizational properties. With regard to unfolded peptides Raman spectroscopy can be used as a technique complementary to infrared (IR) and vibrational circular dichroism (VCD) spectroscopy. This chapter describes how high quality polarized Raman spectra of peptide can be recorded with a Raman microspectrometer and how the structure sensitive amide I band profiles of isotropic and anisotropic Raman scattering can be analyzed in conjunction with the respective IR and VCD profiles to obtain conformational distributions of short unfolded peptides.

UV resonance Raman finds peptide bond-Arg side chain electronic interactions

We measured the UV resonance Raman excitation profiles and Raman depolarization ratios of the arginine (Arg) vibrations of the amino acid monomer as well as Arg in the 21-residue predominantly alanine peptide AAAAA(AAARA)(3)A (AP) between 194 and 218 nm. Excitation within the π → π* peptide bond electronic transitions result in UVRR spectra dominated by amide peptide bond vibrations. The Raman cross sections and excitation profiles indicate that the Arg side chain electronic transitions mix with the AP peptide bond electronic transitions. The Arg Raman bands in AP exhibit Raman excitation profiles similar to those of the amide bands in AP which are conformation specific. These Arg excitation profiles distinctly differ from the Arg monomer. The Raman depolarization ratios of Arg in monomeric solution are quite simple with ρ = 0.33 indicating enhancement by a single electronic transition. In contrast, we see very complex depolarization ratios of Arg in AP that indicate that the Arg residues are resonance enhanced by multiple electronic transitions.

Theoretical spectroscopy of floppy peptides at room temperature. A DFTMD perspective: gas and aqueous phase

Theoretical spectroscopy is mandatory for a precise understanding and assignment of experimental spectra recorded at finite temperature. We review here room temperature DFT-based molecular dynamics simulations for the purpose of interpreting finite temperature infrared spectra of peptides of increasing size and complexity, in terms of temperature-dependent conformational dynamics and flexibility, and vibrational anharmonicities (potential energy surface anharmonicities, vibrational mode couplings and dipole anharmonicities). We take examples from our research projects in order to illustrate the main key-points and strengths of dynamical spectra modeling in that context. The calculations are presented in relation to room temperature gas phase IR-MPD experiments and room temperature liquid phase IR absorption experiments. These illustrations of floppy polypeptides have been chosen in order to convey the following ideas: temperature-dependent spectra modeling is pivotal for a precise understanding of gas phase spectra recorded at room temperature, including conformational dynamics and vibrational anharmonicities; harmonic spectroscopy (as commonly performed in the literature) can be misleading and even erroneous for a proper interpretation of spectra recorded at finite temperature; taking into account vibrational anharmonicities is pivotal for a proper interplay between theory and experiments; amide I-III bands are not necessarily the most relevant fingerprints for unraveling the local structures of peptides and more complex systems; liquid phase simulations have unraveled relationships between the zwitterionic properties of the peptide bonds and infrared signatures. The review presents a state-of-the-art account of the domain and offers perspectives and new developments for future still more challenging applications.

Structures and bonding on a colloidal silver surface of the various length carboxyl terminal fragments of bombesin

Raman (RS) and surface-enhanced Raman scattering spectra (SERS) were measured for various length carboxyl terminal fragments (X-14 of amino acid sequence) of bombesin ( BN): BN13-14, BN12-14, BN11-14, BN10-14, BN9-14, and BN8-14 in silver colloidal solutions. Density functional theory (DFT) calculations of Raman wavenumbers and intensities with extended basis sets (B3LYP/6-31++G**) were performed with the aim of providing the definitive band allocations to the normal coordinates. The proposed band assignment is consistent with the assignment for similar compounds reported in the literature. The nonadsorbed and adsorbed molecular structures were deducted by detailed spectral analysis of the RS and SERS spectra, respectively. This analysis also allowed us to propose the particular surface geometry and orientation of these peptides on silver surface, and their specific interaction with the surface. For example, a SERS spectrum of BN8-14 indicates that the interaction of a thioether atom and Trp8 with the silver surface is favorable and may dictate the orientation and conformation of adsorbed peptide. One of the most prominent and common features in all of the fragments' SERS spectra is a approximately 692 cm (-1) band due to nu(C-S) accompanied by two or three bands of different C-S conformers for all, except BN8-14, which suggests that all of the above-mentioned compounds adsorb on the silver surface through the thioether atom and that the attachment of Trp8 produces limitation in a number of possible C-S conformers adopted on this surface. Our results also show clearly that His12 and CO do not interact with the colloid surface, which supports our earlier results.

UV resonance raman investigation of electronic transitions in alpha-helical and polyproline II-like conformations

UV resonance Raman (UVRR) excitation profiles and Raman depolarization ratios were measured for a 21-residue predominantly alanine peptide, AAAAA(AAARA) 3A (AP), excited between 194 and 218 nm. Excitation within the pi-->pi* electronic transitions of the amide group results in UVRR spectra dominated by amide vibrations. The Raman cross sections and excitation profiles provide information about the nature of the electronic transitions of the alpha-helix and polyproline II (PPII)-like peptide conformations. AP is known to be predominantly alpha-helical at low temperatures and to take on a PPII helix-like conformation at high temperatures. The PPII-like and alpha-helix conformations show distinctly different Raman excitation profiles. The PPII-like conformation cross sections are approximately twice those of the alpha-helix. This is due to hypochromism that results from excitonic interactions between the NV 1 transition of one amide group with higher energy electronic transitions of other amide groups, which decreases the alpha-helical NV 1 (pi-->pi*) oscillator strengths. Excitation profiles of the alpha-helix and PPII-like conformations indicate that the highest signal-to-noise Raman spectra of alpha-helix and PPII-like conformations are obtained at excitation wavelengths of 194 and 198 nm, respectively. We also see evidence of at least two electronic transitions underlying the Raman excitation profiles of both the alpha-helical and the PPII-like conformations. In addition to the well-known approximately 190 nm pi-->pi* transitions, the Raman excitation profiles and Raman depolarization ratio measurements show features between 205-207 nm, which in the alpha-helix likely results from the parallel excitonic component. The PPII-like helix appears to also undergo excitonic splitting of its pi-->pi* transition which leads to a 207 nm feature.

Dependence of amide vibrations on hydrogen bonding

The effect of hydrogen bonding on the amide group vibrational spectra has traditionally been rationalized by invoking a resonance model where hydrogen bonding impacts the amide functional group by stabilizing its [(-)O-C=NH (+)] structure over the [O=C-NH] structure. However, Triggs and Valentini's UV-Raman study of solvation and hydrogen bonding effects on epsilon-caprolactum, N, N-dimethylacetamide (DMA), and N-methylacetamide (NMA) ( Triggs, N. E.; Valentini, J. J. J. Phys. Chem. 1992, 96, 6922-6931) casts doubt on the validity of this model by demonstrating that, contrary to the resonance model prediction, carbonyl hydrogen bonding does not impact the AmII' frequency of DMA. In this study, we utilize density functional theory (DFT) calculations to examine the impact of hydrogen bonding on the C=O and N-H functional groups of NMA, which is typically used as a simple model of the peptide bond. Our calculations indicate that, as expected, the hydrogen bonding frequency dependence of the AmI vibration predominantly derives from the C=O group, whereas the hydrogen bonding frequency dependence of the AmII vibration primarily derives from N-H hydrogen bonding. In contrast, the hydrogen bonding dependence of the conformation-sensitive AmIII band derives equally from both C=O and N-H groups and thus, is equally responsive to hydrogen bonding at the C=O or N-H site. Our work shows that a clear understanding of the normal mode composition of the amide vibrations is crucial for an accurate interpretation of the hydrogen bonding dependence of amide vibrational frequencies.

Time-Domain Calculations of the Infrared and Polarized Raman Spectra of Tetraalanine in Aqueous Solution

The IR and polarized (isotropic and anisotropic) Raman spectra are calculated for the amide I band of tetraalanine ((Ala)4) in aqueous solution by using a time-domain computational method, which includes the effects of the diagonal frequency modulations (of individual peptide groups), the off-diagonal (interpeptide) vibrational couplings, and structural dynamics. It is shown that the calculated band profiles, especially the existence of a large negative noncoincidence effect (i.e., large frequency separations between the IR, isotropic Raman, and anisotropic Raman bands, with the isotropic Raman being higher in frequency), are in reasonable agreement with the experimental results. This negative noncoincidence effect derives from two conditions: the positive coupling between the amide I vibrations of peptide groups and the angle larger than 90 degrees between the transition dipoles of the coupled vibrations. This result means that the dynamically changing structures mainly in the polyproline II and beta-type conformations containing some repeated interconversions obtained from the molecular dynamics calculation are consistent with the existence of a large negative noncoincidence effect, as far as the structures satisfy the above two conditions. It is also shown that the electric fields from solvent water molecules induce larger frequency shifts than those of intrachain interactions, with rapid underdamped oscillatory modulations ( approximately 43 fs) due to the librational motions of water molecules that give rise to motional narrowing effect on the spectra. The reason for the difference from the behavior seen for the O-H stretching mode of liquid water is discussed. The time-domain analysis of the mode identity shows that the system proceeds halfway to complete mode mixing with a similar time scale ( approximately 60 fs), suggesting the importance of the nonadiabatic effect, which is included in a natural way in the present computational method.

Vibrational Spectral Simulation for Peptides of Mixed Secondary Structure: Method Comparisons with the Trpzip Model Hairpin

Infrared absorption and vibrational circular dichroism (IR and VCD) spectra of model fragments of TrpZip-style beta-hairpin structures are simulated using density functional theory (DFT) methods to estimate the influence of fragment size, end effects, conformational irregularities, peptide side chains, and solvent. Different fragmentation schemes, computing the strands and turn segments separately, were tested by varying the sizes of each and their respective overlaps. For suitably overlapping fragments, atomic property tensors were found to be reliably transferable, as tested by their ability to generate simulated spectra in good agreement with results from ab initio DFT computations for the entire peptide. This fragment approach significantly reduces computational times and opens up a wider range of systems that can be studied with a DFT-based approach as compared to previous methods based on uniform repeating sequences. However, vacuum calculations do not adequately represent the frequency dispersion of solvated molecules, and thus, some alternate strategies for solvation correction are explored for improving the simulation accuracy. Unlike for regular periodic secondary structure, the solvent significantly impacts the spectral shapes of hairpins, due to the different degrees of hydration of individual amide groups, which can be exposed to or shielded from water due to external vs internal hydrogen bonding. This is amplified by the shielding of selected amides from the solvent due to bulky side chains. The peptide plus solvent was structurally modeled with molecular dynamics methods, and then an electrostatic field-based parametrization correction was added to the force field and intensity tensors to compensate for the solvent dipolar field. The effect of the shielding and subsequent reordering of modes has a larger impact on VCD than IR band shapes.

How to Model Solvation of Peptides? Insights from a Quantum-mechanical and Molecular Dynamics Study ofN-Methylacetamide. 1. Geometries, Infrared, and Ultraviolet Spectra in Water

This paper represents the first part of a study of solvation in peptides using quantum-mechanical and classical approaches. In this study, the peptide is modeled as its simplest analogue, namely, N-methyl-acetamide, and the effects of the solvent (here water, and in the second part of the study, water and acetone) are introduced at three different levels, e.g., through a continuum description, using solute-solvent clusters, and using the same clusters embedded in an external continuum. In turn, the solute-solvent clusters have been obtained in two alternative ways, either by using QM optimization procedures or extracting a proper set of structures from MD simulations. In this part of the study, geometries, IR, and UV spectra are calculated in terms of the different solvation models, and the results are analyzed and compared to get insights about different aspects of solvation involving dynamic and static effects on one hand and bulk or specific interactions on the other hand.

Solvent Effects on IR and VCD Spectra of Helical Peptides: DFT-Based Static Spectral Simulations with Explicit Water

Simulations of IR and VCD spectra are carried out for model alpha-helical, 3(10)-helical, and 3(1)-helical (polyProII-like) oligopeptides, with up to 21 amide groups, and including explicit consideration of effects of directly hydrogen-bonded solvent (water). Parameters used were obtained from ab initio density functional theory (DFT) computations of force field, atomic polar and axial tensors for oligopeptides of 5 to 7 amides, whose structures were constrained in (phi,psi) to target the secondary structure type but otherwise fully optimized. By comparison with experimental data as well as with calculations for identical but isolated (gas phase) peptides, the computed effects of an inner shell of aqueous solvent on the vibrational spectra of helical oligopeptides are illustrated. The interaction with solvent causes significant frequency shifts of the amide bands, but only minor changes in the characteristic IR intensity distributions and splittings and the VCD band shapes. Better agreement with experimental band shapes is achieved for the alpha-helical amide I' (N-deuterated) VCD by inclusion of explicit solvent in the calculations. Some improvements are also observed in theoretical VCD predictions for 13C labeled alpha-helical peptides when solvated models are used. However, the qualitative isotopic splitting patterns are preserved and just shifted in frequency due to consistent, solvent independent interamide coupling constants. The critical match of experiment and theory for relative positions of transitions in peptides with specifically separated 13C=O labels, including and neglecting solvent, confirms the stability of the coupling interactions. Despite these solvation effects, the calculated VCD band shape of the amide I mode is shown to be a reliable conformational probe, since it remains basically insensitive to frequency shifts caused by environment. Thus theoretical VCD simulations, even vacuum calculations, are shown to provide useful spectral predictions for solution-phase peptides.