蛋白质共进化分析研究进展

一些对蛋白质活性很重要的残基在进化过程中是高度保守的,另有一些残基通过共进化来维持蛋白质结构和功能上的稳定。由于共进化残基分析可在未知蛋白质结构时,仅依据序列推断残基间的相互作用,因此在蛋白质结构和功能预测上具有重要的研究意义。当前分析共进化残基的方法主要有基于相关系数的方法、基于微扰理论的方法、参数检验法等。然而,由于存在蛋白质系统进化的背景干扰,目前共进化残基分析的精度仍有待进一步提高。本文概述了蛋白质共进化分析的方法及其研究进展,并对其发展趋势进行了预测。     

残基突变对冷休克蛋白热稳定性影响的计算分析研究(英文)

残基突变是提高蛋白质热稳定性最直接有效的方式。在本文中,我们选取一对冷休克蛋白质作为研究对象,其中一个来自嗜温的Bacillus subtilis(Bs-CspB),另一个来自嗜热的Bacillus caldolyticus(Bc-Csp),这两个蛋白质在序列和结构上具有高度的相似性,但两者的耐热能力却相差很大。我们利用全原子模型计算残基突变前后蛋白质的自由能和氨基酸之间相互作用能的变化,分析残基突变对冷休克蛋白热稳定性的影响。通过对比两个蛋白质对应位置上残基的能量,我们成功鉴别出对Bc-Csp的高热稳定性有突出贡献的残基。我们计算了这些残基突变前后,该残基的静电相互作用和范德华相互作用的变化,以分析该残基对Bc-Csp高热稳定性的主要贡献。同时,我们分析了离子键对蛋白质热稳定性的贡献。我们的计算结果和实验结果吻合得很好,关键在于利用该方法可以详细地说明残基突变影响蛋白质热稳定性的根本原因。本文为研究残基突变对蛋白质热稳定性的影响提供了一种计算思路和方法,并有助于设计具有高耐热能力的蛋白质。     

距离约束的HMC采样算法在蛋白质结构预测中的运用

蛋白质结构预测中,采样是指在构象空间中生成具有最小自由能的状态。传统的采样方法是对自由度直接赋值。这种方法在处理较少的残基时能取得好的效果。但是对于包含100个残基以上的蛋白质结构,由于构象空间的急剧增长,难以得到理想的结构。本文引入深度学习中的HMC(Hybrid Monte Carlo)采样方法,以概率分布为依据对蛋白质的自由度进行采样,能够对包含100、200甚至更多个残基的蛋白质结构进行采样。并且,在采样的过程中加入残基间的距离约束,使得一个结构中,相对于Rosetta的ab initio最多有75%(平均40%)的残基对得到优化,满足距离约束。     

化学交联在蛋白质结构研究中的应用

在交联剂分子两端各有一个相同或不同的活性基团,它们能与蛋白质侧链上的氨基、巯基、羟基等形成共价交联.利用交联反应,可以测定寡聚蛋白质的亚基数量、研究蛋白质的高级结构、测量氨基酸残基间的距离及研究蛋白质间的相互作用.     

p53蛋白质Gly249和Ser249替换型三维结构的分子动力学研究

利用p53蛋白质核心区晶体结构作分子动力学发现,除了生化方面的稳定性之外,该区还具有分子力学上的高度稳定性.在此基础上作的R249残基替换分子动力学研究显示,p53蛋白质核心区249位点精氨酸被其他残基替换后能引起p53蛋白质核心区L2、L3结构域间的密切联系趋于松散,正常的空间构象发生改变并使整个核心区结构稳定性受到破坏.这一研究从三维结构变化上,直观地解释了R249残基替换造成的p53蛋白质免疫和生化特性改变的结构机制.     

PrPC与蛋白质X结合位点可及性的初步探讨

根据PDB提供的PrPC的原子坐标,利用MSMS程序,对PrPC氨基酸残基溶剂可及表面积进行了计算和分析.结果表明:(1) PrPC氨基酸残基可及性具有总体一致性特点;(2) PrPC蛋白质序列中非保守残基与种间屏障有一定关系;(3) 在PrPC向PrP Sc转变过程中,由于蛋白质X的结合,PrPC可能会出现一定的构象变化,这种变化利于PrPC向PrPSc发生转变.     

蛋白质-核酸复合物界面氨基酸与核苷酸偏好性分析

