基于HP模型的大麻素受体CB2折叠研究

Cannabinoid receptor Type 2(简称CB2)是大麻素受体的一种亚型,因为其无中枢神经副作用,不会产生成瘾性及耐受性,显示出了非常好的开发前景和潜在的应用价值。其作为免疫调节剂、神经保护剂和抗癌药等将具有巨大市场价值。目前,CB2蛋白的空间结构还未被测定出来,对于CB2的折叠问题研究也开展的较少,为了研究大麻素受体亚型蛋白CB2的折叠问题以及方便更多的研究人员对CB2空间结构和相关药理特性的研究,本文提出了一种基于HP模型的折叠求解方法。通过使用回溯机制和蒙特卡罗方法对此优化问题进行求解,算法可有效的在全局范围内进行寻找最优解,避免了掉入局部最优问题。实验结果表明,本文方法获取的CB2蛋白空间构象具有较低的能量值,折叠情况较好。     

蛋白质的折叠

重点介绍了蛋白质折叠的热力学控制学说和动力学控制学说,简单介绍了几种蛋白质折叠模型并分析了多肽链在体内进行快速折叠的原因。     

错误折叠蛋白质的聚集效应及其对策

错误折叠蛋白质大量聚集将引发蛋白质构象紊乱症(protein conformational disorder,PCD)。目前的研究发现,蛋白质聚集过程中的中间体(纤维前体,pre—fibfillar)导致细胞膜结构受损,从而诱发细胞凋亡。根据这一原理设计出的抗纤维前体抗体和干扰肽可以作为PCD治疗的一般性方法。此外,该就消除错误折叠蛋白质聚集的研究方向作了展望。     

蛋白质折叠机理的理论研究

简要综述了近年来蛋白质折叠机理的理论研究。首先回顾了蛋白质折叠理论的发展历程,然后对折叠中间体的研究现状作了较详细的介绍。同时,对折叠机理理论研究中的几种理论模型和模拟算法作了细致评述,分析了其现状和存在的问题。最后,总结和讨论了折叠机理理论研究的现存问题及研究热点,并展望了该领域研究的发展趋势。     

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

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

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

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

基于HNP三态模型及相对熵方法的蛋白质折叠研究

把20种氨基酸简化为3类:疏水氨基酸(hydrophobic,H)、亲水氨基酸(hydrophilic,P)及中性氨基酸(neutral,N),每个氨基酸简化为一个点,用其C!原子来代替.采用非格点模型,以相对熵作为优化函数,进行蛋白质三维结构预测.为了与基于相对熵方法的蛋白质设计工作进行统一,采用了新的接触强度函数.选用蛋白质数据库中的天然蛋白质作为测试靶蛋白,结果表明,采用该模型和方法取得了较好的结果,预测结构相对于天然结构的均方根偏差在0.30～0.70nm之间.该工作为基于相对熵及HNP模型的蛋白质设计研究打下了基础.     

蛋白质折叠速率预测研究进展

蛋白质折叠速率预测是当今生物物理学最具挑战性的课题之一。近年来,该领域的研究取得了很大的进展,提出了许多经验参数,例如:接触序、长程序、总接触距离、链拓扑参数、二级结构含量、有效长度、螺旋参数、n-阶接触距离等。这些参数都和蛋白质的折叠速率有很好的相关性,基于这些参数的各种预测方法所得到的预测结果也与实验数据较好地吻合。     

蛋白质工程和蛋白质折叠

蛋白质一级结构决定着高级结构。蛋白质肽链在适宜条件下会自动卷曲形成其相应的高级结构,即自动发生蛋白质折叠,其自动发生的原因和过程仍不十分清楚,但是随着蛋白质工程的日益兴起,这些与折叠有关的问题也愈显重要,就此已有文章进行过讨论[1,2]。反之,如把新兴的蛋白质工程手段(尤其是基因定点诱变技术)应用来研究这些折叠问题,必将推动蛋白质折叠的研究。本文将就蛋白质折叠与蛋白质工程相互影响的一些例子进行讨论。     

α－苦瓜子蛋白晶体结构中的次级相互作用的研究Ⅱ疏水相互作用的分析

α－苦瓜子蛋白结构模型中共有疏水残基１４１个,占践基总数的５７％。其中１２８个疏水残基有相互作用,占疏水残基的９１％。疏水残基之间的相互作用与其是否在α－螺旋或β片层内没有多大关系。本文着重分析了α－螺旋间和α螺旋与β结构之间的疏水相互作用;讨论了α５和大β片层中疏水残基与亲本残基的分布特点,这一特点对该蛋白分子最后折叠成特定的空间结构有重要影响。类似的情况在同源蛋白天花粉蛋白中也是存在的。     

α－苦瓜子蛋白晶体结构中的次级相互作用的研究Ⅱ疏水相互作用的分析

α－苦瓜子蛋白结构模型中共有疏水残基１４１个,占践基总数的５７％。其中１２８个疏水残基有相互作用,占疏水残基的９１％。疏水残基之间的相互作用与其是否在α－螺旋或β片层内没有多大关系。本文着重分析了α－螺旋间和α螺旋与β结构之间的疏水相互作用;讨论了α５和大β片层中疏水残基与亲本残基的分布特点,这一特点对该蛋白分子最后折叠成特定的空间结构有重要影响。类似的情况在同源蛋白天花粉蛋白中也是存在的。     

