A propensity scale for type II polyproline helices (PPII): aromatic amino acids in proline-rich sequences strongly disfavor PPII due to proline-aromatic interactions

Type II polyproline helices (PPII) are a fundamental secondary structure of proteins, common in globular and nonglobular regions and important in cellular signaling. We developed a propensity scale for PPII using a host-guest system with sequence Ac-GPPXPPGY-NH(2), where X represents any amino acid. We found that proline has the highest PPII propensity, but most other amino acids display significant PPII propensities. The PPII propensity of leucine was the highest of all propensities of non-proline residues. Alanine and residues with linear side chains displayed the next highest PPII propensities. Three classes of residues displayed lower PPII propensities: β-branched amino acids (Thr, Val, and Ile), short amino acids with polar side chains (Asn, protonated Asp, Ser, Thr, and Cys), and aromatic amino acids (Phe, Tyr, and Trp). tert-Leucine particularly disfavored PPII. The basis of the low PPII propensities of aromatic amino acids in this context was significant cis-trans isomerism, with proline-rich peptides containing aromatic residues exhibiting 45-60% cis amide bonds, due to Pro-cis-Pro-aromatic and aromatic-cis-Pro amide bonds.     

Foldability and chameleon propensity of fold‐switching protein sequences

It has been shown recently by Porter & Looger that a significant number of proteins exist that can form more than one stable fold. This note examines the sequences of these fold‐switching proteins by (a) calculating their foldability scores recently introduced by the present author and (b) comparing the propensity of chameleon sequences in fold switchers and in non fold switchers. It has been found that the average foldability score of the fold switchers indicates weaker foldability. As for the propensity of chameleon sequences of length 5 to 7 it was found, somewhat surprisingly, that there is only a very small difference between the fold switchers and the non fold switchers. Furthermore, when looking at amino acid propensities, for several amino acids there was even an opposing trend in the deviation of their propensities from the overall amino acid propensities.     

Local and nonlocal interactions in globular proteins and mechanisms of alcohol denaturation.

How important are helical propensities in determining the conformations of globular proteins? Using the two-dimensional lattice model and two monomer types, H (hydrophobic) and P (polar), we explore both nonlocal interactions, through an HH contact energy, as developed in earlier work, and local interactions, through a helix energy, σ. By computer enumeration, the partition functions for short chains are obtained without approximation for the full range of both types of energy. When nonlocal interactions dominate, some sequences undergo coil-globule collapse to a unique native structure. When local interactions dominate, all sequences undergo helix–coil transitions. For two different conformational properties, the closest correspondence between the lattice model and proteins in the Protein Data Bank is obtained if the model local interactions are made small compared to the HH contact interaction, suggesting that helical propensities may be only weak determinants of globular protein structures in water. For some HP sequences, varying σ/ leads to additional sharp transitions (sometimes several) and to “conformational switching” between unique conformations. This behavior resembles the transitions of globular proteins in water to helical states in alcohols. In particular, comparison with experiments shows that whereas urea as a denaturant is best modeled as weakening both local and nonlocal interactions, trifluoroethanol is best modeled as mainly weakening HH interactions and slightly enhancing local helical interactions.     

A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix.

In this work, we have analyzed the relative importance of secondary versus tertiary interactions in stabilizing and guiding protein folding. For this purpose, we have designed four different mutants to replace the alpha-helix of the GB1 domain by a sequence with strong beta-hairpin propensity in isolation. In particular, we have chosen the sequence of the second beta-hairpin of the GB1 domain, which populates the native conformation in aqueous solution to a significant extent. The resulting protein has roughly 30 % of its sequence duplicated and maintains the 3D-structure of the wild-type protein, but with lower stability (up to -5 kcal/mol). The loss of intrinsic helix stability accounts for about 80 % of the decrease in free energy, illustrating the importance of local interactions in protein stability. Interestingly enough, all the mutant proteins, included the one with the duplicated beta-hairpin sequence, fold with similar rates as the GB1 domain. Essentially, it is the nature of the rate-limiting step in the folding reaction that determines whether a particular interaction will speed up, or not, the folding rates. While local contacts are important in determining protein stability, residues involved in tertiary contacts in combination with the topology of the native fold, seem to be responsible for the specificity of protein structures. Proteins with non-native secondary structure tendencies can adopt stable folds and be as efficient in folding as those proteins with native-like propensities.     

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.     

Does Lack of Secondary Structure Imply Intrinsic Disorder in Proteins? A Sequence Analysis

