Statistical analysis of the distribution of amino acid residues among helical and non-helical regions in globular proteins

Intrinsically disordered proteins: Chronology of a discovery

Intrinsic disorder is a new reality that appears to penetrate every corner of modern protein science. It is difficult to imagine that only 20 years ago the situation was completely different, and almost nobody had heard about 'structure-less' but functional proteins. As a matter of fact, for many at that time, this idea was completely heretical when viewed in light of the then dominating lock-and-key model describing the protein structure-function relationship, where a unique amino acid sequence defines a unique crystal-like 3D structure that serves as a prerequisite for a unique function of a protein. It seems like the entire field of protein intrinsic disorder has magically emerged at the turn of the century due to a revelation to a small group of researchers. Although this may very well be true, literature shows that the first observations contradicting the lock-and-key view of protein functionality started to appear almost immediately after this model was proposed. The goal of this article is to provide a brief chronology (though admittedly a subjective one) of the events in the field of protein science that eventually culminated in the discovery of the protein intrinsic disorder phenomenon. The entire process represents a good example of the "dwarf standing on the shoulders of giants" (Latin: nanos gigantum humeris insidentes) metaphor, where the truth is discovered by building on previous discoveries.

The Future of Protein Secondary Structure Prediction Was Invented by Oleg Ptitsyn

When Oleg Ptitsyn and his group published the first secondary structure prediction for a protein sequence, they started a research field that is still active today. Oleg Ptitsyn combined fundamental rules of physics with human understanding of protein structures. Most followers in this field, however, use machine learning methods and aim at the highest (average) percentage correctly predicted residues in a set of proteins that were not used to train the prediction method. We show that one single method is unlikely to predict the secondary structure of all protein sequences, with the exception, perhaps, of future deep learning methods based on very large neural networks, and we suggest that some concepts pioneered by Oleg Ptitsyn and his group in the 70s of the previous century likely are today’s best way forward in the protein secondary structure prediction field.

Comparison of Peptide Ion Conformers Arising from Non-Helical and Helical Peptides Using Ion Mobility Spectrometry and Gas-Phase Hydrogen/Deuterium Exchange

The dominant gas-phase conformer of [M+3H]³⁺ ions of the model peptide acetyl-PSSSSKSSSSKSSSSKSSSSK has been examined with ion mobility spectrometry (IMS), gas-phase hydrogen deuterium exchange (HDX), and mass spectrometry (MS) techniques. The [M+3H]³⁺ peptide ions are observed predominantly as a relatively compact conformer type. Upon subjecting these ions to electron transfer dissociation (ETD), the level of protection for each amino acid residue in the peptide sequence is assessed. The overall per-residue deuterium uptake is observed to be relatively more efficient for the neutral residues than for the model peptide acetyl-PAAAAKAAAAKAAAAKAAAAK. In comparison, the N-terminal and C-terminal regions of the serine peptide show greater relative protection compared with interior residues. Molecular dynamics (MD) simulations have been used to generate candidate structures for collision cross section and HDX reactivity matching. Hydrogen accessibility scoring (HAS) for select structural candidates from MD simulations has been used to suggest conformer types that could contribute to the observed HDX patterns. The results are discussed with respect to recent studies employing extensive MD simulations of gas-phase structure establishment of a peptide system. Graphical Abstract ᅟ.

Application of expert networks for predicting proteins secondary structure

The present study utilizes expert neural networks for the prediction of proteins secondary structure. We use three independent networks, one for each structure (alpha, beta and coil) as the first-level processing unit; decision upon the chosen structure for each residue is carried out by a second-level, post-processing unit, which utilizes the Chou and Fasman frequency values Falpha and Fbeta in order to strengthen and/or deplete the probability of the specific structure under investigation. The highest prediction case was 76%. Our method requires primitive computational means and a relatively small training set, while still been comparable to previous work. It is not meant to be an alternative to the determination of secondary structure by means of free energy minimization, integration of dynamic equations of motion or crystallography, which are expensive, time-consuming and complicated, but to provide additional constrains, which might be considered and incorporated into larger computing setups in order to reduce the initial search space for the above methods.

