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

LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana

BACKGROUND

LEA (late embryogenesis abundant) proteins have first been described about 25 years ago as accumulating late in plant seed development. They were later found in vegetative plant tissues following environmental stress and also in desiccation tolerant bacteria and invertebrates. Although they are widely assumed to play crucial roles in cellular dehydration tolerance, their physiological and biochemical functions are largely unknown.

RESULTS

We present a genome-wide analysis of LEA proteins and their encoding genes in Arabidopsis thaliana. We identified 51 LEA protein encoding genes in the Arabidopsis genome that could be classified into nine distinct groups. Expression studies were performed on all genes at different developmental stages, in different plant organs and under different stress and hormone treatments using quantitative RT-PCR. We found evidence of expression for all 51 genes. There was only little overlap between genes expressed in vegetative tissues and in seeds and expression levels were generally higher in seeds. Most genes encoding LEA proteins had abscisic acid response (ABRE) and/or low temperature response (LTRE) elements in their promoters and many genes containing the respective promoter elements were induced by abscisic acid, cold or drought. We also found that 33% of all Arabidopsis LEA protein encoding genes are arranged in tandem repeats and that 43% are part of homeologous pairs. The majority of LEA proteins were predicted to be highly hydrophilic and natively unstructured, but some were predicted to be folded.

CONCLUSION

The analyses indicate a wide range of sequence diversity, intracellular localizations, and expression patterns. The high fraction of retained duplicate genes and the inferred functional diversification indicate that they confer an evolutionary advantage for an organism under varying stressful environmental conditions. This comprehensive analysis will be an important starting point for future efforts to elucidate the functional role of these enigmatic proteins.

Time-Domain Calculations of the Infrared and Polarized Raman Spectra of Tetraalanine in Aqueous Solution

The IR and polarized (isotropic and anisotropic) Raman spectra are calculated for the amide I band of tetraalanine ((Ala)4) in aqueous solution by using a time-domain computational method, which includes the effects of the diagonal frequency modulations (of individual peptide groups), the off-diagonal (interpeptide) vibrational couplings, and structural dynamics. It is shown that the calculated band profiles, especially the existence of a large negative noncoincidence effect (i.e., large frequency separations between the IR, isotropic Raman, and anisotropic Raman bands, with the isotropic Raman being higher in frequency), are in reasonable agreement with the experimental results. This negative noncoincidence effect derives from two conditions: the positive coupling between the amide I vibrations of peptide groups and the angle larger than 90 degrees between the transition dipoles of the coupled vibrations. This result means that the dynamically changing structures mainly in the polyproline II and beta-type conformations containing some repeated interconversions obtained from the molecular dynamics calculation are consistent with the existence of a large negative noncoincidence effect, as far as the structures satisfy the above two conditions. It is also shown that the electric fields from solvent water molecules induce larger frequency shifts than those of intrachain interactions, with rapid underdamped oscillatory modulations ( approximately 43 fs) due to the librational motions of water molecules that give rise to motional narrowing effect on the spectra. The reason for the difference from the behavior seen for the O-H stretching mode of liquid water is discussed. The time-domain analysis of the mode identity shows that the system proceeds halfway to complete mode mixing with a similar time scale ( approximately 60 fs), suggesting the importance of the nonadiabatic effect, which is included in a natural way in the present computational method.

Secondary-structure analysis of denatured proteins by vacuum-ultraviolet circular dichroism spectroscopy

To elucidate the structure of denatured proteins, we measured the vacuum-ultraviolet circular dichroism (VUVCD) spectra from 260 to 172 nm of three proteins (metmyoglobin, staphylococcal nuclease, and thioredoxin) in the native and the acid-, cold-, and heat-denatured states, using a synchrotron-radiation VUVCD spectrophotometer. The circular dichroism spectra of proteins fully unfolded by guanidine hydrochloride (GdnHCl) were also measured down to 197 nm for comparison. These denatured proteins exhibited characteristic VUVCD spectra that reflected a considerable amount of residual secondary structures. The contents of alpha-helices, beta-strands, turns, poly-L-proline type II (PPII), and unordered structures were estimated for each denatured state of the three proteins using the SELCON3 program with Protein Data Bank data and the VUVCD spectra of 31 reference proteins reported in our previous study. Based on these contents, the characteristics of the four types of denaturation were discussed for each protein. In all types of denaturation, a decrease in alpha-helices was accompanied by increases in beta-strands, PPII, and unordered structures. About 20% beta-strands were present even in the proteins fully unfolded by GdnHCl in which beta-sheets should be broken. From these results, we propose that denatured proteins constitute an ensemble of residual alpha-helices and beta-sheets, partly unfolded (or distorted) alpha-helices and beta-strands, PPII, and unordered structures.

UV Raman spatially resolved melting dynamics of isotopically labeled polyalanyl peptide: slow alpha-helix melting follows 3(10)-helices and pi-bulges premelting

We used UV resonance Raman (UVRR) to examine the spatial dependence of the T-jump secondary structure relaxation of an isotopically labeled 21-residue mainly Ala peptide, AdP. The AdP penultimate Ala residues were perdeuterated, leaving the central residues hydrogenated, to allow separate monitoring of melting of the middle versus the end peptide bonds. For 5 to 30 degrees C T-jumps, the central peptide bonds show a approximately 2-fold slower relaxation time (189 +/- 31 ns) than do the exterior peptide bonds (97 +/- 15 ns). In contrast, for a 20 to 40 degrees C T-jump, the central peptide bond relaxation appears to be faster (56 +/- 6 ns) than that of the penultimate peptide bonds (131 +/- 46 ns). We show that, if the data are modeled as a two-state transition, we find that only exterior peptide bonds show anti-Arrhenius folding behavior; the middle peptide bonds show both normal Arrhenius-like folding and unfolding. This anti-Arrhenius behavior results from the involvement of pi-bulges/helices and 3(10)-helix states in the melting. The unusual temperature dependence of the (un)folding rates of the interior and exterior peptide bonds is due to the different relative (un)folding rates of 3(10)-helices, alpha-helices, and pi-bulges/helices. Pure alpha-helix unfolding rates are approximately 12-fold slower (approximately 1 micros) than that of pi-bulges and 3(10)-helices. In addition, we also find that the alpha-helix is most stable at the AdP N-terminus where eight consecutive Ala occur, whereas the three hydrophilic Arg located in the middle and at the C-terminus destabilize the alpha-helix in these regions and induce defects such as pi-bulges and 3(10)-helices.

