Theoretical prediction of the basic helix types in alpha,beta-hybrid peptides

This study provides a complete overview on all possible helical- folding patterns, their stabilities, and their detailed molecular structure in the novel foldamer class of alpha,beta-hybrid peptides on the basis of ab initio molecular orbital (MO) theory. The results indicate a considerable intrinsic potential of backbone folding. As found for other peptide foldamers, representatives of mixed or beta-helices are most stable in more apolar media, whereas polar environments favor the helices with the hydrogen bonds pointing in only one direction. The theoretical results confirm the hydrogen-bonding patterns found in the first experimental studies on these hybrid peptides. Selecting special backbone substitution patterns, the secondary structure potential of the alpha,beta-hybrid peptides could be of great importance for a rational peptide and protein design.     

Helix propensities of the amino acids measured in alanine-based peptides without helix-stabilizing side-chain interactions.

Helix propensities of the amino acids have been measured in alanine-based peptides in the absence of helix-stabilizing side-chain interactions. Fifty-eight peptides have been studied. A modified form of the Lifson-Roig theory for the helix-coil transition, which includes helix capping (Doig AJ, Chakrabartty A, Klingler TM, Baldwin RL, 1994, Biochemistry 33:3396-3403), was used to analyze the results. Substitutions were made at various positions of homologous helical peptides. Helix-capping interactions were found to contribute to helix stability, even when the substitution site was not at the end of the peptide. Analysis of our data with the original Lifson-Roig theory, which neglects capping effects, does not produce as good a fit to the experimental data as does analysis with the modified Lifson-Roig theory. At 0 degrees C, Ala is a strong helix former, Leu and Arg are helix-indifferent, and all other amino acids are helix breakers of varying severity. Because Ala has a small side chain that cannot interact significantly with other side chains, helix formation by Ala is stabilized predominantly by the backbone ("peptide H-bonds"). The implication for protein folding is that formation of peptide H-bonds can largely offset the unfavorable entropy change caused by fixing the peptide backbone. The helix propensities of most amino acids oppose folding; consequently, the majority of isolated helices derived from proteins are unstable, unless specific side-chain interactions stabilize them.     

Helix nucleation kinetics from molecular simulations in explicit solvent

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.     

Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein.

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Explicit-chain model of native-state hydrogen exchange: implications for event ordering and cooperativity in protein folding

Native-state hydrogen exchange experiments on several proteins have revealed partially unfolded conformations with diverse stabilities. These equilibrium observations have been used to support kinetic arguments that folding proceeds via a sequential "pathway." This interpretative logic is evaluated here by analyzing the relationship between thermodynamic behavior and folding kinetics in a class of simplified lattice protein models. The chain models studied have varying degrees of cooperative interplay (coupling) between local helical conformational preference and favorable nonlocal interactions. When model cooperativity is high, as native conditions are weakened, "isotherms" of free energy of exchange for residues belonging to the same helix merge together before global unfolding. The point of merger depends on the model energetic favorability of the helix. This trend is similar to the corresponding experimental observations. Kinetically, we find that the ordering of helix formation in the very last stage of native core assembly tends to follow the stabilities of their converged isotherms. In a majority (but not all) of folding trajectories, the final assembly of helices that are thermodynamically more stable against exchange precedes that of helices that are less stable against exchange. These model features are in partial agreement with common experimental interpretations. However, the model results also underscore the ensemble nature of the folding process: the kinetics of helix formation is not a discrete, strictly "all-or-none" process as that envisioned by certain non-explicit-chain models. Helices generally undergo many cycles of partial formation and dissolution before their conformations are fixed in the final assembly stage of folding, a kinetic stage that takes up only approximately 2% of the average folding time in the present model; and the ordering of the helices' final assembly in some trajectories can be different from the dominant ordering stipulated by the exchange isotherms.     

Proline residues in transmembrane alpha helices affect the folding of bacteriorhodopsin