A theoretical model of Aquifex pyrophilus flagellin: implications for its thermostability

Aquifex pyrophilus is a flagellated hyperthermophilic eubacterial species that grows optimally at 85 degrees C. The thermostable A. pyrophilus flagellar filament is primarily composed of a single protein called flagellin (FlaA). The N- and C-terminal sequence regions of FlaA are important for self-assembly and share high sequence similarity with mesophilic bacterial flagellins. We have developed a predictive 3D-structure of FlaA, using the published structure of mesophilic Salmonella typhimurium flagellin (FliC) as a template and analyzed it with respect to possible determinants of thermostability. A sequence comparison of FlaA and FliC revealed a +7.0% increase in FlaA hydrophobic residues, a +0.6% increase in charged residues and a corresponding decrease of -6.0% in polar residues. The FlaA N- and C-termini also have higher proportions of hydrophobic and charged residues at the expense of polar residues and higher non-polar surface areas. Thus, a predominant stabilizing factor in FlaA appears to be increased hydrophobicity, which often confers greater rigidity to proteins. Fewer intramolecular ion pairs were observed in FlaA than FliC, although an increase in the positive charge potential of the FlaA D0 and D1 domains was also observed; increased intermolecular salt bridges may also contribute to the thermal stability of the oligomeric flagellar fiber. [Figure: see text].

Amino Acid Propensities Revisited

Statistical analysis of amino acid patterns in approximately 160,000 alpha-helices in experimentally determined structures revealed di-, tri-, and tetrapeptides, whose frequencies deviate most from the statistical model. Importantly, some sequences were never found in alpha- helices. This fact was detected initially with tripeptides, where nearly 1% of the possible sequences were never seen in the helical segments. For tetrapeptides, this effect is very strong and significant; almost 43% of the possible sequences never appear in alpha-helices. It is possible that there are some steric and energetic restrictions that do not allow these tetrameric amino acid sequences to form alpha-helical structure.

Thermodynamics Of α‐Helix Formation

The alpha-helix was the first proposed and experimentally confirmed secondary structure. The elegant simplicity of the alpha-helical structure, stabilized by hydrogen bonding between the backbone carbonyl oxygen and the peptide amide four residues away, has captivated the scientific community. In proteins, alpha-helices are also stabilized by the so-called capping interactions that occur at both the C- and the N-termini of the helix. This chapter provides a brief historical overview of the thermodynamic studies of the energetics of helix formation, and reviews recent progress in our understanding of the thermodynamics of helix formation.

Elucidation of factors responsible for enhanced thermal stability of proteins: a structural genomics based study

Understanding the molecular basis for the enhanced stability of proteins from thermophiles has been hindered by a lack of structural data for homologous pairs of proteins from thermophiles and mesophiles. To overcome this difficulty, complete genome sequences from 9 thermophilic and 21 mesophilic bacterial genomes were aligned with protein sequences with known structures from the protein data bank. Sequences with high homology to proteins with known structures were chosen for further analysis. High quality models of these chosen sequences were obtained using homology modeling. The current study is based on a data set of models of 900 mesophilic and 300 thermophilic protein single chains and also includes 178 templates of known structure. Structural comparisons of models of homologous proteins allowed several factors responsible for enhanced thermostability to be identified. Several statistically significant, specific amino acid substitutions that occur going from mesophiles to thermophiles are identified. Most of these are at solvent-exposed sites. Salt bridges occur significantly more often in thermophiles. The additional salt bridges in thermophiles are almost exclusively in solvent-exposed regions, and 35% are in the same element of secondary structure. Helices in thermophiles are stabilized by intrahelical salt bridges and by an increase in negative charge at the N-terminus. There is an approximate decrease of 1% in the overall loop content and a corresponding increase in helical content in thermophiles. Previously overlooked cation-pi interactions, estimated to be twice as strong as ion-pairs, are significantly enriched in thermophiles. At buried sites, statistically significant hydrophobic amino acid substitutions are typically consistent with decreased side chain conformational entropy.

A helix propensity scale based on experimental studies of peptides and proteins

