Amino acid sequence predicts folding rate for middle-size two-state proteins

The significant correlation between protein folding rates and the sequence-predicted secondary structure suggests that folding rates are largely determined by the amino acid sequence. Here, we present a method for predicting the folding rates of proteins from sequences using the intrinsic properties of amino acids, which does not require any information on secondary structure prediction and structural topology. The contribution of residue to the folding rate is expressed by the residue's Omega value. For a given residue, its Omega depends on the amino acid properties (amino acid rigidity and dislike of amino acid for secondary structures). Our investigation achieves 82% correlation with folding rates determined experimentally for simple, two-state proteins studied until the present, suggesting that the amino acid sequence of a protein is an important determinant of the protein-folding rate and mechanism.     

Direct correlation between proteins' folding rates and their amino acid compositions: an ab initio folding rate prediction

Discovering the mechanism of protein folding, in molecular biology, is a great challenge. A key step to this end is to find factors that correlate with protein folding rates. Over the past few years, many empirical parameters, such as contact order, long-range order, total contact distance, secondary structure contents, have been developed to reflect the correlation between folding rates and protein tertiary or secondary structures. However, the correlation between proteins' folding rates and their amino acid compositions has not been explored. In the present work, we examined systematically the correlation between proteins' folding rates and their amino acid compositions for two-state and multistate folders and found that different amino acids contributed differently to the folding progress. The relation between the amino acids' molecular weight and degeneracy and the folding rates was examined, and the role of hydrophobicity in the protein folding process was also inspected. As a consequence, a new indicator called composition index was derived, which takes no structure factors into account and is merely determined by the amino acid composition of a protein. Such an indicator is found to be highly correlated with the protein's folding rate (r > 0.7). From the results of this work, three points of concluding remarks are evident. (1) Two-state folders and multistate folders have different rate-determining amino acids. (2) The main determining information of a protein's folding rate is largely reflected in its amino acid composition. (3) Composition index may be the best predictor for an ab initio protein folding rate prediction directly from protein sequence from the standpoint of practical application.     

Effect of protein backbone folding on the stability of protein-ligand complexes

The role played by the degree of folding of protein backbones in explaining the binding energetics of protein-ligand interactions has been studied. We analyzed the protein/peptide interactions in the RNase-S system in which amino acids at two positions of the peptide S have been mutated. The global degree of folding of the protein S correlates in a significant way with the free energy and enthalpy of the protein-peptide interactions. A much better correlation is found with the local contribution to the degree of folding of one amino acid residue: Thr36. This residue is shown to have a destabilizing interaction with Lys41, which interacts directly with peptide S. Another system, consisting of the interactions of small organic molecules with HIV-1 protease was also studied. In this case, the global change in the degree of folding of the protease backbone does not explain the binding energetics of protein-ligand interactions. However, a significant correlation is observed between the free energy of binding and the contribution of two amino acid residues in the HVI-1 protease: Gly49 and Ile66. In general, it was observed that the changes in the degree of folding are not restricted to the binding site of the protein chain but are distributed along the whole protein backbone. This study provides a basis for further consideration of the degree of folding as a parameter for empirical structural parametrizations of the binding energetics of protein folding and binding.     

Protein structure, neighbor effect, and a new index of amino acid dissimilarities.

Amino acids interact with each other, especially with neighboring amino acids, to generate protein structures. We studied the pattern of association and repulsion of amino acids based on 24,748 protein-coding genes from human, 11,321 from mouse, and 15,028 from Escherichia coli, and documented the pattern of neighbor preference of amino acids. All amino acids have different preferences for neighbors. We have also analyzed 7,342 proteins with known secondary structure and estimated the propensity of the 20 amino acids occurring in three of the major secondary structures, i.e., helices, sheets, and turns. Much of the neighbor preference can be explained by the propensity of the amino acids in forming different secondary structures, but there are also a number of intriguing association and repulsion patterns. The similarity in neighbor preference among amino acids is significantly correlated with the number of amino acid substitutions in both mitochondrial and nuclear genes, with amino acids having similar sets of neighbors replacing each other more frequently than those having very different sets of neighbors. This similarity in neighbor preference is incorporated into a new index of amino acid dissimilarities that can predict nonsynonymous codon substitutions better than the two existing indices of amino acid dissimilarities, i.e., Grantham's and Miyata's distances.     

