一个识别蛋白质折叠模式的SVM分类器

蛋白质折叠模式识别是一种分析蛋白质结构的重要方法。以序列相似性较低的蛋白质为训练集,提取蛋白质序列信息频数及疏水性等信息作为折叠类型特征,从SCOP数据库中已分类蛋白质构建1 393种折叠模式的数据集,采用SVM预测蛋白质1 393种折叠模式。封闭测试准确率达99.612 2%,基于SCOP的开放测试准确率达79.632 9%。基于另一个权威测试集的开放测试折叠准确率达64.705 9%,SCOP类准确率达76.470 6%,可以有效地对蛋白质折叠模式进行预测,从而为蛋白质从头预测提供参考。     

多酚化合物对蛋白质构象疾病的抑制作用研究进展

多酚化合物是广泛存在于中草药或植物性食物中的生物活性分子。大量研究表明多酚化合物对蛋白质构象疾病具有防治作用，如可通过疏水作用抑制蛋白质或多肽形成淀粉样纤维或解聚淀粉样纤维，或者通过其抗氧化性降低淀粉样纤维对细胞的氧化损伤。本文综述了近年来中草药或食物来源的天然多酚化合物体外抑制蛋白质或多肽错误折叠和淀粉样纤维形成，以及在动物疾病模型体内的抑制效果及作用机理的研究进展，以期为该类化合物预防和减轻蛋白质构象疾病的应用提供基础参考。     

蚂蚁群落优化算法在蛋白质折叠二维亲-疏水格点模型中的应用

氨基酸的亲疏水格点模型是研究蛋白质折叠的一种重要的简化模型,其优化问题是一个非确定型的多项式问题。采用蚂蚁群落优化算法对这一问题进行了研究,对测试数据的计算结果表明,在一定规模下,此算法能够有效地获得亲-疏水格点模型的最优解,其效率优于传统的Monte Carlo仿真等方法。     

模拟蛋白质折叠过程的新算法研究

蛋白质折叠过程模拟是当前蛋白质研究领域的一个难点问题。针对这一问题,提出了一个描述蛋白质折叠过程的算法-拟蛇算法,并且从分子振荡和分子动力学理论两个方面来证明该算法的核心函数是可行和正确的。经过实验总结出所有蛋白质空间结构都可以通过两种类型函数构造出来,提出了描述蛋白质折叠过程模型。与其它蛋白质折叠过程模拟算法的实验结果比较表明,拟蛇算法所构造的空间结构能量值最小、相似度最好。进而说明拟蛇算法和蛋白质折叠过程模型在描述蛋白质折叠过程方面具有明显优势。     

基于序列疏水值震荡的折叠速率预测

蛋白质折叠速率的正确预测对理解蛋白质的折叠机理非常重要。本文从伪氨基酸组成的方法出发,提出利用序列疏水值震荡的方法来提取蛋白质氨基酸的序列顺序信息,建立线性回归模型进行折叠速率预测。该方法不需要蛋白质的任何二级结构、三级结构信息或结构类信息,可直接从序列对蛋白质折叠速率进行预测。对含有62个蛋白质的数据集,经过Jack．knife交互检验验证,相关系数达到0．804,表示折叠速率预测值与实验值有很好的相关性,说明了氨基酸序列信息对蛋白质折叠速率影响重要。同其他方法相比,本文的方法具有计算简单,输入参数少等特点。     

分子伴侣及其在蛋白质折叠中的作用研究进展

蛋白质折叠是一个复杂的、动态的过程,蛋白质的折叠不是自发的,需要其他物质的帮助.了解分子伴侣在蛋白质折叠过程中的的作用,有助于进一步研究蛋白质折叠机制.本文介绍了分子伴侣及其分类,重点综述了各类分子伴侣在蛋白质折叠中的机制,并提出了研究分子伴侣在蛋白质折叠中的作用的重要意义.     

基于结构比较的蛋白质模建系统及其评估：Ⅲ.灵敏的蛋…

本文按二级结构、疏水性和侧链氢链把蛋白质的局部结构环境划分为６４类,按环境依赖的残基替代频数表构成残基与环境的兼容性分数表,可用来评估一个蛋白质的整体折叠正确与否和检测蛋白质折叠的局部错误。     

未折叠蛋白质的普遍初始热力学亚稳态