The average globular protein contains 30% alpha-helix, the most common type of secondary structure. Some amino acids occur more frequently in alpha-helices than others; this tendency is known as helix propensity. Here we derive a helix propensity scale for solvent-exposed residues in the middle positions of alpha-helices. The scale is based on measurements of helix propensity in 11 systems, including both proteins and peptides. Alanine has the highest helix propensity, and, excluding proline, glycine has the lowest, approximately 1 kcal/mol less favorable than alanine. Based on our analysis, the helix propensities of the amino acids are as follows (kcal/mol): Ala = 0, Leu = 0.21, Arg = 0.21, Met = 0.24, Lys = 0.26, Gln = 0.39, Glu = 0.40, Ile = 0.41, Trp = 0.49, Ser = 0.50, Tyr = 0. 53, Phe = 0.54, Val = 0.61, His = 0.61, Asn = 0.65, Thr = 0.66, Cys = 0.68, Asp = 0.69, and Gly = 1.

Tn5 transposase mutants that alter DNA binding specificity

Tn5 transposase (Tnp) binds to Tn5 and IS50 end inverted repeats, the outside end (OE) and the inside end (IE), to initiate transposition. We report the isolation of four Tnp mutants (YH41, TP47, EK54 and EV54) that increase the OE-mediated transposition frequency and enhance the binding affinity of Tnp for OE DNA. In addition, two of the Tnp mutants (TP47 and EK54) appear to be change-of-specificity mutants, since they alter the recognition of OE versus IE relative to the wild-type Tnp. EK54 enhances OE recognition but decreases IE recognition. TP47 enhances both OE and IE recognition but with a much greater enhancement for IE than for OE. This change-of-specificity effect of TP47 is observed only when TP47 Tnp is synthesized in cis to the DNA that contains the ends. We propose that Lys54 makes a favorable interaction with an OE-specific nucleotide pair(s), while Pro47 may cause a more favorable interaction with an IE-specific nucleotide pair(s) than it does with the corresponding OE-specific nucleotide pair(s). A model to explain the preference of TP47 Tnp for the IE in cis but not in trans is proposed.

Conformation of thymosin beta 9 in water/fluoroalcohol solution determined by NMR spectroscopy

The conformation of thymosin beta 9 in solution of 40% (v/v) 1,1,1,3,3,3-hexafluoro-2-propanol-d2 in water has been investigated by two-dimensional 1H-nmr spectroscopy. Under this condition thymosin beta 9 adopts an ordered structure. The determination of the conformation of the peptide was based on a set of 304 approximate interproton distance constraints derived from nuclear Overhauser enhancement measurements. The conformation of thymosin beta 9 includes two helical regions from residues 4 to 27 and 32 to 41. The two helices are separated by a poorly defined loop region between amino acids 28 and 31; the N-terminus of thymosin beta 9 shows random-coil structure only.

The alpha-helix as seen from the protein tertiary structure: a 3-D structural classification

Helices are selected from globular protein structures defined at high resolution by X-ray analysis. We cluster alpha-helices in two ways: according to their position in the tertiary structure by considering patterns of solvent inaccessible residues and according to the arc of the solvent inaccessible face. For each class of helices we have defined propensities for amino acids at each position; these can be used to calculate templates for recognition of a member of that class. The analysis provides a basis for the prediction of alpha-helices and estimation of their approximate position in a protein tertiary structure. It also provides an approach to estimating the probability of finding amino acid sequences as helices in solution and in a folded protein, thus indicating those helices that might be involved in nucleation of protein folding.

Molecular dynamics simulations of a winter flounder "antifreeze" polypeptide in aqueous solution

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute-solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge.

Use of conditional probabilities for determining relationships between amino acid sequence and protein secondary structure

