Discrimination between peptide 3(10)- and alpha-helices. Theoretical analysis of the impact of alpha-methyl substitution on experimental spectra

Increased accuracy of vibrational circular dichroism calculations for isotopically labeled helical peptides

Vibrational circular dichroism (VCD) spectra have been computed with qualitatively correct sign patterns for α-helical peptides using various methods, ranging from empirical models to ab initio quantum mechanical computations. However, some details, such as deuteration effects and isotope substitution shifts and sign patterns for the resultant amide I' band shape, have remained a predictive challenge. Fully optimized computations for a 25-residue Ala-rich peptide, including implicit solvent corrections and explicit side chains that experimentally stabilize these model helical peptides in water, have been carried out using density functional theory (DFT). These fully minimized structures show minor changes in the (ϕ,ψ) torsions at the termini and yield an extra negative band to the low energy side of the characteristic amide I' couplet VCD, in agreement with experiments. Additionally, these calculations give the right sign and relative intensity patterns, as compared to experimental results, for several 13C=O substituted variants. The differences from previously reported computations that used ideal helical structures and vacuum conditions imply that inclusion of distorted termini and solvent effects can have an impact on the final detailed spectral patterns. Inclusion of side chains in these calculations had very little effect on the computed amide I' IR and VCD. Tests of constrained geometries, varying dielectric, and different functionals indicate that each can affect the band shapes, particularly for the 12C=O components, but these aspects do not fully explain the difference from previous spectral simulations. Inclusion of long-range amide coupling, as obtained from DFT computation of the full structure, or transfer of parameters from a somewhat longer peptide model, rather than shorter model, seems to be more important for the final detailed band shape under isotopic substitution. However, these corrections can also induce other changes, suggesting that previously reported, limited calculations may have been qualitatively useful due to a balance of errors. This may also explain the success of simple empirical IR models.

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

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

In Situ Spectroscopic Detection of Large-Scale Reorientations of Transmembrane Helices During Influenza A M2 Channel Opening

Viroporins are small ion channels in membranes of enveloped viruses that play key roles during viral life cycles. To use viroporins as drug targets against viral infection requires in-depth mechanistic understanding and, with that, methods that enable investigations under in situ conditions. Here, we apply surface-enhanced infrared absorption (SEIRA) spectroscopy to Influenza A M2 reconstituted within a solid-supported membrane, to shed light on the mechanics of its viroporin function. M2 is a paradigm of pH-activated proton channels and controls the proton flux into the viral interior during viral infection. We use SEIRA to track the large-scale reorientation of M2's transmembrane α-helices in situ during pH-activated channel opening. We quantify this event as a helical tilt from 26° to 40° by correlating the experimental results with solid-state nuclear magnetic resonance-informed computational spectroscopy. This mechanical motion is impeded upon addition of the inhibitor rimantadine, giving a direct spectroscopic marker to test antiviral activity. The presented approach provides a spectroscopic tool to quantify large-scale structural changes and to track the function and inhibition of the growing number of viroporins from pathogenic viruses in future studies.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Optically Active Vibrational Spectroscopy of α-Aminoisobutyric Acid Foldamers in Organic Solvents and Phospholipid Bilayers

Helical α-aminoisobutyric acid (Aib) foldamers show great potential as devices for the communication of conformational information across phospholipid bilayers, but determining their conformation in bilayers remains a challenge. In the present study, Raman, Raman optical activity (ROA), infrared (IR) and vibrational circular dichroism (VCD) spectroscopies have been used to analyze the conformational preferences of Aib foldamers in solution and when interacting with bilayers. A 3₁₀ -helix marker band at 1665-1668 cm^-1 in Raman spectra was used to show that net helical content increased strongly with oligomer length. ROA and VCD spectra of chiral Aib foldamers provided the chiroptical signature for both left- and right-handed 3₁₀ -helices in organic solvents, with VCD establishing that foldamer screw-sense was preserved when the foldamers became embedded within bilayers. However, the population distribution between different secondary structures was perturbed by the chiral phospholipid. These studies indicate that ROA and VCD spectroscopies are valuable tools for the study of biomimetic structures, such as artificial signal transduction molecules, in phospholipid bilayers.