探索和理解蛋白质折叠问题一直是分子生物学、结构生物学和生物物理学的终极挑战.未折叠的蛋白质应该存在一种普遍初始热力学亚稳态,否则无法解释蛋白质是如何在剧烈的热振动干扰下完成快速精确折叠的.本文通过分析水溶液环境和蛋白质折叠的相关性,揭示了一种由水分子屏蔽效应引起的未折叠蛋白质的普遍初始热力学亚稳态,该亚稳态的存在是水溶液环境中水分子的物理性质决定,并赋予未折叠蛋白质抵抗热扰动和避免错误折叠的能力.我们通过研究已发表的实验数据和建立分子模型,找到了该初始热力学亚稳态存在的相关证据,并推测了该亚稳态导致蛋白质精确折叠的相关物理学机制.     

蛋白质折叠速率数据集的构建及分析

近年来,随着高精度的蛋白质折叠速率实验数据的不断积累,使得从蛋白质折叠速率角度研究蛋白质折叠机制的理论工作者,迎来了前所未有的机遇和挑战。然而,却有约100多个蛋白质的折叠速率实验数据散落在2个数据库和若干文献中。为了方便今后的理论工作分析,作者将这些散落数据汇集整理出来,构建了一个包含109个非冗余单体野生型蛋白质的折叠速率数据集,称为PFRD109(protein folding rate dataset 109)。PFRD109所包含的109个蛋白质中,有69个二态蛋白和40个多态蛋白,折叠速率从10-4到106s-1,跨度为10个数量级。链长最短的为16 aa,最长为390 aa,二态蛋白平均长度为78 aa,多态蛋白平均长度为137 aa。当前,生物信息学对蛋白质折叠速率的研究,主要集中于寻找与折叠速率和折叠动力学相关的各种生化参数或拓扑参数,进而实现对蛋白质折叠速率和蛋白质折叠动力学类型的预测。因此,本文还针对PFRD109数据集,就这两个方面进行了一些参数的统计分析。     

电导法研究表面活性剂与蛋白质的相互作用

较系统地研究了溶液中离子型表面活性剂与蛋白质相互作用时电导率的变化,并根据实验现象,得出表面活性剂与蛋白质作用的两种模式－疏水作用模式和电荷作用模式。表面活性剂采用疏水作用模式与蛋白质结合时,蛋白质的二硫键逐一断裂,三级结构逐渐打开,电导率曲线出现一些小“平台”,采用电荷作用模式吸附蛋白质时,首先形成疏水复合体,产生白色浑浊,随离子型表面活性剂浓度的增加,疏水复合体转变成亲水复合体,白色浑浊     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

NMR characterization of a peptide model provides evidence for significant structure in the unfolded state of the villin headpiece helical subdomain

The villin headpiece subdomain (HP36) is the smallest naturally occurring protein that folds cooperatively. The protein folds on a microsecond time scale. Its small size and very rapid folding have made it a popular target for biophysical studies of protein folding. Temperature-dependent one-dimensional (1D) NMR studies of the full-length protein together with CD and 1D NMR studies of the 21-residue peptide fragment (HP21) derived from HP36 have shown that there is significant structure in the unfolded state of HP36 and have demonstrated that HP21 is a good model of these interactions. Here, we characterized the model peptide HP21 in detail by two-dimensional NMR. Strongly upfield shifted C(alpha) protons, the magnitude of the 3J(NH,alpha) coupling constants, and the pattern of backbone-backbone and backbone-side chain NOEs indicate that the ensemble of structures populated by HP21 contains alpha-helical structure and native as well as non-native hydrophobic contacts. The hydrogen-bonded secondary structure inferred from the NOEs is, however, not sufficient to confer significant protection against amide H-D exchange. These studies indicate that there is significant secondary structure and hydrophobic clustering in the unfolded state of HP36. The implications for the folding of HP36 are discussed.     

Detailed folding structures of M-lycotoxin-Hc1a and its mutageneses using 2D HP model

Although the hydrophobic-polar (HP) model was proposed a decade ago, it applies almost to no real-case study because of its intense computation. In this study, a 2D HP model was applied to study the folding structures of M-lycotoxin-Hc1a, an antimicrobial peptide, in order to get full pictures of its numerous folding structures. The normalised hydrophobicity index was used to convert M-lycotoxin-Hc1a and its six mutageneses into HP sequences, and then the 2D HP model was used to compute all the possible folding structures (3²⁴ = 282,429,536,481), and finally the normalised hydrophobicity index was used to distinguish the native state. The results showed that M-lycotoxin-Hc1a had 6 and 138 folding structures at their native state with the minimal energy of ? 13 at pH 2 and pH 7 when glycine served as hydrophobic amino acid. When glycine serves as polar amino acid, M-lycotoxin-Hc1a had 12 and 54 folding structures at their native state with the minimal energy of ? 12 and ? 13 at pH 2 and pH 7, respectively. This study advanced the knowledge on how to apply the HP model to real-life study, and how the mutageneses influenced the folding structures of M-lycotoxin-Hc1a, their native states and minimal energy at different pH levels.     

