Analysis of the code relating sequence to conformation in globular proteins. Development of a stereochemical alphabet on the basis of intra-residue information

1. The relation of primary sequence to all residue backbone conformations was explored to test out starting conformations for protein folding. 2. Information theory was used to obtain measures of information which quantitate the role of each residue in determining its own conformation; i.e. intra-residue information. 3. The information measures are plotted as a function of varphi, psi peptide-backbone angles and varphi, psi contour maps obtained for each of the 20 amino acids. These show characteristic differences between residues. 4. To find practical ways of relating sequence to varphi, psi angles, several types of stereochemical alphabet were investigated. The value of these was tested by using them to predict the varphi, psi angles of nine different proteins. 5. A difference plot was constructed to show regions of the sequence that require little or no information extra to the intra-residue information in order to predict a correct conformation. These regions are suggested to be candidates for nucleating sites in the protein.     

Analysis of the code relating sequence to conformation in globular proteins. The distribution of residue pairs in turns and kinks in the backbone chain

1. The residue pair is considered as the fundamental unit which differentiates alpha-helix, beta-pleated sheet and the various turns and kink structures of the protein backbone. 2. The HPLG alphabet (Robson & Pain, 1974) is used to group pairs of residues, giving 16 possible conformational pairs, all of which are found with differing frequencies in the nine proteins examined. 3. The frequencies of occurrence of the 16 different types of turn or kink are analysed in relation to the constituent amino acids. Those containing the L or G conformation are of low frequency and are grouped for purposes of this analysis. 4. The distribution of amino acids within all the conformational pairs is non-random, with distinct preferences shown by certain residues. 5. All pairs containing an L or G conformation require the presence of a glycine or a proton-donor side chain. 6. The results are discussed in terms of the determination of these ;random' structures by local interactions.     

Analysis of the code relating sequence to conformation in globular proteins. An informational analysis of the role of the residue in determining the conformation of its neighbours in the primary sequence

1. The effect exerted by a residue on the conformation of neighbouring residues was analysed by using data from nine globular proteins of known sequence and conformation. 2. An information measure was used which estimated the role of a residue in influencing neighbouring conformations and also its tendency to influence the lengths of runs of residues in that conformation. This measure was estimated for each residue in all conformations defined by domains on the varphi, psi diagram. 3. Plots of the information measure yielded an intercept, which was a measure of intra-residue information for a residue. The slope was a measure of the statistical co-operativity or tendency of the residue to influence the occurrence of its neighbours in runs of a particular conformation. Both parameters are a function of the residue type. Statistical co-operativity is found in the alpha(1)-helical (H(1)) and beta-pleated-sheet (P(2)) conformations and, to a lesser extent, in their distorted variants H(2) and P(1). 4. The directional nature of these influences for H(1) and P(2) conformations is illustrated by plots of the information measure against the distance m from the residue, for m=-10 to +10. 5. The results for statistical co-operativity are discussed in relation to theories of helix-coil and pleated-sheet-coil transitions. The value of the information-theory-derived parameters in obtaining s parameters for the Zimm & Bragg (1959) equations is illustrated. 6. Directional effects are discussed with particular relation to mechanisms of the termination of helices and the involvement of the alpha(II) conformation and also to discontinuities in pleated-sheet conformations.     

A code relating sequence to structure in nucleic acids

Nucleic acids are elucidated in configuration space. An algorithm relating sequence to stability in A and B helical secondary structures, is stated to incorporate NMR conformational and optical melting data. This made possible a classification of elementary sequences in terms of configuration forces driving between A and B states, a finding useful in prediction of structural behavior of different sequences of DNA, RNA and their hybrids.     

Analysis of the relationship between side-chain conformation and secondary structure in globular proteins

The relationship between the preferred side-chain dihedral angles and the secondary structure of a residue was examined. The structures of 61 proteins solved to a resolution of 2.0 A (1 A = 0.1 nm) or better were analysed using a relational database to store the information. The strongest feature observed was that the chi 1 distribution for most side-chains in an alpha-helix showed an absence of the g- conformation and a shift towards the t conformation when compared to the non-alpha/beta structures. The exceptions to this tendency were for short polar side-chains that form hydrogen bonds with the main-chain which prefer g+. Shifts in the chi 1 preferences for residues in the beta-sheet were observed. Other side-chain dihedral angles (chi 2, chi 3, chi 4) were found to be influenced by the main-chain. This paper presents more accurate distributions for the side-chain dihedral angles which were obtained from the increased number of proteins determined to high resolution. The means and standard deviations for chi 1 and chi 2 angles are presented for all residues according to the secondary structure of the main-chain. The means and standard deviations are given for the most popular conformations for side-chains in which chi 3 and chi 4 rotations affect the position of C atoms.     