The conditional probability, P(sigma/x), is a statement of the probability that the value of sigma will be found given the prior information that a value of x has been observed. Here sigma represents any one of the secondary structure types, alpha, beta, tau, and rho for helix, sheet, turn, and random, respectively, and x represents a sequence attribute, including, but not limited to: (1) hydropathy; (2) hydrophobic moments assuming helix and sheet; (3) Richardson and Richardson helical N-cap and C-cap values; (4) Chou-Fasman conformational parameters for helix, P alpha, for sheet, P beta, and for turn, P tau; and (5) Garnier, Osguthorpe, and Robson (GOR) information values for helix, I alpha, for sheet, I beta, for turn, I tau, and for random structure, I rho. Plots of P(sigma/x) vs. x are demonstrated to provide information about the correlation between structure and attribute, sigma and x. The separations between different P(sigma/x) vs. x curves indicate the capacity of a given attribute to discriminate between different secondary structural types and permit comparison of different attributes. P(alpha/x), P(beta/x), P(tau/x) and P(rho/x) vs. x plots show that the most useful attributes for discriminating helix are, in order: hydrophobic moment assuming helix greater than P alpha much greater than N-cap greater than C-cap approximately I alpha approximately I tau. The information value for turns, I tau, was found to discriminate helix better than turns. Discrimination for sheet was found to be in the following order: I beta much greater than P beta approximately hydropathy greater than I rho approximately hydrophobic moment assuming sheet. Three attributes, at their low values, were found to give significant discrimination for the absence of helix: I alpha approximately P alpha approximately hydrophobic moment assuming helix. Also, three other attributes were found to indicate the absence of sheet: P beta much greater than I rho approximately hydropathy. Indications of the absence of sigma could be as useful for some applications as the indication of the presence of sigma.

Does helix dipole have any role in binding metal ions in protein structures?

Positions of metal-binding residues with respect to helical terminii in protein structures have been analyzed in order to determine if the location of these ligands is influenced by the helix dipole. Most ligands do not show any preference for the amino- or carboxy-terminus of a helix. For steric reasons, peptide ligands can be located only at the C-terminus. The availability of a second ligand residue closely placed along the sequence may be of more importance, rather than the electrostatic interaction involving helix dipole, in cases where ligands are found near the C-terminus. The location of heme-binding histidine residues at the C-terminus may be due to the steric requirements of the heme group and also the intrahelical hydrogen bond that the residues can form at this position. Considerations based on such geometrical features and not just the helix dipole may help us to understand the observed distribution of charged residues along alpha-helices and the favorable role these amino acids have on folding of isolated helices.

Physical reasons for secondary structure stability: alpha-helices in short peptides

It was recently found that some short peptides (including C- and S-peptide fragments of RNase A) can have considerable helicity in solution, which was considered to be surprising. Does the observed helicity require a new explanation, or is it consistent with previous understanding? In this work we show that this helicity is consistent with the physical theory of secondary structure based on an extension of the conventional Zimm-Bragg model. Without any special modifications, this theory explains reasonably well almost all the experimentally observed dependencies of helicity on pH, temperature, and amino acid replacements. We conclude that the observed "general level" of helicity of C- and S-peptides (5-30% at room temperature and 10-50% near 0 degrees C) is "normal" for short peptides consisting mainly of helix-forming and helix-indifferent residues. The helicity is modified by a multitude of weak specific side chain interactions, many of which are taken into account by the present theory; some discrepancies between the theory and experiment can be explained by weak side-chain-side chain interactions that were neglected. A reasonable coincidence of the theory with experiment suggests that it had been used to investigate the role of local interactions in the formation of alpha-helical "embryos" in unfolded protein chains.

The general mitochondrial matrix processing protease from rat liver: structural characterization of the catalytic subunit

A critical step in the import of nuclear-encoded precursor proteins into mitochondria involves proteolytic cleavage of their amino-terminal leader peptides by processing proteases found in the mitochondrial matrix. We report here the characterization of the general matrix processing protease from rat liver mitochondria. The final enzyme preparation consisted of two polypeptides, a catalytically active 55-kDa subunit and a 52-kDa one. To deduce the complete primary structure of the 55-kDa subunit, we first sequenced its mature amino terminus and several tryptic peptides derived from the pure protein. Next, using mixed oligonucleotide primers that had sequences based on two of these peptides, we synthesized a partial cDNA probe by selective amplification of liver RNA with the polymerase chain reaction. The amplified probe was then used to obtain a nearly full-length clone from a rat liver cDNA library. This cDNA codes for 508 amino acid residues, including 16 residues of an amino-terminal leader peptide, the cleavage site of which is located two polypeptide bonds downstream from an arginine residue. The mature portion has a predicted molecular mass of 55.2 kDa; it shows 36% identity with the mitochondrial processing peptidases of Saccharomyces cerevisiae and Neurospora crassa. A conserved structural feature is a putative, negatively charged alpha-helix, located in the amino-terminal half of the subunit; this element might be important for the recognition of positively charged leader peptides characteristic of mitochondrial precursor proteins.

