Criteria that discriminate between native proteins and incorrectly folded models

Various theoretical concepts, such as free energy potentials, electrostatic interaction potentials, atomic packing, solvent-exposed surface, and surface charge distribution, were tested for their ability to discriminate between native proteins and misfolded protein models. Misfolded models were constructed by introducing incorrect side chains onto polypeptide backbones: side chains of the alpha-helical hemerythrin were modeled on the beta-sheeted backbone of immunoglobulin VL domain, whereas those of the VL domain were similarly modeled on the hemerythrin backbone. CONGEN, a conformational space sampling program, was used to construct the side chains, in contrast to the previous work, where incorrect side chains were modeled in all trans conformations. Capability of the conformational search procedure to reproduce native conformations was gauged first by rebuilding (the correct) side chains in hemerythrin and the VL domain: constructs with r.m.s. differences from the x-ray side chains 2.2-2.4 A were produced, and many calculated conformations matched the native ones quite well. Incorrectly folded models were then constructed by the same conformational protocol applied to incorrect amino acid sequences. All CONGEN constructs, both correctly and incorrectly folded, were characterized by exceptionally small molecular surfaces and low potential energies. Surface charge density, atomic packing, and Coulomb formula-based electrostatic interactions of the misfolded structures and the correctly folded proteins were similar, and therefore of little interest for diagnosing incorrect folds. The following criteria clearly favored the native structures over the misfolded ones: 1) solvent-exposed side-chain nonpolar surface, 2) number of buried ionizable groups, and 3) empirical free energy functions that incorporate solvent effects.     

Free energy landscapes of two model peptides: α-helical and β-hairpin peptides explored with Brownian dynamics simulation

We applied an atomistic Brownian dynamics (BD) simulation with multiple time step method for the folding simulation of a 13-mer α-helical peptide and a 12-mer β-hairpin peptide, giving successful folding simulations. In this model, the driving energy contribution towards folding came from both electrostatic and van der Waals interactions for the α-helical peptide and from van der Waals interactions for the β-hairpin peptide. Although, many non-native structures having the same or lower energy than that of native structure were observed, the folded states formed the most populated cluster when the structures obtained by the BD simulations were subjected to the cluster analysis based on distance-based root mean square deviation of side-chains between different structures. This result indicates that we can predict the native structures from conformations sampled by BD simulation.     

Solution structure of the tissue-type plasminogen activator kringle 2 domain complexed to 6-aminohexanoic acid an antifibrinolytic drug.

The solution structure of a recombinant tissue-type plasminogen activator kringle 2 domain, complexed with the antifibrinolytic drug 6-aminohexanoic acid (6-AHA) was determined via 1H nuclear magnetic resonance spectroscopy and dynamical simulated annealing calculations. The structure determination is based on 610 intramolecular kringle 2 and 14 intermolecular kringle 2-6-AHA interproton distance restraints, as well as on 82 torsion angle restraints. Three sets of simulated annealing structures were computed from three different classes of starting structures: (1) random conformations devoid of disulfide bridges; (2) random conformations that contain correct disulfide bonds; and (3) a folded conformation modeled after the homologous prothrombin kringle 1 X-ray crystallographic structure. All three sets of structures are well defined, with averaged atomic root-mean-square deviations between individual structures and mean set structures of 0.77, 0.99 and 0.70 A for backbone atoms, and 1.36, 1.55 and 1.41 A for all atoms, respectively. Kringle 2 is an oblate ellipsoid with overall dimensions of approximately 34 A x 30 A x 17 A. It exhibits a compact globular conformation characterized by a number of turns and loop elements as well as by one right-handed alpha-helix and five (1 extended and 4 rudimentary) antiparallel beta-sheets. The extended beta-sheet exhibits a right-handed twist. Close van der Waals' contacts between the Cys22-Cys63 and Cys51-Cys75 disulfide bridges and the central hydrophobic core composed of the Trp25, Leu46, His48a and Trp62 side-chains are among the distinguishing features of the kringle 2 fold. The binding site for 6-AHA appears as a rather exposed cleft with a negatively charged locus defined by the Asp55 and Asp57 side-chains, and with an aromatic pocket structured by the Tyr36, Trp62, His64 and Trp72 side-chains. The Trp62 and His64 rings line the back surface of the pocket, while the Tyr36 and Trp72 rings confine it from two sides. The Trp62 and Trp72 indole rings conform a V-shaped groove. The methyl groups of Val35 also contribute lipophilic character to the ligand-interacting surface. It is suggested that the positively charged side-chains of Lys34 and, potentially, Arg69 may favor interactions with the carboxylate group of the ligand. The Trp25 and Tyr74 aromatic rings, although conserved elements of the binding site structure, seem not to undergo direct contacts with the ligand.     