Calmodulin signaling: analysis and prediction of a disorder-dependent molecular recognition

Calmodulin (CaM) signaling involves important, wide spread eukaryotic protein-protein interactions. The solved structures of CaM associated with several of its binding targets, the distinctive binding mechanism of CaM, and the significant trypsin sensitivity of the binding targets combine to indicate that the process of association likely involves coupled binding and folding for both CaM and its binding targets. Here, we use bioinformatics approaches to test the hypothesis that CaM-binding targets are intrinsically disordered. We developed a predictor of CaM-binding regions and estimated its performance. Per residue accuracy of this predictor reached 81%, which, in combination with a high recall/precision balance at the binding region level, suggests high predictability of CaM-binding partners. An analysis of putative CaM-binding proteins in yeast and human strongly indicates that their molecular functions are related to those of intrinsically disordered proteins. These findings add to the growing list of examples in which intrinsically disordered protein regions are indicated to provide the basis for cell signaling and regulation.

What determines the spectrum of protein native state structures?

We present a brief summary of the key factors underlying protein structure, as developed in the investigations of Pauling, Ramachandran, and Rose. We then outline a simplified physical model of proteins that focuses on geometry and symmetry. Although this model superficially appears unrelated to the detailed chemical descriptions commonly applied to proteins, we show that it captures the essential elements of the chemistry and provides a unified framework for understanding the common characteristics of folded proteins. We suggest that the spectrum of protein native state structures is determined by geometry and symmetry and the role of the sequence is to choose its native state structure from this predetermined menu.

Influence of local and residual structures on the scaling behavior and dimensions of unfolded proteins

Although recent spectroscopic studies of chemically denatured proteins hint at significant nonrandom residual structure, the results of extensive small angle X-ray scattering studies suggest random coil behavior, calling for a coherent understanding of these seemingly contradicting observations. Here, we report the results of a Monte Carlo study of the effects of two types of local structures, alpha helix and Polyproline II (PPII) helix, on the dimensions of random coil polyalanine chains viewed as a model of highly denatured proteins. We find that although Flory's power law scaling, long regarded as a signature of random coil behavior, holds for chains containing up to 90% alpha or PPII helix, the absolute magnitude of the chain dimensions is sensitive to helix content. As residual alpha helix content increases, the chain contracts until it reaches a minimum radius at approximately 70% helix, after which the chain dimensions expand rapidly. With an alpha helix content of approximately 20%, corresponding to the Ramachandran probability of being in the helical basin, experimentally observed radii of gyration are recovered. Experimental radii are similarly recovered at an alpha helix content of approximately 87%, providing an explanation for the previously puzzling experimental finding that the dimensions of the highly helical methanol-induced unfolded state are experimentally indistinguishable from those of the helix-poor urea-unfolded state. In contrast, the radius of gyration increases monotonically with increasing PPII content, and is always more expanded than the dimensions observed experimentally. These results suggest that PPII is unlikely the sole, dominant preferred conformation for unfolded proteins.

Site-Specific Hydrogen-Bonding Interaction between N-Acetylproline Amide and Protic Solvent Molecules: Comparisons of IR and VCD Measurements with MD Simulations

The effects of solute-solvent interactions on solution structures of small peptides have been paid a great deal of attention. To study the effect of hydrogen-bonding interactions on peptide solution structures, we measured the amide I IR and VCD spectra of N-acetylproline amide (AP) in various protic solvents, i.e., D2O, MeOD, EtOD, and PrOD, and directly compared them with theoretically simulated ones. The numbers of protic solvent molecules hydrogen-bonded to the two peptide bonds in the AP were quantitatively determined by carrying out the molecular dynamics (MD) simulations and then compared with the spectral analyses of the experimentally measured amide I bands. The two peptides in the AP have different propensities of forming H-bonds with protic solvent molecules, and the H-bond population distribution is found to be strongly site-specific and solvent-dependent. However, it is found that adoption of the polyproline II (PII) conformation by AP in protic solvents does not strongly depend on the hydrogen bond network-forming ability of protic solvents nor on the solvent polarity. We present a brief discussion on the validity as well as limitation of the currently available force field parameters used for the present MD simulation study.

Secondary structure provides a template for the folding of nearby polypeptides

Although protein structures are primarily encoded by their sequences, they are also critically dependent on environmental factors such as solvents and interactions with other molecules. Here we investigate how the folding-energy landscape of a short peptide is altered by interactions with another peptide, by performing atomistic replica-exchange molecular dynamics simulations of polyalanines in various environments. We analyzed the free-energy landscapes of Ala7 and Ala8 in isolation, near an alpha-helix template, and near a beta-strand template. The isolated Ala7 and Ala8 at 270 K were mainly in polyproline II helix conformations and in equilibrium between the alpha-helix and polyproline II helix, respectively, in harmony with the experiment. Interestingly, we found remarkably strong secondary-structure "templating"; namely, the alpha-helix template enhanced alpha-helix conformation and the beta-strand template induced beta-strand conformation in the simulated Ala8. The alpha-helix template lowered the nearby dielectric constant, which strengthened hydrogen bonds in the simulated Ala8, leading to alpha-helix stabilization. The beta-strand template provided hydrogen bond positions to the simulated Ala8, sharply inducing beta-strand structure. With or without templates, the energy landscape of Ala8 is always funnel-like and centered at the alpha-helix conformation, whereas entropic contribution disfavors the alpha-helix, leading to subtle competition. Secondary-structure templating may play a critical role in protein conformation dynamics in the cellular environment.

