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.     

Simulated folding in polypeptides of diversified molecular tacticity: implications for protein folding and de novo design

Stereochemistry could be a powerful variable for conformational tune up of polypeptides for de novo design. It may be also useful probe of possible role of interamide energetics in selection and stabilization of conformation. The homopolypeptides Ac-Xxx30-NHMe, with Xxx = Ala, Val, and Leu, of diversified stereochemical structure are generated by simulated racemization with a modified GROMOS-96 force field. The polypeptides, and other systematic stereochemical variants, are folded by simulated annealing with another modified GROMOS-96 force field under the dielectric constant values 1, 4, and 10. The resultant 15,000 molecular folds of isotactic (poly-L-chiral), syndiotactic (alternating L,D-chiral), and heterotactic (random-L,D-chiral) stereochemical structure, belonging to three polypeptide series, achieved under three different folding conditions, are assessed statistically for structure-to-energy-to-conformation relationship. The results suggest that interamide electrostatics could be a major factor in secondary-structure selection in polypeptides while main-chain stereochemistry could dictate molecular packing and therefore the relative magnitude of hydrogen-bond and Lennard-Jones (LJ) contributions in conformational energy. A method for computational design of heterotactic molecular folds in polypeptide structure has been developed, and the first road map for a chiral tune up of polypeptide structure based on stereochemical engineering has been laid down. Broad implications for protein structure, folding, and de novo design are briefly discussed.     

Conformational strain in the hydrophobic core and its implications for protein folding and design

We have designed de novo 13 divergent spectrin SH3 core sequences to determine their folding properties. Kinetic analysis of the variants with stability similar to that of the wild type protein shows accelerated unfolding and refolding rates compatible with a preferential stabilization of the transition state. This is most likely caused by conformational strain in the native state, as deletion of a methyl group (Ile-->Val) leads to deceleration in unfolding and increased stability (up to 2 kcal x mol(-1)). Several of these Ile-->Val mutants have negative phi(-U) values, indicating that some noncanonical phi(-U) values might result from conformational strain. Thus, producing a stable protein does not necessarily mean that the design process has been entirely successful. Strained interactions could have been introduced, and a reduction in the buried volume could result in a large increase in stability and a reduction in unfolding rates.     

Finding the lowest free energy conformation of a protein is an NP-hard problem: Proof and implications

The protein folding problem and the notion of NP-completeness and NP-hardness are discussed. A lattice model is suggested to capture the essece of protein folding. For this model we present a proof that finding the lowest free energy conformation belongs to the class of NP-hard problems. The implications of the proof are discussed and we suggest that the natural folding process cannot be considered as a search for the global free energy minimum. However, we suggest an explanation as to why, for many proteins, the native functional conformation maycoincide with the lowest free energy conformation.     

The implications of gene heterozygosity for protein folding and protein turnover

The offspring of closely related parents often suffer from inbreeding depression, sometimes resulting in a slower growth rate for inbred offspring relative to non-inbred offspring. Previous research has shown that some of the slower growth rate of inbred organisms can be attributed to the inbred organisms’ increased levels of protein turnover. This paper attempts to show that the higher levels of protein turnover among inbred organisms can be attributed to accumulations of misfolded and aggregated proteins that require degradation by the inbred organisms’ protein quality control systems. The accumulation of misfolded and aggregated proteins within inbred organisms are the result of more negative free energies of folding for proteins encoded at homozygous gene loci and higher concentrations of potentially aggregating non-native protein species within the cell. The theory presented here makes several quantitative predictions that suggest a connection between protein misfolding/aggregation and polyploidy that can be tested by future research.     

Extracting information on folding from the amino acid sequence: accurate predictions for protein regions with preferred conformation in the absence of tertiary interactions.

