Significant Residue Features Revealed By Eigenvalue Decomposition Analysis Of Blosum Matrices

Here we attempt to characterize protein evolution by residue features which dominate residue substitution in homologous proteins. Evolutionary information contained in residue substitution matrix is abstracted with the method of eigenvalue decomposition. Top eigenvectors in the eigenvalue spectrums are analyzed as function of the level of similarity, i.e. sequence identity (SI) between homologous proteins. It is found that hydrophobicity and volume are two significant residue features conserved in protein evolution. There is a transition point at SI approximate to 45%. Residue hydrophobicity is a feature governing residue substitution as SI >= 45%. Whereas below this SI level, residue volume is a dominant feature. (C) 2007 Elsevier B.V. All rights reserved.

A scheme for multiple sequences alignment optimization-an improvement based on family representative mechanics features

As a basic tool of modern biology, sequence alignment can provide us useful information in fold, function, and active site of protein. For many cases, the increased quality of sequence alignment means a better performance. The motivation of present work is to increase ability of the existing scoring scheme/algorithm by considering residue–residue correlations better. Based on a coarse-grained approach, the hydrophobic force between each pair of residues is written out from protein sequence. It results in the construction of an intramolecular hydrophobic force network that describes the whole residue–residue interactions of each protein molecule, and characterizes protein's biological properties in the hydrophobic aspect. A former work has suggested that such network can characterize the top weighted feature regarding hydrophobicity. Moreover, for each homologous protein of a family, the corresponding network shares some common and representative family characters that eventually govern the conservation of biological properties during protein evolution. In present work, we score such family representative characters of a protein by the deviation of its intramolecular hydrophobic force network from that of background. Such score can assist the existing scoring schemes/algorithms, and boost up the ability of multiple sequences alignment, e.g. achieving a prominent increase (50%) in searching the structurally alike residue segments at a low identity level. As the theoretical basis is different, the present scheme can assist most existing algorithms, and improve their efficiency remarkably.

从基因组角度探讨果蝇蛋白质的进化

蛋白质进化的驱动力是分子进化中最基本的问题之一，研究者们早已在群体 遗传学研究中做出了许多模型和推测。基因组时代的到来使得证实这些经典理论 成为可能果蝇12 个基因组测序计划产生了大量新的数据，通过分析这些新的数 据，可以获得对果蝇蛋白质进化的更深入认识，也可以为分子进化研究的许多经 典理论提供验证。 我们利用黑腹果蝇种组的六个基因组，结合比较基因组学和功能基因组学的 研究手段，对黑腹果蝇种组中蛋白质的进化速率与许多可能影响蛋白质进化的功 能基因组学参数进行了统计学分析，从全基因组水平上探讨了果蝇蛋白质的进 化。 研究结果表明，表达特异性是决定果蝇蛋白质进化的最显著因素，广泛表达 的持家基因具有更低的进化速率，而表达量尽管与蛋白质进化呈负相关，但是这 个效应并不强，表达对蛋白质进化的影响与密码子偏好性和基因的功能多样性无 关，可能是由翻译后选择导致的。基因的功能重要性和基因的功能密度都对蛋白 质的进化有显著的影响，二者之间存在一定的联系，但是基因的功能密度对蛋白 质进化的影响要大于基因的功能重要性。蛋白质长度受到表达量选择，高表达的 基因往往具有更短的蛋白长度，因此，高表达蛋白的蛋白质长度与蛋白质进化速 率不相关，但是蛋白质长度在非高表达基因中对蛋白质进化存在显著的影响，其 机制可以部分的由基因的功能密度来解释，更确切的机制还需进一步研究。第一 个内含子的长度与蛋白质的进化速率显著相关，可能是由于第一个内含子中存在 较多的调控序列导致。基因重复的固定机制主要是高复杂度基因发生重复后的亚 功能化，而功能补偿效应和剂量平衡效应并不重要。基因的重组率在全基因组水 平上与蛋白质进化速率显著负相关，可能是因为基因组中大多数蛋白并不参与到 正选择的过程当中所致。 本研究系统的分析了影响果蝇蛋白质进化速率的各种功能基因组学参数的 效应，除了重新分析了过往研究中提到的基因表达水平，蛋白质相互作用，蛋白质长度这些参数在果蝇蛋白质进化中的影响，本工作还分析了许多新的参数，包 括表达特异度，内含子长度，基因功能重要性，基因功能密度，基因重复，基因 重组等。我们的工作不仅进一步证实了许多在进化基因组学研究中被广泛接受的 理论，还获得了一些与过去研究不同的结果，例如，果蝇蛋白质长度与蛋白质进 化速率呈负相关，第一个内含子在蛋白质进化中有很重要的作用等。

H5N1 avian influenza re-emergence of Lake Communication Qinghai: phylogenetic and antigenic analyses of the newly isolated viruses and roles of migratory birds in virus circulation

