A simple model for polyproline II structure in unfolded states of alanine-based peptides

The striking similarity between observed circular dichroism spectra of nonprolyl homopolymers and that of regular left-handed polyproline II (P_II) helices prompted Tiffany and Krimm to propose in 1968 that unordered peptides and unfolded proteins are built of P_II segments linked by sharp bends. A large body of experimental evidence, accumulated over the past three decades, provides compelling evidence in support of the original hypothesis of Tiffany and Krimm. Of particular interest are the recent experiments of Shi et al. who find significant P_II structure in a short unfolded alanine-based peptide. What is the physical basis for P_II helices in peptide and protein unfolded states? The widely accepted view is that favorable chain-solvent hydrogen bonds lead to a preference for dynamical fluctuations about noncooperative P_II helices in water. Is this preference simply a consequence of hydrogen bonding or is it a manifestation of a more general trend for unfolded states which are appropriately viewed as chains in a good solvent? The prevalence of closely packed interiors in folded proteins suggests that under conditions that favor folding, water—which is a better solvent for itself than for any polypeptide chain—expels the chain from its midst, thereby maximizing chain packing. Implicit in this view is a complementary idea: under conditions that favor unfolding, chain-solvent interactions are preferred and in a so-called good solvent, chain packing density is minimized. In this work we show that minimization of chain packing density leads to preferred fluctuations for short polyalanyl chains around canonical, noncooperative P_II-like conformations. Minimization of chain packing is modeled using a purely repulsive soft-core potential between polypeptide atoms. Details of chain-solvent interactions are ignored. Remarkably, the simple model captures the essential physics behind the preference of short unfolded alanine-based peptides for P_II helices. Our results are based on a detailed analysis of the potential energy landscape which determines the system''s structural and thermodynamic preferences. We use the inherent structure formalism of Stillinger and Weber, according to which the energy landscape is partitioned into basins of attraction around local minima. We find that the landscape for the experimentally studied seven-residue alanine-based peptide is dominated by fluctuations about two noncooperative structures: the left-handed polyproline II helix and its symmetry mate.     

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.     

A comprehensive library of blocked dipeptides reveals intrinsic backbone conformational propensities of unfolded proteins

Despite prolonged scientific efforts to elucidate the intrinsic peptide backbone preferences of amino-acids based on understanding of intermolecular forces, many open questions remain, particularly concerning neighboring peptide interaction effects on the backbone conformational distribution of short peptides and unfolded proteins. Here, we show that spectroscopic studies of a complete library of 400 dipeptides reveal that, irrespective of side-chain properties, the backbone conformation distribution is narrow and they adopt polyproline II and β-strand, indicating the importance of backbone peptide solvation and electronic effects. By directly comparing the dipeptide circular dichroism and NMR results with those of unfolded proteins, the comprehensive dipeptides form a complete set of structural motifs of unfolded proteins. We thus anticipate that the present dipeptide library with spectroscopic data can serve as a useful database for understanding the nature of unfolded protein structures and for further refinements of molecular mechanical parameters.     

A tetrapeptide-based method for polyproline II-type secondary structure prediction

We describe a new method for polyproline II-type (PPII) secondary structure prediction based on tetrapeptide conformation properties using data obtained from all globular proteins in the Protein Data Bank (PDB). This is the first method for PPII prediction with a relatively high level of accuracy (approximately 60%). Our method uses only frequencies of different conformations among oligopeptides without any additional parameters. We also attempted to predict alpha-helices and beta-strands using the same approach. We find that the application of our method reveals interrelation between sequence and structure even for very short oligopeptides (tetrapeptides).     

Left-handed polyproline II helix formation is (very) locally driven