Intrinsically disordered proteins (IDPs)/protein regions (IDPRs) lack unique three-dimensional structure at the level of secondary and/or tertiary structure and are represented as an ensemble of interchanging conformations. To investigate the role of presence/absence of secondary structures in promoting intrinsic disorder in proteins, a comparative sequence analysis of IDPs, IDPRs and proteins with minimal secondary structures (less than 5%) is required. A sequence analysis reveals proteins with minimal secondary structure content have high mean net positive charge, low mean net hydrophobicity and low sequence complexity. Interestingly, analysis of the relative local electrostatic interactions reveal that an increase in the relative repulsive interactions between amino acids separated by three or four residues lead to either loss of secondary structure or intrinsic disorder. IDPRs show increase in both local negative-negative and positive-positive repulsive interactions. While IDPs show a marked increase in the local negative-negative interactions, proteins with minimal secondary structure depict an increase in the local positive-positive interactions. IDPs and IDPRs are enriched in D, E and Q residues, while proteins with minimal secondary structure are depleted of these residues. Proteins with minimal secondary structures have higher content of G and C, while IDPs and IDPRs are depleted of these residues. These results confirm that proteins with minimal secondary structure have a distinctly different propensity for charge, hydrophobicity, specific amino acids and local electrostatic interactions as compared to IDPs/IDPRs. Thus we conclude that lack of secondary structure may be a necessary but not a sufficient condition for intrinsic disorder in proteins.     

Folding type-specific secondary structure propensities of amino acids,derived from α-helical, β-sheet, α/β, and α+β proteins of known structures

Folding type-specific secondary structure propensities of 20 naturally occurring amino acids have been derived from α-helical, β-sheet, α/β, and α+β proteins of known structures. These data show that each residue type of amino acids has intrinsic propensities in different regions of secondary structures for different folding types of proteins. Each of the folding types shows markedly different rank ordering, indicating folding type-specific effects on the secondary structure propensities of amino acids. Rigorous statistical tests have been made to validate the folding type-specific effects. It should be noted that α and β proteins have relatively small α-helices and β-strands forming propensities respectively compared with those of α+β and α/β proteins. This may suggest that, with more complex architectures than α and β proteins, α+β and α/β proteins require larger propensities to distinguish from interacting α-helices and β-strands. Our finding of folding type-specific secondary structure propensities suggests that sequence space accessible to each folding type may have differing features. Differing sequence space features might be constrained by topological requirement for each of the folding types. Almost all strong β-sheet forming residues are hydrophobic in character regardless of folding types, thus suggesting the hydrophobicities of side chains as a key determinant of β-sheet structures. In contrast, conformational entropy of side chains is a major determinant of the helical propensities of amino acids, although other interactions such as hydrophobicities and charged interactions cannot be neglected. These results will be helpful to protein design, class-based secondary structure prediction, and protein folding. © 1998 John Wiley & Sons, Inc. Biopoly 45: 35–49, 1998     

Amino-acid substitutions at the fully exposed P1 site of bovine pancreatic trypsin inhibitor affect its stability

It is widely accepted that solvent-exposed sites in proteins play only a negligible role in determining protein energetics. In this paper we show that amino acid substitutions at the fully exposed Lys15 in bovine pancreatic trypsin inhibitor (BPTI) influenced the CD- and DSC-monitored stability: The T(den) difference between the least (P1 Trp) and the most stable (P1 His) mutant is 11.2 degrees C at pH 2.0. The DeltaH(den) versus T(den) plot for all the variants at three pH values (2.0, 2.5, 3.0) is linear (DeltaC(p,den) = 0.41 kcal* mole(-1) * K(-1); 1 cal = 4.18 J) leading to a DeltaG(den) difference of 2.1 kcal*mole(-1). Thermal denaturation of the variants monitored by CD signal at pH 2.0 in the presence of 6 M GdmCl again showed differences in their stability, albeit somewhat smaller (DeltaT(den) =7.1 degrees C). Selective reduction of the Cys14-Cys 38 disulfide bond, which is located in the vicinity of the P1 position did not eliminate the stability differences. A correlation analysis of the P1 stability with different properties of amino acids suggests that two mechanisms may be responsible for the observed stability differences: the reverse hydrophobic effect and amino acid propensities to occur in nonoptimal dihedral angles adopted by the P1 position. The former effect operates at the denatured state level and causes a drop in protein stability for hydrophobic side chains, due to their decreased exposure upon denaturation. The latter factor influences the native state energetics and results from intrinsic properties of amino acids in a way similar to those observed for secondary structure propensities. In conclusion, our results suggest that the protein-stability-derived secondary structure propensity scales should be taken with more caution.     

Conformational properties of a peptide model for unfolded alpha-helices