Proline residues occur frequently in transmembrane alpha helices, which contrasts with their behaviour as helix-breakers in water-soluble proteins. The three membrane-embedded proline residues of bacteriorhodopsin have been replaced individually by alanine and glycine to give P50A, or P50G on helix B, P91A, or P91G on helix C, and P186A or P186G on helix F, and the effect on the protein folding kinetics has been investigated. The rate-limiting apoprotein folding step, which results in formation of a seven transmembrane, alpha helical state, was slower than wild-type protein for the Pro50 and Pro91 mutants, regardless of whether they were mutated to Ala or Gly. These proline residues give rise to several inter-helix contacts, which are therefore important in folding to the seven transmembrane helix state. No evidence for cis-trans isomerisations of the peptidyl prolyl bonds was found during this rate-limiting apoprotein folding step. Mutations of all three membrane-embedded proline residues affected the subsequent retinal binding and final folding to bacteriorhodopsin, suggesting that these proline residues contribute to formation of the retinal binding pocket within the helix bundle, again via helix/helix interactions. These results point to proline residues in transmembrane alpha helices being important in the folding of integral membrane proteins. The helix/helix interactions and hydrogen bonds that arise from the presence of proline residues in transmembrane alpha helices can affect the formation of transmembrane alpha helix bundles as well as cofactor binding pockets.     

Peptide mimics of the M13 coat protein transmembrane segment. Retention of helix-helix interaction motifs

Sequence-specific noncovalent helix-helix interactions between transmembrane (TM) segments in proteins are investigated by incorporating selected TM sequences into synthetic peptides using the construct CKKK-TM-KKK. The peptides are of suitable hydrophobicity for spontaneous membrane insertion, whereas formation of an N-terminal S-S bond can bring pairs of TM helices into proximity and promote their parallel orientation. Using the propensity of the protein to undergo thermally induced alpha-helix --> beta-sheet transitions as a parameter for helix stability, we compared the wild type and mutant (V29A and V31A) bacteriophage M13 coat proteins with their corresponding TM peptide constructs (M13 residues 24-42). Our results demonstrated that the relevant helix-helix tertiary contacts found in the intact proteins persist in the peptide mimics. Molecular dynamics simulations support the tight "two in-two out" dimerization motif for V31A consistent with mutagenesis data. The overall results reinforce the notion of TM segments as autonomous folding domains and suggest that the generic peptide construct provides a viable reductionist system for membrane protein structural and computational analysis.     

Preformed secondary structure drives the association reaction of GCN4-p1, a model coiled-coil system

The structure of the transition state for the rate-limiting step in the folding and association of the homodimeric coiled-coil peptide GCN4-p1, was probed by mutational analysis. A series of quadruple amino acid replacements that spanned the helix propensity scale were made at the four external f positions in the heptad repeat. Equilibrium and kinetic circular dichroism studies demonstrate that both the stability and the unfolding and refolding rate constants vary with helix propensity but also reflect interactions of the altered side-chains with their local environments. Pairwise replacements and fragment studies show that the two C-terminal heptads are the likely source of the nucleating helices. Helix-helix recognition between preformed elements of secondary structure plays an important role in this fundamental folding reaction.     

QHELIX: A Computational Tool for the Improved Measurement of Inter-Helical Angles in Proteins

Knowledge about the assembled structures of the secondary elements in proteins is essential to understanding protein folding and functionality. In particular, the analysis of helix geometry is required to study helix packing with the rest of the protein and formation of super secondary structures, such as, coiled coils and helix bundles, formed by packing of two or more helices. Here we present an improved computational method, QHELIX, for the calculation of the orientation angles between helices. Since a large number of helices are known to be in curved shapes, an appropriate definition of helical axes is a prerequisite for calculating the orientation angle between helices. The present method provides a quantitative measure on the irregularity of helical shape, resulting in discriminating irregular-shaped helices from helices with an ideal geometry in a large-scale analysis of helix geometry. It is also capable of straightforwardly assigning the direction of orientation angles in a consistent way. These improvements will find applications in finding a new insight on the assembly of protein secondary structure. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.     

Computer simulations of membrane protein folding: structure and dynamics

A lattice model of membrane proteins with a composite energy function is proposed to study their folding dynamics and native structures using Monte Carlo simulations. This model successfully predicts the seven helix bundle structure of sensory rhodopsin I by practicing a three-stage folding. Folding dynamics of a transmembrane segment into a helix is further investigated by varying the cooperativity in the formation of alpha helices for both random folding and assisted folding. The chain length dependence of the folding time of a hydrophobic segment to a helical state is studied for both free and anchored chains. An unusual length dependence in the folding time of anchored chains is observed.     