Left-handed helix conformation of poly-L-proline II type in globular proteins. Statistics of incidence and a role of sequence]

Regions of left-handed polyproline II type conformation in globular proteins were studied throughout the PDB bank. The length and sequence of corresponding fragments were analyzed. It was found that a lot of tetrapeptides (from combinatorial possible ones) show the tendency to be included in the left-handed helices. Much more tetrapeptides do not occur in this structure type.     

Theory of reversible denaturation of globular proteins.

A theoretical method is developed by which the character of the process of protein denaturation (e.g., whether or not it is of the all-or-none type) can be discussed in terms of conformation of native proteins and the forces stabilizing it. An important role is played by a quantity S(H): entropy of a protein molecule in solution in the conformational states with a given value of enthalpy H. It is demonstrated that the all-or-none type denaturation of proteins is a rather direct consequence of the globularity and specificity of the native conformations. Denaturations with significant intermediate states are discussed. Denaturations induced by added denaturants are also discussed.     

Structural and sequence characteristics of long alpha helices in globular proteins.

Elucidation of the detailed structural features and sequence requirements for alpha helices of various lengths could be very important in understanding secondary structure formation in proteins and, hence, in the protein folding mechanism. An algorithm to characterize the geometry of an alpha helix from its C(alpha) coordinates has been developed and used to analyze the structures of long alpha helices (number of residues > or = 25) found in globular proteins, the crystal structure coordinates of which are available from the Brookhaven Protein Data Bank. All long alpha helices can be unambiguously characterized as belonging to one of three classes: linear, curved, or kinked, with a majority being curved. Analysis of the sequences of these helices reveals that the long alpha helices have unique sequence characteristics that distinguish them from the short alpha helices in globular proteins. The distribution and statistical propensities of individual amino acids to occur in long alpha helices are different from those found in short alpha helices, with amino acids having longer side chains and/or having a greater number of functional groups occurring more frequently in these helices. The sequences of the long alpha helices can be correlated with their gross structural features, i.e., whether they are curved, linear, or kinked, and in case of the curved helices, with their curvature.     

Correlation of sequence and tertiary structure in globular proteins

The x-ray crystal structures of 19 selected proteins are examined empirically for correlations between the amino acid sequence and long-range, tertiary conformation. There is clear evidence for preferential associations of certain types of amino acids, particularly among the hydrophobic aliphatic, aromatic, and cysteine residues. However, the likelihoods of forming these residue-pair contacts are all less than 12%, so packing and geometric requirements must often take precedent over energetic considerations. The prediction of long-range contacts is not substantially improved by taking into account the sequentially previous residues. The analysis of atom–atom contacts shows a similar lack of predictive ability, but the results show that a good approximation to the interresidue energy function must include different types of interactions at two or three different sites on some amino acids. Backbone–backbone long-range interactions are relatively rare and nonspecific, whereas some “polar” side chains form hydrogen bonds from the polar groups while occasionally forming hydrophobic contacts with the remainder of the chain.     

Theory for the folding and stability of globular proteins

Using lattice statistical mechanics, we develop theory to account for the folding of a heteropolymer molecule such as a protein to the globular and soluble state. Folding is assumed to be driven by the association of solvophobic monomers to avoid solvent and opposed by the chain configurational entropy. Theory predicts a phase transition as a function of temperature or solvent character. Molecules that are too short or too long or that have too few solvophobic residues are predicted not to fold. Globular molecules should have a largely solvophobic core, but there is an entropic tendency for some residues to be "out of place", particularly in small molecules. For long chains, molecules comprised of globular domains are predicted to be thermodynamically more stable than spherical molecules. The number of accessible conformations in the globular state is calculated to be an exceedingly small fraction of the number available to the random coil. Previous estimates of this number, which have motivated kinetic theories of folding, err by many tens of orders of magnitude.     

Maps of preferential conformation of dipeptides in the structural regions of globular proteins

According to X-ray crystallographic analysis of 46 globular proteins the probability of dipeptides frequency in alpha-helical beta-sheet and random coil conformations has been studied on the basis of which the maps of preferentially conformational state of dipeptide in different elements of secondary structure of proteins have been obtained.     