The left-handed polyproline II helix (PPII) is believed to be the preferred conformation for proline-rich regions of sequence in proteins. Such regions have been postulated to be protein-protein interaction domains. The formation of this structure is studied here using simple Monte Carlo computer simulations employing the hard sphere potential. It is found that polyproline sequences adopt only the PPII structure in the simulations. Non-proline, non-glycine residues inserted as guests into polyproline host peptides are conformationally restricted by the following proline residues and tend to be part of the PPII helix. It is found through insertion of two alanine residues into polyproline that the PPII structure is not propagated through more than one non-proline residue. This finding calls into question the hypothesis that proline-rich regions will preferentially adopt this structure since many such sequences are comprised of less than 50% proline residues. Proteins 33:218–226, 1998. © 1998 Wiley-Liss, Inc.     

Exploring the impact of polyproline II (PII) conformational bias on the binding of peptides to the SEM-5 SH3 domain

The left-handed polyproline II helical structure (P(II)) is observed to be a dominant conformation in the disordered states of protein and small polypeptide chains, even when no prolines are present in the sequence. Recently, in work by Ferreon and Hilser, the energetics associated with Ala and Gly substitutions at a surface exposed proline site were determined calorimetrically by measuring the binding energetics of Sos peptide variants to the C-terminal Src Homology 3 domain of SEM-5. The results were interpreted as a significant conformational bias toward the bound conformation (i.e., P(II)), even when the ligand is unbound. That study was not able to determine, however, whether the conformational bias of the peptides could be explained in terms other than that of a P(II) preference. Here, we test, using a computer algorithm based on the hard sphere collision (HSC) model, the notion of whether a bias in the unbound states of the peptide ligands is specific for the P(II) conformation, or if a bias to any other region of (phi, psi) space can also result in the same observed binding energetics. The results of these computer simulations indicate that, of the regions of (phi, psi) modeled for bias in the small peptides, only the bias to the P(II) conformation, and at rates of bias similar to the experimentally observed rates, quantitatively reproduced the experimental binding energetics.     

Unfolded state of polyalanine is a segmented polyproline II helix

Definition of the unfolded state of proteins is essential for understanding their stability and folding on biological timescales. Here, we find that under near physiological conditions the configurational ensemble of the unfolded state of the simplest protein structure, polyalanine alpha-helix, cannot be described by the commonly used Flory random coil model, in which configurational probabilities are derived from conformational preferences of individual residues. We utilize novel effectively ergodic sampling algorithms in the presence of explicit aqueous solvation, and observe water-mediated formation of polyproline II helical (P(II)) structure in the natively unfolded state of polyalanine, and its facilitation of alpha-helix formation in longer peptides. The segmented P(II) helical coil preorganizes the unfolded state ensemble for folding pathway entry by reducing the conformational space available to the diffusive search. Thus, as much as half of the folding search in polyalanine is accomplished by preorganization of the unfolded state.     

Native-like mean structure in the unfolded ensemble of small proteins

The nature of the unfolded state plays a great role in our understanding of proteins. However, accurately studying the unfolded state with computer simulation is difficult, due to its complexity and the great deal of sampling required. Using a supercluster of over 10,000 processors we have performed close to 800 micros of molecular dynamics simulation in atomistic detail of the folded and unfolded states of three polypeptides from a range of structural classes: the all-alpha villin headpiece molecule, the beta hairpin tryptophan zipper, and a designed alpha-beta zinc finger mimic. A comparison between the folded and the unfolded ensembles reveals that, even though virtually none of the individual members of the unfolded ensemble exhibits native-like features, the mean unfolded structure (averaged over the entire unfolded ensemble) has a native-like geometry. This suggests several novel implications for protein folding and structure prediction as well as new interpretations for experiments which find structure in ensemble-averaged measurements.     

Distinct circular dichroism spectroscopic signatures of polyproline II and unordered secondary structures: Applications in secondary structure analyses