Proline and Glycine Control Protein Self-Organization into Elastomeric or Amyloid Fibrils

Elastin provides extensible tissues, including arteries and skin, with the propensity for elastic recoil, whereas amyloid fibrils are associated with tissue-degenerative diseases, such as Alzheimer's. Although both elastin-like and amyloid-like materials result from the self-organization of proteins into fibrils, the molecular basis of their differing physical properties is poorly understood. Using molecular simulations of monomeric and aggregated states, we demonstrate that elastin-like and amyloid-like peptides are separable on the basis of backbone hydration and peptide-peptide hydrogen bonding. The analysis of diverse sequences, including those of elastin, amyloids, spider silks, wheat gluten, and insect resilin, reveals a threshold in proline and glycine composition above which amyloid formation is impeded and elastomeric properties become apparent. The predictive capacity of this threshold is confirmed by the self-assembly of recombinant peptides into either amyloid or elastin-like fibrils. Our findings support a unified model of protein aggregation in which hydration and conformational disorder are fundamental requirements for elastomeric function.

A backbone-based theory of protein folding

Under physiological conditions, a protein undergoes a spontaneous disorder order transition called "folding." The protein polymer is highly flexible when unfolded but adopts its unique native, three-dimensional structure when folded. Current experimental knowledge comes primarily from thermodynamic measurements in solution or the structures of individual molecules, elucidated by either x-ray crystallography or NMR spectroscopy. From the former, we know the enthalpy, entropy, and free energy differences between the folded and unfolded forms of hundreds of proteins under a variety of solvent/cosolvent conditions. From the latter, we know the structures of approximately 35,000 proteins, which are built on scaffolds of hydrogen-bonded structural elements, alpha-helix and beta-sheet. Anfinsen showed that the amino acid sequence alone is sufficient to determine a protein's structure, but the molecular mechanism responsible for self-assembly remains an open question, probably the most fundamental open question in biochemistry. This perspective is a hybrid: partly review, partly proposal. First, we summarize key ideas regarding protein folding developed over the past half-century and culminating in the current mindset. In this view, the energetics of side-chain interactions dominate the folding process, driving the chain to self-organize under folding conditions. Next, having taken stock, we propose an alternative model that inverts the prevailing side-chain/backbone paradigm. Here, the energetics of backbone hydrogen bonds dominate the folding process, with preorganization in the unfolded state. Then, under folding conditions, the resultant fold is selected from a limited repertoire of structural possibilities, each corresponding to a distinct hydrogen-bonded arrangement of alpha-helices and/or strands of beta-sheet.

Conformational Preferences of Non-Prolyl and Prolyl Residues

The conformational study on Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in the backbone conformational preference and the cis-trans isomerization for the non-prolyl and prolyl residues in the gas phase and in the solutions (chloroform and water). For the alanine and proline dipeptides, with the increase of solvent polarity, the populations of the conformation tC with an intramolecular C(7) hydrogen bond significantly decrease, and those of the polyproline II-like conformation tF and the alpha-helical conformation tA increase, which is in good agreement with the results from circular dichroism and NMR experiments. For both the dipeptides, as the solvent polarity increases, the relative free energy of the cis conformer to the trans conformer decreases and the rotational barrier to the cis-trans isomerization increases. It is found that the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about the non-prolyl and prolyl peptide bonds in both the gas phase and the solutions. The pertinent distance d(N...H-N(NHMe)) can successfully describe the increase in the rotational barriers for the non-prolyl and prolyl trans-cis isomerization as the solvent polarity increases and the higher barriers for the non-prolyl residue than for the prolyl residue, as seen in experimental and calculated results. By analysis of the contributions to rotational barriers, the cis-trans isomerization for the non-prolyl and prolyl peptide bonds is proven to be entirely enthalpy driven in the gas phase and in the solutions. The calculated cis populations and rotational barriers to the cis-trans isomerization for both the dipeptides in chloroform and/or water accord with the experimental values.

Toward an accurate theoretical framework for describing ensembles for proteins under strongly denaturing conditions

Our focus is on an appropriate theoretical framework for describing highly denatured proteins. In high concentrations of denaturants, proteins behave like polymers in a good solvent and ensembles for denatured proteins can be modeled by ignoring all interactions except excluded volume (EV) effects. To assay conformational preferences of highly denatured proteins, we quantify a variety of properties for EV-limit ensembles of 23 two-state proteins. We find that modeled denatured proteins can be best described as follows. Average shapes are consistent with prolate ellipsoids. Ensembles are characterized by large correlated fluctuations. Sequence-specific conformational preferences are restricted to local length scales that span five to nine residues. Beyond local length scales, chain properties follow well-defined power laws that are expected for generic polymers in the EV limit. The average available volume is filled inefficiently, and cavities of all sizes are found within the interiors of denatured proteins. All properties characterized from simulated ensembles match predictions from rigorous field theories. We use our results to resolve between conflicting proposals for structure in ensembles for highly denatured states.

Analysis of molecular recognition features (MoRFs)

