Theoretical study of the Raman optical activity spectra of 3(10)-helical polypeptides

Raman optical activity (ROA) spectroscopy is a promising analytical method for studying the structure and conformation of polypeptides and proteins in solution. However, the structural information obtained from such vibrational spectra is only indirect and theoretical studies are often necessary to identify how the structure determines the observed spectra. One particular target is the identification and discrimination of different helical secondary structure elements. Herein, a theoretical investigation of the ROA spectra of a series of 3(10)-helical polypeptides is presented. In particular, the effect of chain length, C(α)-substitution pattern, the introduction of larger aliphatic side chains, and the variation of their conformation on the ROA spectra is studied. To extract general principles from these calculations, the positions, intensities, and shapes of the ROA bands are analyzed in terms of localized modes, which makes it possible to identify possible ROA signatures of 3(10) -helical structures, but also provides fundamental insight into the generation of ROA signals in complex polypeptides. Finally, the calculated spectra can be compared to the previously reported ROA spectrum of a specifically designed 3(10) -helical heptapeptide. This allows most of the features in the experimental spectrum to be assigned.     

Identifying Protein β‐Turns with Vibrational Raman Optical Activity

β‐turns belong to the most important secondary structure elements in proteins. On the basis of density functional calculations, vibrational Raman optical activity signatures of different types of β‐turns are established and compared as well as related to other signatures proposed in the literature earlier. Our findings indicate that there are much more characteristic ROA signals of β‐turns than have been hitherto suggested. These suggested signatures are, however, found to be valid for the most important type of β‐turns. Moreover, we compare the influence of different amino acid side chains on these signatures and investigate the discrimination of β‐turns from other secondary structure elements, namely α‐ and 3₁₀‐helices.     

Demonstration of the ring conformation in polyproline by the Raman optical activity

Raman and Raman optical activity (ROA) spectra of poly-L-proline were recorded in a wide frequency range and analyzed with respect to the proline side chain conformation. The analysis was based on comparison to ab initio simulations of spectral band positions and intensities. The presence of two conformer states of the five-member ring was found, approximately equally populated in the polypeptide. Additionally, Raman and ROA spectral shapes indicated that the peptide adopts the polyproline II helical conformation, in both aqueous and TFE solutions. The helix, however, is perturbed by fluctuations, which affects the vibrational coupling among amino acid residues and broadens the ROA bands. Contributions of the side and main peptide chains to the polyproline ROA intensities have comparable magnitudes. Thus understanding of the origins of both signals is important for determination of the peptide structure by ROA.     

Charge‐Induced Unzipping of Isolated Proteins to a Defined Secondary Structure

Here we present a combined experimental and theoretical study on the secondary structure of isolated proteins as a function of charge state. In infrared spectra of the proteins ubiquitin and cytochrome c, amide I (C=O stretch) and amide II (N–H bend) bands can be found at positions that are typical for condensed‐phase proteins. For high charge states a new band appears, substantially red‐shifted from the amide II band observed at lower charge states. The observations are interpreted in terms of Coulomb‐driven transitions in secondary structures from mostly helical to extended C₅‐type hydrogen‐bonded structures. Support for this interpretation comes from simple energy considerations as well as from quantum chemical calculations on model peptides. This transition in secondary structure is most likely universal for isolated proteins that occur in mass spectrometric experiments.     

MOVIPAC: Vibrational spectroscopy with a robust meta‐program for massively parallel standard and inverse calculations

We present the software package MO VI PAC for calculations of vibrational spectra, namely infrared, Raman, and Raman Optical Activity (ROA) spectra, in a massively parallelized fashion. MO VI PAC unites the latest versions of the programs SNF and AKIRA alongside with a range of helpful add‐ons to analyze and interpret the data obtained in the calculations. With its efficient parallelization and meta‐program design, MO VI PAC focuses in particular on the calculation of vibrational spectra of very large molecules containing on the order of a hundred atoms. For this purpose, it also offers different subsystem approaches such as Mode‐ and Intensity‐Tracking to selectively calculate specific features of the full spectrum. Furthermore, an approximation to the entire spectrum can be obtained using the Cartesian Tensor Transfer Method. We illustrate these capabilities using the example of a large π‐helix consisting of 20 (S)‐alanine residues. In particular, we investigate the ROA spectrum of this structure and compare it to the spectra of α‐ and 3₁₀‐helical analogs. © 2012 Wiley Periodicals, Inc.     