Circular dichroism (CD) spectroscopy is a valuable method for defining canonical secondary structure contents of proteins based on empirically‐defined spectroscopic signatures derived from proteins with known three‐dimensional structures. Many proteins identified as being “Intrinsically Disordered Proteins” have a significant amount of their structure that is neither sheet, helix, nor turn; this type of structure is often classified by CD as “other”, “random coil”, “unordered”, or “disordered”. However the “other” category can also include polyproline II (PPII)‐type structures, whose spectral properties have not been well‐distinguished from those of unordered structures. In this study, synchrotron radiation circular dichroism spectroscopy was used to investigate the spectral properties of collagen and polyproline, which both contain PPII‐type structures. Their native spectra were compared as representatives of PPII structures. In addition, their spectra before and after treatment with various conditions to produce unfolded or denatured structures were also compared, with the aim of defining the differences between CD spectra of PPII and disordered structures. We conclude that the spectral features of collagen are more appropriate than those of polyproline for use as the representative spectrum for PPII structures present in typical amino acid‐containing proteins, and that the single most characteristic spectroscopic feature distinguishing a PPII structure from a disordered structure is the presence of a positive peak around 220nm in the former but not in the latter. These spectra are now available for inclusion in new reference data sets used for CD analyses of the secondary structures of soluble proteins.     

The synaptic acetylcholinesterase tetramer assembles around a polyproline II helix

Functional localization of acetylcholinesterase (AChE) in vertebrate muscle and brain depends on interaction of the tryptophan amphiphilic tetramerization (WAT) sequence, at the C-terminus of its major splice variant (T), with a proline-rich attachment domain (PRAD), of the anchoring proteins, collagenous (ColQ) and proline-rich membrane anchor. The crystal structure of the WAT/PRAD complex reveals a novel supercoil structure in which four parallel WAT chains form a left-handed superhelix around an antiparallel left-handed PRAD helix resembling polyproline II. The WAT coiled coils possess a WWW motif making repetitive hydrophobic stacking and hydrogen-bond interactions with the PRAD. The WAT chains are related by an approximately 4-fold screw axis around the PRAD. Each WAT makes similar but unique interactions, consistent with an asymmetric pattern of disulfide linkages between the AChE tetramer subunits and ColQ. The P59Q mutation in ColQ, which causes congenital endplate AChE deficiency, and is located within the PRAD, disrupts crucial WAT-WAT and WAT-PRAD interactions. A model is proposed for the synaptic AChE(T) tetramer.     

Evolutionary conservation of the polyproline II conformation surrounding intrinsically disordered phosphorylation sites

Intrinsically disordered (ID) proteins function in the absence of a unique stable structure and appear to challenge the classic structure-function paradigm. The extent to which ID proteins take advantage of subtle conformational biases to perform functions, and whether signals for such mechanism can be identified in proteome-wide studies is not well understood. Of particular interest is the polyproline II (PII) conformation, suggested to be highly populated in unfolded proteins. We experimentally determine a complete calorimetric propensity scale for the PII conformation. Projection of the scale into representative eukaryotic proteomes reveals significant PII bias in regions coding for ID proteins. Importantly, enrichment of PII in ID proteins, or protein segments, is also captured by other PII scales, indicating that this enrichment is robustly encoded and universally detectable regardless of the method of PII propensity determination. Gene ontology (GO) terms obtained using our PII scale and other scales demonstrate a consensus for molecular functions performed by high PII proteins across the proteome. Perhaps the most striking result of the GO analysis is conserved enrichment (P < 10⁻⁸) of phosphorylation sites in high PII regions found by all PII scales. Subsequent conformational analysis reveals a phosphorylation-dependent modulation of PII, suggestive of a conserved “tunability” within these regions. In summary, the application of an experimentally determined polyproline II (PII) propensity scale to proteome-wide sequence analysis and gene ontology reveals an enrichment of PII bias near disordered phosphorylation sites that is conserved throughout eukaryotes.     

Spectroscopic evidence revealing polyproline II structure in hydrophobic, putatively elastomeric sequences encoded by specific exons of human tropoelastin