Distribution of charged residues stabilizes individual helices in myoglobins

The effect of the distribution of charged residues on stability of alpha helices in isolated peptides and in globular proteins exemplified by myoglobins from 62 different species is discussed. A highly simplified set of rules is used to account for the interaction of charged groups with the dipole of an alpha helix. Only the position and sign of a charge with respect to the center of the helix and its ability to participate in intrahelical salt bridges determine its effect. These rules lead to a linear correlation between the helicity in variant C-peptide helices from RNAse and the extent to which the charge distribution opposes the helix dipole. Of the sample of 496 helices in the myoglobins studied, 456 exhibit arrangements of charges which oppose the effective dipole moment of the helix according to this calculation. A number of variants occur which leave the backbone moment of helices A-D unchanged, or even add to it. However no such variants exist in the sequences of helices E-H. We suggest that the E, F, G and H helices in myoglobins which show the strongest reversal of the helix dipole participate in the structures of early intermediates in folding of the chain. Stable helix structures should be more likely to occur in these isolated sequences also, and introduction of charge alterations in helices E to H should affect the initial refolding rate of mutant myoglobins.

Further studies of the helix dipole model: effects of a free alpha-NH3+ or alpha-COO- group on helix stability

Interactions between the alpha-helix peptide dipoles and charged groups close to the ends of the helix were found to be an important determinant of alpha-helix stability in a previous study. The charge on the N-terminal residue of the C-peptide from ribonuclease A was varied chiefly by changing the alpha-NH2 blocking group, and the correlation of helix stability with N-terminal charge was demonstrated. An alternative explanation for some of those results is that the succinyl and acetyl blocking groups stabilize the helix by hydrogen bonding to an unsatisfied main-chain NH group. The helix dipole model is tested here with peptides that contain either a free alpha-NH3+ or alpha-COO- group, and no other charged groups that would titrate with similar pKa's. This model predicts that alpha-NH3+ and alpha-COO- groups are helix-destabilizing and that the destabilizing interactions are electrostatic in origin. The hydrogen bonding model predicts that alpha-NH3+ and alpha-COO- groups are not themselves helix-destabilizing, but that an acetyl or amide blocking group at the N- or C-terminus, respectively, stabilizes the helix by hydrogen bonding to an unsatisfied main-chain NH or CO group. The results are as follows: (1) Removal of the charge from alpha-NH3+ and alpha-COO- groups by pH titration stabilizes an alpha-helix. (2) The increase in helix stability on pH titration of these groups is close to the increase produced by adding an acetyl or amide blocking group.(ABSTRACT TRUNCATED AT 250 WORDS)

Amino acid preferences for specific locations at the ends of alpha helices

A definition based on alpha-carbon positions and a sample of 215 alpha helices from 45 different globular protein structures were used to tabulate amino acid preferences for 16 individual positions relative to the helix ends. The interface residue, which is half in and half out of the helix, is called the N-cap or C-cap, whichever is appropriate. The results confirm earlier observations, such as asymmetrical charge distributions in the first and last helical turn, but several new, sharp preferences are found as well. The most striking of these are a 3.5:1 preference for Asn at the N-cap position, and a preference of 2.6:1 for Pro at N-cap + 1. The C-cap position is overwhelmingly dominated by Gly, which ends 34 percent of the helices. Hydrophobic residues peak at positions N-cap + 4 and C-cap - 4.

Secondary structure predictions and medium range interactions

