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.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

Effects of i and i+3 residue identity on cis-trans isomerism of the aromatic(i+1)-prolyl(i+2) amide bond: implications for type VI beta-turn formation

Cis-trans isomerization of amide bonds plays critical roles in protein molecular recognition, protein folding, protein misfolding, and disease. Aromatic-proline sequences are particularly prone to exhibit cis amide bonds. The roles of residues adjacent to a tyrosine-proline residue pair on cis-trans isomerism were examined. A short series of peptides XYPZ was synthesized and cis-trans isomerism was analyzed. Based on these initial studies, a series of peptides XYPN, X = all 20 canonical amino acids, was synthesized and analyzed by NMR for i residue effects on cis-trans isomerization. The following effects were observed: (a) aromatic residues immediately preceding Tyr-Pro disfavor cis amide bonds, with K(trans/cis)= 5.7-8.0, W > Y > F; (b) proline residues preceding Tyr-Pro lead to multiple species, exhibiting cis-trans isomerization of either or both X-Pro amide bonds; and (c) other residues exhibit similar values of K(trans/cis) (= 2.9-4.2), with Thr and protonated His exhibiting the highest fraction cis. beta-Branched and short polar residues were somewhat more favorable in stabilizing the cis conformation. Phosphorylation of serine at the i position modestly increases the stability of the cis conformer. In addition, the effect of the i+3 residue was examined in a limited series of peptides TYPZ. NMR data indicated that aromatic residues, Pro, Asn, Ala, and Val at the i+3 residue all favor cis amide bonds, with aromatic residues and Asn favoring more compact phi at Tyr(cis) and Ala and Pro favoring more extended phi at Tyr(cis). D-Alanine at the i+3 position particularly disfavors cis amide bonds.     

Amino acid propensities for the collagen triple-helix

Determination of the tendencies of amino acids to form alpha-helical and beta-sheet structures has been important in clarifying stabilizing interactions, protein design, and the protein folding problem. In this study, we have determined for the first time a complete scale of amino acid propensities for another important protein motif: the collagen triple-helix conformation with its Gly-X-Y repeating sequence. Guest triplets of the form Gly-X-Hyp and Gly-Pro-Y are used to quantitate the conformational propensities of all 20 amino acids for the X and Y positions in the context of a (Gly-Pro-Hyp)(8) host peptide. The rankings for both the X and Y positions show the highly stabilizing nature of imino acids and the destabilizing effects of Gly and aromatic residues. Many residues show differing propensities in the X versus Y position, related to the nonequivalence of these positions in terms of interchain interactions and solvent exposure. The propensity of amino acids to adopt a polyproline II-like conformation plays a role in their triple-helix rankings, as shown by a moderate correlation of triple-helix propensity with frequency of occurrence in polyproline II-like regions. The high propensity of ionizable residues in the X position suggests the importance of interchain hydrogen bonding directly or through water to backbone carbonyls or hydroxyprolines. The low propensity of side chains with branching at the C(delta) in the Y position supports models suggesting these groups block solvent access to backbone C=O groups. These data provide a first step in defining sequence-dependent variations in local triple-helix stability and binding, and are important for a general understanding of side chain interactions in all proteins.     

固有无序蛋白与蛋白质相互作用位点残基特征分析