Several proteomic studies in the last decade revealed that many proteins are either completely disordered or possess long structurally flexible regions. Many such regions were shown to be of functional importance, often allowing a protein to interact with a large number of diverse partners. Parallel to these findings, during the last five years structural bioinformatics has produced an explosion of results regarding protein-protein interactions and their importance for cell signaling. We studied the occurrence of relatively short (10-70 residues), loosely structured protein regions within longer, largely disordered sequences that were characterized as bound to larger proteins. We call these regions molecular recognition features (MoRFs, also known as molecular recognition elements, MoREs). Interestingly, upon binding to their partner(s), MoRFs undergo disorder-to-order transitions. Thus, in our interpretation, MoRFs represent a class of disordered region that exhibits molecular recognition and binding functions. This work extends previous research showing the importance of flexibility and disorder for molecular recognition. We describe the development of a database of MoRFs derived from the RCSB Protein Data Bank and present preliminary results of bioinformatics analyses of these sequences. Based on the structure adopted upon binding, at least three basic types of MoRFs are found: alpha-MoRFs, beta-MoRFs, and iota-MoRFs, which form alpha-helices, beta-strands, and irregular secondary structure when bound, respectively. Our data suggest that functionally significant residual structure can exist in MoRF regions prior to the actual binding event. The contribution of intrinsic protein disorder to the nature and function of MoRFs has also been addressed. The results of this study will advance the understanding of protein-protein interactions and help towards the future development of useful protein-protein binding site predictors.

Evidence for polyproline II helical structure in short polyglutamine tracts

Nine neurodegenerative diseases, including Huntington's disease, are associated with the aggregation of proteins containing expanded polyglutamine sequences. The end result of polyglutamine aggregation is a beta-sheet-rich deposit. There exists evidence that an important intermediate in the aggregation process involves intramolecular beta-hairpin structures. However, little is known about the starting state, monomeric polyglutamine. Most experimental studies of monomeric polyglutamine have concluded that the backbone is completely disordered. However, such studies are hampered by the inherent tendency for polyglutamine to aggregate. A recent computational study suggested that the glutamine residues in polyglutamine tracts have a significant propensity to adopt the left-handed polyproline II (P(II)) helical conformation. In this work, we use NMR spectroscopy to demonstrate that glutamine residues possess a high propensity to adopt the P(II) conformation. We present circular dichroism spectra that indicate the presence of significant amounts of P(II) helical structure in short glutamine tracts. These data demonstrate that the propensity to adopt the P(II) structure is retained for glutamine repeats of up to at least 15 residues. Although other structures, such as alpha-helices and beta-sheets, become possible at greater lengths, our data indicate that glutamine residues in monomeric polyglutamine have a significant propensity to adopt the P(II) structure, although not necessarily in long contiguous helical stretches. We note that we have no evidence to suggest that the observed P(II) helical structure is a precursor to polyglutamine aggregation. Nonetheless, increased understanding of monomeric polyglutamine structures will aid our understanding of the aggregation process.

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by Förster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.

Polyproline II conformation is one of many local conformational states and is not an overall conformation of unfolded peptides and proteins

The alanine-based peptide Ac-XX(A)7OO-NH2, referred to as XAO (where X, A, and O denote diaminobutyric acid, alanine, and ornithine, respectively), has recently been proposed to possess a well defined polyproline II (P(II)) conformation at low temperatures. Based on the results of extensive NMR and CD investigations combined with theoretical calculations, reported here, we present evidence that, on the contrary, this peptide does not have any significant amount of organized P(II) structure but exists in an ensemble of conformations with a distorted bend in the N- and C-terminal regions. The conformational ensemble was obtained by molecular dynamics/simulated annealing calculations using the amber suite of programs with time-averaged distance and dihedral-angle restraints obtained from rotating-frame nuclear Overhauser effect (ROE) volumes and vicinal coupling constants 3J(HN Eta alpha), respectively. The computed ensemble-averaged radius of gyration Rg (7.4 +/- 1.0) A is in excellent agreement with that measured by small-angle x-ray scattering (SAXS) whereas, if the XAO peptide were in the P(II) conformation, Rg would be 11.6 A. Depending on the pH, peptide concentration, and temperature, the CD spectra of XAO do or do not possess the maximum with positive ellipticity in the 217-nm region, which is characteristic of the P(II) structure, reflecting a shifting conformational equilibrium rather than an all-or-none transition. The "P(II) conformation" should, therefore, be considered as one of the accessible conformational states of individual amino acid residues in peptides and proteins rather than as a structure of most of the chain in the early stage of folding.

The interplay of sequence and stereochemistry in defining conformation in proteins and polypeptides

Sequential specification of conformation in proteins and polypeptides is a triangular interplay involving the system of linked peptides, the sequences in side chains, and water as solvent. Stereochemistry in side chain linkages is obviously important in the interaction between all of the players, but no specification of its explicit role, if any, in linking sequence with conformation has been made. Flory and coworkers made a puzzling observation in 1967 that, when mutated from poly-L to alternating-L,D stereochemical structure, polypeptides will suffer a reduction in overall dimension or characteristic ratio by an astonishing factor of 10 and to a value even lower than that predicted for free rotation (Miller, W. G.; Brant, D. A.; Flory, P. J. J Mol Biol 1967, 23, 67-80). Enquiring into this longstanding puzzle, Durani and coworkers found that the stereochemical modification will also abolish conformational sensitivity in polypeptide structure to solvent, because electrostatic interactions in the system of linked peptides are transformed from a condition of mutual conflict to one of harmony (Ramakrishnan, V.; Ranbhor, R.; Kumar, A.; Durani, S. J Phys Chem B 2006, 110, 9314-9323). Thus, poly-L stereochemistry could be the fulcrum linking sequences with phi,psis in protein and polypeptide structures, via dielectric arbitrations in a conflicting type of interpeptide electrostatics, in agreement with the electrostatic screening model of Avbelj and Moult (Avbelj, F.; Moult, J. Biochemistry 1995, 34, 755-764).