Molecular cartography of globular proteins with application to antigenic sites

A method, molecular cartography, is introduced as a way to quantitate the topographic structure of a protein surface. The method is applied to the problem of antigenic determinants, and it is used to examine local and global topography of reported antigenic regions on the surface of myoglobin and lysozyme. In nine antigenic sites taken from the literature and studied in detail, no local property was found in sites that was not also found in remaining regions of the surface. However, a strong correlation was found between antigenic sites and regions of the surface that are globally exposed. This finding suggests that global exposure of the protein surface may play a primary role in determining the antigenic structure of the protein. Molecular cartography may be useful in other instances of protein–protein interactions such as those between proteolytic enzymes and their substrates.     

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.     

Structural changes and fluctuations of proteins. II. Analysis of the denaturation of globular proteins.

The statistical thermodynamic model of protein structure proposed in paper I is developed with special attention to the hydrophobic interaction. Calorimetric measurements of the thermal denaturation of five globular proteins, ribonuclease A, lysozyme, alpha-chymotrypsin, cytochrome c, and myoglobin, are quantitatively analyzed using the model. The thermodynamic parameters obtained by the least squares method reflect the global, average properties of proteins and are in good agreement with the expected values estimated from experimental and theoretical studies for model peptides. The average bond energy epsilon is well related to the tertiary structure of each protein. However, the difference in the parameters between different proteins is not observed for the cooperative energy ZJ and the chain entropy alpha. The individuality of a protein as far as its structural stability is concerned, is mainly reflected by the parameter gamma specifying the hydrophobic nature of a protein. The model is further applied in the analysis of several aspects of the structural stability of globular proteins. Denaturation induced by denaturants, salts, and pH are also explained by the model in a unified manner.     

Local sequence patterns of hydrophobicity and solvent accessibility in soluble globular proteins

We examined the variation in the solvent accessibility and hydrophobicity of the amino acids along the sequences of 58 soluble globular proteins with known tertiary structure. We found that there is a significant tendency for the accessibilities to run in clusters along the sequence but that the hydrophobicities are distributed without such nonrandom clusters. Theseresults suggest severe limitations on the power of sequence analysis tools that use average hydrophobicity scores of overlapping subsequences to predict accessibility.     

The prediction of the conformation of membrane proteins from the sequence of amino acids.

The methods of Chou & Fasman [Biochemistry (1974) 13, 211-222, 222-245] and of Lim [J. Mol. Biol. (1974)88, 857-872, 873-894] for predicting secondary structure from amino acid sequence have been applied to five predominantly helical membrane-associated peptides. The predictions from the method of Lim (1974a,b) are consistent with the experimental observations, whereas those from Chou & Fasman (1974a,b), although not inconsistent with alpha-helix, favour a beta-structure for several very hydrophobic regions. The results may be rationalized in terms of the effect of the solvent on the conformation of a polypeptide.     

Disulfide bonds and the stability of globular proteins.

An understanding of the forces that contribute to stability is pivotal in solving the protein-folding problem. Classical theory suggests that disulfide bonds stabilize proteins by reducing the entropy of the denatured state. More recent theories have attempted to expand this idea, suggesting that in addition to configurational entropic effects, enthalpic and native-state effects occur and cannot be neglected. Experimental thermodynamic evidence is examined from two sources: (1) the disruption of naturally occurring disulfides, and (2) the insertion of novel disulfides. The data confirm that enthalpic and native-state effects are often significant. The experimental changes in free energy are compared to those predicted by different theories. The differences between theory and experiment are large near 300 K and do not lend support to any of the current theories regarding the stabilization of proteins by disulfide bonds. This observation is a result of not only deficiencies in the theoretical models but also from difficulties in determining the effects of disulfide bonds on protein stability against the backdrop of numerous subtle stabilizing factors (in both the native and denatured states), which they may also affect.     

Prediction of the structural class of globular proteins from their amino acid sequence

Two-dimensional representation of consequence of 32 proteins with known three-dimensional structure has been obtained on 20 X 20 matrix of the distribution of amino acid pairs (nearest neighbours). Prediction algorithm of the structural class of globular proteins has been worked out on the basis of the comparison of 20 X 20 matrix of the distribution of amino acid pairs for the proteins of different structural classes. The accuracy of structural class predictions of 32 proteins has been carried out (all the proteins are taken from numerous ones used to obtain the algorithm).