Assembly reflects evolution of protein complexes

A homomer is formed by self-interacting copies of a protein unit. This is functionally important, as in allostery, and structurally crucial because mis-assembly of homomers is implicated in disease. Homomers are widespread, with 50-70% of proteins with a known quaternary state assembling into such structures. Despite their prevalence, their role in the evolution of cellular machinery and the potential for their use in the design of new molecular machines, little is known about the mechanisms that drive formation of homomers at the level of evolution and assembly in the cell. Here we present an analysis of over 5,000 unique atomic structures and show that the quaternary structure of homomers is conserved in over 70% of protein pairs sharing as little as 30% sequence identity. Where quaternary structure is not conserved among the members of a protein family, a detailed investigation revealed well-defined evolutionary pathways by which proteins transit between different quaternary structure types. Furthermore, we show by perturbing subunit interfaces within complexes and by mass spectrometry analysis, that the (dis)assembly pathway mimics the evolutionary pathway. These data represent a molecular analogy to Haeckel's evolutionary paradigm of embryonic development, where an intermediate in the assembly of a complex represents a form that appeared in its own evolutionary history. Our model of self-assembly allows reliable prediction of evolution and assembly of a complex solely from its crystal structure.     

基于置换距离度量的蛋白质多序列比对算法性能评估

 蛋白质多序列比对是一种重要的生物信息学工具,在生物的进化分析以及蛋白质的结构预测方面有着重要的应用。各种比对算法在这个领域都取得了很大的成功,但是每种算法都有其固有的缺陷。提出置换距离法,对当前流行的几种蛋白质多序列比对算法进行对比评价。由于置换距离法仅关注于不同蛋白质间进化距离的相对次序,而不考虑这些进化距离之间的细微差异,因而得到的评价结论更具有鲁棒性。另外,采用最长公共子序法度量置换距离可以比较准确的反映不同置换之间的差异性。基于该算法,对Dialign, Tcoffee, ClustalW和Muscle多序列比对算法进行了性能评估。     

Plakalbumin, alpha 1-antitrypsin, antithrombin and the mechanism of inflammatory thrombosis

An old puzzle in protein biochemistry concerns the ready conversion of ovalbumin, by proteolysis, to the much more stable derivative, plakalbumin. Ovalbumin is now known to belong to the serpin superfamily, most of which are serine proteinase inhibitors. We report here studies of two such members of the family, the human plasma proteins alpha 1-antitrypsin and antithrombin, and show that they undergo a similar change in stability on selective proteolysis. This change, which is accompanied by a loss of inhibitory activity, can best be considered as an irreversible molecular transition from a native stressed (S) conformation, to a more ordered relaxed (R) form. The maintenance of the native S conformation, and hence the maintenance of inhibitory activity, is critically dependent on the integrity of an exposed loop of polypeptide. We propose that the susceptibility of this peptide loop to proteolytic cleavage gives it an incidental role as a physiological switch which allows the inactivation of individual inhibitors by specific proteolysis. The vulnerability of this exposed loop in each inhibitor also explains the pathological action of a number of venoms and toxins. In particular, the demonstration here of the cleavage of antithrombin, by leukocyte elastase, explains an observed change in blood coagulation that accompanies severe inflammation and which can result in fatal thrombosis.     

Ancestral lysozymes reconstructed, neutrality tested, and thermostability linked to hydrocarbon packing

The controversy surrounding the idea that neutral mutations dominate protein evolution is attributable in part to the inadequacy of the tools available to evolutionary investigators. With a few exceptions, most investigations into the force driving protein evolution have relied on indirect criteria for distinguishing neutral and non-neutral variants. To investigate a particular pathway of molecular evolution, we have reconstructed by site-directed mutagenesis likely ancestral variants of the lysozymes of modern game birds (order Galliformes), tested their activity and thermostability and determined their three-dimensional structure. We focused on amino acids at three positions that are occupied in all known game birds either by the triplet Thr 40, Ile 55, Ser 91, or by the triplet Ser 40, Val 55, Thr 91. We have synthesized proteins representing intermediates along the possible three-step evolutionary pathways between these triplets. Although all of these are active and stable, none of these intermediates is found in known lysozymes. A comparison of the structures and thermostabilities of the variants reveals a linear correlation between the side-chain volume of the triplet and the thermostability of the protein. Each pathway connecting the two extant triplet sequences includes a variant with a thermostability outside the range of the extant proteins. This observation is consistent with a non-neutral evolutionary pathway. The existence of variants that are more stable than the extant proteins suggests that selection for maximum thermostability may not have been an important factor in the evolution of this enzyme.     