Elastin, the protein responsible for tissue elasticity, is contained in arterial walls, lungs, and skin. Given the cassette like organization of the human tropoelastin gene, giving rise to alternating exons encoding for crosslink domains and elastomeric domains, it is tempting to suggest that polypeptides encoded by different exons could adopt structures independent of the other exons. The results obtained with the polypeptide sequences encoded by exons 3, 7, and 30 of human tropoelastin are described. It is shown that these hydrophobic exons may partly assume the polyproline II (PPII) structure, as found by circular dichroism studies in aqueous solution. Classical Raman spectroscopy evidences a specific sharp band at 1314 cm(-1), which is assigned to the PPII structure adopted by these exons in the solid state. As these sequences are among those putatively responsible for elastomeric properties, these findings are of particular interest in relation to the current models of the elasticity of elastin.     

Direct characterization of the folded, unfolded and urea-denatured states of the C-terminal domain of the ribosomal protein L9

The stability of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) is strongly dependent upon pH. Below pH 4.2, the folded and unfolded states are both populated significantly. Their interconversion is slow on the NMR chemical shift time-scale and separate, well-resolved resonances from each state are observed. This allows the hydrodynamic properties of both states to be studied under identical conditions by using pulse field gradient NMR experiments. Hydrodynamic radii of the folded, unfolded and urea denatured protein molecules at pD 3.8 have been derived. The acid-denatured protein has a significantly smaller hydrodynamic radius, 28.2A, compared to that of the urea-denatured protein, which is 33.6A at pD 3.8. Far-UV CD spectra show that there is more residual secondary structure retained in the acid-denatured ensemble than in the urea-denatured one. ANS binding experiments and analysis of the CD data show that this acid-denatured species is not a molten globule state. Diffusion measurements of CTL9 were conducted over the pD range from 2.1 to 7.0. The hydrodynamic radii of both the folded and the acid-unfolded protein start to increase below pD 4, with the radius of hydration of the acid-unfolded state increasing from 25.1A at pD 4.2 to 33.5A at pD 2.1. The hydrodynamic radius of the urea-denatured protein is much less sensitive to pH. The unfolded protein at pD 2.1, no urea, has almost the same hydrodynamic radius as the urea-denatured protein at pD 3.8. The CD spectra, however, show significant differences in residual secondary structure, and the acid-denatured state contains more structure.     

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.     

Polyproline II helix is the preferred conformation for unfolded polyalanine in water

Does aqueous solvent discriminate among peptide conformers? To address this question, we computed the solvation free energy of a blocked, 12‐residue polyalanyl‐peptide in explicit water and analyzed its solvent structure. The peptide was modeled in each of 4 conformers: α‐helix, antiparallel β‐strand, parallel β‐strand, and polyproline II helix (P_II). Monte Carlo simulations in the canonical ensemble were performed at 300 K using the CHARMM 22 forcefield with TIP3P water. The simulations indicate that the solvation free energy of P_II is favored over that of other conformers for reasons that defy conventional explanation. Specifically, in these 4 conformers, an almost perfect correlation is found between a residue's solvent‐accessible surface area and the volume of its first solvent shell, but neither quantity is correlated with the observed differences in solvation free energy. Instead, solvation free energy tracks with the interaction energy between the peptide and its first‐shell water. An additional, previously unrecognized contribution involves the conformation‐dependent perturbation of first‐shell solvent organization. Unlike P_II, β‐strands induce formation of entropically disfavored peptide:water bridges that order vicinal water in a manner reminiscent of the hydrophobic effect. The use of explicit water allows us to capture and characterize these dynamic water bridges that form and dissolve during our simulations. Proteins 2004. © 2004 Wiley‐Liss, Inc.     

Cooperative interactions and a non-native buried Trp in the unfolded state of an SH3 domain

