A new method of research on molecular evolution of proteinase superfamily

The molecular evolutionary tree, also known as a phylogenetic tree, of the serine proteinase superfamily was constructed by means of structural alignment. Three-dimensional structures of proteins were aligned by the SSAP program of Orengo and Taylor to obtain evolutionary distances. The resulting evolutionary tree provides a topology graph that can reflect the evolution of structure and function of homology proteinase. Moreover, study on evolution of the serine proteinase superfamily can lead to better understanding of the relationship and evolutionary difference among proteins of the superfamily, and is of significance to protein engineering, molecular design and protein structure prediction. Structure alignment is one of the useful methods of research on molecular evolution of protein.     

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.     

基于混沌游戏图形表示的细菌必需基因与非必需基因的Hurst指数特征对比分析

必需基因是维持生物体生命活动必不可少的基因,基因特征的选取对于必需基因的理论预测识别方法至关重要.作者采用分形理论对DEG必需基因数据库中的所有细菌物种进行了基于Hurst指数的必需基因与非必需基因特征对比分析,旨在提取必需基因区别于非必需基因的显著特征.结果表明,绝大部分必需基因与非必需基因序列具有长程相关性,所分析细菌物种的必需基因和非必需基因的Hurst指数均较好地遵从Gamma分布规律,且两类基因的Hurst指数具有显著性差异.     

The proteins of linked genes evolve at similar rates

Much more variation in the rate of protein evolution occurs than is expected by chance. But why some proteins evolve rapidly but others slowly is poorly resolved. It was proposed, for example, that essential genes might evolve slower than dispensable ones, but this is not the case; and despite earlier claims, rates of evolution do not correlate with amino-acid composition. A few patterns have been found: proteins involved in antagonistic co-evolution (for example, immune genes, parasite antigens and reproductive conflict genes) tend to be rapidly evolving, and there is a correlation between the rate of protein evolution and the mutation rate of the gene. Here we report a new highly statistically significant predictor of a protein's rate of evolution, and show that linked genes have similar rates of protein evolution. There is also a weaker similarity of rates of silent site evolution (see ref. 13), which appears to be, in part, a consequence of the similarity in rates of protein evolution. The similarity in rates of protein evolution is not a consequence of underlying mutational patterns. A pronounced negative correlation between the rate of protein evolution and a covariant of the recombination rate indicates that rates of protein evolution possibly reflect, in part, the local strength of stabilizing selection.     

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.     

A new way of enhancing the thermostability of proteases

Thermostability of enzyme can be enhanced by single amino acid substitutions. Recent advances in genetic engineering have made it possible to create novel proteins in a predictable manner where structural information for the protein is available. This 'protein engineering' has already been used to enhance enzyme thermostability, but it is usually not clear which amino acid substitutions should be made. We consider that the following approach should be helpful in engineering proteins with enhanced thermostability: highly conserved residues should be left unchanged; the sequences of known mesophilic and thermophilic proteins should be used to suggest the kinds of changes likely to increase thermostability; and substitutions should be made that increase internal hydrophobicity and that stabilize helices for strong internal packing. We describe here the use of this approach to alter the thermostability of the thermostable neutral protease of Bacillus stearothermophilus, the sequence of which is known. Surprisingly we find that a single mutation that decreases thermostability can require two mutations that increase stability to compensate for it. The effects on stability are not additive, suggesting cooperativity.     

同源建模方法预测蛋白质突变结构的适用性分析

通过从Protein Data Bank(PDB)结构数据库中提取单氨基酸突变的晶体结构,构建了一组无冗余的测试数据集,对目前应用最广泛的两款同源建模预测软件(SWISS-MODEL和MODELLER)进行了测试分析,发现它们对蛋白质的整体结构预测效果良好,均方根偏差小于0.5埃(RMSD0.5),但在突变导致结构显著变化(RMSD1.5)的情况下却均不能得到准确结果.分类统计显示,发生在蛋白质结构内部和极性氨基酸之间的突变结构变化小,两款软件预测效果较好(RMSD1.0).突变导致结构显著变化的可能性不高(5%),但它对蛋白质功能的影响不可忽视,因此应用同源建模方法对于蛋白质突变的模拟并不完全适用,还需要开发新方法来提高准确性.     

The complete DNA sequence of yeast chromosome III.

The entire DNA sequence of chromosome III of the yeast Saccharomyces cerevisiae has been determined. This is the first complete sequence analysis of an entire chromosome from any organism. The 315-kilobase sequence reveals 182 open reading frames for proteins longer than 100 amino acids, of which 37 correspond to known genes and 29 more show some similarity to sequences in databases. Of 55 new open reading frames analysed by gene disruption, three are essential genes; of 42 non-essential genes that were tested, 14 show some discernible effect on phenotype and the remaining 28 have no overt function.     

New yeast actin-like gene required late in the cell cycle.