Sequence periodicity and secondary structure propensity in model proteins

We explore the question of whether local effects (originating from the amino acids intrinsic secondary structure propensities) or nonlocal effects (reflecting the sequence of amino acids as a whole) play a larger role in determining the fold of globular proteins. Earlier circular dichroism studies have shown that the pattern of polar, non polar amino acids (nonlocal effect) dominates over the amino acid intrinsic propensity (local effect) in determining the secondary structure of oligomeric peptides. In this article, we present a coarse grained computational model that allows us to quantitatively estimate the role of local and nonlocal factors in determining both the secondary and tertiary structure of small, globular proteins. The amino acid intrinsic secondary structure propensity is modeled by a dihedral potential term. This dihedral potential is parametrized to match with experimental measurements of secondary structure propensity. Similarly, the magnitude of the attraction between hydrophobic residues is parametrized to match the experimental transfer free energies of hydrophobic amino acids. Under these parametrization conditions, we systematically explore the degree of frustration a given polar, non polar pattern can tolerate when the secondary structure intrinsic propensities are in opposition to it. When the parameters are in the biophysically relevant range, we observe that the fold of small, globular proteins is determined by the pattern of polar, non polar amino acids regardless of their instrinsic secondary structure propensities. Our simulations shed new light on previous observations that tertiary interactions are more influential in determining protein structure than secondary structure propensity. The fact that this can be inferred using a simple polymer model that lacks most of the biochemical details points to the fundamental importance of binary patterning in governing folding.     

Relationship between protein folding kinetics and amino acid properties

The successful prediction of protein-folding rates based on the sequence-predicted secondary structure suggests that the folding rates might be predicted from sequence alone. To pursue this question, we directly predict the folding rates from amino acid sequences, which do not require any information on secondary or tertiary structure. Our work achieves 88% correlation with folding rates determined experimentally for proteins of all folding types and peptide, suggesting that almost all of the information needed to specify a protein's folding kinetics and mechanism is comprised within its amino acid sequence. The influence of residue on folding rate is related to amino acid properties. Hydrophobic character of amino acids may be an important determinant of folding kinetics, whereas other properties, size, flexibility, polarity and isoelectric point, of amino acids have contributed little to the folding rate constant.     

Molecular basis of co-operativity in protein folding. III. Structural identification of cooperative folding units and folding intermediates.

The hierarchical partition function formalism for protein folding developed earlier has been extended through the use of three-dimensional polar and apolar contact plots. For each amino acid residue in the protein, these plots indicate the apolar and polar surfaces that are buried from the solvent, the identity of all amino acid residues that contribute to this shielding, and the magnitude of their contributions. These contact plots are then used to examine the distribution of the free energy of stabilization throughout the protein molecule. Analysis of these data allows identification of co-operative folding units and their hierarchical levels, and the identification of partially folded intermediates with a significant probability of being populated. The overall folding/unfolding thermodynamics of 12 globular proteins, for which crystallographic and experimental thermodynamics are available, is predicted within error. An energetic classification of partially folded intermediates is presented and the results compared to those cases for which structural and thermodynamic experimental information is available. Four different types of partially folded states and their structural energies are considered. (1) Local intermediates, in which only a local region of the protein loses secondary and tertiary interactions, while the rest of the protein remains intact. (2) Global intermediates, corresponding to the standard molten globule definition, in which significant secondary structure is maintained but native-like tertiary structure contacts are disrupted. (3) Extended intermediates characterized by the existence of secondary structure elements (e.g. alpha-helices) exposed to solvent. (4) Folding intermediates in proteins with two structural domains. The structure and energetics of folding intermediates of apo-myoglobin, alpha-lactalbumin, phosphoglycerate kinase and arabinose-binding protein are considered in detail.     