Models of protein folding often hypothesize that the first step is local secondary structure formation. The assumption is that unfolded polypeptide chains possess an intrinsic propensity to form these local secondary structures. On the basis of this idea, it is tempting to model the local conformational properties of unfolded proteins using well-established residue secondary structure propensities, in particular, alpha-helix forming propensities. We have used spectroscopic methods to investigate the conformational behavior of a host-guest series of peptides designed to model unfolded alpha-helices. A suitable peptide model for unfolded alpha-helices was determined from studies of the length dependence of the conformational properties of alanine-based peptides. The chosen host peptide possessed a small, detectable, alpha-helix content. Substituting various representative guest residues into the central position of the host peptide at times changed the conformational behavior dramatically, and often in ways that could not be predicted from known alpha-helix forming propensities. The data presented can be used to rationalize some of these propensities. However, it is clear that secondary structure propensities cannot be used to predict the local conformational properties of unfolded proteins.     

A helix propensity scale based on experimental studies of peptides and proteins.

The average globular protein contains 30% alpha-helix, the most common type of secondary structure. Some amino acids occur more frequently in alpha-helices than others; this tendency is known as helix propensity. Here we derive a helix propensity scale for solvent-exposed residues in the middle positions of alpha-helices. The scale is based on measurements of helix propensity in 11 systems, including both proteins and peptides. Alanine has the highest helix propensity, and, excluding proline, glycine has the lowest, approximately 1 kcal/mol less favorable than alanine. Based on our analysis, the helix propensities of the amino acids are as follows (kcal/mol): Ala = 0, Leu = 0.21, Arg = 0.21, Met = 0.24, Lys = 0.26, Gln = 0.39, Glu = 0.40, Ile = 0.41, Trp = 0.49, Ser = 0.50, Tyr = 0. 53, Phe = 0.54, Val = 0.61, His = 0.61, Asn = 0.65, Thr = 0.66, Cys = 0.68, Asp = 0.69, and Gly = 1.     

Environment-specific amino acid substitution tables: tertiary templates and prediction of protein folds.

The local environment of an amino acid in a folded protein determines the acceptability of mutations at that position. In order to characterize and quantify these structural constraints, we have made a comparative analysis of families of homologous proteins. Residues in each structure are classified according to amino acid type, secondary structure, accessibility of the side chain, and existence of hydrogen bonds from the side chains. Analysis of the pattern of observed substitutions as a function of local environment shows that there are distinct patterns, especially for buried polar residues. The substitution data tables are available on diskette with Protein Science. Given the fold of a protein, one is able to predict sequences compatible with the fold (profiles or templates) and potentially to discriminate between a correctly folded and misfolded protein. Conversely, analysis of residue variation across a family of aligned sequences in terms of substitution profiles can allow prediction of secondary structure or tertiary environment.     

Tertiary structural constraints on protein evolutionary diversity: templates, key residues and structure prediction

The pattern of residue substitution in divergently evolving families of globular proteins is highly variable. At each position in a fold there are constraints on the identities of amino acids from both the three-dimensional structure and the function of the protein. To characterize and quantify the structural constraints, we have made a comparative analysis of families of homologous globular proteins. Residues are classified according to amino acid type, secondary structure, accessibility of the sidechain, and existence of hydrogen bonds from sidechain to other sidechains or peptide carbonyl or amide functions. There are distinct patterns of substitution especially where residues are both solvent inaccessible and hydrogen bonded through their sidechains. The patterns of residue substitution can be used to construct templates or to identify 'key' residues if one or more structures are known. Conversely, analysis of conversation and substitution across a large family of aligned sequences in terms of substitution profiles can allow prediction of tertiary environment or indicate a functional role. Similar analyses can be used to test the validity of putative structures if several homologous sequences are available.     

Dihedral angle and secondary structure database of short amino acid fragments

Dihedral angles of amino acids are of considerable importance in protein tertiary structure prediction as they define the backbone of a protein and hence almost define the protein's entire conformation. Most ab initio protein structure prediction methods predict the secondary structure of a protein before predicting the tertiary structure because three-dimensional fold consists of repeating units of secondary structures. Hence, both dihedral angles and secondary structures are important in tertiary structure prediction of proteins. Here we describe a database called DASSD (Dihedral Angle and Secondary Structure Database of Short Amino acid Fragments) that contains dihedral angle values and secondary structure details of short amino acid fragments of lengths 1, 3 and 5. Information stored in this database was extracted from a set of 5,227 non-redundant high resolution (less than 2-angstroms) protein structures. In total, DASSD stores details for about 733,000 fragments. This database finds application in the development of ab initio protein structure prediction methods using fragment libraries and fragment assembly techniques. It is also useful in protein secondary structure prediction.