本文对固有无序蛋白(IDPs)与其他蛋白质相互作用位点残基特征进行了研究.首先在数据库中选出满足条件的109条IDPs蛋白质链及与其他配体蛋白形成的299个IDPs-蛋白质复合物,然后提取复合物中作为相互作用位点的IDPs-蛋白质残基.这109条IDPs链中共含有50 031个氨基酸残基,其中处于作用位点的残基有4 822个.通过分析发现,20种氨基酸在形成IDPs-蛋白质相互作用位点残基时具有不同的倾向性,根据形成作用位点残基的倾向性,20种氨基酸可分成三大类:倾向型氨基酸(ILE、LEU、ARG、PHE、TYR、MET、TRP)、中间型氨基酸(GLN、GLU、THR、LYS、VAL、ASP、HIS)、非倾向型氨基酸(PRO、SER、GLY、ALA、ASN、CYS).研究结果还进一步表明,不同氨基酸在有序区域与无序区域形成IDPs-蛋白质作用位点残基的倾向性不同.其中,氨基酸TRP、LEU、ILE、CYS在有序和无序区域形成作用位点残基的差异性尤为明显,而氨基酸GLU、PHE、HIS、ALA则基本没有多大差别.对IDPs-蛋白质相互作用位点残基理化特征进行分析发现:疏水性强、侧链净电荷量较少、极性较小、溶剂可及性表面积较大、侧链体积较大、极化率较大的氨基酸比较倾向于形成作用位点残基.主成分分析结果显示,残基的极化率、侧链体积和溶剂可及表面积对作用位点残基影响最大.     

Host-guest scale of left-handed polyproline II helix formation

Despite the clear importance of the left-handed polyproline II (PPII) helical conformation in many physiologically important processes as well as its potential significance in protein unfolded states, little is known about the physical determinants of this conformation. We present here a scale of relative PPII helix-forming propensities measured for all residues, except tyrosine and tryptophan, in a proline-based host peptide system. Proline has the highest measured propensity in this system, a result of strong steric interactions that occur between adjacent prolyl rings. The other measured propensities are consistent with backbone solvation being an important component in PPII helix formation. Side chain to backbone hydrogen bonding may also play a role in stabilizing this conformation. The PPII helix-forming propensity scale will prove useful in future studies of the conformational properties of proline-rich sequences as well as provide insights into the prevalence of PPII helices in protein unfolded states.     

Hydrogen bonds between short polar side chains and peptide backbone: prevalence in proteins and effects on helix-forming propensities

A survey of 322 proteins showed that the short polar (SP) side chains of four residues, Thr, Ser, Asp, and Asn, have a very strong tendency to form hydrogen bonds with neighboring backbone amides. Specifically, 32% of Thr, 29% of Ser, 26% of Asp, and 19% of Asn engage in such hydrogen bonds. When an SP residue caps the N terminal of a helix, the contribution to helix stability by a hydrogen bond with the amide of the N3 or N2 residue is well established. When an SP residue is in the middle of a helix, the side chain is unlikely to form hydrogen bonds with neighboring backbone amides for steric and geometric reasons. In essence the SP side chain competes with the backbone carbonyl for the same hydrogen-bonding partner (i.e., the backbone amide) and thus SP residues tend to break backbone carbonyl-amide hydrogen bonds. The proposition that this is the origin for the low propensities of SP residues in the middle of alpha helices (relative to those of nonpolar residues) was tested. The combined effects of restricting side-chain rotamer conformations (documented by Creamer and Rose, Proc Acad Sci USA, 1992;89:5937-5941; Proteins, 1994;19:85-97) and excluding side- chain to backbone hydrogen bonds by the helix were quantitatively analyzed. These were found to correlate strongly with four experimentally determined scales of helix-forming propensities. The correlation coefficients ranged from 0.72 to 0.87, which are comparable to those found for nonpolar residues (for which only the loss of side-chain conformational entropy needs to be considered).     

Polyproline helices in protein structures: A statistical survey

A statistical survey of polyproline II (PPII) helices extracted from protein crystal structures is here reported. The average hydrophobicity of these helices is intermediate between those displayed by beta-strands and coil regions and is similar to that of alpha-helices. PPII helices with amphipathic properties have been identified and classified. Amino acid propensities for PPII helices derived in this study differ significantly from those previously reported. They show a little albeit significant correlation with propensities for alpha-helices whereas they are fully non-correlated to propensities for beta-sheets. Finally, PPII propensities have been correlated with amino acid frequencies in structural proteins, such as collagen and extensins.     