Actin, a major cytoskeletal component of all eukaryotic cells, is one of the most highly conserved proteins. It is involved in various cellular processes such as motility, cytoplasmic streaming, chromosome segregation and cytokinesis. The actin from the yeast Saccharomyces cerevisiae, encoded by the essential ACT1 gene, is 89% identical to mouse cytoplasmic actin and is involved in the organization and polarized growth of the cell surface. We report here the characterization of ACT2, a previously undescribed yeast split gene encoding a putative protein (391 amino acids, relative molecular mass (Mr) 44,073) that is 47% identical to yeast actin. The requirement of the ACT2 gene for vegetative growth of yeast cells and the existence of related genes in other eukaryotes indicate an important and conserved role for these actin-like proteins. Superimposition of the Act2 polypeptide onto the three-dimensional structure of known actins reveals that most of the divergence occurred in loops involved in actin polymerization, DNase I and myosin binding, leaving the core domain mainly unaffected. To our knowledge, the Act2 protein from S. cerevisiae is the first highly divergent actin molecule described. Structural and physiological data suggest that the Act2 protein might have an important role in cytoskeletal reorganization during the cell cycle.     

Metabolic network analysis of the causes and evolution of enzyme dispensability in yeast

Under laboratory conditions 80% of yeast genes seem not to be essential for viability. This raises the question of what the mechanistic basis for dispensability is, and whether it is the result of selection for buffering or an incidental side product. Here we analyse these issues using an in silico flux model of the yeast metabolic network. The model correctly predicts the knockout fitness effects in 88% of the genes studied and in vivo fluxes. Dispensable genes might be important, but under conditions not yet examined in the laboratory. Our model indicates that this is the dominant explanation for apparent dispensability, accounting for 37-68% of dispensable genes, whereas 15-28% of them are compensated by a duplicate, and only 4-17% are buffered by metabolic network flux reorganization. For over one-half of those not important under nutrient-rich conditions, we can predict conditions when they will be important. As expected, such condition-specific genes have a more restricted phylogenetic distribution. Gene duplicates catalysing the same reaction are not more common for indispensable reactions, suggesting that the reason for their retention is not to provide compensation. Instead their presence is better explained by selection for high enzymatic flux.     

A yeast prion provides a mechanism for genetic variation and phenotypic diversity

A major enigma in evolutionary biology is that new forms or functions often require the concerted effects of several independent genetic changes. It is unclear how such changes might accumulate when they are likely to be deleterious individually and be lost by selective pressure. The Saccharomyces cerevisiae prion [PSI+] is an epigenetic modifier of the fidelity of translation termination, but its impact on yeast biology has been unclear. Here we show that [PSI+] provides the means to uncover hidden genetic variation and produce new heritable phenotypes. Moreover, in each of the seven genetic backgrounds tested, the constellation of phenotypes produced was unique. We propose that the epigenetic and metastable nature of [PSI+] inheritance allows yeast cells to exploit pre-existing genetic variation to thrive in fluctuating environments. Further, the capacity of [PSI+] to convert previously neutral genetic variation to a non-neutral state may facilitate the evolution of new traits.     

Direct estimation of per nucleotide and genomic deleterious mutation rates in Drosophila

Spontaneous mutations are the source of genetic variation required for evolutionary change, and are therefore important for many aspects of evolutionary biology. For example, the divergence between taxa at neutrally evolving sites in the genome is proportional to the per nucleotide mutation rate, u (ref. 1), and this can be used to date speciation events by assuming a molecular clock. The overall rate of occurrence of deleterious mutations in the genome each generation (U) appears in theories of nucleotide divergence and polymorphism, the evolution of sex and recombination, and the evolutionary consequences of inbreeding. However, estimates of U based on changes in allozymes or DNA sequences and fitness traits are discordant. Here we directly estimate u in Drosophila melanogaster by scanning 20 million bases of DNA from three sets of mutation accumulation lines by using denaturing high-performance liquid chromatography. From 37 mutation events that we detected, we obtained a mean estimate for u of 8.4 x 10(-9) per generation. Moreover, we detected significant heterogeneity in u among the three mutation-accumulation-line genotypes. By multiplying u by an estimate of the fraction of mutations that are deleterious in natural populations of Drosophila, we estimate that U is 1.2 per diploid genome. This high rate suggests that selection against deleterious mutations may have a key role in explaining patterns of genetic variation in the genome, and help to maintain recombination and sexual reproduction.     

Interference among deleterious mutations favours sex and recombination in finite populations