A recently developed procedure to predict backbone structure from the amino acid sequence [Rooman, M., Kocher, J. P., & Wodak, S. (1991) J. Mol. Biol, 221, 961-979] is fine tuned to identify protein segments, of length 5-15 residues, that adopt well-defined conformations in the absence of tertiary interactions. These segments are obtained by requiring that their predicted lowest energy structures have a sizable energy gap relative to other computed conformations. Applying this procedure to 69 proteins of known structure, we find that regions with largest energy gaps--those having highly preferred conformations--are also the most accurately predicted ones. On the basis of previous findings that such regions correlate well with sites that become structured early during folding, our approach provides the means of identifying such sites in proteins without prior knowledge of the tertiary structure. Furthermore, when predictions are performed so as to ignore the influence of residues flanking each segment along the sequence, a situation akin to excising the considered peptide from the rest of the chain, they offer the possibility of identifying protein segments liable to adopt well-defined conformations on their own. The described approach should have useful applications in experimental and theoretical investigations of protein folding and stability, and aid in designing peptide drugs and vaccines.     

Gel electrophoresis in studies of protein conformation and folding

Electrophoresis through polyacrylamide gels is a useful method for distinguishing conformational states of proteins and analyzing the thermodynamic and kinetic properties of transitions between conformations. Although the relationship between protein conformation and electrophoretic mobility is quite complex, relative mobilities provide qualitative estimates of compactness. Conformational states which interconvert slowly on the time scale of the electrophoretic separation can often be resolved, and the rates of interconversion can be estimated. If the transitions are more rapid, then the electrophoretic mobility represents the equilibrium distribution of conformations. Protein unfolding transitions induced by urea are readily studied using slab gels containing a gradient of urea concentration perpendicular to the direction of electrophoresis. Protein applied across the top of such a gel migrates in the presence of continuously varying urea concentrations, and a profile of the unfolding transition is generated directly. Transitions induced by other agents could be studied using analogous gradient gels. Electrophoretic methods are especially suited for studying small quantities of protein, and complex mixtures, since the different components can be separated during the electrophoresis.     

Residue contacts in protein structures and implications for protein folding

The preferential association of amino acid side groups with specific side chain atoms are examined in 44 known protein structures. The resulting association potentials among residue side groups are used to detect structural homology in proteins displaying little or no homology in their primary sequences. Suggestions are also made regarding the nature of the protein folding process. They are based on statistical observations that delineate the extent of short and long range interactions and that display side group bias in association with other side chain atoms on their N-terminal side.     

Snapshots of protein folding problem: implications of folding and misfolding studies

Deciphering the code that determines the three-dimensional structure of proteins and the ability to predict the final folded form of a protein is still elusive to molecular biophysists. In the case of several proteins a similar tertiary structure is not accompanied by any significant sequence similarity. The question now remains whether a code beyond the genetic code that describes the arrangement of the amino acid within a three dimensional protein structure. The available data undoubtedly demonstrates that the redundancy of this code must be tremendous. Several techniques such as nuclear magnetic resonance spectroscopy and laser detection techniques, coupled with fast initiation of the folding reaction, can now probe the folding events in milliseconds or even faster and provide highly relevant information. The thermodynamic analysis of the folding process and of kinetic intermediates opens whole new avenue of understanding. Breaking the protein folding code would enable scientists to look at a gene whose function is unknown and predict the three-dimensional structure of the protein it encodes. This would give them a very good idea of what the gene does. In this review we hope to bring together the information available about protein folding with particular emphasis on folding intermediate(s). Additionally, the practical consequences of the solution of the protein folding problem in medicine and biotechnology are also discussed.     

Specific non-native hydrophobic interactions in a hidden folding intermediate: implications for protein folding

Structures of intermediates and transition states in protein folding are usually characterized by amide hydrogen exchange and protein engineering methods and interpreted on the basis of the assumption that they have native-like conformations. We were able to stabilize and determine the high-resolution structure of a partially unfolded intermediate that exists after the rate-limiting step of a four-helix bundle protein, Rd-apocyt b(562), by multidimensional NMR methods. The intermediate has partial native-like secondary structure and backbone topology, consistent with our earlier native state hydrogen exchange results. However, non-native hydrophobic interactions exist throughout the structure. These and other results in the literature suggest that non-native hydrophobic interactions may occur generally in partially folded states. This can alter the interpretation of mutational protein engineering results in terms of native-like side chain interactions. In addition, since the intermediate exists after the rate-limiting step and Rd-apocyt b(562) folds very rapidly (k(f) approximately 10(4) s(-1)), these results suggest that non-native hydrophobic interactions, in the absence of topological misfolding, are repaired too rapidly to slow folding and cause the accumulation of folding intermediates. More generally, these results illustrate an approach for determining the high-resolution structure of folding intermediates.     