4,4'-Dianilino-1,1'-binaphthyl-5,5'-sulfonate, a novel molecule having chaperone-like activity

4,4'-Dianilino-1,1'-binaphthyl-5,5'-sulfonate (bis-ANS) and 1-anilinonaphthalene-8-sulfonate (ANS) are hydrophobic probes that are widely used in protein folding studies, using their capacity to bind to hydrophobic regions of partially unfolded proteins and in turn leading to an increase in fluorescence. Here we reveal a novel chaperone-like activity for bis-ANS, which acted as a highly effective inhibitor for the thermal- or chemical-induced aggregation of alcohol dehydrogenase, insulin or the whole cell extract of Escherichia coli, with ANS showing a much weaker effect. The studies to elucidate the mechanism underlying this activity show that bis-ANS is able to form stable soluble aggregates with the denaturing proteins and dramatically increase its fluorescence intensity upon incubation with aggregation-prone proteins. Moreover, we found that bis-ANS is able to prevent the heat inactivation of citrate synthase. These observations suggest that bis-ANS is able to block the exposed hydrophobic surfaces to suppress protein aggregation, acting in a way similar to what small heat shock proteins (one sub-class of molecular chaperones) do. The data presented here, together with the report that bis-ANS was able to suppress the amyloid formation of the prion peptide [J. Biol. Chem. 279 (2004) 5346], suggest that this molecule may be used as a potential protein stabilizer in addition to its current application as a hydrophobic probe.     

Initiation of chromosomal DNA replication which is stimulated without oversupply of DnaA protein in Escherichia coli

Summary The temperature-sensitive dnaA46 mutation in Escherichia coli can be phenotypically suppressed at 42° C by oversupply of GroELS proteins, and the suppressed cells grow extremely slowly at 30° C. We found that the phenotype of dnaA46 showing this cold sensitivity was dominant over the phenotype of dnaA ⁺, and could not be rescued by introduction of oriC-independent replication systems. These results suggest that the cold sensitivity was not caused by a simple defect in replication. When a growing culture of a dnaA46 strain with a GroELS-overproducing plasmid was shifted from 42° to 30° C in the presence of chloramphenicol, the chromosomal DNA replicated excessively. Initiation of replication occurred at the site of oriC repeatedly four or five times during a 4 h incubation period without concomitant protein synthesis, indicating an excessive capacity for initiation. Such overreplication did not take place at 42° C in the suppressed dnaA46 strain, or at either temperature in GroELS-oversupplied dnaA ⁺ cells. No significant difference was detected between the cellular content of DnaA protein in suppressed cells where the initiation capacity was abnormally high, and that in wild-type cells in which the initiation capacity was normal. Thus, DnaA protein might function in vivo through some phase control mechanism for initiation, apart from a simple regulation by its total amount. A possible mechanism is proposed based on the participation of GroELS proteins in protein folding.A preliminary account of this work was presented at the Annual Meeting of the Molecular Biology Society of Japan in 1989.     

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain–side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein’s structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein–protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at . The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation. Figure Fuzzy-oil-drop model applied to identify the ligation site in lysozyme complexed with N-acetylglucosamine (PDB ID:1LMQ) in form of hydrophobicity deficiency (ΔH) profile and three-dimensional distribution of on protein surface     

Funnel-like organization in sequence space determines the distributions of protein stability and folding rate preferred by evolution

To understand the physical and evolutionary determinants of protein folding, we map out the complete organization of thermodynamic and kinetic properties for protein sequences that share the same fold. The exhaustive nature of our study necessitates using simplified models of protein folding. We obtain a stability map and a folding rate map in sequence space. Comparison of the two maps reveals a common organizational principle: optimality decreases more or less uniformly with distance from the optimal sequence in the sequence space. This gives a funnel-shaped optimality surface. Evolutionary dynamics of a sequence population on these two maps reveal how the simple organization of sequence space affects the distributions of stability and folding rate preferred by evolution.     

Role of Hydrophobic Core on the Thermal Stability of Proteins—Molecular Dynamics Simulations on a Single Point Mutant of Sso7d

Abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Enumerating Designing Sequences in the HP Model

The hydrophobic/polar HP model on the square lattice has been widely used toinvestigate basics of protein folding. In the cases where all designing sequences (sequences with unique ground states) were enumerated without restrictions on the number of contacts, the upper limit on the chain length N has been 18–20 because of the rapid exponential growth of thenumbers of conformations and sequences. We show how a few optimizations push this limit by about 5 units. Based on these calculations, we study the statistical distribution of hydrophobicity along designing sequences. We find that the average number of hydrophobic and polar clumps along the chains is larger for designing sequences than for random ones, which is in agreement with earlier findings for N 18 and with results for real enzymes. We also show that this deviation from randomness disappears if the calculations are restricted to maximally compact structures.