Monitoring the Orientational Changes of Alamethicin during Incorporation into Bilayer Lipid Membranes

Antimicrobial peptides (AMPs) are the first line of defense after contact of an infectious invader, for example, bacterium or virus, with a host and an integral part of the innate immune system of humans. Their broad spectrum of biological functions ranges from cell membrane disruption over facilitation of chemotaxis to interaction with membrane-bound or intracellular receptors, thus providing novel strategies to overcome bacterial resistances. Especially, the clarification of the mechanisms and dynamics of AMP incorporation into bacterial membranes is of high interest, and different mechanistic models are still under discussion. In this work, we studied the incorporation of the peptaibol alamethicin (ALM) into tethered bilayer lipid membranes on electrodes in combination with surface-enhanced infrared absorption (SEIRA) spectroscopy. This approach allows monitoring the spontaneous and potential-induced ion channel formation of ALM in situ. The complex incorporation kinetics revealed a multistep mechanism that points to peptide-peptide interactions prior to penetrating the membrane and adopting the transmembrane configuration. On the basis of the anisotropy of the backbone amide I and II infrared absorptions determined by density functional theory calculations, we employed a mathematical model to evaluate ALM reorientations monitored by SEIRA spectroscopy. Accordingly, ALM was found to adopt inclination angles of ca. 69°-78° and 21° in its interfacially adsorbed and transmembrane incorporated states, respectively. These orientations can be stabilized efficiently by the dipolar interaction with lipid head groups or by the application of a potential gradient. The presented potential-controlled mechanistic study suggests an N-terminal integration of ALM into membranes as monomers or parallel oligomers to form ion channels composed of parallel-oriented helices, whereas antiparallel oligomers are barred from intrusion.

Effective representation of amide III, II, I, and A modes on local vibrational modes: Analysis of ab initio quantum calculation results

The Hamiltonian matrix for the first excited vibrational states of a protein can be effectively represented by local vibrational modes constituting amide III, II, I, and A modes to simulate various vibrational spectra. Methods for obtaining the Hamiltonian matrix from ab initio quantum calculation results are discussed, where the methods consist of three steps: selection of local vibrational mode coordinates, calculation of a reduced Hessian matrix, and extraction of the Hamiltonian matrix from the Hessian matrix. We introduce several methods for each step. The methods were assessed based on the density functional theory calculation results of 24 oligopeptides with four different peptide lengths and six different secondary structures. The completeness of a Hamiltonian matrix represented in the reduced local mode space is improved by adopting a specific atom group for each amide mode and reducing the effect of ignored local modes. The calculation results are also compared to previous models using C=O stretching vibration and transition dipole couplings. We found that local electric transition dipole moments of the amide modes are mainly bound on the local peptide planes. Their direction and magnitude are well conserved except amide A modes, which show large variation. Contrary to amide I modes, the vibrational coupling constants of amide III, II, and A modes obtained by analysis of a dipeptide are not transferable to oligopeptides with the same secondary conformation because coupling constants are affected by the surrounding atomic environment.

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.

Comparison of vibrational circular dichroism instruments: development of a new dispersive VCD

A dispersive vibrational circular dichroism (VCD) instrument has been designed and optimized for the measurement of mid-infrared (MIR) bands such as the amide I and amide II vibrational modes of peptides and proteins. The major design considerations were to construct a compact VCD instrument for biological molecules, to increase signal-to-noise (S/N) ratio, to simultaneously collect and digitize the sample transmission and polarization modulation signals, and to digitally ratio them to yield a VCD spectrum. These were realized by assembling new components using design factors adapted from previous VCD instruments. A collection of spectra for peptides and proteins having different dominant secondary structures (alpha-helix, beta-sheet, and random coil) measured for identical samples under the same conditions showed that the new instrument had substantially improved S/N as compared with our previous dispersive VCD instrument. These instruments both provide protein VCD for the amide I that are comparable to or somewhat better than those measurable with commercial Fourier transform (FT) VCD instruments if just the amide I band in the spectra is obtained at modest resolution (8 cm(-1)) with the same total data collection time on each type of instrument.