Structural mapping of the ClpB ATPases of Plasmodium falciparum: Targeting protein folding and secretion for antimalarial drug design

Caseinolytic chaperones and proteases (Clp) belong to the AAA+ protein superfamily and are part of the protein quality control machinery in cells. The eukaryotic parasite Plasmodium falciparum, the causative agent of malaria, has evolved an elaborate network of Clp proteins including two distinct ClpB ATPases. ClpB1 and ClpB2 are involved in different aspects of parasitic proteostasis. ClpB1 is present in the apicoplast, a parasite-specific and plastid-like organelle hosting various metabolic pathways necessary for parasite growth. ClpB2 localizes to the parasitophorous vacuole membrane where it drives protein export as core subunit of a parasite-derived protein secretion complex, the Plasmodium Translocon of Exported proteins (PTEX); this process is central to parasite virulence and survival in the human host. The functional associations of these two chaperones with parasite-specific metabolism and protein secretion make them prime drug targets. ClpB proteins function as unfoldases and disaggregases and share a common architecture consisting of four domains—a variable N-terminal domain that binds different protein substrates, followed by two highly conserved catalytic ATPase domains, and a C-terminal domain. Here, we report and compare the first crystal structures of the N terminal domains of ClpB1 and ClpB2 from Plasmodium and analyze their molecular surfaces. Solution scattering analysis of the N domain of ClpB2 shows that the average solution conformation is similar to the crystalline structure. These structures represent the first step towards the characterization of these two malarial chaperones and the reconstitution of the entire PTEX to aid structure-based design of novel anti-malarial drugs.     

The extracellular regions of PSMA and the transferrin receptor contain an aminopeptidase domain: implications for drug design

The Aeromonas proteolytica aminopeptidase (AMP), Pseudomonas sp. (RS-16) carboxypeptidase G2 (CPG2), and Streptomyces griseus aminopeptidase (SGAP) are zinc dependent proteolytic enzymes with cocatalytic zinc ion centers and a conserved aminopeptidase fold. A BLAST search with the sequence of the solved AMP structure indicated that a similar domain could be found in prostate-specific membrane antigen (PSMA) and the transferrin receptor (TfR). When the PSMA or TfR sequence was input into the THREADER program, the top structural matches were SGAP and AMP confirming that these are structurally conserved domains. Optimal sequence alignment of PSMA and TfR using the known three-dimensional structures of AMP, CPG2, and SGAP shows that the critical amino acids involved in forming the catalytic pocket are conserved in PSMA but absent in the TfR. The specificity pocket in AMP is formed from four aromatic side chains and the equivalent region in CPG2/PSMA has a changed sequence pattern. Since CPG2 and PSMA are folate hydrolases, the changed specificity pocket leaves space to accommodate the large pteroate moiety of folic acid. In contrast, no enzyme function has been ascribed to the TfR.     

The conformation of a nascent polypeptide inside the ribosome tunnel affects protein targeting and protein folding

In this report, we describe insights into the function of the ribosome tunnel that were obtained through an analysis of an unusual 25 residue N‐terminal motif (EspP_1‐25) associated with the signal peptide of the Escherichia coli EspP protein. It was previously shown that EspP_1‐25 inhibits signal peptide recognition by the signal recognition particle, and we now show that fusion of EspP_1‐25 to a cytoplasmic protein causes it to aggregate. We obtained two lines of evidence that both of these effects are attributable to the conformation of EspP_1‐25 inside the ribosome tunnel. First, we found that mutations in EspP_1‐25 that abolished its effects on protein targeting and protein folding altered the cross‐linking of short nascent chains to ribosomal components. Second, we found that a mutation in L22 that distorts the tunnel mimicked the effects of the EspP_1‐25 mutations on protein biogenesis. Our results provide evidence that the conformation of a polypeptide inside the ribosome tunnel can influence protein folding under physiological conditions and suggest that ribosomal mutations might increase the solubility of at least some aggregation‐prone proteins produced in E. coli.     