高低危型HPV E7的分析以及58型E7结构与结合位点预测

 对18种HPV E7蛋白进行生物信息学分析,同时对HPV58 E7进行结构以及可结合位点的预测.从NCBI数据库中获取蛋白序列,通过ExPASy Protparam进行生物信息分析,Clustal W进行多序列比对,MEGA软件构建进化树,利用Zhang lab QUARK以及Q-SiteFinder能量依赖的检测对HPV58 E7进行结构预测与结合位点预测.得到了E7生物信息,多序列比对数据,构建成进化树,获得HPV58 E7蛋白结构,并且找到pRb与E7结合的三维位点.不同类型HPV E7蛋白存在进化差异,得到的HPV58型结构以及与pRb的可以结合位点,为抑制HPV58 E7蛋白功能的相关实验提供了理论依据.     

Hormone binding globulins undergo serpin conformational change in inflammation

A surprising recent finding is that thyroxine binding globulin (TBG) and cortisol binding globulin (CBG), are members of the serine protease inhibitor (serpin) superfamily. Apparently evolution has completely adapted the serpin structure for its new role in these proteins as a transport agent, as there is no evidence of any retained protease inhibitory activity. This drastic change in function raises the question as to why such a complex molecular framework has been selected for the relatively simple task of hormone transport? To function as inhibitors the serpins have a native stressed (S) conformation that makes them vulnerable to proteolytic cleavage, the cleavage being accompanied by an irreversible transition to a stable relaxed (R) form. We demonstrate here that TBG and CBG have retained the stressed native structure typical of the inhibitor members of the family and we provide evidence that the S-R transition has been adapted to allow altered hormone delivery at inflammatory sites.     

Familial dementia caused by polymerization of mutant neuroserpin.

Aberrant protein processing with tissue deposition is associated with many common neurodegenerative disorders; however, the complex interplay of genetic and environmental factors has made it difficult to decipher the sequence of events linking protein aggregation with clinical disease. Substantial progress has been made toward understanding the pathophysiology of prototypical conformational diseases and protein polymerization in the superfamily of serine proteinase inhibitors (serpins). Here we describe a new disease, familial encephalopathy with neuroserpin inclusion bodies, characterized clinically as an autosomal dominantly inherited dementia, histologically by unique neuronal inclusion bodies and biochemically by polymers of the neuron-specific serpin, neuroserpin. We report the cosegregation of point mutations in the neuroserpin gene (PI12) with the disease in two families. The significance of one mutation, S49P, is evident from its homology to a previously described serpin mutations, whereas that of the other, S52R, is predicted by modelling of the serpin template. Our findings provide a molecular mechanism for a familial dementia and imply that inhibitors of protein polymerization may be effective therapies for this disorder and perhaps for other more common neurodegenerative diseases.     

分子进化的几个普适性质全国理论生物物理研讨会(武汉,1999)上的报告

第２ 至第６ 节为全国理论生物物理研讨会（武汉,１９９９）上的报告,论述了分子进化中的５ 个普适性质：分子手性,遗传密码,核苷关联的短程为主性及其进化相关性,以最大信息原理概括突变和选择的进化机制,由二维结构序列决定的蛋白质框架结构．第一节为作者在开幕式上致词的一部分,讨论了理论生物学在生命科学中的地位与作用,特应听众要求,发表于此．     

Protein dispensability and rate of evolution.