Several authors have proposed that predictions of protein secondary structure derived from statistical information about the known structures can be improved when information about neighboring residues participating in short and medium range interactions is included. A substantial improvement shown here indicates that current methods of including this information are not more successful than methods that do not. Evaluations of the Chou and Fasman method (Adv. Enzymol. 47 (1978) 45-148), that does not include information about interactions (except in averaging), have shown it to be about 49% correct for three states (helix, beta-sheet and undefined). In comparison, the method of Garnier et al. (J. Mol. Biol. 120 (1978) 97-120), that explicitly includes information about neighboring residues, has an accuracy of 57% residues correct for three states. However, we have obtained an 8% improvement for predictions of secondary structure based on the algorithm by Chou and Fasman. The improvements are obtained by eliminating many rules and by choosing the best decision constants for structure assignments. The simplified method described here is 57% correct for three states using preference values calculated in 1978.

Amino acid distribution in protein secondary structures

The compositional distribution of the twenty amino acids was examined for particular positions within secondary structures (alpha-helices, beta-strands, and turns) taken from a 44-protein sample. Correlation coefficients calculated between positional composition of the amino acids and various of their physico-chemical characteristics indicated considerable asymmetry in the properties of the residues comprising regions within and adjacent to secondary structures, modes of helix formation, physical parameters most sensitive to the buriedness of residues in beta-strands, and possible improvements in the accuracy of secondary structure prediction methodologies.

Similarities of protein topologies: evolutionary divergence, functional convergence or principles of folding?

(A)Evolutionary similarities of protein structuresTwo decades have passed from the time that the three dimensional structure of the first globular protein, sperm whale myoglobin, was decoded (Kendrewet al.1960). Its structure, which now looks so simple and habitual, then seemed to be unusually complicated. The decoding of the subsequent proteins, lysozyme (Blakeet al.1965), ribonuclease (Kartha, Bello & Harker, 1967), chymotrypsin (Matthewset al.1967), carboxypeptidase (Lipscombet al.1969) redoubled the feeling of amazement and even of some confusion before the extremely complicated, intricate and, above all, absolutely unlike protein structures. Some consolation against this background was the evident and far-reaching similarity between the three-dimensional structures of myoglobin and hemoglobin subunits (Perutz, Kendrew & Watson, 1965) and an analogous similarity between the structures of chymotrypsin and other serine proteases, elastase (Shotton & Watson, 1970) and trypsin (Stroud, Kay & Dickerson, 1972). However this similarity was easily explained by the far-reaching homology between the primary structures of myoglobin and hemoglobin and between the primary structures of serine proteases.

A re-examination of the conformation of secretin in water

An examination of a series of peptides corresponding to partial sequences of secretin and the application of empirical conformational parameters to the sequence has reopened the question about the position of the short helical stretch in the folded chain. In earlier studies, this was tentatively placed near the N-terminus, while newly accumulated evidence points to the C-terminal area. The possible role of ion pairs in the stabilization of the folding was investigated by the preparation and examination of synthetic analogues of the C-terminal tricosapeptide part of secretin. The carboxyl groups of glutamic acid (in position 9) and of aspartic acid (in position 15) were replaced by carboxamides. The 9-glutamine and the 15-asparagine analogues show a significant decrease in helical character. This loss of "structure" is even more pronounced in the 9-glutamine-15-asparagine tricosapeptide. Thus, ion-pair formation is indeed implicated as one of the forces which stabilize the folded conformation of the chain. Possible correlations between biological activity and secretin-like architecture were studied on several smooth muscle preparations.

Analysis of code relating sequences to conformation in globular prtoeins. Theory and application of expected information

1. An information theory analysis of the folding of a globular protein is proposed. 2. The folding is seen as a transfer of information between two messages, the primary sequence and the biologically active conformation. 3. It is shown how the information transferred was estimated by inspection of proteins of known primary sequence and conformation. 4. In this estimation, concerted use of subjective (Bayesian) probabilities leads to a more robust approach which can be employed whether the number of proteins of known sequence and conformation is large or small. 5. Further, it is demonstrated that the problem then becomes a very simple algebraic formulation for information estimates. 6. Finally, it is shown how this process of information theory analysis can be reversed to predict the conformation of a protein by using its primary sequence and the above information estimates obtained from other proteins. 7. The present paper provides the theoretical basis for the derivation and application of a stereochemical alphabet (Robson & Pain, 1974a,c), and for an investigation of the effects of residues on the conformations of their neighbours (Robson & Pain, 1974b).

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.