FTIR spectroscopic studies on aggregation process of the beta-amyloid 11-28 fragment and its variants

Aggregation of Abeta peptides is a seminal event in Alzheimer's disease. Detailed understanding of the Abeta assembly process would facilitate the targeting and design of fibrillogenesis inhibitors. Here, conformational studies using FTIR spectroscopy are presented. As a model peptide, the 11-28 fragment of Abeta was used. This model peptide is known to contain the core region responsible for Abeta aggregation. The structural behavior of the peptide during aggregation provoked by the addition of water to Abeta(11-28) solution in hexafluoroisopropanol was compared with the properties of its variants corresponding to natural, clinically relevant mutants at positions 21-23 (A21G, E22K, E22G, E22Q and D23N). The results showed that the aggregation of the peptides proceeds via a helical intermediate, and it is possible that the formation of alpha-helical structures is preceded by creation of 3(10)-helix/3(10)-turn structures.

Alanine polypeptide structural fingerprints at room temperature: what can be gained from non-harmonic Car-Parrinello molecular dynamics simulations

Structural infrared fingerprints of neutral gas phase alanine peptides of increasing size and complexity (dipeptide, octapeptide, and beta-strand peptide) are characterized through DFT-based Car-Parrinello molecular dynamics simulations. Harmonic and nonharmonic vibrational signatures are calculated from the time correlation of the dipole moment of the gas phase peptide in a direct way (without any approximation) respectively from low temperature (20 K) and room temperature (300 K) molecular dynamics. Our main purpose is to answer the two following questions: (i) Is the direct inclusion of temperature for the calculation of infrared spectra mandatory for the comprehension of the vibrational signatures experimentally recorded at room temperature? (ii) To what extent is the amide I, II, and III domain sensitive enough to the local structure of the peptides, to provide vibrational signatures that can be definitely used to assess the peptide conformation at 300 K?

Conformational changes in cryptophane having C1-symmetry studied by vibrational circular dichroism

Vibrational circular dichroism (VCD) measurements and density functional theory (DFT) calculations were used to obtain the absolute configuration of a cryptophane molecule having C1-symmetry (labeled cryptophane-H). This molecule exhibits chiroptical properties different from those published for cryptophane-A having D3-symmetry [Brotin et al. J. Am. Chem. Soc. 2006, 128, 5533-5540]. In particular, we have shown that the conformation of the aliphatic linkers is very dependent on the solvent used and its ability to enter (CDCl3 solution) or not (C2D2Cl4 solution) in the cryptophane cavity. Calculations performed at the DFT (B3PW91/6-31G*) level establish, besides the absolute configuration, the preferential anti and gauche conformations of the aliphatic linkers of the chloroform@cryptophane-H complex and the empty cryptophane-H molecule, respectively. Polarimetric measurements performed in several solvents reflect also the change of conformation of the bridges upon guest encapsulation.

A density functional theory study of vibrational coupling in the amide I band of beta-sheet models

We report the first molecular orbital/density functional theory (DFT) calculations on the vibrational frequencies involved in the amide I band of completely geometrically optimized models for beta-sheet peptides based upon (up to 16) glycine residues. These calculations use the B3LYP/D95** level of DFT. The primary means of vibrational coupling occurs through H bond, rather than through space, interactions, which is consistent with a previous report on alpha-helical polyalanines and H-bonding chains of both formamides and 4-pyridones. We decoupled the C=O stretching vibrations using selected 14C substitutions to probe the coupling mechanism and to determine "natural" frequencies for individual 14C=Os. The intermolecular H-bonding interactions affect the geometries of the amide groups. Those near the center of H-bonding chains have long C=O bonds. The C=O bond lengths correlate with these "natural" frequencies, The frequencies obtained from the DFT calculations are generally more coupled, and the most intense are more red shifted than those calculated by transition dipole coupling (TDC). TDC inverts the order of the shifted frequencies compared to DFT in several cases.

Vibrational-exciton couplings for the amide I, II, III, and A modes of peptides