If protein evolution is due in large part to slightly deleterious amino acid substitutions, then the rate of evolution should be greater in proteins that contribute less to individual fitness. The rationale for this prediction is that relatively dispensable proteins should be subject to weaker purifying selection, and should therefore accumulate mildly deleterious substitutions more rapidly. Although this argument was presented over twenty years ago, and is fundamental to many applications of evolutionary theory, the prediction has proved difficult to confirm. In fact, a recent study showed that essential mouse genes do not evolve more slowly than non-essential ones. Thus, although a variety of factors influencing the rate of protein evolution have been supported by extensive sequence analysis, the relationship between protein dispensability and evolutionary rate has remained unconfirmed. Here we use the results from a highly parallel growth assay of single gene deletions in yeast to assess protein dispensability, which we relate to evolutionary rate estimates that are based on comparisons of sequences drawn from twenty-one fully annotated genomes. Our analysis reveals a highly significant relationship between protein dispensability and evolutionary rate, and explains why this relationship is not detectable by categorical comparison of essential versus non-essential proteins. The relationship is highly conserved, so that protein dispensability in yeast is also predictive of evolutionary rate in a nematode worm.     

Bioinformatical study on the proteomics and evolution of SARS-CoV

A novel coronavirus has been identified as the causative agent of the severe acute respiratory syndrome (SARS). For all the SARS-CoV associated proteins derivated from the SARS-CoV genome, the physiochemical properties such as the molecular weight, isoelectric point and extinction coefficient of each protein were calculated. The transmembrane segments and subeellular localization (SubLoeation) prediction and conserved protein motifs search against database were employed to analyze the function of SARS-CoV proteins. Also, the homology protein sequence alignment and evolutionary distance matrix calculation between SARS-CoV associated proteins and the corresponding proteins of other coronaviruses were employed to identify the classification and phylogenetic relationship between SARS-CoV and other coronaviruses. The results showed that SARS-CoV is a novel coronavirus which is different from any of the three previously known groups of coronviruses, but it is closer to BoCoV and MHV than to other coronaviruses. This study is in aid of experimental determination of SARS-CoV proteomics and the development of antiviral vaccine.     

菠菜SoSWEET蛋白家族基因鉴定与表达分析

以现有菠菜基因组信息为基础,通过生物信息学方法筛选鉴定16个菠菜SoSWEET蛋白家族成员,命名为SoSWEET1~SoSWEET16.氨基酸残基数量在648~1 140之间,分子质量在54 070.32~95 868.64 u之间,理论等电点（pI）在5.06~5.19之间.亚细胞定位预测有6个SoSWEET蛋白定位于细胞膜,5个SoSWEET蛋白定位于内质网,5个SoSWEET蛋白定位于细胞膜、内质网.系统进化分析将菠菜SoSWEET蛋白家族分成4个亚族,在此基础上对基因结构、保守基序、顺式作用调控元件等进行分析.共鉴定了10个高度保守基序,其中所有菠菜SoSWEET蛋白都包含基序1,2和4,是构成菠菜SoSWEET蛋白中最高度保守的部分.所有菠菜SoSWEET蛋白家族成员都含MtN3_slv和PQ-loop superfamily结构域.大多数SoSWEET蛋白家族的基因含有5个内含子.顺式作用元件预测结果表明,菠菜SoSWEET基因启动子上包含光响应、生长发育、植物激素响应和逆境胁迫响应等顺式作用元件.组织表达分析表明,所有SoSWEET基因在根、茎、叶和叶柄中都有表达,霜霉病胁迫处理后16个基因表现出不同响应变化.本研究为后续深入研究菠菜SoSWEET蛋白家族成员的功能提供了重要参考.     

Structure of Sindbis virus core protein reveals a chymotrypsin-like serine proteinase and the organization of the virion

