Approaching a complete classification of protein secondary structure

A complete classification of types of the protein secondary structure is developed on the basis of computer analysis of the crystallographic structural data deposited in the protein Data Bank. The majority of amino acid residues fall into five conformation types. A conclusion is drawn that the number of sequence variants of torsion angles phi, psi in globular proteins is limited and is essentially less than the number of possible amino acid sequences for this chain length. Along with alpha-helix and beta-structure, the distribution analysis assigning every maximum of distribution of amino acid conformations on Ramachandran map to a certain type of the secondary structure exposed a third type of the secondary structure that was previously neglected. This type of the structure is extended left-handed helical conformation, designated as mobile (M-) conformation. A full set of M-conformation fragments that seems to play a major role in protein globule dynamics has been obtained, a small radius of correlation for the polypeptide chain in M-conformation is demonstrated. It explains a prevalence of short segments of mobile conformation revealed in globular proteins. For secondary structure types, the frequency of occurrence of amino acid residues has been computed.     

Third type of secondary structure: noncooperative mobile conformation. Protein Data Bank analysis

Analysis of 68 proteins from Protein Data Bank disclosed a new widely spread type of the secondary structure that is designated as mobile (M-) conformation. Helical parameters of M-conformation are close to the poly-L-proline II type helix. Its occurrence in globular proteins approximates that of the beta-sheet. The angles corresponding to the position of the M-conformation maximum in distribution of amino acid residues on a conformational map are phi: -65 degrees, psi: 140 degrees. Unique features and high occurrence in proteins make it possible to distinguish the M-conformation as an independent third type of the secondary structure in globular proteins, that should be included in the present classification.     

Streptokinase is a flexible multi-domain protein

The structure of streptokinase in solution has been studied by dynamic light scattering, small-angle X-ray scattering and circular dichroism spectroscopy. The Stokes' radius and radius of gyration of the protein monomer are 3.58 nm and 4.03 nm, respectively. The maximum intraparticle distance of the molecule is 14 nm. More than half of the amino acids of the molecule are organized in regular secondary structures. The X-ray scattering curve, the results from dynamic light scattering, and the finding that at least 50% of the amino acid residues are organized in regularly folded secondary structures are consistent with the following structural model. Streptokinase consists of four compact, separately folded, domains linked by mobile segments of the protein chain. The molecule exhibits the conformation of a flexible string-of-beads in solution.     

Conformational preference functions for predicting helices in membrane proteins

A suite of FORTRAN programs, PREF, is described for calculating preference functions from the data base of known protein structures and for comparing smoothed profiles of sequence-dependent preferences in proteins of unknown structure. Amino acid preferences for a secondary structure are considered as functions of a sequence environment. Sequence environment of amino acid residue in a protein is defined as an average over some physical, chemical, or statistical property of its primary structure neighbors. The frequency distribution of sequence environments in the data base of soluble protein structures is approximately normal for each amino acid type of known secondary conformation. An analytical expression for the dependence of preferences on sequence environment is obtained after each frequency distribution is replaced by corresponding Gaussian function. The preference for the α-helical conformation increases for each amino acid type with the increase of sequence environment of buried solvent-accessible surface areas. We show that a set of preference functions based on buried surface area is useful for predicting folding motifs in α-class proteins and in integral membrane proteins. The prediction accuracy for helical residues is 79% for 5 integral membrane proteins and 74% for 11 α-class soluble proteins. Most residues found in transmembrane segments of membrane proteins with known α-helical structure are predicted to be indeed in the helical conformation because of very high middle helix preferences. Both extramembrane and transmembrane helices in the photosynthetic reaction center M and L subunits are correctly predicted. We point out in the discussion that our method of conformational preference functions can identify what physical properties of the amino acids are important in the formation of particular secondary structure elements. © 1993 John Wiley & Sons, Inc.     

The disordered mobile loop of GroES folds into a defined beta-hairpin upon binding GroEL

The GroES mobile loop is a stretch of approximately 16 amino acids that exhibits a high degree of flexible disorder in the free protein. This loop is responsible for the interaction between GroES and GroEL, and it undergoes a folding transition upon binding to GroEL. Results derived from a combination of transferred nuclear Overhauser effect NMR experiments and molecular dynamics simulations indicate that the mobile loop adopts a beta-hairpin structure with a Type I, G1 Bulge turn. This structure is distinct from the conformation of the loop in the co-crystal of GroES with GroEL-ADP but identical to the conformation of the bacteriophage-panned "strongly binding peptide" in the co-crystal with GroEL. Analysis of sequence conservation suggests that sequences of the mobile loop and strongly binding peptide were selected for the ability to adopt this hairpin conformation.     