Simulations as analytical tools to understand protein aggregation and predict amyloid conformation

Computational tools are increasingly being applied to solve the protein aggregation problem, providing insight into amyloid structures and aggregation mechanisms. The paradigm of Abeta amyloid structure elucidation provides an example of an innovative experimental design and endeavor, echoing the computational testing of possible molecular associations, all reflected in the current Ma-Nussinov-Tycko model of the Abeta amyloid. Simulations have shown that dimer formation can lock some misfolded conformations, and catalyze the shift of the equilibrium away from the native state. In most cases, a stable amyloid seed requires at least two-layered beta-sheets with properly registered side-chains. Under kinetic control, the final protein aggregations are the outcome of maximizing the van der Waals interactions between side chains and backbone hydrogen bonds.     

Secondary structure length as a determinant of folding rate of proteins with two- and three-state kinetics

We present a simple method for determining the folding rates of two- and three-state proteins from the number of residues in their secondary structures (secondary structure length). The method is based on the hypothesis that two- and three-state foldings share a common pattern. Three-state proteins first condense into metastable intermediates, subsequent forming of alpha-helices, turns, and beta-sheets at slow rate-limiting step. The folding rate of such proteins anticorrelate with the length of these beta-secondary structures. It is also assumed that in two-state folding, rapidly folded alpha-helices and turns may facilitate formation of fleeting unobservable intermediates and thus show two-state behavior. There is an inverse relationship between the folding rate and the length of beta-sheets and loops. Our study achieves 94.0 and 88.1% correlations with folding rates determined experimentally for 21 three- and 38 two-state proteins, respectively, suggesting that protein-folding rates are determined by the secondary structure length. The kinetic kinds are selected on the basis of a competitive formation of hydrophobic collapse and alpha-structure in early intermediates.     

Conformation of peptide fragments of proteins in aqueous solution: implications for initiation of protein folding

Applications of sensitive new technologies, in particular, two-dimensional NMR spectroscopy, have allowed detection of folded structures in short peptide fragments of proteins in aqueous solution under conditions where native proteins fold. These structures are in rapid dynamic exchange with unfolded states. These observations provide evidence in support of models for protein folding which postulate localized regions of folded structure as initiation sites for the folding process. Since these initiation processes are expected to be rapid, such models are consistent with kinetic evidence that the rate-determining steps of protein folding occur late in the process and probably involve rearrangement of incorrectly folded intermediates.     

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.     

Construction of side-chains in homology modelling. Application to the C-terminal lobe of rhizopuspepsin

A detailed and rule-based side-chain modelling procedure for globular proteins is presented. It uses the conformational information contained in a homologous (template) structure as a starting point and includes recipes for atom placement and for checking and improving the atomic positions. The scheme does not rely on intuitive judgements or visual examination of the model during construction or refinement. It comprises four stages; the first three are relatively simple and the fourth is more complex. In the first stage, initial conformations for as many atoms as possible are transferred from the template structure based on the application of trends reported previously. Second, these trends are used to correct poor van der Waals overlaps. Third, the remaining side-chains atoms (those for which no information is contained in the template) are placed by evaluating their rigid rotation, van der Waals surfaces. The fourth stage consists of a hierarchial series of conformational checks. They involve the evaluation of individual residue energies in the absence and presence of the rest of the protein relative to statistical trends observed in the template structure, the comparison of hydrogen-bonding patterns and side-chain accessibilities in the model and template and brief energy minimization followed by an evaluation of the rigid rotation potential energy surfaces of each side-chain. The checks pinpoint "incorrectly" modelled side-chains, suggest conformational changes and provide a means for determining the portions of the model that are likely to be correct and those likely to be in error. The procedure developed in the paper is tested by modelling the side-chains of the C-terminal lobe of the aspartyl proteinase rhizopuspepsin, using the rhizopuspepsin backbone and the homologous protein, penicillopepsin, as a template for the side-chains. The resultant model was compared to the high-resolution X-ray structure of rhizopuspepsin. Using penicillopepsin data only (stage I), 58% of the chi 1 dihedrals and 44% of the chi 2 dihedrals were modelled correctly. Once poor van der Waals overlaps had been corrected and all of the atoms had been placed (stages II and III), 86% of the chi 1 dihedrals and 75% of the chi 2 dihedrals were correct. After the refinement had been completed (stage IV), 92% of the chi 1 dihedrals and 81% of the chi 2 dihedrals were correctly positioned.(ABSTRACT TRUNCATED AT 400 WORDS)     