Vibrational spectroscopic characteristics of secondary structure polypeptides in liquid water: Constrained MD simulation studies

Using the constrained MD simulation method in combination with quantum chemistry calculation, Hessian matrix reconstruction, and fragmentation approximation methods, we established a computational scheme for numerical simulations of amide I IR absorption, vibrational circular dichroism (VCD), and 2D IR photon echo spectra of peptides in solution. Six different secondary structure peptides, i.e., alpha-helix, 3(10)-helix, pi-helix, antiparallel and parallel beta-sheets, and polyproline II (P(II)), are considered, and the vibrational characteristic features in their linear and nonlinear spectra in the amide I band region are discussed. Isotope-labeling effects on IR and VCD spectra are notable only for alpha- and pi-helical peptides due to the strong vibrational couplings between two nearest neighboring amide I local oscillators. The amplitudes of difference 2D IR spectra are shown to be strongly dependent on both the extent of mode delocalization and the relative orientation of local mode transition dipoles determined by secondary structure.

Hydrogen-bonded turns in proteins: the case for a recount

Beta-turns are sites at which proteins change their overall chain direction, and they occur with high frequency in globular proteins. The Protein Data Bank has many instances of conformations that resemble beta-turns but lack the characteristic N-H(i) --> O=C(i - 3) hydrogen bond of an authentic beta-turn. Here, we identify potential hydrogen-bonded beta-turns in the coil library, a Web-accessible database utility comprised of all residues not in repetitive secondary structure, neither alpha-helix nor beta-sheet (http://www.roselab.jhu.edu/coil). In particular, candidate turns were identified as four-residue segments satisfying highly relaxed geometric criteria but lacking a strictly defined hydrogen bond. Such candidates were then subjected to a minimization protocol to determine whether slight changes in torsion angles are sufficient to shift the conformation into reference-quality geometry without deviating significantly from the original structure. This approach of applying constrained minimization to known structures reveals a substantial population of previously unidentified, stringently defined, hydrogen-bonded beta-turns. In particular, 33% of coil library residues were classified as beta-turns prior to minimization. After minimization, 45% of such residues could be classified as beta-turns, with another 8% in 3(10) helixes (which closely resemble type III beta-turns). Of the remaining coil library residues, 37% have backbone dihedral angles in left-handed polyproline II structure.

Building native protein conformation from highly approximate backbone torsion angles

Reconstructing a protein in three dimensions from its backbone torsion angles is an ongoing challenge because minor inaccuracies in these angles produce major errors in the structure. As a familiar example, a small change in an elbow angle causes a large displacement at the end of your arm, the longer the arm, the larger the displacement. Even accurate knowledge of the backbone torsions and Psi is insufficient, owing to the small, but cumulative, deviations from ideality in backbone planarity, which, if ignored, also lead to major errors in the structure. Against this background, we conducted a computational experiment to assess whether protein conformation can be determined from highly approximate backbone torsion angles, the kind of information that is now obtained readily from NMR. Specifically, backbone torsion angles were taken from proteins of known structure and mapped into 60 degrees x 60 degrees grid squares, called mesostates. Side-chain atoms beyond the beta -carbon were discarded. A mesostate representation of the protein backbone was then used to extract likely candidates from a fragment library of mesostate pentamers, followed by Monte Carlo-based fragment-assembly simulations to identify stable conformations compatible with the given mesostate sequence. Only three simple energy terms were used to gauge stability: molecular compaction, soft-sphere repulsion, and hydrogen bonding. For the six representative proteins described here, stable conformers can be partitioned into a remarkably small number of topologically distinct clusters. Among these, the native topology is found with high frequency and can be identified as the cluster with the most favorable energy.

Sterics and solvation winnow accessible conformational space for unfolded proteins

The magnitude of protein conformational space is over-estimated by the traditional random-coil model, in which local steric restrictions arise exclusively from interactions between adjacent chain neighbors. Using a five-state model, we assessed the extent to which steric hindrance and hydrogen bond satisfaction, energetically significant factors, impose additional conformational restrictions on polypeptide chains, beyond adjacent residues. Steric hindrance is repulsive: the distance of closest approach between any two atoms cannot be less than the sum of their van der Waals radii. Hydrogen bond satisfaction is attractive: polar backbone atoms must form hydrogen bonds, either intramolecularly or to solvent water. To gauge the impact of these two factors on the magnitude of conformational space, we systematically enumerated and classified the disfavored conformations that restrict short polyalanyl backbone chains. Applying such restrictions to longer chains, we derived a scaling law to estimate conformational restriction as a function of chain length. Disfavored conformations predicted by the model were tested against experimentally determined structures in the coil library, a non-helix, non-strand subset of the PDB. These disfavored conformations are usually absent from the coil library, and exceptions can be uniformly rationalized.

Dynamical binding of proline-rich peptides to their recognition domains

Recognition of proline-rich sequences plays an important role for the assembly of multi-protein complexes during the course of eukaryotic signal transduction and is mediated by a set of protein folds that share characteristic features. For many complex systems containing proline-rich sequences, multiple binding modes have been found by theoretical and/or experimental studies. In this review, we discuss the different binding modes as well as the correlated dynamics of the peptides and their recognition domains, and some implications to their biological functions. Furthermore, we give an outlook of the systems in the context of systems biology.

Unusual compactness of a polyproline type II structure

Polyproline type II (PPII) helix has emerged recently as the dominant paradigm for describing the conformation of unfolded polypeptides. However, most experimental observables used to characterize unfolded proteins typically provide only short-range, sequence-local structural information that is both time- and ensemble-averaged, giving limited detail about the long-range structure of the chain. Here, we report a study of a long-range property: the radius of gyration of an alanine-based peptide, Ace-(diaminobutyric acid)2-(Ala)7-(ornithine)2-NH2. This molecule has previously been studied as a model for the unfolded state of proteins under folding conditions and is believed to adopt a PPII fold based on short-range techniques such as NMR and CD. By using synchrotron radiation and small-angle x-ray scattering, we have determined the radius of gyration of this peptide to be 7.4 +/- 0.5 angstroms, which is significantly less than the value expected from an ideal PPII helix in solution (13.1 angstroms). To further study this contradiction, we have used molecular dynamics simulations using six variants of the AMBER force field and the GROMOS 53A6 force field. However, in all cases, the simulated ensembles underestimate the PPII content while overestimating the experimental radius of gyration. The conformational model that we propose, based on our small angle x-ray scattering results and what is known about this molecule from before, is that of a very flexible, fluctuating structure that on the level of individual residues explores a wide basin around the ideal PPII geometry but is never, or only rarely, in the ideal extended PPII helical conformation.

The structure of "unstructured" regions in peptides and proteins: role of the polyproline II helix in protein folding and recognition

Classical descriptions of the three-dimensional shapes of proteins usually invoke three main structures: alpha-helix, beta-sheet, and beta-turn. More recently, the polyproline II (PPII) structure has been implicated in diverse biological activities including signal transduction, transcription, cell motility, and immune response. Concurrently, evidence is accumulating that PPII structure has a significant role in the unfolded states of proteins. In this article, we connect the structural properties of PPII helices to their roles in protein recognition and protein unfolded states. The properties unique to the PPII conformation are linked to the exploitation of this structure for the molecular recognition of proteins, using peptide ligands of the Src homology 3 (SH3) domain as an example. The evidence supporting a role for PPII conformation in protein-unfolded states is also presented in the context of the forces that may stabilize the PPII conformation in unfolded polypeptides.

Are proteins made from a limited parts list?

Understanding the process of protein folding has been recognized as an important challenge for >70 years. It is, quintessentially, a thermodynamic problem and, arguably, thermodynamics is our most powerful discipline for understanding biological systems. Yet, despite all this, we still lack predictive understanding of protein folding. Is something missing from this picture?

A novel method reveals that solvent water favors polyproline II over beta-strand conformation in peptides and unfolded proteins: conditional hydrophobic accessible surface area (CHASA)

In aqueous solution, the ensemble of conformations sampled by peptides and unfolded proteins is largely determined by their interaction with water. It has been a long-standing goal to capture these solute-water energetics accurately and efficiently in calculations. Historically, accessible surface area (ASA) has been used to estimate these energies, but this method breaks down when applied to amphipathic peptides and proteins. Here we introduce a novel method in which hydrophobic ASA is determined after first positioning water oxygens in hydrogen-bonded orientations proximate to all accessible peptide/protein backbone N and O atoms. This conditional hydrophobic accessible surface area is termed CHASA. The CHASA method was validated by predicting the polyproline-II (P(II)) and beta-strand conformational preferences of non-proline residues in the coil library (i.e., non-alpha-helix, non-beta-strand, non-beta-turn library derived from X-ray elucidated structures). Further, the method successfully rationalizes the previously unexplained solvation energies in polyalanyl peptides and compares favorably with published experimentally determined P(II) residue propensities. We dedicate this paper to Frederic M. Richards.

Measuring unfolding of proteins in the presence of denaturant using fluorescence correlation spectroscopy

IFABP is a small (15 kDa) protein consisting mostly of antiparallel beta-strands that surround a large cavity into which ligands bind. We have previously used FCS to show that the native protein, labeled with fluorescein, exhibits dynamic fluctuation with a relaxation time of 35 micros. Here we report the use of FCS to study the unfolding of the protein induced by guanidine hydrochloride. Although the application of this technique to measure diffusion coefficients and molecular dynamics is straightforward, the FCS results need to be corrected for both viscosity and refractive index changes as the guanidine hydrochloride concentration increases. We present here a detailed study of the effects of viscosity and refractive index of guanidine hydrochloride solutions to calibrate FCS data. After correction, the increase in the diffusion time of IFABP corresponds well with the unfolding transition monitored by far ultraviolet circular dichroism. We also show that the magnitude of the 35 micros phase, reflecting the conformational fluctuation in the native state, decreases sharply as the concentration of denaturant increases and the protein unfolds. Although FCS experiments indicate that the unfolded state at pH 2 is rather compact and native-like, the radius in the presence of guanidine hydrochloride falls well within the range expected for a random coil.

Conformational prerequisites for formation of amyloid fibrils from histones

We demonstrate that bovine core histones are natively unfolded proteins in solutions with low ionic strength due to their high net positive charge at pH 7.5. Using a variety of biophysical techniques we characterized their conformation as a function of pH and ionic strength, as well as correlating the conformation with aggregation and amyloid fibril formation. Tertiary structure was absent under all conditions except at pH 7.5 and high ionic strength. The addition of trifluoroethanol or high ionic strength induced significant alpha-helical secondary structure at pH 7.5. At low pH and high salt concentration, small-angle X-ray scattering and SEC HPLC indicate the histones are present as a hexadecamer of globular subunits. The secondary structure at low pH was independent of the ionic strength or presence of TFE, as judged by FTIR. The data indicate that histones are able to adopt five different relatively stable conformations; this conformational variability probably reflects, in part, their intrinsically disordered structure. Under most of the conditions studied the histones formed amyloid fibrils with typical morphology as seen by electron microscopy. In contrast to most aggregation/amyloidogenic systems, the kinetics of fibrillation showed an inverse dependence on histone concentration; we attribute this to partitioning to a faster pathway leading to non-fibrillar self-associated aggregates at higher protein concentrations. The rate of fibril formation was maximal at low pH, and decreased to zero by pH 10. The kinetics of fibrillation were very dependent on the ionic strength, increasing with increasing salt concentration, and showing marked dependence on the nature of the ions; interestingly Gdn.HCl increased the rate of fibrillation, although much less than NaCl. Different ions also differentially affected the rate of nucleation and the rate of fibril elongation.

The polyproline II conformation in short alanine peptides is noncooperative

The finding that short alanine peptides possess a high fraction of polyproline II (PII) structure (Phi=-75 degrees, Psi=+145 degrees ) at low temperature has broad implications for unfolded states of proteins. An important question concerns whether or not this structure is locally determined or cooperative. We have monitored the conformation of alanine in a series of model peptides AcGGAnGGNH2 (n=1-3) over a temperature range from -10 degrees C to +80 degrees C. Use of 15N-labeled alanine substitutions makes it possible to measure 3JalphaN coupling constants accurately over the full temperature range. Based on a 1D next-neighbor model, the cooperative parameter sigma of PII nucleation is evaluated from the coupling constant data. The finding that sigma is close to unity (1 +/- 0.2) indicates a noncooperative role for alanine in PII structure formation, consistent with statistical surveys of the Protein Data Bank that suggest that most PII structure occurs in isolated residues. Lack of cooperativity in these models implies that hydration effects that influence PII conformation in water are highly localized. Using a nuclear Overhauser effect ratio strategy to define the alanine Psi angle, we estimate that, at 40 degrees C, the time-averaged alanine conformation (Phi=-80 degrees, Psi=+170 degrees ) deviates from canonical PII structure, indicating that PII melts at high temperature. Thus, the high-temperature state of short alanine peptides seems to be an unfolded ensemble with higher distribution in the extended beta structure basin, but not a coil.

Unified perspective on proteins: a physics approach

We study a physical system which, while devoid of the complexity one usually associates with proteins, nevertheless displays a remarkable array of proteinlike properties. The constructive hypothesis that this striking resemblance is not accidental not only leads to a unified framework for understanding protein folding, amyloid formation, and protein interactions but also has implications for natural selection.

Revisiting the Ramachandran plot: hard-sphere repulsion, electrostatics, and H-bonding in the alpha-helix

What determines the shape of the allowed regions in the Ramachandran plot? Although Ramachandran explained these regions in terms of 1-4 hard-sphere repulsions, there are discrepancies with the data where, in particular, the alphaR, alphaL, and beta-strand regions are diagonal. The alphaR-region also varies along the alpha-helix where it is constrained at the center and the amino terminus but diffuse at the carboxyl terminus. By analyzing a high-resolution database of protein structures, we find that certain 1-4 hard-sphere repulsions in the standard steric map of Ramachandran do not affect the statistical distributions. By ignoring these steric clashes (NH(i+1) and O(i-1)C), we identify a revised set of steric clashes (CbetaO, O(i-1)N(i+1), CbetaN(i+1), O(i-1)Cbeta, and O(i-1)O) that produce a better match with the data. We also find that the strictly forbidden region in the Ramachandran plot is excluded by multiple steric clashes, whereas the outlier region is excluded by only one significant steric clash. However, steric clashes alone do not account for the diagonal regions. Using electrostatics to analyze the conformational dependence of specific interatomic interactions, we find that the diagonal shape of the alphaR and alphaL-regions also depends on the optimization of the NH(i+1) and O(i-1)C interactions, and the diagonal beta-strand region is due to the alignment of the CO and NH dipoles. Finally, we reproduce the variation of the Ramachandran plot along the alpha-helix in a simple model that uses only H-bonding constraints. This allows us to rationalize the difference between the amino terminus and the carboxyl terminus of the alpha-helix in terms of backbone entropy.

Reassessing random-coil statistics in unfolded proteins

The Gaussian-distributed random coil has been the dominant model for denatured proteins since the 1950s, and it has long been interpreted to mean that proteins are featureless, statistical coils in 6 M guanidinium chloride. Here, we demonstrate that random-coil statistics are not a unique signature of featureless polymers. The random-coil model does predict the experimentally determined coil dimensions of denatured proteins successfully. Yet, other equally convincing experiments have shown that denatured proteins are biased toward specific conformations, in apparent conflict with the random-coil model. We seek to resolve this paradox by introducing a contrived counterexample in which largely native protein ensembles nevertheless exhibit random-coil characteristics. Specifically, proteins of known structure were used to generate disordered conformers by varying backbone torsion angles at random for approximately 8% of the residues; the remaining approximately 92% of the residues remained fixed in their native conformation. Ensembles of these disordered structures were generated for 33 proteins by using a torsion-angle Monte Carlo algorithm with hard-sphere sterics; bulk statistics were then calculated for each ensemble. Despite this extreme degree of imposed internal structure, these ensembles have end-to-end distances and mean radii of gyration that agree well with random-coil expectations in all but two cases.

Origin of the neighboring residue effect on peptide backbone conformation

Unfolded peptides in water have some residual structure that may be important in the folding process, and the nature of the residual structure is currently of much interest. There is a neighboring residue effect on backbone conformation, discovered in 1997 from measurements of (3)J(HN alpha) coupling constants. The neighboring residue effect appears also in the "coil library" of Protein Data Bank structures of residues not in alpha-helix and not in beta-structure. When a neighboring residue (i - 1 or i + 1) belongs to class L (aromatic and beta-branched amino acids, FHITVWY) rather than class S (all others, G and P excluded), then the backbone angle of residue i is more negative for essentially all amino acids. Calculated values of peptide solvation (electrostatic solvation free energy, ESF) predict basic properties of the neighboring residue effect. We show that L amino acids reduce the solvation of neighboring peptide groups more than S amino acids. When tripeptides from the coil library are excised to allow solvation, the central residues have more negative values of but less negative values of with L than with S neighbors. The coil library values of <(3)J(HN alpha)>, which vary strikingly among the amino acids, are correlated with the neighboring residue effect seen by ESF. Moreover, values for the "blocking effect" of side chains on the hydrogen exchange rates of peptide NH protons are correlated with ESF values.

Preferred peptide backbone conformations in the unfolded state revealed by the structure analysis of alanine-based (AXA) tripeptides in aqueous solution

We have combined Fourier transform IR, polarized Raman spectroscopy, and vibrational CD measurements of the amide I' band profile of alanyl-X-alanine tripeptides in (2)H(2)O to obtain the dihedral angles of their central amino acid residue. X represents glycine, valine, methionine, histidine, serine, proline, lysine, leucine, tryptophan, tyrosine, and phenylalanine. The experimental data were analyzed by means of a recently developed algorithm, which exploits the excitonic coupling between the amide modes of the two peptide groups. The results were checked by measuring the respective electronic CD spectra. The investigated peptides can be sorted into three classes. Valine, phenylalanine, tryptophan, histidine, and serine predominantly adopt an extended beta-strand conformation. Cationic lysine and proline prefer a polyproline II-like structure. Alanine, methionine, glycine, and leucine populate these two conformations with comparable probability. Our results are in variance with the prediction of the random-coil model, but supportive of Flory's isolated-pair hypothesis. We combined the obtained structural propensities of the investigated residues and similar information about other residues in the literature (i.e., glutamate, aspartate, isoleucine, and glutamine) to predict possible conformations of the monomeric amyloid beta peptide A beta(1-42) in aqueous solution, which reproduces results from most recent spectroscopic studies. Thus, it is demonstrated that the unfolded state of peptides can be understood in terms of the intrinsic structural propensities of their amino acid residues.

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: alpha-helix, antiparallel beta-strand, parallel beta-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), beta-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.