Raman Optical Activity (ROA) as a New Tool to Elucidate the Helical Structure of Poly(phenylacetylene)s

Poly(phenylacetylene)s are a family of helical polymers constituted by conjugated double bonds. Raman spectra of these polymers show a structural fingerprint of the polyene backbone which, in combination with its helical orientation, makes them good candidates to be studied by Raman optical activity (ROA). Four different well‐known poly(phenylacetylene)s adopting different scaffolds and ten different helical senses have been prepared. Raman and ROA spectra were recorded and allowed to establish ROA‐spectrum/helical‐sense relationships: a left/right‐handed orientation of the polyene backbone (M_helix/P_helix) produces a triplet of positive/negative ROA bands. Raman and ROA spectra of each polymer exhibited the same profile, and the sign of the ROA spectrum was opposite to the lowest‐energy electronic circular dichroism (ECD) band, indicating a resonance effect. Resonance ROA appears then as an indicator of the helical sense of poly(phenylacetylene)s, especially for those with an extra Cotton band in the ECD spectrum, where a wrong helical sense is assigned based on ECD, while ROA alerts of this misassignment.     

Raman Optical Activity of a Purely σ‐Bonded Helical Chromophore: (−)‐(M)‐σ‐[4]Helicene

The recent synthesis of enantiomerically pure (−)‐(M)‐σ‐[4]helicene has provided an archetype helical model system for vibrational optical activity, comparable to what π‐helicenes represent in the field of electronic optical activity. We present the first measurements and the first calculations of the Raman optical activity (ROA) of this interesting molecule. Observed and calculated ROA is large throughout the vibrational spectrum, in agreement with expectations, and spectacular effects, with Δ values close to 0.5%, occur in the 900‐cm⁻¹ region. Agreement between the experimental spectrum and the theoretical one, calculated with density‐functional theory for the vibrational part and Hartree‐Fock linear response theory for the molecular electronic tensors, is excellent, clearly the best that has been achieved to date in the field. This allows us to place confidence in the results of an analysis of Raman and ROA scattering generation in the molecule, obtained by a newly developed graphical procedure for extracting this kind of information from ab initio calculations. One finds that relative contributions made by carbon and hydrogen atoms can be comparable in size, but can also vary considerably, even between closely lying vibrations, and that, for most vibrations, the generation of ROA difference intensity is distributed rather differently than that of Raman intensity over the shape of the molecule. The sign of the ROA is, for the set of vibrations in the 900‐cm⁻¹ region, which we analyze in detail, determined by coupling terms between the two halves of the molecule, while Raman intensity is primarily generated within the two fragments, with coupling terms between them only adding to or substracting from it.     

Can Raman Optical Activity Separate Axial from Local Chirality? A Theoretical Study of Helical Deca‐Alanine

Motivated by experimental work on the distinction of protein secondary structure motifs by Raman optical activity (ROA) spectroscopy, we demonstrate using density functional theory that axial chirality in structures with different local chirality can be filtered out by ROA spectroscopy. To this purpose, two diastereomers of right-handed helical deca-alanine, the (all-S) and the (R,S,R,S,R,S,R,S,R,S) form, are compared. Furthermore, we suggest to interpret calculated ROA spectra of large molecules in terms of vibrational bands rather than individual peaks. This is due to the non-homogeneous effect of the harmonic approximation as well as of the chosen electronic structure method onto the vibrational frequencies, which in a dense region of many vibrations will strongly determine the shape of the spectrum. In addition, the calculated ROA spectrum of (all-S)-deca-alanine is compared to the experimental spectrum of poly-(L)-alanine in solution.     

Influence of the Helical Backbone in the Behavior of a Simple Model for Dimeric Coiled‐Coil Proteins