Sindbis virus consists of a nucleocapsid core surrounded by a lipid membrane through which penetrate 80 glycoprotein trimers. The structure of the core protein comprising the coat surrounding the genomic RNA has been determined. The polypeptide fold from residue 114 to residue 264 is homologous to that of chymotrypsin-like serine proteinases with catalytic residues His 141, Asp 163 and Ser 215 of the core protein positioned as in other serine proteinases. The C-terminal tryptophan remains in the P1 substrate site subsequent to the autocatalytic cis cleavage of the capsid protein, thus rendering the proteinase inactive. Model building of the Sindbis core protein dimer shows that the nucleocapsid is likely to have T = 4 quasisymmetry.     

水稻S5-ORF3蛋白质的进化分析

生殖隔离在形成新物种和物种同一性的保持中有重要作用,因而在水稻优化育种中有重要的意义。S5位点是一个已克隆的控制水稻生殖隔离的基因,它产生的S5-ORF3,ORF4和ORF5蛋白质共同调控着indica-japonica杂交后代的育性,特别是S5-ORF3蛋白质在打破水稻亚种indica-japonica之间的生殖隔离与促进物种间的基因交流中发挥重要作用。针对S5-ORF3的进化分析对于研究它的功能和起源很重要,但目前尚无报道。本文基于序列和进化分析,提出S5-ORF3为一种定位在内质网中的新的HSP70家族蛋白,但缺乏HSP70的C端的多肽结合结构域,推测S5-ORF3可能通过影响alpha-淀粉酶的合成来作用于内质网压力;找到了ORF3的19条同源蛋白,并用极大似然算法构建了可靠的进化树,发现水稻的S5-ORF3与节节麦的luminal-binding蛋白进化关系最为密切;作了置信度100%的S5-ORF3蛋白质的三维结构预测,预测了居中的配体绑定位点;发现了水稻的S5-ORF3与其他HSP70蛋白相比独特的基序,为进一步研究S5-ORF3的作用机理和演化历史提供了线索和数据支持。     

The phylogenetic history of immunodeficiency viruses

Knowledge of the phylogenetic history of the human immunodeficiency viruses (HIV-1 and HIV-2) is important for our understanding of the epidemiology of AIDS, the disease caused by these viruses. Reconstruction of the evolutionary tree is hampered, however, by two problems. One is the high variation in nucleotide sequence between the known HIV isolates which can create formidable difficulties in identifying homologous genomic sites that may be used in a molecular phylogenetic reconstruction. Another impediment has been the lack of unequivocal time calibration points: there is only a sparse 'fossil record' for HIV and limited historical epidemiological data. We have largely overcome these difficulties by: (1) a thorough optimal-sequence alignment analysis; (2) the inclusion of sequences of an early (1976) HIV-1 isolate, a recent (1986) HIV-2 isolate and two simian immunodeficiency viruses (SIV) along with five other HIV-1 isolates; and (3) the reconstruction of a minimum-length evolutionary tree based on the envelope-gene variable positions. We conclude that HIV-1 may have evolved from its common ancestor with HIV-2 as recently as 40 years ago.     

Modest stabilization by most hydrogen-bonded side-chain interactions in membrane proteins

Understanding the energetics of molecular interactions is fundamental to all of the central quests of structural biology including structure prediction and design, mapping evolutionary pathways, learning how mutations cause disease, drug design, and relating structure to function. Hydrogen-bonding is widely regarded as an important force in a membrane environment because of the low dielectric constant of membranes and a lack of competition from water. Indeed, polar residue substitutions are the most common disease-causing mutations in membrane proteins. Because of limited structural information and technical challenges, however, there have been few quantitative tests of hydrogen-bond strength in the context of large membrane proteins. Here we show, by using a double-mutant cycle analysis, that the average contribution of eight interhelical side-chain hydrogen-bonding interactions throughout bacteriorhodopsin is only 0.6 kcal mol(-1). In agreement with these experiments, we find that 4% of polar atoms in the non-polar core regions of membrane proteins have no hydrogen-bond partner and the lengths of buried hydrogen bonds in soluble proteins and membrane protein transmembrane regions are statistically identical. Our results indicate that most hydrogen-bond interactions in membrane proteins are only modestly stabilizing. Weak hydrogen-bonding should be reflected in considerations of membrane protein folding, dynamics, design, evolution and function.     