Sex and recombination are widespread, but explaining these phenomena has been one of the most difficult problems in evolutionary biology. Recombination is advantageous when different individuals in a population carry different advantageous alleles. By bringing together advantageous alleles onto the same chromosome, recombination speeds up the process of adaptation and opposes the fixation of harmful mutations by means of Muller's ratchet. Nevertheless, adaptive substitutions favour sex and recombination only if the rate of adaptive mutation is high, and Muller's ratchet operates only in small or asexual populations. Here, by tracking the fate of modifier alleles that alter the frequency of sex and recombination, we show that background selection against deleterious mutant alleles provides a stochastic advantage to sex and recombination that increases with population size. The advantage arises because, with low levels of recombination, selection at other loci severely reduces the effective population size and genetic variance in fitness at a focal locus (the Hill-Robertson effect), making a population less able to respond to selection and to rid itself of deleterious mutations. Sex and recombination reveal the hidden genetic variance in fitness by combining chromosomes of intermediate fitness to create chromosomes that are relatively free of (or are loaded with) deleterious mutations. This increase in genetic variance within finite populations improves the response to selection and generates a substantial advantage to sex and recombination that is fairly insensitive to the form of epistatic interactions between deleterious alleles. The mechanism supported by our results offers a robust and broadly applicable explanation for the evolutionary advantage of recombination and can explain the spread of costly sex.     

Experimental niche evolution alters the strength of the diversity–productivity relationship

The relationship between biodiversity and ecosystem functioning (BEF) has become a cornerstone of community and ecosystem ecology and an essential criterion for making decisions in conservation biology and policy planning. It has recently been proposed that evolutionary history should influence the BEF relationship because it determines species traits and, thus, species’ ability to exploit resources. Here we test this hypothesis by combining experimental evolution with a BEF experiment. We isolated 20 bacterial strains from a marine environment and evolved each to be generalists or specialists. We then tested the effect of evolutionary history on the strength of the BEF relationship with assemblages of 1 to 20 species constructed from the specialists, generalists and ancestors. Assemblages of generalists were more productive on average because of their superior ability to exploit the environmental heterogeneity. The slope of the BEF relationship was, however, stronger for the specialist assemblages because of enhanced niche complementarity. These results show how the BEF relationship depends critically on the legacy of past evolutionary events.     

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.     

Hsp104 is a highly conserved protein with two essential nucleotide-binding sites

Most eukaryotic cells produce proteins with relative molecular masses in the range of 100,000 to 110,000 after exposure to high temperatures. These proteins have been studied only in yeast and mammalian cells. In Saccharomyces cerevisiae, heat-shock protein hsp104 is vital for tolerance to heat, ethanol and other stresses. The mammalian hsp110 protein is nucleolar and redistributes with growth state, nutritional conditions and heat shock. The relationships between hsp110, hsp104 and the high molecular mass heat-shock proteins of other organisms were unknown. We report here that hsp104 is a member of the highly conserved ClpA/ClpB protein family first identified in Escherichia coli and that additional heat-inducible members of this family are present in Schizosaccharomyces pombe and in mammals. Mutagenesis of two putative nucleotide-binding sites in hsp104 indicates that both are essential for function in thermotolerance.     

A Reliable Neighbor-Based Method for Identifying Essential Proteins by Integrating Gene Expressions,Orthology,and Subcellular Localization Information

Essential proteins are those necessary for the survival or reproduction of species and discovering such essential proteins is fundamental for understanding the minimal requirements for cellular life, which is also meaningful to the disease study and drug design. With the development of high-throughput techniques, a large number of Protein-Protein Interactions(PPIs) can be used to identify essential proteins at the network level. Up to now, though a series of network-based computational methods have been proposed, it is still a challenge to improve the prediction precision as the high false positives in PPI networks. In this paper, we propose a new method GOS to identify essential proteins by integrating the Gene expressions, Orthology, and Subcellular localization information.The gene expressions and subcellular localization information are used to determine whether a neighbor in the PPI network is reliable. Only reliable neighbors are considered when we analyze the topological characteristics of a protein in a PPI network. We also analyze the orthologous attributes of each protein to reflect its conservative features, and use a random walk model to integrate a protein's topological characteristics and its orthology. The experimental results on the yeast PPI network show that the proposed method GOS outperforms the ten existing methods DC, BC, CC, SC, EC, IC, NC, Pe C, ION, and CSC.     

Haploidy or diploidy: which is better?

Although the evolutionary advantages of sexual reproduction have been extensively discussed, much less attention has been paid to haploid and diploid phases of the sexual life cycle. The relative lengths of these phases differ greatly in various taxa, including as extremes those with one or the other phase reduced to a single cell. Here we consider the efficiency of elimination of deleterious mutations as an evolutionary force and compare the mutation loads under haploid and diploid selection, Ln and L2n. With truncation-like selection, partial dominance, and heterozygous effect of a mutation less than about 1/4 its hemizygous effect, L2n less than Ln; otherwise L2n greater than Ln. The difference becomes important when the genomic deleterious mutation rate exceeds about 1 per genome. This suggests that the mutation rate, degree of dominance and mode of selection can be important in life-cycle evolution.