Using computer simulations as a tool for thought experiments, we investigate the influence of the helical backbone geometry in the association process and the final structures of a simple model which mimics parallel, two‐stranded coiled‐coil proteins. We define three types of helices: two of them have straight helical axes and 3.5 or 3.6 residues per helical turn; the third type presents a coiled helical axis, according to the canonical scheme defined by Crick. By using a Monte Carlo simulation algorithm, we observe that the three models exhibit different transition temperatures for the formation of the dimeric structure from two independent peptides, and a different behavior concerning the appearance of out‐of‐register structures. The energy minimized dimer structures present strong deviations from the correct association for straight helices with 3.6 residues/turn, especially for long peptides, deviations which are absent for the other two types when only the burial of hydrophobic residues is considered. A careful analysis of the energies for the out‐of‐register configurations and the contact maps reveals also differences between dimers resulting from the model with Crick parameterization and with 3.5 residues/turn. The results presented in this paper may be relevant for the design of simple models which use rigid α‐helices built from predicted elements of secondary structure.

Top views of the helical models used in this work.     

α‐Aminoxy Oligopeptides: Synthesis,Secondary Structure,and Cytotoxicity of a New Class of Anticancer Foldamers

α‐Aminoxy peptides are peptidomimetic foldamers with high proteolytic and conformational stability. To gain an improved synthetic access to α‐aminoxy oligopeptides we used a straightforward combination of solution‐ and solid‐phase‐supported methods and obtained oligomers that showed a remarkable anticancer activity against a panel of cancer cell lines. We solved the first X‐ray crystal structure of an α‐aminoxy peptide with multiple turns around the helical axis. The crystal structure revealed a right‐handed 2₈‐helical conformation with precisely two residues per turn and a helical pitch of 5.8 Å. By 2D ROESY experiments, molecular dynamics simulations, and CD spectroscopy we were able to identify the 2₈‐helix as the predominant conformation in organic solvents. In aqueous solution, the α‐aminoxy peptides exist in the 2₈‐helical conformation at acidic pH, but exhibit remarkable changes in the secondary structure with increasing pH. The most cytotoxic α‐aminoxy peptides have an increased propensity to take up a 2₈‐helical conformation in the presence of a model membrane. This indicates a correlation between the 2₈‐helical conformation and the membranolytic activity observed in mode of action studies, thereby providing novel insights in the folding properties and the biological activity of α‐aminoxy peptides.     

The Influence of the Amino Acid Side Chains on the Raman Optical Activity Spectra of Proteins

The Raman optical activity (ROA) spectra of proteins show distinct patterns arising from the secondary structure. It is generally believed that the spectral contributions of the side-chains largely cancel out because of their flexibility and the occurrence of many side-chains with different conformations. Yet, the influence of the side-chains on the ROA patterns assigned to different secondary structures is unknown. Here, the first systematic study of the influence of all amino acid side-chains on the ROA patterns is presented based on density functional theory (DFT) calculations of an extensive collection of peptide models that include many different side-chain and secondary structure conformations. It was shown that the contributions of the side-chains to a large extent average out with conformational flexibility. However, specific side-chain conformations can have significant contributions to the ROA patterns. It was also shown that α-helical structure is very sensitive to both the exact backbone conformation and the side-chain conformation. Side-chains with χ₁≈−60° generate ROA patterns alike those in experiment. Aromatic side-chains strongly influence the amide III ROA patterns. Because of the huge structural sensitivity of ROA, the spectral patterns of proteins arise from extensive conformational averaging of both the backbone and the side-chains. The averaging results in the fine spectral details and relative intensity differences observed in experimental spectra.     

Construction of a Triangle‐Shaped Trimer and a Tetrahedron Using an α‐Helix‐Inserted Circular Permutant of Cytochrome c555

Highly‐ordered protein structures have gained interest for future uses for biomaterials. Herein, we constructed a building block protein (BBP) by the circular permutation of the hyperthermostable Aquifex aeolicus cytochrome (cyt) c₅₅₅, and assembled BBP into a triangle‐shaped trimer and a tetrahedron. The angle of the intermolecular interactions of BBP was controlled by cleaving the domain‐swapping hinge loop of cyt c₅₅₅ and connecting the original N‐ and C‐terminal α‐helices with an α‐helical linker. We obtained BBP oligomers up to ≈40 mers, with a relatively large amount of trimers. According to the X‐ray crystallographic analysis of the BBP trimer, the N‐terminal region of one BBP molecule interacted intermolecularly with the C‐terminal region of another BBP molecule, resulting in a triangle‐shaped structure with an edge length of 68 Å. Additionally, four trimers assembled into a unique tetrahedron in the crystal. These results demonstrate that the circular permutation connecting the original N‐ and C‐terminal α‐helices with an α‐helical linker may be useful for constructing organized protein structures.     