A bacterial dynamin-like protein

Dynamins form a superfamily of large mechano-chemical GTPases that includes the classical dynamins and dynamin-like proteins (DLPs). They are found throughout the Eukarya, functioning in core cellular processes such as endocytosis and organelle division. Many bacteria are predicted by sequence to possess large GTPases with the same multidomain architecture that is found in DLPs. Mechanistic dissection of dynamin family members has been impeded by a lack of high-resolution structural data currently restricted to the GTPase and pleckstrin homology domains, and the dynamin-related human guanylate-binding protein. Here we present the crystal structure of a cyanobacterial DLP in both nucleotide-free and GDP-associated conformation. The bacterial DLP shows dynamin-like qualities, such as helical self-assembly and tubulation of a lipid bilayer. In vivo, it localizes to the membrane in a manner reminiscent of FZL, a chloroplast-specific dynamin-related protein with which it shares sequence similarity. Our results provide structural and mechanistic insight that may be relevant across the dynamin superfamily. Concurrently, we show compelling similarity between a cyanobacterial and chloroplast DLP that, given the endosymbiotic ancestry of chloroplasts, questions the evolutionary origins of dynamins.     

Evolution of increased complexity in a molecular machine

Many cellular processes are carried out by molecular 'machines'-assemblies of multiple differentiated proteins that physically interact to execute biological functions. Despite much speculation, strong evidence of the mechanisms by which these assemblies evolved is lacking. Here we use ancestral gene resurrection and manipulative genetic experiments to determine how the complexity of an essential molecular machine--the hexameric transmembrane ring of the eukaryotic V-ATPase proton pump--increased hundreds of millions of years ago. We show that the ring of Fungi, which is composed of three paralogous proteins, evolved from a more ancient two-paralogue complex because of a gene duplication that was followed by loss in each daughter copy of specific interfaces by which it interacts with other ring proteins. These losses were complementary, so both copies became obligate components with restricted spatial roles in the complex. Reintroducing a single historical mutation from each paralogue lineage into the resurrected ancestral proteins is sufficient to recapitulate their asymmetric degeneration and trigger the requirement for the more elaborate three-component ring. Our experiments show that increased complexity in an essential molecular machine evolved because of simple, high-probability evolutionary processes, without the apparent evolution of novel functions. They point to a plausible mechanism for the evolution of complexity in other multi-paralogue protein complexes.     

Cloning and characterization of a new actin gene from Oryza sativa L.

Using Rho family member osRACD as bait, a new member of actin gene family -Act was isolated from Oryza sativa by yeast two-hybrid system. The full-length cDNA was cloned with 5' RACE technology, which contains an open reading frame of 1134 bp with a predicted protein of 377 amino acids. Sequence alignment revealed 96% to 81.8% identities with some known actin proteins in plants. The method of bioinformatics was used to analyze the protein modification sites, structure and evolution of the gene. Southern blot analysis showed that Act is a single-copy gene in the genome. The result of RT-PCR showed it is ubiquitously expressed in root, shoot, callus and panicle in a temporal fashion. The relationship between Rho family and actin family in evolution and function was also studied.     

Insulin-stimulated MAP-2 kinase phosphorylates and activates ribosomal protein S6 kinase II

Ribosomal protein S6 is a component of the eukaryotic 40S ribosomal subunit that becomes phosphorylated on multiple serine residues in response to a variety of mitogens, including insulin, growth factors, and transforming proteins of many oncogenic viruses. Recently, an activated S6 kinase (S6 K II) has been purified to homogeneity from Xenopus eggs, and characterized immunologically and at the molecular level. Purified S6 K II can be deactivated in vitro by incubation with either protein phosphatase 1 or protein phosphatase 2A. Reactivation and phosphorylation of S6 K II occurs in vitro with an insulin-stimulated microtubule-associated protein-2 (MAP-2) protein kinase which is itself a phosphoprotein that can be deactivated by protein phosphatase 2A. These studies suggest that a step in insulin signalling involves sequential activation by phosphorylation of at least two serine/threonine protein kinases.