Alpha-synuclein multistate folding thermodynamics: implications for protein misfolding and aggregation

Alpha-synuclein aggregation has been tightly linked with the pathogenesis of Parkinson's disease and other neurodegenerative disorders. Despite the protein's putative function in presynaptic vesicle regulation, the roles of lipid binding in modulating alpha-synuclein conformations and the aggregation process remain to be fully understood. This study focuses on a detailed thermodynamic characterization of monomeric alpha-synuclein folding in the presence of SDS, a well-studied lipid mimetic. Far-UV CD spectroscopy was employed for detection of conformational transitions induced by SDS, temperature, and pH. The data we present here clearly demonstrate the multistate nature of alpha-synuclein folding, which involves two predominantly alpha-helical partially folded thermodynamic intermediates that we designate as F (most folded) and I (intermediately folded) states. Likely structures of these alpha-synuclein conformational states are also discussed. These partially folded forms can exist in the presence of either monomeric or micellar forms of SDS, which suggests that alpha-synuclein has an intrinsic propensity for adopting multiple alpha-helical structures even in the absence of micelle or membrane binding, a feature that may have implications for its biological activity and toxicity. Additionally, we discuss the relation between alpha-synuclein three-state folding and its aggregation, within the context of isothermal titration calorimetry and transmission electron microscopy measurements of SDS-initiated oligomer formation.     

Distributions in protein conformation space: implications for structure prediction and entropy

By considering how polymer structures are distributed in conformation space, we show that it is possible to quantify the difficulty of structural prediction and to provide a measure of progress for prediction calculations. The critical issue is the probability that a conformation is found within a specified distance of another conformer. We address this question by constructing a cumulative distribution function (CDF) for the average probability from observations about its limiting behavior at small displacements and numerical simulations of polyalanine chains. We can use the CDF to estimate the likelihood that a structure prediction is better than random chance. For example, the chance of randomly predicting the native backbone structure of a 150-amino-acid protein to low resolution, say within 6 A, is 10(-14). A high-resolution structural prediction, say to 2 A, is immensely more difficult (10(-57)). With additional assumptions, the CDF yields the conformational entropy of protein folding from native-state coordinate variance. Or, using values of the conformational entropy change on folding, we can estimate the native state's conformational span. For example, for a 150-mer protein, equilibrium alpha-carbon displacements in the native ensemble would be 0.3-0.5 A based on T Delta S of 1.42 kcal/(mol residue).     

Membrane-bound states of alpha-lactalbumin: implications for the protein stability and conformation.

alpha-Lactalbumin (alpha-LA) associates with dimyristoylphosphatidylcholine (DMPC) or egg lecithin (EPC) liposomes. Thermal denaturation of isolated DMPC or EPC alpha-LA complexes was dependent on the metal bound state of the protein. The intrinsic fluorescence of thermally denatured DMPC-alpha-LA was sensitive to two thermal transitions: the Tc of the lipid vesicles, and the denaturation of the protein. Quenching experiments suggested that tryptophan accessibility increased upon protein-DMPC association, in contrast with earlier suggestions that the limited emission red shift upon association with the liposome was due to partial insertion of tryptophan into the apolar phase of the bilayer (Hanssens I et al., 1985, Biochim Biophys Acta 817:154-166). On the other hand, above the protein transition (70 degrees C), the spectral blue shifts and reduced accessibility to quencher suggested that tryptophan interacts significantly with the apolar phase of either DMPC and EPC. At pH 2, where the protein inserts into the bilayer rapidly, the isolated DMPC-alpha-LA complex showed a distinct fluorescence thermal transition between 40 and 60 degrees C, consistent with a partially inserted form that possesses some degree of tertiary structure and unfolds cooperatively. This result is significant in light of earlier findings of increased helicity for the acid form, i.e., molten globule state of the protein (Hanssens I et al., 1985, Biochim Biophys Acta 817:154-166). These results suggest a model where a limited expansion of conformation occurs upon association with the membrane at neutral pH and physiological temperatures, with a concomitant increase in the exposure of tryptophan to external quenchers; i.e., the current data do not support a model where an apolar, tryptophan-containing surface is covered by the lipid phase of the bilayer.     

