Infrared spectroscopy of the alanine dipeptide analog in liquid water with DFT-MD. Direct evidence for P(II)/beta conformations

The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins

Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.

Ion-specific binding of cations to the carboxylate and of anions to the amide of alanylalanine

Studies of ion-specific effects on oligopeptides have aided our understanding of Hofmeister effects on proteins, yet the use of different model peptides and different experimental sensitivities have led to conflicting conclusions. To resolve these controversies, we study a small model peptide, L-Alanyl-L-alanine (2Ala), carrying all fundamental chemical protein motifs: C-terminus, amide bond, and N-terminus. We elucidate the effect of GdmCl, LiCl, KCl, KI, and KSCN by combining dielectric relaxation, nuclear magnetic resonance (¹H-NMR), and (two-dimensional) infrared spectroscopy. Our dielectric results show that all ions reduce the rotational mobility of 2Ala, yet the magnitude of the reduction is larger for denaturing cations than for anions. The NMR chemical shifts of the amide group are particularly sensitive to denaturing anions, indicative of anion-amide interactions. Infrared experiments reveal that LiCl alters the spectral homogeneity and dynamics of the carboxylate, but not the amide group. Interaction of LiCl with the negatively charged pole of 2Ala, the COO^- group, can explain the marked cationic effect on dipolar rotation, while interaction of anions between the poles, at the amide, only weakly perturbs dipolar dynamics. As such, our results provide a unifying view on ions' preferential interaction sites at 2Ala and help rationalize Hofmeister effects on proteins.

Effect of the Hydration Shell on the Carbonyl Vibration in the Ala-Leu-Ala-Leu Peptide

The vibrational spectrum of the Ala-Leu-Ala-Leu peptide in solution, computed from first-principles simulations, shows a prominent band in the amide I region that is assigned to stretching of carbonyl groups. Close inspection reveals combined but slightly different contributions by the three carbonyl groups of the peptide. The shift in their exact vibrational signature is in agreement with the different probabilities of these groups to form hydrogen bonds with the solvent. The central carbonyl group has a hydrogen bond probability intermediate to the other two groups due to interchanges between different hydrogen-bonded states. Analysis of the interaction energies of individual water molecules with that group shows that shifts in its frequency are directly related to the interactions with the water molecules in the first hydration shell. The interaction strength is well correlated with the hydrogen bond distance and hydrogen bond angle, though there is no perfect match, allowing geometrical criteria for hydrogen bonds to be used as long as the sampling is sufficient to consider averages. The hydrogen bond state of a carbonyl group can therefore serve as an indicator of the solvent's effect on the vibrational frequency.

Anharmonic Infrared Spectroscopy through the Fourier Transform of Time Correlation Function Formalism in ONETEP

Density functional theory molecular dynamics (DFT-MD) provides an efficient framework for accurately computing several types of spectra. The major benefit of DFT-MD approaches lies in the ability to naturally take into account the effects of temperature and anharmonicity, without having to introduce any ad hoc or a posteriori corrections. Consequently, computational spectroscopy based on DFT-MD approaches plays a pivotal role in the understanding and assignment of experimental peaks and bands at finite temperature, particularly in the case of floppy molecules. Linear-scaling DFT methods can be used to study large and complex systems, such as peptides, DNA strands, amorphous solids, and molecules in solution. Here, we present the implementation of DFT-MD IR spectroscopy in the ONETEP linear-scaling code. In addition, two methods for partitioning the dipole moment within the ONETEP framework are presented. Dipole moment partitioning allows us to compute spectra of molecules in solution, which fully include the effects of the solvent, while at the same time removing the solvent contribution from the spectra.

Can far-IR action spectroscopy combined with BOMD simulations be conformation selective?

The combination of conformation selective far-IR/UV double resonance spectroscopy with Born–Oppenheimer molecular dynamics (BOMD) simulations is presented here for the structural characterization of the Ac-Phe-Pro-NH₂ peptide in the far-infrared spectral domain, i.e. for radiation below 800 cm⁻¹.

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.

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.

Modeling the infrared and circular dichroism spectroscopy of a bridged cyclic diamide

Density functional theory based molecular dynamics simulations are used to study the structure, infrared (IR) spectroscopy, circular dichroism (CD) spectroscopy, and coupling between the amide I vibrations of a bridged cyclic diamide in the gas phase and in aqueous solution. IR spectra computed via the dipole moment time correlation function show a large red-shift of 30 cm(-1) in the amide I vibration in solution compared to the gas phase, and are in good agreement with experiment. Conformationally averaged CD spectra computed using the CIS(D) method are highly sensitive to the structures used, and structures sampled in the aqueous phase simulation are required to obtain qualitatively correct CD spectra. Analysis of the coupling between the amide I modes shows that in the aqueous phase there is an increased localization of the vibrations on the individual peptide groups and a reduction in the mode coupling parameter compared to the gas phase. Overall, the results illustrate the significance of incorporating molecular dynamics in the simulation of IR and CD spectra.