Mechanism of Cis-Inhibition of PolyQ Fibrillation by PolyP: PPII Oligomers and the Hydrophobic Effect

PolyQ peptides teeter between polyproline II (PPII) and β-sheet conformations. In tandem polyQ-polyP peptides, the polyP segment tips the balance toward PPII, increasing the threshold number of Gln residues needed for fibrillation. To investigate the mechanism of cis-inhibition by flanking polyP segments on polyQ fibrillation, we examined short polyQ, polyP, and tandem polyQ-polyP peptides. These polyQ peptides have only three glutamines and cannot form β-sheet fibrils. We demonstrate that polyQ-polyP peptides form small, soluble oligomers at high concentrations (as shown by size exclusion chromatography and diffusion coefficient measurements) with PPII structure (as shown by circular dichroism spectroscopy and ³J_HN-Cα constants of Gln residues from constant time correlation spectroscopy NMR). Nuclear Overhauser effect spectroscopy and molecular modeling suggest that self-association of these peptides occurs as a result of both hydrophobic and steric effects. Pro side chains present three methylenes to solvent, favoring self-association of polyP through the hydrophobic effect. Gln side chains, with two methylene groups, can adopt a conformation similar to that of Pro side chains, also permitting self-association through the hydrophobic effect. Furthermore, steric clashes between Gln and Pro side chains to the C-terminal side of the polyQ segment favor adoption of the PPII-like structure in the polyQ segment. The conformational adaptability of the polyQ segment permits the cis-inhibitory effect of polyP segments on fibrillation by the polyQ segments in proteins such as huntingtin.     

Salmon calcitonin and amyloid beta: two peptides with amyloidogenic capacity adopt different conformational manifolds in their unfolded states

The molecular conformations of salmon calcitonin in aqueous solution have been investigated by exploiting the different influences of excitonic coupling on the amide I band profile in the isotropic and anisotropic Raman, FTIR, and vibrational circular dichroism spectra of a polypeptide. The N-terminal loop, caused by a disulfide bridge between cysteines at positions 1 and 7, was modeled by performing a conformational search by molecular mechanics calculations. The remaining part of the peptide chain was modeled as a mixture of three sequences containing different fractions of residues adopting poly-l-proline II (PPII), extended beta-strand, and alpha-helix-like conformations. This yielded an excellent reproduction of the experimentally observed amide I' band profiles. A comparison with recent data on the beta-amyloid fragment Abeta(1)(-)(28) revealed a lower PPII content and more conformational heterogeneity for calcitonin. Thus, our results underscore the notion that individual structural propensities of amino acid residues give rise to structural differences between the unfolded states of even long peptide chains, at variance with expectations based on a random or statistical coil model.     

Intramolecular surface contacts contain information about protein-protein interface regions

MOTIVATION: Some amino acids clearly show preferences over others in protein-protein interfaces. These preferences, or so-called interface propensities can be used for a priori interface prediction. We investigated whether the prediction accuracy could be improved by considering not single but pairs of residues in an interface. Here we present the first systematic analysis of intramolecular surface contacts in interface prediction. RESULTS: We show that preferences do exist for contacts within and around an interface region within one molecule: specific pairs of amino acids are more often occurring than others. Using intramolecular contact propensities in a blind test, higher average scores were assigned to interface residues than to non-interface residues. This effect persisted as small but significant when the contact propensities were corrected to eliminate the influence of single amino acid interface propensity. This indicates that intramolecular contact propensities may replace interface propensities in protein-protein interface prediction. AVAILABILITY: The source code is available on request from the authors.     

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.     

Conformation-Determining role for the N-acetyl group in the O-glycosidic linkage, α-GalNAc-Thr