The couplings between all amide fundamentals and their overtones and combination vibrational states are calculated. Combined with the level energies reported previously (Hayashi, T.; Zhuang, W.; Mukamel, S. J. Phys. Chem. A 2005, 109, 9747), we obtain a complete effective vibrational Hamiltonian for the entire amide system. Couplings between neighboring peptide units are obtained using the anharmonic vibrational Hamiltonian of glycine dipeptide (GLDP) at the BPW91/6-31G(d,p) level. Electrostatic couplings between non-neighboring units are calculated by the fourth rank transition multipole coupling (TMC) expansion, including 1/R3 (dipole-dipole), 1/R4 (quadrupole-dipole), and 1/R5 (quadrupole-quadrupole and octapole-dipole) interactions. Exciton delocalization length and its variation with frequency in the various amide bands are calculated. The simulated infrared amide I and II absorptions and CD spectra of 24 residue alpha-helical motifs (SPE3) are in good agreement with experiment.

Through Hydrogen-Bond Vibrational Coupling in Hydrogen-Bonding Chains of 4-Pyridones with Implications for Peptide Amide I Absorptions: Density Functional Theory Compared with Transition Dipole Coupling

We present B3LYP/D95** calculations on the C=O and N-H couplings in H-bonded chains of 4-pyridones. 14C-substitutions are used to decouple various vibrations for purposes of illustration. The coupled C=O vibrations bear analogy to the amide I bands of proteins and peptides. The coupling of the C=O's occurs primarily via the cooperative H-bonds rather than transition dipole coupling (TDC), as demonstrated by the fact that (1) the couplings are greater than previously reported for similar studies on formamides despite the larger distance between the C=O's in the pyridone chains (TDC coupling decreases with distance) and (2) the red shifts (also greater than for formamides) can be attributed to the changes in the geometries (particularly the C=O bond lengths) of the individual 4-pyridones in the H-bonding chains induced by the H-bonds and resulting polarization of the monomers.

Two-Dimensional Circularly Polarized IR Photon Echo Spectroscopy of Polypeptides: Four-Wave-Mixing Optical Activity Measurement

A coherent two-dimensional (2D) optical spectroscopy utilizing circularly polarized (CP) beams, which was shown to be useful in studying molecular chirality in condensed phases, was theoretically proposed recently [Cho et al. J. Chem. Phys. 2003, 119, 7003]. A photon echo (PE) version of 2D optical activity spectroscopy is discussed in this paper. Considering various dipeptide and polypeptide systems, where the amide I local modes constitute the set of basis modes used to describe exciton and biexciton states as linear combinations of those basis modes, we present numerically simulated 2D circularly polarized IR PE spectra. It is shown that this novel spectroscopic method can provide additional information on the angles between the transition magnetic dipole and the transition electric dipole of two different vibrationally excited states, which are highly sensitive to the 3D structure and chirality of a given polypeptide. Also, a hierarchical relation of IR absorption, vibrational circular dichroism, 2D IR PE, and 2D CP-IR PE is discussed to show advantages of 2D optical activity spectroscopy in general.

Doubly resonant three-wave-mixing spectroscopy of a chiral coupled-chromophore system in solution: coherent two-dimensional optical activity spectroscopy

Theoretical descriptions of doubly resonant two-dimensional (2D) sum-frequency-generation (SFG) and difference-frequency-generation (DFG) spectroscopies of coupled-chromophore systems are presented. Despite that each electronic or vibrational chromophore is achiral, the interaction-induced chirality of a coupled multichromophore system in solution can be measured by using the doubly resonant 2D three-wave-mixing (3WM) spectroscopic method. An electronically coupled dimer, where each monomer is modeled as a simple two-level system, can have nonvanishing SFG (or DFG) properties, e.g., susceptibility in frequency domain or nonlinear response function in time domain, if the induced dipole vector of the dimer is not orthogonal to the vector product of the two monomer electronic transition dipole vectors. In order to demonstrate that these 2D 3WM spectroscopic methods can be used to determine the solution structure of a polypeptide, the authors carried out quantum chemistry calculations for an alanine dipeptide and obtained first- and second-order dipole derivatives associated with the amide I vibrational transitions of the dipeptide. It is shown that the numerically simulated 2D IR-IR SFG spectrum is highly sensitive to the dipeptide secondary structure and provides rich information on the one- and two-exciton states. It is believed that the theoretically proposed doubly resonant 2D 3WM spectroscopy, which can be considered to be an optical activity spectroscopy, will be of use in studying both structural and dynamical aspects of coupled multichromophore systems, such as proteins, nucleic acids, nanoparticle aggregates etc.