The presence of residual structure in the unfolded state of the N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) has been investigated using far- and near-UV circular dichroism (CD), fluorescence, and NMR spectroscopy. The unfolded (U(exch)) state of the drkN SH3 domain is significantly populated and exists in equilibrium with the folded (F(exch)) state under non-denaturing conditions near physiological pH. Denaturation experiments have been performed on the drkN SH3 domain in order to monitor the change in ellipticity, fluorescence intensity, and chemical shift between the U(exch) state and chemically or thermally denatured states. Differences between the unfolded and chemically or thermally denatured states highlight specific areas of residual structure in the unfolded state that are cooperatively disrupted upon denaturation. Results provide evidence for cooperative interactions in the unfolded state involving residues of the central beta-sheet, particularly the beta4 strand. Denaturation as well as hydrogen-exchange experiments demonstrate a non-native burial of the Trp ring within this "cooperative" core of the unfolded state. These findings support the presence of non-native hydrophobic clusters, organised by Trp rings, within disordered states.     

Transient structure formation in unfolded acyl-coenzyme A-binding protein observed by site-directed spin labelling

Paramagnetic relaxation has been used to monitor the formation of structure in the folding peptide chain of guanidinium chloride-denatured acyl-coenzyme A-binding protein. The spin label (1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-methyl)methanesulfonate (MTSL) was covalently bound to a single cysteine residue introduced into five different positions in the amino acid sequence. It was shown that the formation of structure in the folding peptide chain at conditions where 95% of the sample is unfolded brings the relaxation probe close to a wide range of residues in the peptide chain, which are not affected in the native folded structure. It is suggested that the experiment is recording the formation of many discrete and transient structures in the polypeptide chain in the preface of protein folding. Analysis of secondary chemical shifts shows a high propensity for alpha-helix formation in the C-terminal part of the polypeptide chain, which forms an alpha-helix in the native structure and a high propensity for turn formation in two regions of the polypeptide that form turns in the native structure. The results contribute to the idea that native-like structural elements form transiently in the unfolded state, and that these may be of importance to the initiation of protein folding.     

Stereoelectronic effects on the transition barrier of polyproline conformational interconversion

There has been growing interest in polyproline type II (PPII) helices since PPII helices have been found in folded and unfolded proteins and involved in a variety of biological activities. Polyproline can also form type I helices (PPI) which are very different from PPII conformation and only exist in certain organic solvents. Recent studies have shown that stereoelectronic effects play a critical role in stabilizing a PPI or PPII helix. Here, we have synthesized a series of host–guest peptides with an electron‐withdrawing substituent at the 4R or 4S position of proline and used a kinetic approach to further explore stereoelectronic effects on the transition barrier of the interconversion between PPI and PPII conformations. Time‐dependent circular dichroism measurements revealed that the rates of PPII → PPI conversion were reduced upon incorporating the hydroxyl‐, fluoro‐, and methoxy‐groups at the 4R position while the rates would be increased if these substituents were at the 4S position. We quantified the changes in transition free energy by comparing their rate constants. (4R,2S)‐4‐Fluoroproline and (4S,2S)‐4‐fluoroproline have the largest effect on the transition energy barrier for PPII → PPI conversion. Our results provide important insights into the role of stereoelectronic effects on the PPII → PPI transition state barrier, which has not been reported in past thermodynamic studies.     

NMR identification of left-handed polyproline type II helices

NMR characteristics of a model left-handed 3(1)-helical peptide are reported in this study. With temperature and sequence corrections on the predicted random coil (15)N chemical shifts, a significant (15)N chemical shift deviation is observed for the model 3(1) peptide. The (15)N chemical shift differences also correlate well with the molar ellipticities (at 220 nm) of the CD spectra at different temperatures, indicating that the (15)N chemical shift is a sensitive probe for 3(1)-helices. The average (3)J(HNalpha) and (1)J(CalphaHalpha) values of the model peptide are determined to be 6.5 and 142.6 Hz, respectively, which are consistent with the values calculated from the geometry of 3(1)-helices. With careful measurements of amide (15)N chemical shifts and incorporating temperature and sequence effect corrections, the (15)N chemical shifts can be used together with (3)J(HNalpha) and (1)J(CalphaHalpha) to differentiate 3(1)-helices from random coils with high confidence. Based on the observed NMR characteristics, a strategy is developed for probing left-handed 3(1)-helical structures from other secondary structures.