An ¹H-nmr study of 2-acetamido-2-deoxy-3,4,6-tri-O-acetyl-D-galactopyranose (AcGalNAc) glycosylated Thr-containing tripeptides in Me₂SO-d₆ solution reveals two mutually exclusive intramolecular hydrogen bonds. In Z-Thr(AcGalNAc)-Ala-Ala-OMe, there is an intramolecular hydrogen bond between the Thr amide proton and the sugar N-acetyl carbonyl oxygen. The strength of this hydrogen bond will be dependent on the amino acid residues on the Thr C terminal side to some undetermined distance. In Ac-Thr(AcGalNAc)-Ala-Ala-OMe, a different intramolecular hydrogen bond between the sugar N-acetyl amide proton and the Thr carbonyl oxygen exists. The choice of hydrogen bonds seems dependent on the bulkiness of the residues on the Thr N terminal side. The consequence of such strong hydrogen bonds is a clearly defined orientation of the sugar moiety with respect to the peptide backbone. In the former, the plane of the sugar pyranose ring is roughly oriented perpendicularly to the peptide backbone. The latter orientation is where the plane of the sugar ring is roughly in line with the peptide backbone. In both orientations, the sugar moiety can increase the shielding of the neighboring amino acid residues from the solvent. The idea that the amino acid residues near the glycosylated Thr influence orientation of the sugar moiety with respect to the peptide backbone and in turn possibly hinder peptide backbone flexibility has interesting implications in the conformational as well as the biological role of O-glycoproteins.     

The degradation of glycoproteins with lithium borohydride. Isolation and analysis of O-glycopeptides with reduced C-terminal amino acid residue

By the example of fetuin and a blood-group-specific mucin from porcine stomach, we showed that, under conditions of reductive degradation of glycoproteins with LiBH4-LiOH in 70% aqueous tert-butyl alcohol, the reduction and cleavage of amide bonds occur much faster than the simultaneous beta-elimination of carbohydrate chains O-linked with Ser and Thr residues of the peptide chain. The major degradation products containing the O-linked glycans are the O-glycosylated derivatives of 2-aminopropane-1,3-diol and 2-aminobutane-1,3-diol (the products of reduction of glycosylated Ser and Thr) and the glycopeptides containing 2-4 amino acid residues with reduced C-terminal amino acid. Seventeen homogeneous O-glycopeptides were isolated from the fetuin degradation products by ion-exchange and reversed-phase HPLC. Their structures were determined by MALDI-TOF mass spectrometry and by analyses for amino acids, amino alcohols, and carbohydrates. The application of the reaction for characterization of O-glycans and localization of O-glycosylation sites in O- and N,O-glycoproteins is discussed.     

Role of amino acid hydrophobicity, aromaticity, and molecular volume on IAPP(20-29) amyloid self-assembly

Aromatic amino acids strongly promote cross-β amyloid formation; whether the amyloidogenicity of aromatic residues is due to high hydrophobicity and β-sheet propensity or formation of stabilizing π-π interactions has been debated. To clarify the role of aromatic residues on amyloid formation, the islet amyloid polypeptide 20-29 fragment [IAPP(20-29)], which contains a single aromatic residue (Phe 23), was adopted as a model. The side chain of residue 23 does not self-associate in cross-β fibrils of IAPP(20-29) (Nielsen et al., Angew Chem Int Ed 2009;48:2118-2121), allowing investigation of the amyloidogenicity of aromatic amino acids in a context where direct π-π interactions do not occur. We prepared variants of IAPP(20-29) in which Tyr, Leu, Phe, pentafluorophenylalanine (F5-Phe), Trp, cyclohexylalanine (Cha), α-naphthylalanine (1-Nap), or β-naphthylalanine (2-Nap) (in order of increasing peptide hydrophobicity) were incorporated at position 23 (SNNXGAILSS-NH2), and the kinetic and thermodynamic effects of these mutations on cross-β self-assembly were assessed. The Tyr, Leu, and Trp 23 variants failed to readily self-assemble at concentrations up to 1.5 mM, while the Cha 23 mutant fibrillized with attenuated kinetics and similar thermodynamic stability relative to the wild-type Phe 23 peptide. Conversely, the F5-Phe, 1-Nap, and 2-Nap 23 variants self-assembled at enhanced rates, forming fibrils with greater thermodynamic stability than the wild-type peptide. These results indicate that the high amyloidogenicity of aromatic amino acids is a function of hydrophobicity, β-sheet propensity, and planar geometry and not the ability to form stabilizing or directing π-π bonds.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