Review: conformation and folding of novel beta-structural elements in viral fiber proteins: the triple beta-spiral and triple beta-helix

Apart from alpha-helical coiled coils and the collagen triple helices, fibrous proteins can contain beta-structure in various conformations. Elongated enzymes such as pectate lyase and the bacteriophage P22 tailspike protein contain single-stranded beta-helices. Virus and bacteriophage fibers, which are often trimeric, have been shown to contain novel triple-stranded beta-structures such as the triple beta-spiral and the triple beta-helix. The conformation and folding of viral fibers containing beta-structure are discussed.     

Predicting Secondary Structural Folding Kinetics for Nucleic Acids

We report a new computational approach to the prediction of RNA secondary structure folding kinetics. In this approach, each elementary kinetic step is represented as the transformation between two secondary structures that differ by a helix. Based on the free energy landscape analysis, we identify three types of dominant pathways and the rate constants for the kinetic steps: 1), formation; 2), disruption of a helix stem; and 3), helix formation with concomitant partial melting of a competing (incompatible) helix. The third pathway, termed the tunneling pathway, is the low-barrier dominant pathway for the conversion between two incompatible helices. Comparisons with experimental data indicate that this new method is quite reliable in predicting the kinetics for RNA secondary structural folding and structural rearrangements. The approach presented here may provide a robust first step for further systematic development of a predictive theory for the folding kinetics for large RNAs.     

Position-dependence of stabilizing polar interactions of asparagine in transmembrane helical bundles

Recent studies with model peptides and statistical analyses of the crystal structures of membrane proteins have shown that buried polar interactions contribute significantly to the stabilization of the three-dimensional structures of membrane proteins. Here, we probe how the location of these polar groups along the transmembrane helices affect their free energies of interaction. Asn residues were placed singly and in pairs at three positions within a model transmembrane helix, which had previously been shown to support the formation of trimers in micelles. The model helix was designed to form a transmembrane coiled coil, with Val side chains at the "a" positions of the heptad repeat. Variants of this peptide were prepared in which an Asn residue was introduced at one or more of the "a" positions, and their free energies of association were determined by analytical ultracentrifugation. When placed near the middle of the transmembrane helix, the formation of trimers was stabilized by at least -2.0 kcal/mol per Asn side chain. When the Asn was placed at the interface between the hydrophobic and polar regions of the peptide, the substitution was neither stabilizing nor destabilizing (0.0 +/- 0.5 kcal/mol of monomer). Finally, it has previously been shown that a Val-for-Asn mutation in a water-soluble coiled coil destabilizes the structure by approximately 1.5 kcal/mol of monomer [Acharya, A., et al. (2002) Biochemistry 41, 14122-14131]. Thus, the headgroup region of a micelle appears to have a conformational impact intermediate between that of bulk water and the apolar region of micelle. A similarly large dependence on the location of the polar residues was found in a statistical survey of helical transmembrane proteins. The tendency of different types of residues to be buried in the interiors versus being exposed to lipids was analyzed. Asn and Gln show a very strong tendency to be buried when they are located near the middle of a transmembrane helix. However, when placed near the ends of transmembrane helices, they show little preference for the surface versus the interior of the protein. These data show that Asn side chains within the apolar region of the transmembrane helix provide a significantly larger driving force for association than Asn residues near the apolar/polar interface. Thus, although polar interactions are able to strongly stabilize the folding of membrane proteins, the energetics of association depend on their location within the hydrophobic region of a transmembrane helix.     

Folding propensities of peptide fragments of myoglobin.