Vibrational raman optical activity as a mean for revealing the helicity of oligosilanes: a quantum chemical investigation

Using theoretical simulations based on Hartree-Fock and density-functional theory calculations, the simulated vibrational Raman optical activity spectra of helical conformers of heptasilane are shown to present signatures sensitive to the helicity. These signatures are associated with collective wagging, twisting, and rocking motions. These simulated spectra have been obtained by combining analytical and numerical differentiation procedures to evaluate the geometry derivatives of the optical tensors entering into the expressions of the vibrational Raman optical activity intensities. From an investigation of basis set and electron correlation effects, it is shown that, like for local vibrations, diffuse functions are compulsory for evaluating the vibrational Raman optical activity intensities of collective vibrational motions.     

New Helical Folds in α‐Peptides with Alternating Chirality

In α‐peptides, the 8/10 helix is theoretically predicted to be energetically unstable and has not been experimentally observed so far. Based on our earlier studies on ‘helical induction’ and ‘hybrid helices’, we have adopted the ‘end‐capping’ strategy to induce the 8/10 helix in α‐peptides by using short α/β‐peptides. Thus, α‐peptides containing a regular string of α‐amino acids with alternating chirality were end capped by α/β‐peptides with 11/9‐helical motifs at the termini. Extensive NMR spectroscopy studies of these peptides revealed the presence of a hitherto unknown 8/10‐helical pattern; the H‐bonds in the shorter pseudorings were rather weak. The approach of using short helical motifs to induce new mixed helices in α‐peptides could provide avenues for more versatile design strategies.     

Length-Dependent Collective Vibrational Dynamics in Alpha-Helices

Functions of protein molecules in nature are closely associated with their well-defined three-dimensional structures and dynamics in body fluid. So far, many efforts have been made to reveal the relation of protein structure, dynamics, and function, but the structural origin of protein dynamics, especially for secondary structures as building blocks of conformation transition, is still ambiguous. Here we theoretically uncover the collective vibrations of elastic poly-alanine α-helices and find vibration patterns that are distinctively different over residue numbers ranging from 20 to 80. Contrary to the decreasing vibration magnitude from ends to the middle region for short helices, the vibration magnitude for long helices takes the minimum at approximately 1/5 of helix length from ends but reaches a peak at the center. Further analysis indicates that major vibrational modes of helical structures strongly depend on their lengths, where the twist mode dominates in the vibrations of short helices while the bend mode dominates the long ones analogous to an elastic Euler beam. The helix-coil transition pathway is also affected by the alternation of the first-order mode in helices with different lengths. The dynamic properties of the helical polypeptides are promising to be harnessed for de novo design of protein-based materials and artificial biomolecules in clinical treatments.     

Visualization and characterization of the infrared active amide I vibrations of proteins

To facilitate the analysis of frequency-structure correlations in the amide I vibrational spectroscopy of proteins, we investigate visualization methods and spatial correlation functions that describe delocalized vibrations of proteins and protein secondary structures. To study those vibrational modes revealed in infrared spectroscopy, we characterize frequency-dependent bright states obtained from doorway mode analysis. Our visualization methods pictorially color code amplitude and phase of each oscillator within the structure to reveal spatially varying patterns characteristic of excitations within sheets and helices. Spatial correlation functions in the amplitude and phase of amide I oscillators quantitatively address the extent of delocalization and the alpha helical and beta sheet character of these modes. Specifically, we investigate the vibrations of idealized antiparallel beta sheets and alpha helices and perform case studies on three proteins: concanavalin A, myoglobin, and ubiquitin.     

Engineering a β‐Helical d,l‐Peptide for Folding in Polar Media