COPS--cis/trans peptide bond conformation prediction of amino acids on the basis of secondary structure information

SUMMARY: COPS predicts for all 20 naturally occurring amino acids whether the peptide bond in a protein is in cis or trans conformation. The algorithm is based only on secondary structure information of amino acid triplets without considering the amino acid sequence information. Conformation parameters are derived from solved 3D structures deposited in the PDB and led to propensities based on modified Chou-Fasman parameters. COPS analyses amino acid triplets taking only their respective secondary structure into consideration and upon application of a set of rules utilizing the conformation parameters, the N-terminal peptide bond conformation of the middle residue is predicted. COPS was tested on a random selection of protein datasets. AVAILABILITY: The COPS program and further information are freely available from the FMP website at http://www.fmp-berlin.de/nmr/cops CONTACT: labudde@fmp-berlin.de.     

Microcalorimetric study of polypeptides in the conformation of the left helix of the poly-L-proline II type

Microcalorimetry shows that polypeptide partial heat capacity changes linearly with temperature in the region of M-conformation existence and that there is a peculiarity in the transition region where this structure is exhausted, without cooperative heat absorption.     

Characterization of fibronectin from fetal human plasma in comparison with adult plasma fibronectin

Peptides containing fewer than 50 amino acids show little ordered structure under physiological conditions. In this paper it is shown that in the receptor environment, secondary structure could be induced in small peptides that involves 87% of all the amino acid residues. The statistical methods of Chou and Fasman are used to predict the conformation of 41 peptide hormones or neuromodulators in the proteinaceous environment of the receptor, and four distinct conformational groupings are elucidated. beta-bend, beta-structure and alpha-helical conformation are possible for distinct groups of linear peptides, and disulfide bridge containing peptides show a common beta-bend beta-structure conformation at the receptor. In the predicted receptor conformation, the peptides show hydrophobic and hydrophilic domains that must reflect the distribution of corresponding regions in the ligand-binding site of the receptor. The predicted ligand conformation should allow a more rational approach to interpreting existing structure activity studies and the design of new analogs of pharmacological interest.     

Prediction of protein secondary structure based on residue pairs

The GOR program for predicting protein secondary structure is extended to include triple correlation. A score system for a residue pair to be at certain conformation state is derived from the conditional weight matrix describing amino acid frequencies at each position of a window flanking the pair under the condition for the pair to be at the fixed state. A program using this score system to predict protein secondary structure is established. After training the model with a learning set created from PDB_SELECT, the program is tested with two test sets. As a method using single sequence for predicting secondary structures, the approach achieves a high accuracy near 70%.     

Conformational prediction and circular dichroism studies on ribosomal protein S4 from Escherichia coli.

The conformation of ribosomal protein S4 from Escherichia coli has been studied by circular dichroism (CD) and shown to possess unique conformation free in solution. The near ultraviolet spectrum suggests the existence of unique tertiary structural environment for the aromatic amino acid residues. The far ultraviolet spectrum gives an estimation of its secondary structure which is 32% alpha-helix and 14% beta-structure in reconstitution buffer at 25 degrees C. The conformation of S4 has been predicted from its sequence, and two models are presented here. An attempt is made to correlate these two molecular models with the available physicochemical data concerning the shape, conformation, and possible RNA binding site of protein S4.     

The Genetic Algorithm and the Conformational Search of Polypeptides and Proteins

Abstract

The genetic algorithm is a technique of function optimization derived from the principles of evolutionary theory. We have adapted it to perform conformational search on polypeptides and proteins. The algorithm was first tested on several small polypeptides and the 46 amino acid protein crambin under the AMBER potential energy function. The probable global minimum conformations of the polypeptides were located 90% of the time and a non-native conformation of crambin was located that was 150kcal/mol lower in potential energy than the minimized crystal structure conformation. Next, we used a knowledge-based potential function to predict the structures of melittin, pancreatic polypeptide, and crambin. A 2.31 Å ΔRMS conformation of melittin and a 5.33 Å ΔRMS conformation of pancreatic polypeptide were located by genetic algorithm-based conformational search under the knowledge-based potential function. Although the ΔRMS of pancreatic polypeptide was somewhat high, most of the secondary structure was correct. The secondary structure of crambin was predicted correctly, but the potential failed to promote packing interactions. Finally, we tested the packing aspects of our potential function by attempting to predict the tertiary structure of cytochrome b ₅₆₂ given correct secondary structure as a constraint. The final predicted conformation of cytochrome b ₅₆₂ was an almost completely extended continuous helix which indicated that the knowledge-based potential was useless for tertiary structure prediction. This work serves as a warning against testing potential functions designed for tertiary structure prediction on small proteins.     