Highly pathogenic avian influenza H5N1 virus has swept west across the globe and caused serious debates on the roles of migratory birds in virus circulation since the first large-scale outbreak in migratory birds of Lake Qinghai, 2005. In May 2006, another outbreak struck Lake Qinghai and six novel strains were isolated. To elucidate these QH06 viruses, the six isolates were subjected to whole-genome sequencing. Phylogenetic analyses show that QH06 viruses are derived from the lineages of Lake Qinghai, 2005. Five of the six novel isolates are adjacent to the strain A/Cygnus olor/Croatia/1/05, and the last one is related to the strain A/duck/Novosibirsk/ 02/05, an isolate of the flyway. Antigenic analyses suggest that QH06 and QH05 viruses are similar to each other. These findings implicate that QH06 viruses of Lake Qinghai may travel back via migratory birds, though not ruling out the possibility of local circulation of viruses of Lake Qinghai.

Molecular cloning and comparative analysis of four beta-galactosidase genes from Bifidobacterium bifidum NCIMB41171

Bifidobacterium bifidum NCIMB41171 carries four genes encoding different beta-galactosidases. One of them, named bbgIII, consisted of an open reading frame of 1,935 amino acid (a.a.) residues encoding a protein with a multidomain structure, commonly identified on cell wall bound enzymes, having a signal peptide, a membrane anchor, FIVAR domains, immunoglobulin Ig-like and discoidin-like domains. The other three genes, termed bbgI, bbgII and bbgIV, encoded proteins of 1,291, 689 and 1,052 a.a. residues, respectively, which were most probably intracellularly located. Two cases of protein evolution between strains of the same species were identified when the a.a. sequences of the BbgI and BbgIII were compared with homologous proteins from B. bifidum DSM20215. The homologous proteins were found to be differentiated at the C-terminal a.a. part either due to a single nucleotide insertion or to a whole DNA sequence insertion, respectively. The bbgIV gene was located in a gene organisation surrounded by divergently transcribed genes putatively for sugar transport (galactoside-symporter) and gene regulation (LacI-transcriptional regulator), a structure that was found to be highly conserved in B. longum, B. adolescentis and B. infantis, suggesting optimal organisation shared amongst those species.

The mitochondrial view of Blastocladiella emersonii

The mitochondrial genome of the chytrid Blastocladiella emersonii was sequenced and annotated, revealing the complete set of oxidative phosphorylation genes and tRNAs/rRNAs necessary for the translation process. Phylogenetic reconstructions reinforce the proposal of the new phylum Blastocladiomycota. Evidences of gene duplication due to inserted elements suggest shared susceptibility to gene invasion/exchange between chytrids and zygomycetes. The gene content of B. emersonii is very similar to Allomyces macrogynus but the content of intronic and changeable elements is much lower suggesting a stronger resistance to this kind of exchange. in addition, a total of 401 potential nuclear transcripts encoding mitochondrial proteins were obtained after B. emersonii EST database scanning using Saccharomyces cerevisiae, Homo sapiens and Arabidopsis thaliana data as probes and TargetP tool to find N-terminal mitochondrial signal in translated sequences. (c) 2008 Elsevier B.V. All rights reserved.

Common structural basis for central nervous system drug design

The main theme of this thesis is that there is a common structural basis for drugs acting on the central nervous system (CNS), and that this concept may be used to design new CNS-active drugs which have greater specificity and hence less side-effects. To develop these ideas, the biological basis of how drugs modify CMS neurotransmission is described, and illustrated using dopaminergic pathways. An account is then given of the use of physicochemical concepts in contemporary drug design. The complete conformational analysis of several antipsychotic drugs is used to illustrate some of these techniques in the development of a model for antipsychotic drug action. After reviewing current structure-activity studies in several classes of CNS drugs (antipsychotics, anti-depressants, stimulants, hal1ucinogens, anticonvulsants and analgesics), a hypothesis for a common structural basis of CNS drug action is proposed- This is based on a topographical comparison of the X-ray structures of eight representative CNS-active drugs, and consists of three parts: 1.there is a common structural basis for the activity of many different CNS-active drug classes; 2. an aromatic ring and a nitrogen atom are the primary binding groups whose topographical arrangement is fundamental to the activity of these drug classes; 3. the nature and placement of secondary binding determines different classes of CNS drug activity. A four-Point model for this common structural basis is then defined using 14- CNS-active drug structures that include the original eight used in proposing the hypothesis. The coordinates of this model are: R1 (0. 3.5, 0), R2 (0, -3.5, O), N (4.8. -0.3, 1.4), and R3 (6.3, 1.3, 0), where R1 and R2 represent the point locations of a hydrophobic interaction of the common aromatic ring with a receptor, and R3 locates the receptor point for a hydrogen bond involving the common nitrogen, N. Extended structures were used to define the receptor points R1, R2 and R3, and the complete conformational space of each of the 14 molecules was considered. It is then shoun that the model may be used to predict whether a given structure is likely to show CNS activity: a search over 1,000 entries in the current Merck Index shows a high probability (82%) of CNS activity in compounds fitting the structural model. Analysis of CNS neurotransmitters and neuropeptides shows that these fit the common model well. Based on the available evidence supporting chemical evolution, protein evolution, and the evolution of neurotransmitter functions, it is surmised that the aromatic ring/nitrogen atom pharmacophore proposed in the common model supports the idea of the evolution of CNS receptors and their neurotransmitters, possibly from an aromatic amine or acety1cho1ine acting as a primaeval communicating molecule. The third point in the hypothesis trilogy is then addressed. The extensive conformation-activity analyses that have resulted in well-defined models for five separate CNS drug classes are used to map out the locations of secondary binding groups relative to the common model for anti-psychotics, antidepressants, analgesics, anticholinergics, and anticonvulsants. With this information, and knowledge derived from receptor-binding data, it is postulated that drugs having specified activity could be designed. In order to generate novel structures having a high probability of CNS-activity, a process of drug design is described in which known CNS structures are superimposed topographically using the common model as a template. Atoms regarded as superfluous may be selectively deleted and the required secondary binding groups added in predicted locations to give novel structures. It is concluded that this process provides the basis for the rational design of new lead compounds which could further be optimized for potent and specific CNS activity.