N- and C-capping preferences for all 20 amino acids in alpha-helical peptides.

We have determined the N- and C-capping preferences of all 20 amino acids by substituting residue X in the peptides NH2-XAKAAAAKAAAAKAAGY-CONH2 and in Ac-YGAAKAAAAKAAAAKAX-CO2H. Helix contents were measured by CD spectroscopy to obtain rank orders of capping preferences. The data were further analyzed by our modified Lifson-Roig helix-coil theory, which includes capping parameters (n and c), to find free energies of capping (-RT ln n and -RT ln c), relative to Ala. Results were obtained for charged and uncharged termini and for different charged states of titratable side chains. N-cap preferences varied from Asn (best) to Gln (worst). We find, as expected, that amino acids that can accept hydrogen bonds from otherwise free backbone NH groups, such as Asn, Asp, Ser, Thr, and Cys generally have the highest N-cap preference. Gly and acetyl group are favored, as are negative charges in side chains and at the N-terminus. Our N-cap preference scale agrees well with preferences in proteins. In contrast, we find little variation when changing the identity of the C-cap residue. We find no preference for Gly at the C-cap in contrast to the situation in proteins. Both N-cap and C-cap results for Tyr and Trp are inaccurate because their aromatic groups affect the CD spectrum. The data presented here are of value in rationalizing mutations at capping sites in proteins and in predicting the helix contents of peptides.     

A survey of left-handed polyproline II helices

Left-handed polyproline II helices (PPII) are contiguous elements of protein secondary structure in which the phi and psi angles of constituent residues are restricted to around -75 degrees and 145 degrees, respectively. They are important in structural proteins, in unfolded states and as ligands for signaling proteins. Here, we present a survey of 274 nonhomologous polypeptide chains from proteins of known structure for regions that form these structures. Such regions are rare, but the majority of proteins contain at least one PPII helix. Most PPII helices are shorter than five residues, although the longest found contained 12 amino acids. Proline predominates in PPII, but Gln and positively charged residues are also favored. The basis of Gln's prevalence is its ability to form an i, i + 1 side-chain to main-chain hydrogen bond with the backbone carbonyl oxygen of the proceeding residue; this helps to fix the psi angle of the Gln and the phi and psi of the proceeding residue in PPII conformations and explains why Gln is favored at the first position in a PPII helix. PPII helices are highly solvent exposed, which explains why apolar amino acids are disfavored despite preferring this region of phi/psi space when in isolation. PPII helices have perfect threefold rotational symmetry and within these structures we find significant correlation between the hydrophobicity of residues at i and i + 3; thus, PPII helices in globular proteins can be considered to be amphipathic.     

The effect of the polyproline II (PPII) conformation on the denatured state entropy

Polyproline II (PPII) is reported to be a dominant conformation in the unfolded state of peptides, even when no prolines are present in the sequence. Here we use isothermal titration calorimetry (ITC) to investigate the PPII bias in the unfolded state by studying the binding of the SH3 domain of SEM-5 to variants of its putative PPII peptide ligand, Sos. The experimental system is unique in that it provides direct access to the conformational entropy change of the substituted amino acids. Results indicate that the denatured ensemble can be characterized by at least two thermodynamically distinct states, the PPII conformation and an unfolded state conforming to the previously held idea of the denatured state as a random collection of conformations determined largely by hard-sphere collision. The probability of the PPII conformation in the denatured states for Ala and Gly were found to be significant, approximately 30% and approximately 10%, respectively, resulting in a dramatic reduction in the conformational entropy of folding.