蛋白质-核酸相互作用机制到目前还不是很清楚,尤其是蛋白质与RNA的相互作用。目前,可得到的蛋白质-核酸复合物结构数据不断增多,作者收集了Protein Data Bank数据库中所有的蛋白质-核酸复合物结构数据,对复合物中结合残基和结合核苷酸的偏好性进行了统计分析。发现：1）不同功能的蛋白质-核酸复合物间的结合残基数量存在显著差异;2）在蛋白 质-DNA和蛋白质-RNA复合物界面,碱性氨基酸都是最受欢迎的;3）氨基酸的极性大小及方向在决定它是否与RNA分子进行结合时起到重要的作用,同时发现氨基酸侧链形成的空间位阻会影响氨基酸残基与RNA分子的相互作用;4）随着定义结合残基距离阈值的增大,其氨基酸使用的特异性降低,而受欢迎与不受欢迎的氨基酸种类均没有变化。     

蛋白质的翻译后修饰与肿瘤代谢

蛋白质的翻译后修饰(post-translational modifications, PTMs)可改变蛋白质的稳定性、活性或细胞定位,是其功能调控的重要方式。常见的修饰类型包括赖氨酸残基上的乙酰化、泛素化、甲基化以及丝氨酸/苏氨酸残基上的磷酸化、糖基化修饰等。研究发现蛋白质的PTMs与肿瘤的代谢及生长密切相关,而靶向蛋白质PTMs与肿瘤代谢调控的治疗策略或将具有良好的临床应用前景。本文将对蛋白质的PTMs与肿瘤代谢及生长调控间相互作用的研究进展做一阐述。     

固有无序蛋白与蛋白质相互作用位点残基特征分析

本文对固有无序蛋白(IDPs)与其他蛋白质相互作用位点残基特征进行了研究.首先在数据库中选出满足条件的109条IDPs蛋白质链及与其他配体蛋白形成的299个IDPs-蛋白质复合物,然后提取复合物中作为相互作用位点的IDPs-蛋白质残基.这109条IDPs链中共含有50 031个氨基酸残基,其中处于作用位点的残基有4 822个.通过分析发现,20种氨基酸在形成IDPs-蛋白质相互作用位点残基时具有不同的倾向性,根据形成作用位点残基的倾向性,20种氨基酸可分成三大类:倾向型氨基酸(ILE、LEU、ARG、PHE、TYR、MET、TRP)、中间型氨基酸(GLN、GLU、THR、LYS、VAL、ASP、HIS)、非倾向型氨基酸(PRO、SER、GLY、ALA、ASN、CYS).研究结果还进一步表明,不同氨基酸在有序区域与无序区域形成IDPs-蛋白质作用位点残基的倾向性不同.其中,氨基酸TRP、LEU、ILE、CYS在有序和无序区域形成作用位点残基的差异性尤为明显,而氨基酸GLU、PHE、HIS、ALA则基本没有多大差别.对IDPs-蛋白质相互作用位点残基理化特征进行分析发现:疏水性强、侧链净电荷量较少、极性较小、溶剂可及性表面积较大、侧链体积较大、极化率较大的氨基酸比较倾向于形成作用位点残基.主成分分析结果显示,残基的极化率、侧链体积和溶剂可及表面积对作用位点残基影响最大.     

不同类型氨基酸网络参量与蛋白质折叠的关系

理论和实验研究表明,蛋白质天然拓扑结构对其折叠过程具有重要的影响．采用复杂网络的方法分析蛋白质天然结构的拓扑特征,并探索蛋白质结构特征与折叠速率之间的内在联系．分别构建了蛋白质氨基酸网络、疏水网、亲水网、亲水-疏水网以及相应的长程网络,研究了这些网络的匹配系数(assortativity coefficient)和聚集系数(clustering coefficient)的统计特性．结果表明,除了亲水-疏水网,上述各网络的匹配系数均为正值,并且氨基酸网和疏水网的匹配系数与折叠速率表现出明显的线性正相关,揭示了疏水残基间相互作用的协同性有助于蛋白质的快速折叠．同时,研究发现疏水网的聚集系数与折叠速率有明显的线性负相关关系,这表明疏水残基间三角结构(triangle construction)的形成不利于蛋白质快速折叠．还进一步构建了相应的长程网络,发现序列上间距较远的残基接触对的形成将使蛋白质折叠进程变慢．     

Conformational changes necessary for gene regulation by Tet repressor assayed by reversible disulfide bond formation.