Availability     

A comparative study of the relationship between protein structure and beta-aggregation in globular and intrinsically disordered proteins

A growing number of proteins are being identified that are biologically active though intrinsically disordered, in sharp contrast with the classic notion that proteins require a well-defined globular structure in order to be functional. At the same time recent work showed that aggregation and amyloidosis are initiated in amino acid sequences that have specific physico-chemical properties in terms of secondary structure propensities, hydrophobicity and charge. In intrinsically disordered proteins (IDPs) such sequences would be almost exclusively solvent-exposed and therefore cause serious solubility problems. Further, some IDPs such as the human prion protein, synuclein and Tau protein are related to major protein conformational diseases. However, this scenario contrasts with the large number of unstructured proteins identified, especially in higher eukaryotes, and the fact that the solubility of these proteins is often particularly good. We have used the algorithm TANGO to compare the beta aggregation tendency of a set of globular proteins derived from SCOP and a set of 296 experimentally verified, non-redundant IDPs but also with a set of IDPs predicted by the algorithms DisEMBL and GlobPlot. Our analysis shows that the beta-aggregation propensity of all-alpha, all-beta and mixed alpha/beta globular proteins as well as membrane-associated proteins is fairly similar. This illustrates firstly that globular structures possess an appreciable amount of structural frustration and secondly that beta-aggregation is not determined by hydrophobicity and beta-sheet propensity alone. We also show that globular proteins contain almost three times as much aggregation nucleating regions as IDPs and that the formation of highly structured globular proteins comes at the cost of a higher beta-aggregation propensity because both structure and aggregation obey very similar physico-chemical constraints. Finally, we discuss the fact that although IDPs have a much lower aggregation propensity than globular proteins, this does not necessarily mean that they have a lower potential for amyloidosis.     

Interaction of 18-residue peptides derived from amphipathic helical segments of globular proteins with model membranes

We investigated the interaction of six 18-residue peptides derived from amphipathic helical segments of globular proteins with model membranes. The net charge of the peptides at neutral pH varies from −1 to +6. Circular dichroism spectra indicate that peptides with a high net positive charge tend to fold into a helical conformation in the presence of negatively charged lipid vesicles. In helical conformation, their average hydrophobic moment and hydrophobicity would render them surface-active. The composition of amino acids on the polar face of the helix in the peptides is considerably different. The peptides show variations in their ability to permeabilise zwitterionic and anionic lipid vesicles. Whereas increased net positive charge favours greater permeabilisation, the distribution of charged residues in the polar face also plays a role in determining membrane activity. The distribution of amino acids in the polar face of the helix in the peptides that were investigated do not fall into the canonical classes described. Amphipathic helices, which are part of proteins, with a pattern of amino acid distribution different from those observed in class L, A and others, could help in providing newer insights into peptide-membrane interactions.     

Dihedral angle preferences of amino acid residues forming various non-local interactions in proteins

In theory, a polypeptide chain can adopt a vast number of conformations, each corresponding to a set of backbone rotation angles. Many of these conformations are excluded due to steric overlaps. Ramachandran and coworkers were the first to look into this problem by plotting backbone dihedral angles in a two-dimensional plot. The conformational space in the Ramachandran map is further refined by considering the energetic contributions of various non-bonded interactions. Alternatively, the conformation adopted by a polypeptide chain may also be examined by investigating interactions between the residues. Since the Ramachandran map essentially focuses on local interactions (residues closer in sequence), out of interest, we have analyzed the dihedral angle preferences of residues that make non-local interactions (residues far away in sequence and closer in space) in the folded structures of proteins. The non-local interactions have been grouped into different types such as hydrogen bond, van der Waals interactions between hydrophobic groups, ion pairs (salt bridges), and ππ-stacking interactions. The results show the propensity of amino acid residues in proteins forming local and non-local interactions. Our results point to the vital role of different types of non-local interactions and their effect on dihedral angles in forming secondary and tertiary structural elements to adopt their native fold.     

A statistical analysis of the PPII propensity of amino acid guests in proline-rich peptides