Exploring the sequence patterns in the alpha-helices of proteins

This paper reports an extensive sequence analysis of the alpha-helices of proteins. alpha-Helices were extracted from the Protein Data Bank (PDB) and were divided into groups according to their sizes. It was found that some amino acids had differential propensity values for adopting helical conformation in short, medium and long alpha-helices. Pro and Trp had a significantly higher propensity for helical conformation in short helices than in medium and long helices. Trp was the strongest helix conformer in short helices. Sequence patterns favoring helical conformation were derived from a neighbor-dependent sequence analysis of proteins, which calculated the effect of neighboring amino acid type on the propensity of residues for adopting a particular secondary structure in proteins. This method produced an enhanced statistical significance scale that allowed us to explore the positional preference of amino acids for alpha-helical conformations. It was shown that the amino acid pair preference for alpha-helix had a unique pattern and this pattern was not always predictable by assuming proportional contributions from the individual propensity values of the amino acids. Our analysis also yielded a series of amino acid dyads that showed preference for alpha-helix conformation. The data presented in this study, along with our previous study on loop sequences of proteins, should prove useful for developing potential 'codes' for recognizing sequence patterns that are favorable for specific secondary structural elements in proteins.     

A CD study of the α-helix nucleation hypothesis

One of the dilemmas in predicting the secondary structure of proteins from their amino acid propensity for a given conformation is the presence of all amino acids in all types of secondary structure, regardless of their propensity for that specific structure. One explanation is the nucleation hypothesis that only a few residues with a strong propensity for the secondary structure, such as the α-helix structure, initiates its formation and propagates the structure through indifferent sequences until strong breakers terminate the growth on both ends. Eight 15-mer peptides were studied to examine the α-helix nucleation hypothesis. The nucleation sequence of VAEAK, with high helix propensity, was mixed with an indifferent sequence of TSDSR in all possible permutations. From the percent α-helix structure derived from the CD at 222 nm, it appears that helicity does not propagate through the indifferent sequence. © 1994 John Wiley & Sons, Inc.     

Conformational studies of the human complex-forming glycoprotein, heterogeneous in charge: protein HC

Human complex-forming glycoprotein, heterogeneous in charge, is a single polypeptide chain widely distributed in physiological fluids. The conformation of the protein has been studied with attention to the secondary and tertiary structures. Circular dichroism and predictive methods from the amino acid sequence have been employed for the characterization of the secondary structure. This is composed of 20% alpha-helix, 21% beta-structure, 29% beta-turns, 30% aperiodic conformation, and an average number of residues per helical segment of nine. Titration of the protein indicated the existence of two groups for the tyrosine residues, each of them composed of three and five residues. The four tryptophan residues of the molecule are located in two different polarity microenvironments, according to the fluorescence studies. These observations are corroborated by studying the hydropathic profile of the protein. From this study, three different domains are observed in the protein, one of them being exposed and containing the main part of the unordered structure of the molecule. The chromophore naturally associated with the protein has been resolved in three fluorescent units not dependent on the protein conformation. These bands have been observed centered around 290, 360, and 410 nm, which do not correspond to any described chromophore.     

Clustering of amino acids for protein secondary structure prediction

Simple hidden Markov models are proposed for predicting secondary structure of a protein from its amino acid sequence. Since the length of protein conformation segments varies in a narrow range, we ignore the duration effect of length distribution, and focus on inclusion of short range correlations of residues and of conformation states in the models. Conformation-independent and -dependent amino acid coarse-graining schemes are designed for the models by means of proper mutual information. We compare models of different level of complexity, and establish a practical model with a high prediction accuracy.     

Analyses of the Sequence and Structural Properties Corresponding to Pentapeptide and Large Palindromes in Proteins