Two-Dimensional Infrared Spectral Signatures of 310- and α-Helical Peptides

Two-dimensional infrared (2D IR) spectra of Calpha-alkylated model octapeptides Z-(Aib)8-OtBu, Z-(Aib)5-L-Leu-(Aib)2-OMe, and Z-[L-(alphaMeVal)]8-OtBu have been measured in the amide I region to acquire 2D spectral signatures characteristic of 3(10)- and alpha-helical conformations. Phase-adjusted 2D absorptive spectra recorded with parallel polarizations are dominated by intense diagonal peaks, whereas 2D rephasing spectra obtained at the double-crossed polarization configuration reveal cross-peak patterns that are essential for structure determination. In CDCl3, all three peptides are of the 3(10)-helix conformation and exhibit a doublet cross-peak pattern. In 1,1,1,3,3,3-hexafluoroisopropanol, Z-[L-(alphaMeVal)]8-OtBu undergoes slow acidolysis and 3(10)-to-alpha-helix transition. In the course of this conformational change, its 2D rephasing spectrum evolves from an elongated doublet, characteristic of a distorted 3(10)-helix, to a multiple-peak pattern, after becoming an alpha-helix. The linear IR and 2D absorptive spectra are much less informative in discerning the structural changes. The experimental spectra are compared to simulations based on a vibrational exciton Hamiltonian model. The through-bond and through-space vibrational couplings are modeled by ab initio coupling maps and transition dipole interactions. The local amide I frequency is evaluated by a new approach that takes into account the effects of hydrogen-bond geometry and sites. The static diagonal and off-diagonal disorders are introduced into the Hamiltonian through statistical models to account for conformational fluctuations and inhomogeneous broadening. The sensitivity of cross-peak patterns to different helical conformations and the chain length dependence of the spectral features for short 3(10)- and alpha-helices are discussed.

Infrared and vibrational CD spectra of partially solvated alpha-helices: DFT-based simulations with explicit solvent

Theoretical simulations are used to investigate the effects of aqueous solvent on the vibrational spectra of model alpha-helices, which are only partly exposed to solvent to mimic alpha-helices in proteins. Infrared absorption (IR) and vibrational circular dichroism (VCD) amide I' spectra for 15-amide alanine alpha-helices are simulated using density functional theory (DFT) calculations combined with the property transfer method. The solvent is modeled by explicit water molecules hydrogen bonded to the solvated amide groups. Simulated spectra for two partially solvated model alpha-helices, one corresponding to a more exposed and the other to a more buried structure, are compared to the fully solvated and unsolvated (gas phase) simulations. The dependence of the amide I spectra on the orientation of the partially solvated helix with respect to the solvent and effects of solvation on the amide I' of 13C isotopically substituted alpha-helices are also investigated. The partial exposure to solvent causes significant broadening of the amide I' bands due to differences in the vibrational frequencies of the explicitly solvated and unsolvated amide groups. The different degree of partial solvation is reflected primarily in the frequency shifts of the unsolvated (buried) amide group vibrations. Depending on which side of the alpha-helix is exposed to solvent, the simulated IR band-shapes exhibit significant changes, from broad and relatively featureless to distinctly split into two maxima. The simulated amide I' VCD band-shapes for the partially solvated alpha-helices parallel the broadening of the IR and exhibit more sign variation, but generally preserve the sign pattern characteristic of the alpha-helical structures and are much less dependent on the alpha-helix orientation with respect to the solvent. The simulated amide I' IR spectra for the model peptides with explicitly hydrogen-bonded water are consistent with the experimental data for small alpha-helical proteins at very low temperatures, but overestimate the effects of solvent on the protein spectra at ambient temperatures, where the peptide-water hydrogen bonds are weakened by thermal motion.