We constructed and characterized four Tet repressor (TetR) variants with engineered cysteine residues which can form disulfide bonds and are located in regions where conformational changes during induction by tetracycline (tc) might occur. All TetR mutants show nearly wild-type activities in vivo, and the reduced proteins also show wild-type activities in vitro. Complete and reversible disulfide bond formation was achieved in vitro for all four mutants. The disulfide bond in NC18RC94 immobilizes the DNA reading head with respect to the protein core and prevents operator binding. Formation of this disulfide bond is possible only in the tc-bound, but not in the operator-bound conformation. Thus, these residues must have different conformations when bound to these ligands. The disulfide bonds in DC106PC159' and EC107NC165' immobilize the variable loop between alpha-helices 8 and 9 located near the tc-binding pocket. A faster rate of disulfide formation in the operator-bound conformation and a lack of induction after disulfide formation show that the variable loop is located closer to the protein core in the operator-bound conformation and that a movement is necessary for induction. The disulfide bond in RC195VC199' connects alpha-helices 10 and 10' of the two subunits in the dimer and is only formed in the tc-bound conformation. The oxidized protein shows reduced operator binding. Thus, this bond prevents formation of the operator-bound conformation. The detection of conformational changes in three different regions is the first biochemical evidence for induction-associated global internal movements in TetR.     

Structure-expression relationship of tumor necrosis factor receptor mutants that increase expression

The extracellular domain of the p55 TNF receptor (TNFrED) is an important therapeutic protein for targeting tumor necrosis factor-alpha (TNF-alpha). The expression level of the TNFrED is low for bioproduction, which is presumably associated with the complication of pairing 24 cysteine residues to form correct disulfide bonds. Here we report the application of the yeast display method to study expression of TNFrED, a multimeric receptor. Randomly mutated libraries of TNFrED were screened, and two mutants were identified that express several-fold higher protein levels compared with the wild type while still retaining normal binding affinity for TNF-alpha. The substituted residues responsible for the higher protein expression in both mutants were identified as proline, and both proline residues are adjacent to cysteine residues involved in disulfide bonds. Analysis of the mutant residues revealed that the improved level of expression is due to conformational restriction of the substituted residues to that of the folded state seen in the crystal structures of TNFrED thereby forcing the neighboring cysteine residues into the correct orientation for proper disulfide bond formation.     

Disulfide bond structure of the atrial natriuretic peptide receptor extracellular domain: conserved disulfide bonds among guanylate cyclase-coupled receptors

The disulfide bond structure of the extracellular domain of rat atrial natriuretic peptide (ANP) receptor (NPR-ECD) has been determined by mass spectrometry (MS) and Edman sequencing. Recombinant NPR-ECD expressed in COS-1 cells and purified from the culture medium binds ANP with as high affinity as the natural ANP receptor. Reaction with iodoacetic acid yielded no S-carboxymethylcysteine, indicating that all six Cys residues in NPR-ECD are involved in disulfide bonds. Electrospray ionization MS of NPR-ECD deglycosylated by peptide-N-glycosidase F gave a molecular mass of 48377.5+/-1.6 Da, which was consistent with the presence of three disulfide bonds. Liquid chromatography MS analysis of a lysylendopeptidase digest yielded three cystine-containing fragments with disulfide bonds Cys(60)-Cys(86), Cys(164)-Cys(213) and Cys(423)-Cys(432) based on their observed masses. These bonds were confirmed by Edman sequencing of each of the three fragments. No evidence for an inter-molecular disulfide bond was found. The six Cys residues in NPR-ECD, forming a 1-2, 3-4, 5-6 disulfide pairing pattern, are strictly conserved among A-type natriuretic peptide receptors and are similar in B-type receptors. We found that in other families of guanylate cyclase-coupled receptors, the Cys residues involved in 1-2 and 5-6 disulfide pairs are conserved in nearly all, suggesting an important contribution of these disulfide bonds to the receptor's structure and function.     

Role of individual disulfide bonds in the structural maturation of a low molecular weight glutenin subunit.

Gliadins and glutenins are the major storage proteins that accumulate in wheat endosperm cells during seed development. Although gliadins are mainly monomeric, glutenins consist of very large disulfide-linked polymers made up of high molecular weight and low molecular weight subunits. These polymers are among the largest protein molecules known in nature and are the most important determinants of the viscoelastic properties of gluten. As a first step toward the elucidation of the folding and assembly pathways that lead to glutenin polymer formation, we have exploited an in vitro system composed of wheat germ extract and bean microsomes to examine the role of disulfide bonds in the structural maturation of a low molecular weight glutenin subunit. When conditions allowing the formation of disulfide bonds were established, the in vitro synthesized low molecular weight glutenin subunit was recovered in monomeric form containing intrachain disulfide bonds. Conversely, synthesis under conditions that did not favor the formation of disulfide bonds led to the production of large aggregates from which the polypeptides could not be rescued by the post-translational generation of a more oxidizing environment. These results indicate that disulfide bond formation is essential for the conformational maturation of the low molecular weight glutenin subunit and suggest that early folding steps may play an important role in this process, allowing the timely pairing of critical cysteine residues. To determine which cysteines were important to maintain the protein in monomeric form, we prepared a set of mutants containing selected cysteine to serine substitutions. Our results show that two conserved cysteine residues form a critical disulfide bond that is essential in preventing the exposure of adhesive domains and the consequent formation of aberrant aggregates.     