There has been considerable debate about the intrinsic PPII propensity of amino-acid residues in denatured polypeptides. Experimentally, the propensity scale is based on the behavior of guest amino-acid residues placed in the middle of polyproline hosts. We have used classical molecular dynamics simulations, with state-of-the-art force fields to carry out a comprehensive analysis of the conformational equilibria of the proline-based host oligopeptides with single guests. The tracked structural characteristics include the PPII content, the cis/trans isomerization of the prolyl bonds, the puckering of the pyrrolidine rings of the proline residues, and the secondary structural motifs. We find no evidence for an intrinsic PPII propensity in any of the guest amino acids other than proline. Instead, the PPII content as derived from experiments may be explained in terms of: 1), a local correlation between the dihedral angles of the guest amino acid and the proline residue immediately preceding it; and 2), a nonlocal correlation between the cis/trans states of the peptide bonds. In terms of the latter, we find that the presence of a guest (other than proline, tyrosine, or tryptophan) increases the trans content of most of the prolyl bonds, which results in an effective increase of the peptide PPII content. With respect to the local dihedral correlations, we find that these are well described in terms of the so-called odds-ratio statistic. Expressed in terms of free energy language, the PPII content based on the odds-ratio of the relevant residues correlate well with the experimentally measured PPII content.     

Intrinsic secondary structure propensities of the amino acids,using statistical ϕ–ψ matrices: Comparison with experimental scales

Today there are several different experimental scales for the intrinsic α-helix as well as β-strand, propensities of the 20 amino acids obtained from the thermodynamic analysis of various model systems. These scales do not compare well with those extracted from statistical analysis of three-dimensional structure databases. Possible explanations for this could be the limited size of the databases used, the definitions of intrinsic propensities, or the theoretical approach. Here we report a statistical determination of α-helix and β-strand propensities derived from the analysis of a database of 279 three-dimensional structures. Contrary to what has been generally done, we have considered a particular residue as in α-helix or β-strand conformation by looking only at its dihedral angles (?–ψ matrices). Neither the identity nor the conformation of the surrounding residues in the amino acid sequence has been taken into consideration. Pseudoenergy empirical scales have been calculated from the statistical propensities. These scales agree very well with the experimental ones in relative and absolute terms. Moreover, its correlation with the average of the experimental scales for α-helix or β-strand is as good as the correlations of the individual experimental scales with the average. These results show that by using a large enough database and a proper definition for the secondary structure propensities, it is possible to obtain a scale as good as any of experimental origin. Interestingly the ?–ψ analysis of the Ramachandran plot suggests that the amino acids could have different β-strand propensities in different subregions of the β-strand area. © 1994 Wiley-Liss, Inc.     

Protein secondary structure: entropy, correlations and prediction

MOTIVATION: Is protein secondary structure primarily determined by local interactions between residues closely spaced along the amino acid backbone or by non-local tertiary interactions? To answer this question, we measure the entropy densities of primary and secondary structure sequences, and the local inter-sequence mutual information density. RESULTS: We find that the important inter-sequence interactions are short ranged, that correlations between neighboring amino acids are essentially uninformative and that only one-fourth of the total information needed to determine the secondary structure is available from local inter-sequence correlations. These observations support the view that the majority of most proteins fold via a cooperative process where secondary and tertiary structure form concurrently. Moreover, existing single-sequence secondary structure prediction algorithms are almost optimal, and we should not expect a dramatic improvement in prediction accuracy. AVAILABILITY: Both the data sets and analysis code are freely available from our Web site at http://compbio.berkeley.edu/     

Rapid assessment of contact-dependent secondary structure propensity: relevance to amyloidogenic sequences

We have previously demonstrated that calculation of contact-dependent secondary structure propensity (CSSP) is highly sensitive in detecting non-native beta-strand propensities in the core sequences of known amyloidogenic proteins. Here we describe a CSSP method based on an artificial neural network that rapidly and accurately quantifies the influence of tertiary contacts (TCs) on secondary structure propensity in local regions of protein sequences. The present method exhibited 72% accuracy in predicting the alternate secondary structure adopted by chameleon sequences located in highly disparate TC regions. Analysis of 1930 nonhomologous protein domains reveals that the alpha-helix and the beta-strand largely share the same sequence context, and that tertiary context is a major determinant of the native conformation. Conversely, it appears that the propensity of random coils for either the alpha-helix or the beta-strand is largely invariant to tertiary effects. The present CSSP method successfully reproduced the amyloidogenic character observed in local regions of the human islet amyloid polypeptide (hIAPP). Furthermore, CSSP profiles were strongly correlated (r = 0.76) with the observed mutational effects on the aggregation rate of acylphosphatase. Taken together, these results provide compelling evidence in support of the present CSSP approach as a sensitive probe useful for analysis of full-length proteins and for detection of core sequences that may trigger amyloid fibril formation. The combined speed and simplicity of the CSSP method lends itself to proteome-wide analysis of the amyloidogenic nature of common proteins.