Progressive stabilization of intermediate and transition states in protein folding reactions by introducing surface hydrophobic residues

It can be argued from the principle of solvent exclusion that the introduction of hydrophobic residues onto the surface of a protein will not destabilize the folded state because the nonpolar side chain will be at least as exposed in the unfolded state as it is when the protein chain is folded. A comparison of the folding pathway of wild type and 11 site-directed mutants of CD2.d1 shows this to be true. In fact, owing to partial burial of nonpolar groups as folding proceeds, we find that the rapidly formed intermediate state and, to a greater extent, the transition state are generally stabilized by hydrophobic surface mutations. This effect is slightly moderated in the folded state presumably by the perturbation of van der Waals' contacts and/or local electrostatic interactions that have a greater influence in this fully compact structure. The fact that in all but one case we find that stabilization of the rapidly collapsed intermediate is accompanied by a faster acquisition of the folded state refutes the argument that I states are generally "off pathway" conformations or ensembles that lead to the inhibition of otherwise more rapid folding trajectories.     

Solvent accessibility, protein surfaces, and protein folding.

Studies of the native structures of proteins, together with measurements of the thermodynamic properties of the transition between unfolded and native states, have defined the major components of the forces that stabilize native protein structures. However, the nature of the intermediates in the folding process remains largely hypothetical. It is a fairly widespread and not implausible assumption that the intermediates in the folding of a monomeric protein contain the same kinds of secondary and tertiary structures that appear in the native conformation, and that, although unstable, their lifetimes are prolonged by forces similar to those that stabilize the native structure. We wished to examine what happens if, during the folding of a monomeric protein, regions of secondary structure come together to form an intermediate of reduced instability. We applied calculations of accessible surface area (a measure of hydrophobic stabilization) and parameterized nonbonded energy calculations (measuring the strengths of van der Waals forces) to identify the kinds of stabilizing interactions that might be available to such an intermediate. First, we analyzed the total buried surface area of two types of proteins into contributions from formation of secondary structure alone, interaction of pairs of secondary-structural elements, the formation of the structure alone, interaction of pairs of secondary-structural elements, the formation of the complete secondary structure without the turns, and the complete native structure. The formation of secondary structure alone, without tertiary-structural interactions, buries roughly half the surface that the complete structure does. We then analyzed in more detail the approach of two alpha-helices to form a complex, as an illustrative example of the nature of the interaction between compact structural units which remain fairly rigid during their interaction. Many features of the results are not limited to the interaction of alpha-helices. (The results therefore neither confirm nor refute the hypothesis that alpha-helices are intermediates in the folding proteins). We find that the first forces to be felt upon approach arise from solvent conditions on the relative position and orientation of the two helices as does the close packing which optimizes the van der Waals interactions at shorter distances apart. Therefore there appears to be a range of distances in which hydrophobic interactions could create a nonspecific complex between two helices in which the side chains might have sufficient time to seek the proper interdigitation observed in the native structure, where the two helices are in intimate contact. Indeed, we find that only in the final stages of approach is the native geometry the most stable; in the region in which solvent-exclusion forces predominate, the conformation with helix axes parallel is more stable than the native conformation, in the cases we examined...     

Prediction of protein side-chain conformation by packing optimization

We have developed a rapid and completely automatic method for prediction of protein side-chain conformation, applying the simulated annealing algorithm to optimization of side-chain packing (van der Waals) interactions. The method directly attacks the combinatorial problem of simultaneously predicting many residues' conformation, solving in 8 to 12 hours problems for which the systematic search would require over 10(300) central processing unit years. Over a test set of nine proteins ranging in size from 46 to 323 residues, the program's predictions for side-chain atoms had a root-mean-square (r.m.s.) deviation of 1.77 A overall versus the native structures. More importantly, the predictions for core residues were especially accurate, with an r.m.s. value of 1.25 A overall: 80 to 90% of the large hydrophobic side-chains dominating the internal core were correctly predicted, versus 30 to 40% for most current methods. The predictions' main errors were in surface residues poorly constrained by packing and small residues with greater steric freedom and hydrogen bonding interactions, which were not included in the program's potential function. van der Waals interactions appear to be the supreme determinant of the arrangement of side-chains in the core, enforcing a unique allowed packing that in every case so far examined matches the native structure.     