Hydrophobic Core Formation and Dehydration in Protein Folding Studied by Generalized-Ensemble Simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 Å. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental ϕ-value analysis.     

Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.     

A branch and bound algorithm for the protein folding problem in the HP lattice model

A branch and bound algorithm is proposed for the two-dimensional protein folding problem in the HP lattice model. In this algorithm, the benefit of each possible location of hydrophobic monomers is evaluated and only promising nodes are kept for further branching at each level. The proposed algorithm is compared with other well-known methods for 10 benchmark sequences with lengths ranging from 20 to 100 monomers. The results indicate that our method is a very efficient and promising tool for the protein folding problem.     

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.     

Identification of the PXW sequence as a structural gatekeeper of the headpiece C-terminal subdomain fold

The HeadPiece (HP) domain, present in several F-actin-binding multi-domain proteins, features a well-conserved, solvent-exposed PXWK motif in its C-terminal subdomain. The latter is an autonomously folding subunit comprised of three alpha-helices organised around a hydrophobic core, with the sequence motif preceding the last helix. We report the contributions of each conserved residue in the PXWK motif to human villin HP function and structure, as well as the structural implications of the naturally occurring Pro to Ala mutation in dematin HP. NMR shift perturbation mapping reveals that substitution of each residue by Ala induces only minor, local perturbations in the full villin HP structure. CD spectroscopic thermal analysis, however, shows that the Pro and Trp residues in the PXWK motif afford stabilising interactions. This indicates that, in addition to the residues in the hydrophobic core, the Trp-Pro stacking within the motif contributes to HP stability. This is reinforced by our data on isolated C-terminal HP subdomains where the Pro is also essential for structure formation, since the villin, but not the dematin, C-terminal subdomain is structured. Proper folding can be induced in the dematin C-terminal subdomain by exchanging the Ala for Pro. Conversely, the reverse substitution in the villin C-terminal subdomain leads to loss of structure. Thus, we demonstrate a crucial role for this proline residue in structural stability and folding potential of HP (sub)domains consistent with Pro-Trp stacking as a more general determinant of protein stability.     

Competition between intradomain and interdomain interactions: a buried salt bridge is essential for villin headpiece folding and actin binding

Villin-type headpiece domains are ～70 residue motifs that reside at the C-terminus of a variety of actin-associated proteins. Villin headpiece (HP67) is a commonly used model system for both experimental and computational studies of protein folding. HP67 is made up of two subdomains that form a tightly packed interface. The isolated C-terminal subdomain of HP67 (HP35) is one of the smallest autonomously folding proteins known. The N-terminal subdomain requires the presence of the C-terminal subdomain to fold. In the structure of HP67, a conserved salt bridge connects N- and C-terminal subdomains. This buried salt bridge between residues E39 and K70 is unusual in a small protein domain. We used mutational analysis, monitored by CD and NMR, and functional assays to determine the role of this buried salt bridge. First, the two residues in the salt bridge were replaced with strictly hydrophobic amino acids, E39M/K70M. Second, the two residues in the salt bridge were swapped, E39K/K70E. Any change from the wild-type salt bridge residues results in unfolding of the N-terminal subdomain, even when the mutations were made in a stabilized variant of HP67. The C-terminal subdomain remains folded in all mutants and is stabilized by some of the mutations. Using actin sedimentation assays, we find that a folded N-terminal domain is essential for specific actin binding. Therefore, the buried salt bridge is required for the specific folding of the N-terminal domain which confers actin-binding activity to villin-type headpiece domains, even though the residues required for this specific interaction destabilize the C-terminal subdomain.     

Wang–Landau sampling of the interplay between surface adsorption and folding of HP lattice proteins

Generic features associated with the adsorption of proteins on solid surfaces are reviewed within the framework of the hydrophobic-polar (HP) lattice protein model. The thermodynamic behaviour and structural properties of various HP protein sequences interacting with attractive surfaces have been studied using extensive Wang–Landau sampling with different types of surfaces, each of which attracts either: all monomers, only hydrophobic (H) monomers or only polar (P) monomers, respectively. Consequently, different types of folding behaviour occur for varied surface strengths. Analysis of the combined patterns of various structural observables, e.g. the derivatives of the number of interaction contacts, together with the specific heat, leads to the identification of fundamental categories of folding and transition hierarchies. We also inferred a connection between the transition categories and the relative surface strengths, i.e. the ratios of the surface attractive strengths to the intra-chain attraction among H monomers. Thus, we believe that the folding hierarchies and identification scheme are generic for different HP sequences interacting with attractive surfaces, regardless of the chain length, sequence or surface attraction.