Amino acid interaction preferences in proteins

Understanding the key factors that influence the interaction preferences of amino acids in the folding of proteins have remained a challenge. Here we present a knowledge‐based approach for determining the effective interactions between amino acids based on amino acid type, their secondary structure, and the contact based environment that they find themselves in the native state structure as measured by their number of neighbors. We find that the optimal information is approximately encoded in a 60 × 60 matrix describing the 20 types of amino acids in three distinct secondary structures (helix, beta strand, and loop). We carry out a clustering scheme to understand the similarity between these interactions and to elucidate a nonredundant set. We demonstrate that the inferred energy parameters can be used for assessing the fit of a given sequence into a putative native state structure.     

How the folding rates of two‐ and multistate proteins depend on the amino acid properties

Proteins fold by either two‐state or multistate kinetic mechanism. We observe that amino acids play different roles in different mechanism. Many residues that are easy to form regular secondary structures (α helices, β sheets and turns) can promote the two‐state folding reactions of small proteins. Most of hydrophilic residues can speed up the multistate folding reactions of large proteins. Folding rates of large proteins are equally responsive to the flexibility of partial amino acids. Other properties of amino acids (including volume, polarity, accessible surface, exposure degree, isoelectric point, and phase transfer energy) have contributed little to folding kinetics of the proteins. Cysteine is a special residue, it triggers two‐state folding reaction and but inhibits multistate folding reaction. These findings not only provide a new insight into protein structure prediction, but also could be used to direct the point mutations that can change folding rate. Proteins 2014; 82:2375–2382. © 2014 Wiley Periodicals, Inc.     

Refolding of Myoglobins of Nonhomologous Sequences by the Island Model

The folding process of sea hare myoglobin was simulated by the island model, which does not rely on sequence homologies or statistical inference from database of known structure. Sea hare myoglobin has low sequence homology (28%), but high structural similarity, with sperm whale myoglobin, which was already simulated by the island model. Their structural similarity is shown physiochemically from the distribution of hydrophobic-residue pairs, that is, the key pairs for packing of the secondary structures. Irrelevant to the sequence homology, the secondary structures can be packed into the tertiary structure through the hydrophobic interactions among the amino acid pairs responsible for the local structure formation. The results on the two species of myoglobins indicate that, in contrast to other prediction methods, the island model is applicable to any type of protein without extra information other than the distribution of hydrophobic-residue pairs and the positions of the secondary structures. Consequently the present results provide another verification of the validity of the island model for elucidating the mechanisms of protein folding and predicting protein structures.     

H5N2亚型禽流感病毒血凝素裂解位点碱性氨基酸对应mRNA核苷酸的二级结构分析

以H5N2亚型禽流感病毒毒株血凝素蛋白裂解位点碱性氨基酸为研究对象,对其密码子偏好性和对应mRNA序列的折叠二级结构特点进行研究和分析。旨在探讨裂解位点氨基酸对应mRNA核苷酸片段的二级结构与病毒致病力的关系,希望能对禽流感病毒的研究提供一些基础性信息。将mRNA样本按照序列等步长递增的方法,用RNAstructure 4．1程序预测这些样本的动态延伸折叠二级结构。序列和结构的分析结果：裂解位点的碱性氨基酸对富含腺嘌呤的密码子有强烈偏好;与碱性氨基酸对应的mRNA片段上的核苷酸主要位于折叠二级结构的单链环区,少数位于配对螺旋区。结果表明：裂解位点氨基酸对应的mRNA核苷酸形成发夹端环的大小与其碱性氨基酸的多少具有正相关性。     

Helix, sheet, and polyproline II frequencies and strong nearest neighbor effects in a restricted coil library