The analyses of 3967 representative proteins selected from the Protein Data Bank revealed the presence of 2803 pentapeptide and large palindrome sequences with known secondary structure conformation. These represent 2014 unique palindrome sequences. 60% palindromes are not associated with any regular secondary structure and 28% are in helix conformation, 11% in strand conformation and 1% in the coil conformation. The average solvent accessibility values are in the range between 0–155.28 Ų suggesting that the palindromes in proteins can be either buried, exposed to the solvent or share an intermittent property. The number of residue neighborhood contacts defined by interactions ≤ 3.2 Ǻ is in the range between 0–29 residues. Palindromes of the same length in helix, strand and coil conformation are associated with different amino acid residue preferences at the individual positions. Nearly, 20% palindromes interact with catalytic/active site residues, ligand or metal ions in proteins and may therefore be important for function in the corresponding protein. The average hydrophobicity values for the pentapeptide and large palindromes range between -4.3 to +4.32 and the number of palindromes is almost equally distributed between the negative and positive hydrophobicity values. The palindromes represent 107 different protein families and the hydrolases, transferases, oxidoreductases and lyases contain relatively large number of palindromes.     

A critique of the utility of the prediction of protein secondary structure

The Chou-Fasman predictive algorithm for determining the secondary structure of proteins from the primary sequence is reviewed. Many examples of its use are presented which illustrate its wide applicability, such as predicting (a) regions with the potential for conformational change, (b) sequences which are capable of assuming several conformations in different environments, (c) effects of single amino acid mutations, (d) amino acid replacements in synthesis of peptides to bring about a change in conformation, (e) guide to the synthesis of polypeptides with definitive secondary structure,e.g. signal sequences, (f) conformational homologues from varying sequences and (g) the amino acid requirements for amphiphilicα-helical peptides.     

A designed Zn2+-binding amphiphilic polypeptide: energetic consequences of pi-helicity.

The pi-helix is a secondary structure with 4.4 amino acids per helical turn. Although it was proposed in 1952, no experimental support for its existence was obtained until the mid-1980s. While short peptides are unlikely to assume a marginally stable secondary structure spontaneously, they might do so in the presence of appropriate structural constraints. In this paper, we describe a peptide that is designed to assume a pi-helical conformation when stabilized by cetyltrimethylammonium bromide (CTAB) micelles and Zn(2+). In the designed peptide, lipophilic amino acids are placed such that it would be amphiphilic in the pi-helical, but not in the alpha-helical, conformation. Also, two His residues are incorporated with i, i + 5 spacing, designed to allow binding of Zn(2+) in a pi-helical but not an alpha-helical conformation. The peptide was found to form moderately stable monolayers at the air-water interface, with a collapse pressure that almost doubled when there was Zn(2+) in the subphase. Also, CTAB micelles induced a marked increase in the helicity of the peptide. In 50% TFE, the peptide had a CD spectrum consistent with an alpha-helical structure. The addition of 1 mM Zn(2+) to this solvent caused a saturable decline in ellipticity to approximately half of its original value. The peptide also bound Zn(2+) when it was bound to CTAB micelles, with Zn(2+) again inducing a decrease in ellipticity. The peptide had slightly greater affinity for Zn(2+) in the presence of the CTAB than in a 50% TFE solution (K(d) = 3.1 x 10(-4) M in CTAB and 2.3 x 10(-4) M in TFE). van't Hoff analysis indicated that thermal denaturation of the peptide in 50% TFE containing 1 mM Zn(2+) was associated with both enthalpic and entropic changes that were greater than those in the absence of Zn(2+). These observations are all consistent with the proposal that the peptide assumed a pi-helical conformation in the presence of Zn(2+) and CTAB micelles, and has allowed the stability of this rare conformation to be assessed.     

Evolved cellular automata for protein secondary structure prediction imitate the determinants for folding observed in nature

We demonstrate the first application of cellular automata to the secondary structure predictions of proteins. Cellular automata use localized interactions to simulate global phenomena, which resembles the protein folding problem where individual residues interact locally to define the global protein conformation. The protein's amino acid sequence was input into the cellular automaton and rules for updating states were evolved using a genetic algorithm. An optimized accuracy (Q3) for the RS126 and CB513 dataset of 58.21% and 56.51%, respectively, could be obtained. Thus, the current work demonstrates the applicability of a rather simple algorithm on a problem as complex as protein secondary structure prediction.