Análise da hidrofobicidade na evolução de proteínas

Pós-graduação em Biofísica Molecular - IBILCE

Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

A distinct tymovirus infecting Cassia hoffmannseggii in Brazil

Leaves of Cassia hoffmannseggii, a wild fabaceous species found in the Atlantic Forest, with a severe mosaic symptom were collected in Pernambuco State, Brazil. By transmission electron microscopy, two types of virus particles were found: the first was recognized as particles of a potyvirus, which was later identified as Cowpea aphid-borne mosaic virus; and the second was isometric and present in high concentration. The observation of vesicles at the periphery of chloroplasts suggested a tymovirus infection, which was confirmed by subsequent assays. A serological assay against several tymovirus antisera resulted in positive reaction of this tymo-like virus with an antiserum of Passion fruit yellow mosaic virus. By means of RT-PCR and using degenerated primers for the conserved region of RNA-dependent RNA polymerase (RdRp) gene of tymoviruses, a specific DNA fragment was amplified and sequenced. Based on this sequence, a specific forward primer was synthesized and successfully used to amplify the 3' terminal genome region, containing the partial RdRp gene and the complete coat protein (CP) sequences. The CP was 188 amino acids (aa) long, and the highest CP aa identity was observed with Kennedya yellow mosaic virus (61 %). Based on the current ICTV demarcation criterion, this isolate was considered as a distinct tymovirus and tentatively named as Cassia yellow mosaic-associated virus.

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

A natural polymorphism in β-lactamase is a global suppressor

A M182T substitution was discovered as a second-site suppressor of a missense mutation in TEM-1 β-lactamase. The combination of the M182T substitution with other substitutions in the enzyme indicates the M182T substitution is a global suppressor of missense mutations in β-lactamase. The M182T substitution also is found in natural variants of TEM-1 β-lactamase with altered substrate specificity that have evolved in response to antibiotic therapy. The M182T substitution may have been selected in natural isolates as a suppressor of folding or stability defects resulting from mutations associated with drug resistance. This pathway of protein evolution may occur in other targets of antimicrobial drugs such as the HIV protease.

Why are some proteins structures so common?

Many biological proteins are observed to fold into one of a limited number of structural motifs. By considering the requirements imposed on proteins by their need to fold rapidly, and the ease with which such requirements can be fulfilled as a function of the native structure, we can explain why certain structures are repeatedly observed among proteins with negligible sequence similarity. This work has implications for the understanding of protein sequence structure relationships as well as protein evolution.

Specificity of protein - Nucleic acid interaction and the biochemical evolution

The water soluble carbodiimide mediated condensation of dipeptides of the general form Gly-X was carried out in the presence of mono- and poly-nucleotides. The observed yield of the tetrapeptide was found to be higher for peptide-nucleotide system of higher interaction specificity following mainly the anticodon-amino acid relationship (Basu, H.S. & Podder, S.K., 1981, Ind. J. Biochem. Biophys.,19, 251-253). The yield of the condensation product of L-peptide was more because of its higher interaction specificity. The extent of the racemization during the condensation of Gly-L-Phe, Gly-L-Tyr and Gly-D-Phe was found to be dependent on the specificity of the interaction -the higher the specificity, the lesser the racemization. The product formed was shown to have a catalytic effect on the condensation reaction. These data thus provide a mechanism showing how the specific interaction between amino acids/dipeptides and nucleic acids could lead to the formation of the lsquoprimitiversquo translation machinery.