A central issue in protein folding is the degree to which each residue's backbone conformational preferences stabilize the native state. We have studied the conformational preferences of each amino acid when the amino acid is not constrained to be in a regular secondary structure. In this large but highly restricted coil library, the backbone preferentially adopts dihedral angles consistent with the polyproline II conformation rather than alpha or beta conformations. The preference for the polyproline II conformation is independent of the degree of solvation. In conjunction with a new masking procedure, the frequencies in our coil library accurately recapitulate both helix and sheet frequencies for the amino acids in structured regions, as well as polyproline II propensities. Therefore, structural propensities for alpha-helices and beta-sheets and for polyproline II conformations in unfolded peptides can be rationalized solely by local effects. In addition, these propensities are often strongly affected by both the chemical nature and the conformation of neighboring residues, contrary to the Flory isolated residue hypothesis.     

Beta-sheet proteins with nearly identical structures have different folding intermediates

The folding mechanisms of two proteins in the family of intracellular lipid binding proteins, ileal lipid binding protein (ILBP) and intestinal fatty acid binding protein (IFABP), were examined. The structures of these all-beta-proteins are very similar, with 123 of the 127 amino acids of ILBP having backbone and C(beta) conformations nearly identical to those of 123 of the 131 residues of IFABP. Despite this structural similarity, the sequences of these proteins have diverged, with 23% sequence identity and an additional 16% sequence similarity. The folding process was completely reversible, and no significant concentrations of intermediates were observed by circular dichroism or fluorescence at equilibrium for either protein. ILBP was less stable than IFABP with a midpoint of 2. 9 M urea compared to 4.0 M urea for IFABP. Stopped-flow kinetic studies showed that both the folding and unfolding of these proteins were not monophasic, suggesting that either multiple paths or intermediate states were present during these processes. Proline isomerization is unlikely to be the cause of the multiphasic kinetics. ILBP had an intermediate state with molten globule-like spectral properties, whereas IFABP had an intermediate state with little if any secondary structure during folding and unfolding. Double-jump experiments showed that these intermediates appear to be on the folding path for each protein. The folding mechanisms of these proteins were markedly different, suggesting that the different sequences of these two proteins dictate different paths through the folding landscape to the same final structure.     

ω-芋螺毒素MVIIC的N及C端修饰对折叠及活性影响

 合成了 ω-芋螺毒素 MVIIC的三种 N及 C端修饰肽 ,应用高压液相色谱、CD及生物体内活性实验 ,研究了其 N及 C端修饰对折叠及活性的影响 .结果表明 :MVIIC N端用 Phe及 Ser修饰后降低其线性肽形成正确折叠的比例及结构的稳定性 ,对小鼠的脑室给药活性也相应降低 ;C端酰胺转为电负性羧基端后活性降低 ,CD谱存在显著差异 .     

A detailed consideration of a principal domain of vertebrate fibrinogen and its relatives.

Vertebrate fibrinogen is a complex multidomained protein, the structure of which has been inferred mainly from electron microscopy and amino acid sequence studies. Among its most prominent features are two terminal globules, moieties that are mostly composed of the carboxyl-terminal two-thirds of the beta and gamma chains. Sequences homologous to the latter segments are found in several other animal proteins, always as the carboxyl-terminal contributions. An alignment of 15 amino acid sequences from various fibrinogens and related proteins has been used to make judgments about secondary structure. The nature of amino acids at each position in the alignment was used to distinguish alpha helices and beta structure on the one hand from loops and turns on the other, and the resulting assignments compared with predictions of secondary structure by other methods. Additionally, constraints imposed by the locations of cystines, carbohydrate attachment residues, and proteinase-sensitive points provided further insights into the general organization of the postulated secondary structures. Other ancillary data, including the effects of bound calcium and the locations of labeled or variant residues, were also considered. An intriguing similarity to a portion of the recently reported structure of a calcium-dependent lectin is noted.     

A Reexamination of Correlations of Amino Acids with Particular Secondary Structures

Using the data from Protein Data Bank the correlations of primary and secondary structures of proteins were analyzed. The correlation values of the amino acids and the eight secondary structure types were calculated, where the position of the amino acid and the position in sequence with the particular secondary structure differ at most 25. The diagrams describing these results indicate that correlations are significant at distances between −9 and 10. The results show that the substituents on Cβ or Cγ atoms of amino acid play major role in their preference for particular secondary structure at the same position in the sequence, while the polarity of amino acid has significant influence on α-helices and strands at some distance in the sequence. The diagrams corresponding to polar amino acids are noticeably asymmetric. The diagrams point out the exchangeability of residues in the proteins; the amino acids with similar diagrams have similar local folding requirements. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Alignment of protein sequences using secondary structure: a modified dynamic programming method