Interdomain salt bridges modulate ligand-induced domain motion of the sulfate receptor protein for active transport.

The refined crystal structure of the liganded form of the Salmonella typhimurium sulfate-binding protein, a periplasmic receptor of active transport, is made up of two globular domains bisected by a deep cleft wherein the dehydrated sulfate is completely engulfed and bound by hydrogen bonds and van der Waals' forces. Two salt bridges (between Glu15 and Arg174 and between Asp68 and Arg134) span the cleft opening. To elucidate the role of the inter-domain salt bridges in the ligand-induced domain motion, the acidic residues were changed (singly and together) to their corresponding amide side-chains by site-directed mutagenesis of the recombinant Escherichia coli sulfate-binding protein. Rapid kinetics and equilibrium measurements of sulfate binding to the purified mutant proteins demonstrate that these salt bridges stabilize the closed liganded form of the receptor and modulate the rate of cleft opening. Our results have new implications in understanding the dynamics of many other multidomain proteins that undergo similar large-scale domain motions.     

The posttranslational concept of protein folding: how valid is it?

Two concepts of protein folding are known. One of them, the cotranslational concept, states that a protein folds during the synthesis of the polypeptide chain on the ribosome. According to the other, the posttranslational concept, the protein starts to fold just after the synthesis of its polypeptide chain. This article attempts to show that the posttranslational concept is hardly suited to solve the problem of protein folding. In our opinion, polypeptide chains cannot be represented as random coils. They are stiff chain-like macromolecules rather than flexible ones: the single bond rotational barriers of a polypeptide substantially exceed the accepted standard values; even in strong denaturing conditions, a protein possesses a considerable amount of residual folded structures. We believe that the popular "hierarchical" models for the protein folding mechanism are not realistic because the formation of secondary and tertiary structures of proteins occurs simultaneously and cooperatively. The time for the elongation of a polypeptide chain by one amino acid residue during biosynthesis exceeds considerably the time of the formation of alpha-helices and beta-sheets in proteins as well as the time supposed for the spatial structure formation of a native protein during renaturation. Thus, we believe that the mechanism of protein folding in vivo cannot be clarified by denaturation-renaturation experiments. In our opinion, the phenomenon of protein renaturation is no more than the restoration of native protein conformation (which initially forms cotranslationally) disrupted during denaturation, and thus denaturation-renaturation experiments cannot serve as a model to clarify the mechanism of protein folding.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

A new protein folding algorithm based on hydrophobic compactness: Rigid Unconnected Secondary Structure Iterative Assembly (RUSSIA). II: Applications

The RUSSIA procedure (Rigid Unconnected Secondary Structure Iterative Assembly) produces structural models of cores of small- and medium-sized proteins. Loops are omitted from this treatment and regular secondary structures are reduced to points, the centers of their hydrophobic faces. This methodology relies on the maximum compactness of the hydrophobic residues, as described in detail in Part I. Starting data are the sequence and the predicted limits and natures of regular secondary structures (alpha or beta). Helices are treated as rigid cylinders, whereas beta-strands are collectively taken into account within beta-sheets modeled by helicoid surfaces. Strands are allowed to shift along their mean axis to allow some flexibility and the alpha-helices can be placed on either side of beta-sheets. Numerous initial conformations are produced by discrete rotations of the helices and sheets around the direction going from the center of their hydrophobic face to the global center of the protein. Selection of proposed models is based upon a criterion lying on the minimization of distances separating hydrophobic residues belonging to different regular secondary structures. The procedure is rapid and appears to be robust relative to the quality of starting data (nature and length of regular secondary structures). This dependence of the quality of the model on secondary structure prediction and in particular the beta-sheet topology, is one of the limits of the present algorithm. We present here some results for a set of 12 proteins (alpha, beta and alpha/beta classes) of lengths 40-166 amino acids. The r.m.s. deviations for core models with respect to the native proteins are in the range 1.4-3.7 A.     

2.2 A resolution structure analysis of two refined N-acetylneuraminyl-lactose--wheat germ agglutinin isolectin complexes