β Helices—helices formed by alternating d,l ‐peptides and stabilized by β‐sheet hydrogen bonding—are found naturally in only a handful of highly hydrophobic peptides. This paper explores the scope of β‐helical structure by presenting the first design and biophysical characterization of a hydrophilic d,l ‐peptide, 1 , that forms a β helix in methanol. The design of 1 is based on the β‐hairpin/β helix—a new supersecondary that had been characterized previously only for hydrophobic peptides in nonpolar solvents. Incorporating polar residues in 1 provided solubility in methanol, in which the peptide adopts the expected β‐hairpin/β‐helical structure, as evidenced by CD, analytical ultracentrifugation (AUC), NMR spectroscopy, and NMR‐based structure calculations. Upon titration with water (at constant peptide concentration), the structure in methanol ( 1 m ) transitions cooperatively to an extended conformation ( 1 w ) resembling a cyclic β‐hairpin; observation of an isodichroic point in the solvent‐dependent CD spectra indicates that this transition is a two‐state process. In contrast, neither 1 m nor 1 w show cooperative thermal melting; instead, their structures appear intact at temperatures as high as 65 °C; this observation suggests that steric constraint is dominant in stabilizing these structures. Finally, the ¹H NMR CαH spectroscopic resonances of 1 m are downfield‐shifted with respect to random‐coil values, a hitherto unreported property for β helices that appears to be a general feature of these structures. These results show for the first time that an appropriately designed β‐helical peptide can fold stably in a polar solvent; furthermore, the structural and spectroscopic data reported should prove useful in the future design and characterization of water‐soluble β helices.     

Improved helix and kink characterization in membrane proteins allows evaluation of kink sequence predictors

Although the α-helical secondary structure of proteins is well-defined, the exact causes and structures of helical kinks are not. This is especially important for transmembrane (TM) helices of integral membrane proteins, many of which contain kinks providing functional diversity despite predominantly helical structure. We have developed a Monte Carlo method based algorithm, MC-HELAN, to determine helical axes alongside positions and angles of helical kinks. Analysis of all nonredundant high-resolution α-helical membrane protein structures (842 TM helices from 205 polypeptide chains) revealed kinks in 64% of TM helices, demonstrating that a significantly greater proportion of TM helices are kinked than those indicated by previous analyses. The residue proline is over-represented by a factor >5 if it is two or three residues C-terminal to a bend. Prolines also cause kinks with larger kink angles than other residues. However, only 33% of TM kinks are in proximity to a proline. Machine learning techniques were used to test for sequence-based predictors of kinks. Although kinks are somewhat predicted by sequence, kink formation appears to be driven predominantly by other factors. This study provides an improved view of the prevalence and architecture of kinks in helical membrane proteins and highlights the fundamental inaccuracy of the typical topological depiction of helical membrane proteins as series of ideal helices.     

New energy terms for reduced protein models implemented in an off‐lattice force field

Parameterization and test calculations of a reduced protein model with new energy terms are presented. The new energy terms retain the steric properties and the most significant degrees of freedom of protein side chains in an efficient way using only one to three virtual atoms per amino acid residue. The energy terms are implemented in a force field containing predefined secondary structure elements as constraints, electrostatic interaction terms, and a solvent‐accessible surface area term to include the effect of solvation. In the force field the main‐chain peptide units are modeled as electric dipoles, which have constant directions in α‐helices and β‐sheets and variable conformation‐dependent directions in loops. Protein secondary structures can be readily modeled using these dipole terms. Parameters of the force field were derived using a large set of experimental protein structures and refined by minimizing RMS errors between the experimental structures and structures generated using molecular dynamics simulations. The final average RMS error was 3.7 Å for the main‐chain virtual atoms (C_α atoms) and 4.2 Å for all virtual atoms for a test set of 10 proteins with 58–294 amino acid residues. The force field was further tested with a substantially larger test set of 608 proteins yielding somewhat lower accuracy. The fold recognition capabilities of the force field were also evaluated using a set of 27,814 misfolded decoy structures. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1229–1242, 2001     

Multiple Fused Anthracenes as Helical Polycyclic Aromatic Hydrocarbon Motif for Chiroptical Performance Enhancement

Polycyclic aromatic hydrocarbons consisting of three fused anthracene units were designed as new π‐conjugated compounds having helical structures. These expanded helicenes were synthesized by Pt‐catalyzed cycloisomerization of the corresponding ethynyl‐substituted precursors. The nonplanar and helical structure was confirmed by X‐ray analysis and DFT calculations, and the barrier to helical inversion was estimated to be 34 kJ mol^?1. The enantiomers of the diphenyl derivative were successfully resolved by chiral HPLC. Enantiopure samples showed good chiroptical performance in the CD (|Δ?| 1380 L mol^?1 cm^?1) and CPL (|g_lum| 0.013) spectra, and these values were considerably large for simple organic molecules. The unique chiroptical properties are discussed on the basis of the molecular structure and the electronic state with the aid of time‐dependent DFT calculations.