A method for comparison of protein sequences based on their primary and secondary structure is described. Protein sequences are annotated with predicted secondary structures (using a modified Chou and Fasman method). Two lettered code sequences are generated (Xx, where X is the amino acid and x is its annotated secondary structure). Sequences are compared with a dynamic programming method (STRALIGN) that includes a similarity matrix for both the amino acids and secondary structures. The similarity value for each paired two-lettered code is a linear combination of similarity values for the paired amino acids and their annotated secondary structures. The method has been applied to eight globin proteins (28 pairs) for which the X-ray structure is known. For protein pairs with high primary sequence similarity (greater than 45%), STRALIGN alignment is identical to that obtained by a dynamic programming method using only primary sequence information. However, alignment of protein pairs with lower primary sequence similarity improves significantly with the addition of secondary structure annotation. Alignment of the pair with the least primary sequence similarity of 16% was improved from 0 to 37% 'correct' alignment using this method. In addition, STRALIGN was successfully applied to seven pairs of distantly related cytochrome c proteins, and three pairs of distantly related picornavirus proteins.     

Influence of the hydrophilic face on the folding ability and stability of alpha-helix bundles: relevance to the peptide catalytic activity.

Although not the sole feature responsible, the packing of amino acid side chains in the interior of proteins is known to contribute to protein conformational specificity. While a number of amphipathic peptide sequences with optimized hydrophobic domains has been designed to fold into a desired aggregation state, the contribution of the amino acids located on the hydrophilic side of such peptides to the final packing has not been investigated thoroughly. A set of self-aggregating 18-mer peptides designed previously to adopt a high level of alpha-helical conformation in benign buffer is used here to evaluate the effect of the nature of the amino acids located on the hydrophilic face on the packing of a four alpha-helical bundle. These peptides differ from one another by only one to four amino acid mutations on the hydrophilic face of the helix and share the same hydrophobic core. The secondary and tertiary structures in the presence or absence of denaturants were determined by circular dichroism in the far- and near-UV regions, fluorescence and nuclear magnetic resonance spectroscopy. Significant differences in folding ability, as well as chemical and thermal stabilities, were found between the peptides studied. In particular, surface salt bridges may form which would increase both the stability and extent of the tertiary structure of the peptides. The structural behavior of the peptides may be related to their ability to catalyze the decarboxylation of oxaloacetate, with peptides that have a well-defined tertiary structure acting as true catalysts.     

Determining the 3D structure of human ASC2 protein involved in apoptosis and inflammation

ASC2 structure has been well defined by 1141 NOE experimental restraints. The model consists of five alpha helices. alpha-Helices are connected by short random structure loops. The sole exception is the loop connecting helices 2 and 3, which has a 20-residue length. Folding generally agrees with the folding of recently published death domain structures in which alpha-helix structures have been reported. In spite of structural similarity, amino acid sequence homology with the most similar protein (ASC1) is just 64%. DD, DED, and CASP protein structures present six helices along their sequences; ASC2 presents 5 well-defined helices due to long distance restraints. However, a helical fragment was observed between amino acids 38 and 42 (representing helix 3) in the death domains when constructing the model.     

Incorporating Secondary Features into the General form of Chou's PseAAC for Predicting Protein Structural Class

Protein structure information is very useful for the confirmation of protein function. The protein structural class can provide information for protein 3D structure analysis, causing the conformation of the protein overall folding type plays a significant part in molecular biology. In this paper, we focus on the prediction of protein structural class which was based on new feature representation. We extract features from the Chou-Fasman parameter, amino acid compositions, amino acids hydrophobicity features, polarity information and pair-coupled amino acid composition. The prediction result by the Support vector machine (SVM) classifier shows that our method is better than some others.