Conformational constraints for amyloid fibrillation: the importance of being unfolded

Recent reports give strong support to the idea that amyloid fibril formation and the subsequent development of protein deposition diseases originate from conformational changes in corresponding amyloidogenic proteins. In this review, recent findings are surveyed to illustrate that protein fibrillogenesis requires a partially folded conformation. This amyloidogenic conformation is relatively unfolded, and shares many structural properties with the pre-molten globule state, a partially folded intermediate frequently observed in the early stages of protein folding and under some equilibrium conditions. The inherent flexibility of such an intermediate is essential in allowing the conformational rearrangements necessary to form the core cross-beta structure of the amyloid fibril.

The conformation of tetraalanine in water determined by polarized Raman, FT-IR, and VCD spectroscopy

The present article reports the conformation of cationic tetraalanine in aqueous solution. The determination of the dihedral angles of the two central amino acid residues was achieved by analyzing the amide I' band profile in the respective polarized visible Raman, Fourier transform-IR, and vibrational circular dichroism (VCD) spectra by means of a novel algorithm which utilizes the excitonic coupling between the amide I modes of nearest neighbor and second nearest peptide groups. It is an extension of a recently developed theory (Schweitzer-Stenner, R. Biophys. J., 2002, 83, 523-532). UV electronic circular dichroism (ECD) spectra of the peptides were used to validate the results of the structure analysis. The analyses yielded the dihedral angles (phi(12), psi(12)) = (-70 degrees, 155 degrees ) and (phi(23), psi(23)) = (-80 degrees, 145 degrees ). The obtained values are very close to the Ramachandran coordinates of the polyproline II helix (PPII). The data suggest that this is the conformation predominantly adopted by the peptide at room temperature. This notion was corroborated by the corresponding electronic circular dichroism spectrum. Tetraalanine exhibits a higher propensity for PPII than trialanine for which a 50:50 mixture of polyproline II and an extended beta-strand-like conformation was obtained from recent spectroscopic studies (Eker et al., J. Am. Chem. Soc. 2002, 124, 14330-14341). The temperature dependence of the CD spectra rule out that any cooperativity is involved in the strand if PPII transition. This led to the conclusion that solvent-peptide interactions give rise to the observed PPII stability. Our result can be utilized to understand why the denaturation of helix-forming peptides generally yields a PPII rather than a heterogeneous random conformation.

A Natively Unfolded Toxin Domain Uses Its Receptor as a Folding Template

Natively unfolded proteins range from molten globules to disordered coils. They are abundant in eukaryotic genomes and commonly involved in molecular interactions. The essential N-terminal translocation domains of colicin toxins from Escherichia coli are disordered bacterial proteins that bind at least one protein of the Tol or Ton family. The colicin N translocation domain (ColN-(1-90)), which binds to the C-terminal domain of TolA (TolA-(296-421)), shows a disordered far-UV CD spectrum, no near-UV CD signal, and non-cooperative thermal unfolding. As expected, TolA-(296-421) displays both secondary structure in far-UV CD and tertiary structure in near-UV CD. Furthermore it shows a cooperative unfolding transition at 65 degrees C. CD spectra of the 1:1 complex show both increased secondary structure and colicin N-specific near-UV CD signals. A new cooperative thermal transition at 35 degrees C is followed by the unchanged unfolding behavior of TolA-(296-421). Fluorescence and surface plasmon resonance confirm that the new unfolding transition accompanies dissociation of ColN-(1-90). Hence upon binding the disordered structure of ColN-(1-90) converts to a cooperatively folded domain without altering the TolA-(296-421) structure.