Myoglobin has been studied extensively as a paradigm for protein folding. As part of an ongoing study of potential folding initiation sites in myoglobin, we have synthetized a series of peptides covering the entire sequence of sperm whale myoglobin. We report here on the conformation preferences of a series of peptides that cover the region from the A helix to the FG turn. Structural propensities were determined using circular dichroism and nuclear magnetic resonance spectroscopy in aqueous solution, trifluoroethanol, and methanol. Peptides corresponding to helical regions in the native protein, namely the B, C, D, and E helices, populate the alpha region of (phi, psi) space in water solution but show no measurable helix formation except in the presence of trifluoroethanol. The F-helix sequence has a much lower propensity to populate helical conformations even in TFE. Despite several attempts, we were not successful in synthesizing a peptide corresponding to the A-helix region that was soluble in water. A peptide termed the AB domain was constructed spanning the A- and B-helix sequences. The AB domain is not soluble in water, but shows extensive helix formation throughout the peptide when dissolved in methanol, with a break in the helix at a site close to the A-B helix junction in the intact folded myoglobin protein. With the exception of one local preference for a turn conformation stabilized by hydrophobic interactions, the peptides corresponding to turns in the folded protein do not measurably populate beta-turn conformations in water, and the addition of trifluoroethanol does not enhance the formation of either helical or turn structure. In contrast to the series of peptides described here, either studies of peptides from the GH region of myoglobin show a marked tendency to populate helical structures (H), nascent helical structures (G), or turn conformations (GH peptide) in water solution. This region, together with the A-helix and part of the B-helix, has been shown to participate in an early folding intermediate. The complete analysis of conformational properties of isolated myoglobin peptides supports the hypothesis that spontaneous secondary structure formation in local regions of the polypeptide may play an important role in the initiation of protein folding.     

Free-energy determinants of alpha-helix insertion into lipid bilayers.

A detailed treatment is provided of the various free-energy terms that contribute to the transfer of a polyalanine alpha-helix from the aqueous phase into lipid bilayers. In agreement with previous work, the hydrophobic effect is found to provide the major driving force for helix insertion. However, an opposing effect of comparable magnitude is also identified and is attributed to the large free-energy penalty associated with the desolvation of peptide hydrogen bonds on transfer to the low dielectric environment of the bilayer. Lipid perturbation effects as well as the entropy loss associated with helix immobilization in the bilayer are also evaluated. Two configurations of a membrane-bound 25mer polyalanine helix were found to be lower in free energy than the isolated helix in the aqueous phase. The first corresponds to the case of vertical insertion, in which a helix terminus protrudes from each side of the bilayer. The second minimum is for the case of horizontal insertion, for which the helix is adsorbed upon the surface of the bilayer. The calculated free-energy minima are found to be in good agreement with recent measurements of related systems. Large free-energy barriers resulting from desolvation of unsatisfied hydrogen-bonding groups at the helix termini are obtained for both insertion processes. The barriers for insertion are significantly reduced if the helix termini are assumed to be "capped" through the formation of hydrogen bonds with polar sidechains. For uncapped helices, our results support recently proposed models in which helices are inserted by first adsorbing on the membrane surface and then having one terminus "swing around" so as to penetrate the bilayer.     

Helices in peptoids of alpha- and beta-peptides

Peptoids of alpha- and beta-peptides (alpha- and beta-peptoids) can be obtained by shifting the amino acid side chains from the backbone carbon atoms of the monomer constituents to the peptide nitrogen atoms. They are, therefore, N-substituted poly-glycines and poly-beta-alanines, respectively. Due to the substituted nitrogen atoms, the ability for hydrogen bond formation between peptide bonds gets lost. It may be very interesting to see whether such non-natural oligomers could be regarded as foldamers, which fold into definite backbone conformers. In this paper, we provide a complete overview on helix formation in alpha- and beta-peptoids on the basis of systematic theoretical conformational analyses employing the methods of ab initio molecular orbital (MO) theory. It can be shown that the alpha- and beta-peptoid structures form helical structures with both trans and cis peptide bonds despite the missing hydrogen bonds. Obviously, the conformational properties of the backbone are more important for folding than the possibility of hydrogen bonding. There are close relationships between the helices of alpha-peptoids and poly-glycine and poly-proline helices of alpha-peptides, whereas the helices of beta-peptoids correspond to the well-known helical structures of beta-peptides as, for instance, the 3(1)-helix of beta-peptides with 14-membered hydrogen-bonded rings. Thus, alpha- and beta-peptoids enrich the field of foldamers and may be used as useful tools in peptide and protein design.