Two-dimensional 1H NMR of two chemically modified analogs of the basic pancreatic trypsin inhibitor. Sequence-specific resonance assignments and sequence location of conformation changes relative to the native protein

Two-dimensional nuclear magnetic resonance was used to obtain sequence specific assignments for the 1H NMR spectra of two chemically modified analogs of the basic pancreatic trypsin inhibitor. In one analog the disulfide bond 14-38 was cleaved, in the second derivative the N-terminus was transaminated. From measurements of the chemical shifts and determination of the sequence locations of slowly exchanging backbone amide protons it was found that conformational differences between the native inhibitor and the chemical modifications occur exclusively near the modification sites and that the internal hydrogen bonds are nearly fully preserved. Intriguing conformation differences with respect to the native protein are that for five residues in the transaminated inhibitor and for one residue in the reduced inhibitor multiple local conformers are indicated, and that the four internal water molecules observed in the crystal structure of the native inhibitor appear not to be preserved after reduction of the disulfide bond 14-38.     

Identification of Allosteric Disulfides from Prestress Analysis

Disulfide bonds serve to form physical cross-links between residues in protein structures, thereby stabilizing the protein fold. Apart from this purely structural role, they can also be chemically active, participating in redox reactions, and they may even potentially act as allosteric switches controlling protein functions. Specific types of disulfide bonds have been identified in static protein structures from their distinctive pattern of dihedral bond angles, and the allosteric function of such bonds is purported to be related to the torsional strain they store. Using all-atom molecular-dynamics simulations for ∼700 disulfide bonded proteins, we analyzed the intramolecular mechanical forces in 20 classes of disulfide bonds. We found that two particular classes, the −RHStaple and the −/+RHHook disulfides, are indeed more stressed than other disulfide bonds, but the stress is carried primarily by stretching of the S-S bond and bending of the neighboring bond angles, rather than by dihedral torsion. This stress corresponds to a tension force of magnitude ∼200 pN, which is balanced by repulsive van der Waals interactions between the cysteine Cα atoms. We confirm stretching of the S-S bond to be a general feature of the −RHStaples and the −/+RHHooks by analyzing ∼20,000 static protein structures. Given that forced stretching of S-S bonds is known to accelerate their cleavage, we propose that prestress of allosteric disulfide bonds has the potential to alter the reactivity of a disulfide, thereby allowing us to readily switch between functional states.     

Amino acid sequence of human tumor derived angiogenin

The amino acid sequence and disulfide bond pairing of human tumor derived angiogenin, the first tumor angiogenesis factor to be isolated in pure form from human sources, have been determined by conventional sequencing techniques adapted and applied to nanomole and subnanomole levels of material. Angiogenin, obtained from conditioned media of a human colonic adenocarcinoma cell line, is a single-chain protein consisting of 123 amino acids with the following sequences: less than Glu1-Asp-Asn-Ser-Arg-Tyr-Thr-His- Phe-Leu-Thr-Gln-His-Tyr-Asp15-Ala-Lys-Pro-Gln-Gly-Arg-Asp-Asp- Arg-Tyr-Cys-Glu-Ser-Ile-Met30- Arg-Arg-Arg-Gly-Leu-Thr-Ser-Pro-Cys-Lys-Asp-Ile-Asn-Thr- Phe45-Ile-His-Gly-Asn-Lys-Arg-Ser -Ile-Lys-Ala-Ile-Cys-Glu-Asn-Lys60-Asn-Gly-Asn-Pro-His-Arg-Glu-Asn -Leu-Arg-Ile -Ser-Lys-Ser-Ser75 -Phe-Gln-Val-Thr-Thr-Cys-Lys-Leu-His-Gly-Gly-Ser-Pro-Trp-Pro90-Pro -Cys-Gln-Tyr -Arg-Ala-Thr-Ala -Gly-Phe-Arg-Asn-Val-Val-Val105-Ala-Cys-Glu-Asn-Gly-Leu-Pro-Val- His-Leu-Asp-Gln-Ser-Ile-Phe120-Arg-Arg-Pro123-OH. Three disulfide bonds link the half-cystinyl residues 26-81, 39-92, and 57-107. The sequence is homologous to that of the pancreatic ribonucleases with 35% identity and many of the remaining residues conservatively replaced. Similarities are especially apparent around the major active-site residues His-12, Lys-41, and His-119 of ribonuclease which are conserved as are three of the four disulfide bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