Extracting information on folding from the amino acid sequence: consensus regions with preferred conformation in homologous proteins.

It is investigated whether protein segments predicted to have a well-defined conformational preference in the absence of tertiary interactions are conserved in families of homologous proteins. The prediction method follows the procedures of Rooman, M., Kocher, J.-P., and Wodak, S. (preceding paper in this issue). It uses a knowledge-based force field that incorporates only local interactions along the sequence and identifies segments whose lowest energy structure displays a sizable energy gap relative to other computed conformations. In 13 of the protein families and subfamilies considered that are sufficiently homologous to have similar 3D structures, at least one region is consistently predicted as having the same preferred conformation in virtually all family members. These regions are between 4 and 26 residues long. They are often located at chain ends and correspond primarily to segments of secondary structure heavily involved in interactions with the rest of the protein, suggesting that they could act as nuclei around which other parts of the structure would assemble. Experimental data on early folding intermediates or on protein fragments with appreciable structure in aqueous solution are available for more than half of the protein families. Comparison of our results with these data is quite favorable. They reveal that each of the experimentally identified early formed, or independently stable, substructures harbors at least one of the segments consistently predicted as having a preferred conformation by our procedure. The implications of our findings for the conservation of folding pathways in homologous proteins are discussed.     

The dendritic branch is the preferred integrative unit for protein synthesis-dependent LTP

The late-phase of long-term potentiation (L-LTP), the cellular correlate of long-term memory, induced at some synapses facilitates L-LTP expression at other synapses receiving stimulation too weak to induce L-LTP by itself. Using glutamate uncaging and two-photon imaging, we demonstrate that the efficacy of this facilitation decreases with increasing time between stimulations, increasing distance between stimulated spines and with the spines being on different dendritic branches. Paradoxically, stimulated spines compete for L-LTP expression if stimulated too closely together in time. Furthermore, the facilitation is temporally bidirectional but asymmetric. Additionally, L-LTP formation is itself biased toward occurring on spines within a branch. These data support the Clustered Plasticity Hypothesis, which states that such spatial and temporal limits lead to stable engram formation, preferentially at synapses clustered within dendritic branches rather than dispersed throughout the dendritic arbor. Thus, dendritic branches rather than individual synapses are the primary functional units for long-term memory storage.     

Studies of peptide conformation: Backbone folding of tetrapeptide derivatives in methanol and water

Derivatives of tetrapeptide sequences considered likely to form β-turns were investigated by the study of their proton magnetic resonances in methanol and in water. Differential broadening of N—H resonances by an added nitroxyl was used to indicate the presence of the sequestered N—H proton expected in β-turn conformations. Transfer of magnetic saturation from solvent water protons to N—H protons was also examined. The evidence is consistent with significant contributions by β-turn-like backbones to the conformational averages in methanol of the sequences Gly-L -Pro-D -Val-Gly, D (or L )-Val-L -Pro-Gly-Gly, and Gly-L -Pro-L -Asn-Gly, but not the sequence Gly-D -Ala-L -Val-Gly. It is suggested that a Type I turn, Likely in Gly-L -Pro-L -Asn-Gly derivatives, is characterized by sequestered N—H protons of both the third and fourth residues. For all of the peptide derivatives, save possibly Ac-L -Val-L -Pro-Gly-Gly-NHNH₂, contributions from folded structures in water are not detectable by line-broadening experiments. However, the transfer of saturation experiments may be interpreted as indicating some degree of chain folding in water.