The crystal structures of complexes of isolectins 1 and 2 of wheat germ agglutinin (WGA1 and WGA2) with N-acetylneuraminyl-lactose (NeuNAc-alpha(2-3)-Gal-beta(1-4)-Glc) have been refined on the basis of data in the 8 to 2.2 A resolution range to final crystallographic R-factors of 17.2% and 15.3% (Fo greater than 1 sigma), respectively. Specific binding interactions and water association, as well as changes in conformation and mobility of the structure upon ligand binding, were compared in the two complexes. The temperature factors (B = 16.3 A2 and 18.4 A2) were found to be much lower compared with those of their respective native structures (19 to 22 A2). Residues involved in sugar binding, dimerization and in lattice contacts exhibit the largest decreases in B-value, suggesting that sugar binding reduces the overall mobility of the protein molecules in the crystal lattice. The binding mode of this sialyl-trisaccharide, an important cell receptor analogue, has been compared in the two isolectins. Only one of the two unique binding sites (4 per dimer), located in the subunit/subunit interface, is occupied in the crystals. This site, termed the "primary" binding site, contains one of the five amino acid substitutions that differentiate WGA1 and WGA2. Superposition of the refined models in each of the independent crystallographic environments indicates a close match only of the terminal non-reducing NeuNAc residue (root-mean-square delta r of 0.5 to 0.6 A). The Gal-Glc portion was found to superimpose poorly, lack electron density, and possess high atomic thermal factors. In both complexes NeuNAc is stabilized through contact with six amino acid side-chains (Ser114 and Glu115 of subunit 1 and Ser62, Tyr64, Tyr(His)66 and Tyr73 of subunit 2), involving all NeuNAc ring substituents. Refinement has allowed accurate assessment of the contact distances for four hydrogen bonds, a strong buried non-polar contact with the acetamido CH3 group and a large number of van der Waals' interactions with the three aromatic side-chains. The higher affinity of N-acetylneuraminyl-lactose observed by nuclear magnetic resonance studies for WGA1 can be explained by the more favorable binding interactions that occur when residue 66 is a Tyr. The tyrosyl side-chain provides a larger surface for van der Waals' stacking against the NeuNAc pyranose ring than His66 and a hydrogen bond contact with Gal (C2-OH), not possible in WGA2.(ABSTRACT TRUNCATED AT 400 WORDS)     

Distributions of water around amino acid residues in proteins

The atomic co-ordinates from 16 high-resolution (less than or equal to 1.7 A = 0.1 nm), non-homologous proteins have been used to study the distributions of water molecule sites around the 20 different amino acid residues. The proportion of residues whose main-chain atoms are in contact with water molecules was fairly constant (between 40% and 60%), irrespective of the nature of the side-chain. However, the proportion of residues whose side-chain atoms were in contact with water molecules showed a clear (inverse) correlation with the hydrophobicity of the residue, being as low as 14% for leucine and isoleucine but greater than 80% for asparagine and arginine. Despite the problems in determining accurate water molecule sites from X-ray diffraction data and the complexity of the protein surface, distinct non-random distributions of water molecules were found. These hydration patterns are consistent with the expected stereochemistry of the potential hydrogen-bonding sites on the polar side-chains. The water molecules around apolar side-chains lie predominantly at van der Waals' contact distances, but most of these have a primary, shorter contact with a neighbouring polar atom. Further analysis of these distributions, combined with energy minimization techniques, should lead to improved modelling of protein structures, including their primary shells of hydration.     

Aromatic-aromatic interactions in the zinc finger motif. Analysis of the two-dimensional nuclear magnetic resonance structure of a mutant domain.

The folding and stability of globular proteins are determined by a variety of chemical mechanisms, including hydrogen bonds, salt bridges and the hydrophobic effect. Of particular interest are weakly polar interactions involving aromatic rings, which are proposed to regulate the geometry of closely packed protein interiors. Such interactions reflect the electrostatic contribution of pi-electrons and, unlike van der Waals' interactions and the hydrophobic effect, may, in principle, introduce a directional force in a protein's hydrophobic core. Although the weakly polar hypothesis is supported by a statistical analysis of protein structures, the general importance of such contributions to protein folding and stability is unclear. Here, we show the presence of alternative aromatic-aromatic interactions in the two-dimensional nuclear magnetic resonance structure of a mutant Zn finger. Changes in aromatic packing lead in turn to local and non-local differences between the structures of a wild-type and mutant domain. The results provide insight into the evolution of Zn finger sequences and have implications for understanding how geometric relationships may be chemically encoded in a simple sequence template.     

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.