Evidence for intramolecular disulfide bond shuffling in the folding of mutant human lysozyme

Our previous results using the Saccharomyces cerevisiae secretion system suggest that intramolecular exchange of disulfide bonds occurs in the folding pathway of human lysozyme in vivo (Taniyama, Y., Yamamoto, Y., Kuroki, R., and Kikuchi, M. (1990) J. Biol. Chem. 265, 7570-7575). Here we report on the results of introducing an artificial disulfide bond in mutants with 2 cysteine residues substituting for Ala83 and Asp91. The mutant (C83/91) protein was not detected in the culture medium of the yeast, probably because of incorrect folding. Thereupon, 2 cysteine residues Cys77 and Cys95 were replaced with Ala in the mutant C83/91, because a native disulfide bond Cys77-Cys95 was found not necessary for correct folding in vivo (Taniyama, Y., Yamamoto, Y., Nakao, M., Kikuchi, M., and Ikehara, M. (1988) Biochem. Biophys. Res. Commun. 152, 962-967). The resultant mutant (AC83/91) was secreted as two proteins (AC83/91-a and AC83/91-b) with different specific activities. Amino acid and peptide mapping analyses showed that two glutathiones appeared to be attached to the thiol groups of the cysteine residues introduced into AC83/91-a and that four disulfide bonds including an artificial disulfide bond existed in the AC83/91-b molecule. The presence of cysteine residues modified with glutathione may indicate that the non-native disulfide bond Cys83-Cys91 is not so easily formed as a native disulfide bond. These results suggest that the introduction of Cys83 and Cys91 may act to suppress the process of native disulfide bond formation through disulfide bond interchange in the folding of human lysozyme.     

Restoration of structural stability and ligand binding after removal of the conserved disulfide bond in tear lipocalin

Disulfide bonds play diverse structural and functional roles in proteins. In tear lipocalin (TL), the conserved sole disulfide bond regulates stability and ligand binding. Probing protein structure often involves thiol selective labeling for which removal of the disulfide bonds may be necessary. Loss of the disulfide bond may destabilize the protein so strategies to retain the native state are needed. Several approaches were tested to regain the native conformational state in the disulfide-less protein. These included the addition of trimethylamine N-oxide (TMAO) and the substitution of the Cys residues of disulfide bond with residues that can either form a potential salt bridge or others that can create a hydrophobic interaction. TMAO stabilized the protein relaxed by removal of the disulfide bond. In the disulfide-less mutants of TL, 1.0 M TMAO increased the free energy change (ΔG⁰) significantly from 2.1 to 3.8 kcal/mol. Moderate recovery was observed for the ligand binding tested with NBD-cholesterol. Because the disulfide bond of TL is solvent exposed, the substitution of the disulfide bond with a potential salt bridge or hydrophobic interaction did not stabilize the protein. This approach should work for buried disulfide bonds. However, for proteins with solvent exposed disulfide bonds, the use of TMAO may be an excellent strategy to restore the native conformational states in disulfide-less analogs of the proteins.     

The disulfide bond pattern of catrocollastatin C, a disintegrin-like/cysteine-rich protein isolated from Crotalus atrox venom

The disulfide bond pattern of catrocollastatin-C was determined by N-terminal sequencing and mass spectrometry. The N-terminal disintegrin-like domain is a compact structure including eight disulfide bonds, seven of them in the same pattern as the disintegrin bitistatin. The protein has two extra cysteine residues (XIII and XVI) that form an additional disulfide bond that is characteristically found in the disintegrin-like domains of cellular metalloproteinases (ADAMs) and PIII snake venom Zn-metalloproteinases (SVMPs). The C-terminal cysteine-rich domain of catrocollastatin-C contains five disulfide bonds between nearest-neighbor cysteines and a long range disulfide bridge between CysV and CysX. These results provide structural evidence for a redefinition of the disintegrin-like and cysteine-rich domain boundaries. An evolutionary pathway for ADAMs, PIII, and PII SVMPs based on disulfide bond engineering is also proposed.