Phylogenetic relationships and classification of thiolases and thiolase-like proteins of Mycobacterium tuberculosis and Mycobacterium smegmatis

Thiolases are enzymes involved in lipid metabolism. Thiolases remove the acetyl-CoA moiety from 3-ketoacyl-CoAs in the degradative reaction. They can also catalyze the reverse Claisen condensation reaction, which is the first step of biosynthetic processes such as the biosynthesis of sterols and ketone bodies. In human, six distinct thiolases have been identified. Each of these thiolases is different from the other with respect to sequence, oligomeric state, substrate specificity and subcellular localization. Four sequence fingerprints, identifying catalytic loops of thiolases, have been described. In this study genome searches of two mycobacterial species (Mycobacterium tuberculosis and Mycobacterium smegmatis), were carried out, using the six human thiolase sequences as queries. Eight and thirteen different thiolase sequences were identified in M. tuberculosis and M. smegmatis, respectively. In addition, thiolase-like proteins (one encoded in the Mtb and two in the Msm genome) were found. The purpose of this study is to classify these mostly uncharacterized thiolases and thiolase-like proteins. Several other sequences obtained by searches of genome databases of bacteria, mammals and the parasitic protist family of the Trypanosomatidae were included in the analysis. Thiolase-like proteins were also found in the trypanosomatid genomes, but not in those of mammals. In order to study the phylogenetic relationships at a high confidence level, additional thiolase sequences were included such that a total of 130 thiolases and thiolase-like protein sequences were used for the multiple sequence alignment. The resulting phylogenetic tree identifies 12 classes of sequences, each possessing a characteristic set of sequence fingerprints for the catalytic loops. From this analysis it is now possible to assign the mycobacterial thiolases to corresponding homologues in other kingdoms of life. The results of this bioinformatics analysis also show interesting differences between the distributions of M. tuberculosis and M. smegmatis thiolases over the 12 different classes. (C) 2014 Elsevier Ltd. All rights reserved.

An Amino Acid Substitution Matrix for Protein Conformation Identification

Amino acid substitution matrices play an essential role in protein sequence alignment, a fundamental task in bioinformatics. Most widely used matrices, such as PAM matrices derived from homologous sequences and BLOSUM matrices derived from aligned segments of PROSITE, did not integrate conformation information in their construction. There are a few structure-based matrices, which are derived from limited data of structure alignment. Using databases PDB_SELECT and DSSP, we create a database of sequence-conformation blocks which explicitly represent sequence-structure relationship. Members in a block are identical in conformation and are highly similar in sequence. From this block database, we derive a conformation-specific amino acid substitution matrix CBSM60. The matrix shows an improved performance in conformational segment search and homolog detection.

A Computational Module Assembled from Different Protease Family Motifs Identifies PI PLC from Bacillus cereus as a Putative Prolyl Peptidase with a Serine Protease Scaffold

Proteolytic enzymes have evolved several mechanisms to cleave peptide bonds. These distinct types have been systematically categorized in the MEROPS database. While a BLAST search on these proteases identifies homologous proteins, sequence alignment methods often fail to identify relationships arising from convergent evolution, exon shuffling, and modular reuse of catalytic units. We have previously established a computational method to detect functions in proteins based on the spatial and electrostatic properties of the catalytic residues (CLASP). CLASP identified a promiscuous serine protease scaffold in alkaline phosphatases (AP) and a scaffold recognizing a beta-lactam (imipenem) in a cold-active Vibrio AP. Subsequently, we defined a methodology to quantify promiscuous activities in a wide range of proteins. Here, we assemble a module which encapsulates the multifarious motifs used by protease families listed in the MEROPS database. Since APs and proteases are an integral component of outer membrane vesicles (OMV), we sought to query other OMV proteins, like phospholipase C (PLC), using this search module. Our analysis indicated that phosphoinositide-specific PLC from Bacillus cereus is a serine protease. This was validated by protease assays, mass spectrometry and by inhibition of the native phospholipase activity of PI-PLC by the well-known serine protease inhibitor AEBSF (IC50 = 0.018 mM). Edman degradation analysis linked the specificity of the protease activity to a proline in the amino terminal, suggesting that the PI-PLC is a prolyl peptidase. Thus, we propose a computational method of extending protein families based on the spatial and electrostatic congruence of active site residues.

Hardware paralelo reconfigurável para identificação de alinhamentos de sequências de DNA.

Amostras de DNA são encontradas em fragmentos, obtidos em vestígios de uma cena de crime, ou coletados de amostras de cabelo ou sangue, para testes genéticos ou de paternidade. Para identificar se esse fragmento pertence ou não a uma sequência de DNA, é necessário compará-los com uma sequência determinada, que pode estar armazenada em um banco de dados para, por exemplo, apontar um suspeito. Para tal, é preciso uma ferramenta eficiente para realizar o alinhamento da sequência de DNA encontrada com a armazenada no banco de dados. O alinhamento de sequências de DNA, em inglês DNA matching, é o campo da bioinformática que tenta entender a relação entre as sequências genéticas e suas relações funcionais e parentais. Essa tarefa é frequentemente realizada através de softwares que varrem clusters de base de dados, demandando alto poder computacional, o que encarece o custo de um projeto de alinhamento de sequências de DNA. Esta dissertação apresenta uma arquitetura de hardware paralela, para o algoritmo BLAST, que permite o alinhamento de um par de sequências de DNA. O algoritmo BLAST é um método heurístico e atualmente é o mais rápido. A estratégia do BLAST é dividir as sequências originais em subsequências menores de tamanho w. Após realizar as comparações nessas pequenas subsequências, as etapas do BLAST analisam apenas as subsequências que forem idênticas. Com isso, o algoritmo diminui o número de testes e combinações necessárias para realizar o alinhamento. Para cada sequência idêntica há três etapas, a serem realizadas pelo algoritmo: semeadura, extensão e avaliação. A solução proposta se inspira nas características do algoritmo para implementar um hardware totalmente paralelo e com pipeline entre as etapas básicas do BLAST. A arquitetura de hardware proposta foi implementada em FPGA e os resultados obtidos mostram a comparação entre área ocupada, número de ciclos e máxima frequência de operação permitida, em função dos parâmetros de alinhamento. O resultado é uma arquitetura de hardware em lógica reconfigurável, escalável, eficiente e de baixo custo, capaz de alinhar pares de sequências utilizando o algoritmo BLAST.

Purification and cloning of cysteine-rich proteins from Trimeresurus jerdonii and Naja atra venoms

Three 26 kDa proteins, named as TJ-CRVP, NA-CRVP1 and NA-CRVP2, were isolated from the venoms of Trimeresurus jerdonii and Naja atra, respectively. The N-terminal sequences of TJ-CRVP and NA-CRVPs were determined. These components were devoid of the enzymatic activities tested, such as phospholipase A(2), arginine esterase, proteolysis, L-amino acid oxidase, 5' nucleotidase, acetylcholinesterase. Furthermore, these three components did not have the following biological activities: coagulant and anticoagulant activities, lethal activity, myotoxicity, hemorrhagic activity, platelet aggregation and platelet aggregation-inhibiting activities. These proteins are named as cysteine-rich venom protein (CRVP) because their sequences showed high level of similarity with mammalian cysteine-rich secretory protein (CRISP) family. Recently, some CRISP-like proteins were also isolated from several different snake venoms, including Agkistrodon blomhoffi, Trimeresurus flavoviridis, Lanticauda semifascita and king cobra. We presumed that CRVP might be a common component in snake venoms. Of particular interest, phylogenetic analysis and sequence alignment showed that NA-CRVP1 and ophanin, both from elapid snakes, share higher similarity with CRVPs from Viperidae snakes. (C) 2003 Elsevier Ltd. All rights reserved.

A new protein structure of P-II class snake venom metalloproteinases: it comprises metalloproteinase and disintegrin domains

A new metalloproteinase-disintegrin, named Jerdonitin, was purified from Trimeresurus jerdonii venom with a molecular weight of 36 kDa on SDS-PAGE. It dose-dependently inhibited ADP-induced human platelet aggregation with IC50 of 120 nM. cDNA cloning and sequencing revealed that Jerdonitin belonged to the class II of snake venom metalloproteinases (SVMPs) (P-II class). Different from other P-II class SVMPs, metalloproteinase and disintegrin domains of its natural protein were not separated, confirmed by internal peptide sequencing. Compared to other P-II class SVMPs, Jerdonitin has two additional cysteines (Cys219 and Cys238) located in the spacer domain and disintegrin domain, respectively. They probably form a disulfide bond and therefore the metalloproteinase and disintegrin domains cannot be separated by posttranslationally processing. In summary, comparison of the amino acid sequences of Jerdonitin with those of other P-II class SVMPs by sequence alignment and phylogenetic analysis, in conjunction with natural protein structure data, suggested that it was a new type of P-II class SVMPs. (C) 2003 Elsevier Inc. All rights reserved.

Molecular characterization of a weak fibrinogen-clotting enzyme from Trimeresurus jerdonii venom

A fibrinogen-clotting enzyme designed as jerdonobin-II was isolated from the venom of Trimeresurus jerdonii. It differed in molecular weight and N-terminal sequence with the previously isolated jerdonobin, a thrombin-like enzyme from the same venom. The enzyme consists of a single polypeptide chain with molecular weights of 30,000 and 32,000 under non-reducing and reducing conditions, respectively. Jerdonobin-II showed weak fibrinogen clotting activity and its activity unit on fibrinogen was calculated to be less than one unit using human thrombin as standard. The precursor protein sequence of jerodonobin-II was deduced from cloned cDNA sequence. The sequence shows high similarity (identity = 89%) to TSV-PA, a specific plasminogen activator from venom of T stejnegeri. Despite of the sequence similarity, jerdonobin-II was found devoid of plasminogen activating effect. Sequence alignment analysis suggested that the replacement of Lys(239) in TSV-PA to Gln(239) in jerdonobin-II might play an important role on their plasminogen activating activity difference. (C) 2005 Elsevier Ltd. All rights reserved.

Identification of triosephosphate isomerase (TIM) genes from Microcystis (Cyanobacteria : Chroococcales) and theoretical analyses of the potential of cyanobacterial TIM as a target for designing specific inhibitors

The genes encoding triosephosphate isomerase (TIM) in three species of Microcystis (M. aeruginosa, M. viridis and M. wesenbergii) were investigated. Reverse transcriptase-polymerase chain reaction indicated that they were transcribed in the cells. Analyses showed that their DNA and deduced amino acid sequences were highly conserved between all the three species, only a single nonsynonymous substitution was seen at position 31, from an Asp in M. aeruginosa and M. viridis to Glu in M. wesenbergii. Sequence alignment of these with 12 other known cyanobacterial TIM sequences showed that all the cyanobacterial TIMs had a very high level of amino acid identity (over 50% between each two). Comparison of the cyanobacterial TIMs with other reported TIMs (from diverse lineages of the three Domains) showed that they possessed common active-site residues and sequence motifs. All cyanobacterial TIMs have two common cysteine residues (Cys127 and Cys176), and the Cys176 is almost cyanobacteria-specific with only one exception in Streptomyces coelicolor. Both secondary structure alignment and comparative modelling of Synechocystis sp. TIM showed that Cys176 was located at the hinge region of the flexible loop-6 and might therefore be critical to the movement of TIM's loop-6, which is important to the function of the enzyme. Thus, the cyanobacterial TIM-specific Cys176 may be a potential site for the discovery of suitable drugs against cyanobacteria, and such drugs may have utility in controlling water blooms due to cyanobacteria.

The three dimensional structure and biological domains of the p53 C-terminus.

It was expected that there are a coil (289 similar to 325) and two a helix (alpha(1)368 similar to 373, alpha(2)381 similar to 388) structures in p53 protein C-terminal region based on its mRNA secondary structure template and Chou-Fasman's protein secondary structure principle of prediction. The result was conformed by the other four methods of protein secondary structure prediction that are based on the multiple sequence alignment (accuracy = 73.20%). Combine with the 31 amino acids crystal structure of the oligomerization, the three dimensional conformation of p53 C-terminal 108 residues was built using the SGI INDIGO(2) computer. This structure further expounds the relationship among those biological function domains of p53 C- terminus at three-dimensional level.

Cloning, expression and subcellular distribution of a Rana grylio virus late gene encoding ERV1 homologue

An essential for respiration and viability (ERV1) homologue, 88R, was cloned and characterized from Rana grylio virus (RGV). Database searches found its homologues in all sequenced iridoviruses, and sequence alignment revealed a highly conserved motif shared by all ERV1 family proteins: Cys-X-X-Cys. RT-PCR and western blot analysis revealed that 88R begins to transcribe and translate at 6 h postinfection (p.i.) and remains detectable at 48 h p.i. during RGV infection course. Furthermore, using drug inhibition analysis by a de novo protein synthesis inhibitor and a viral DNA replication inhibitor, RGV 88R was classified as a late (L) viral gene during the in vitro infection. 88R-EGFP fusion protein was observed in both the cytoplasm and nucleus of pEGFP-N3-88R transfected EPC cells. Although result of immunofluorescence is similar, 88R protein was not detected in viromatrix. Moreover, function of RGV 88R on virus replication were evaluated by RNAi assay. Nevertheless, effect of knockdown of RGV 88R expression on virus replication was not detected in cultured fish cell lines. Collectively, current data indicate that RGV 88R was a late gene of iridovirus encoding protein that distributed both the cytoplasm and nucleus.

Phylogeography of the freshwater catfish Hemibagrus guttatus (Siluriformes, Bagridae): Implications for South China biogeography and influence of sea-level changes

National Natural Science Foundation of China (NSFC) ; [2007CB411600]; [30530120]

Identification and characterization of a novel envelope protein in Rana grylio virus

Viral envelope proteins have been proposed to play significant roles in virus infection and assembly. In this study, an envelope protein gene, 53R, was cloned and characterized from Rana grylio virus (RGV), a member of the family Iridoviridae. Database searches found its homologues in all sequenced iricloviruses, and sequence alignment revealed several conserved structural features shared by virus capsid or envelope proteins: a myristoylation site, two predicted transmembrane domains and two invariant cysteine residues. Subsequently, RT-PCR and Western blot detection revealed that the transcripts encoding RGV 53R and the protein itself appeared late during infection of fathead minnow cells and that their appearance was blocked by viral DNA replication inhibitor, indicating that RGV 53R is a late expression gene. Moreover, immunofluorescence localization found an association of 53R with virus factories in RGV-infected cells, and this association was further confirmed by expressing a 53R-GFP fusion protein in pEGFP-N3/53R-transfected cells. Furthermore, detergent extraction and Western blot detection confirmed that RGV 53R was associated with virion membrane. Therefore, the current data suggest that RGV 53R is a novel viral envelope protein and that it may play an important role in virus assembly. This is thought to be the first report on a viral envelope protein that is conserved in all sequenced iridoviruses.

Predicted glycosyl transferase genes located outside the HEP island are required for formation of heterocyst envelope polysaccharide in Anabaena sp strain PCC 7120

During maturation, heterocysts form an envelope layer of polysaccharide, called heterocyst envelope polysaccharide (HEP), whose synthesis depends on a cluster of genes, the HEP island, and on an additional, distant gene, hepB, or a gene immediately downstream from hepB. We show that HEP formation depends upon the predicted glycosyl transferase genes all4160 at a third locus and alr3699, which is adjacent to hepB and is cotranscribed with it. Mutations in the histidine kinase genes hepN and hepK appear to silence the promoter of hepB and incompletely down-regulate all4160.

Mutations stabilize small subunit ribosomal RNA in desiccation-tolerant cyanobacteria Nostoc

The ribosomal RNA molecule is an ideal model for evaluating the stability of a gene product under desiccation stress. We isolated 8 Nostoc strains that had the capacity to withstand desiccation in habitats and sequenced their 16S rRNA genes. The stabilities of 16S rRNAs secondary structures, indicated by free energy change of folding, were compared among Nostoc and other related species. The results suggested that 163 rRNA secondary structures of the desiccation-tolerant Nostoc strains were more stable than that of planktonic Nostocaceae species. The stabilizing mutations were divided into two categories: (1) those causing GC to replace other types of base pairs in stems and (2) those causing extension of stems. By mapping stabilizing mutations onto the Nostoc phylogenetic tree based on 16S rRNA gene, it was shown that most of stabilizing mutations had evolved during adaptive radiation among Nostoc spp. The evolution of 16S rRNA along the Nostoc lineage is suggested to be selectively advantageous under desiccation stress.

Isolation and characterization of Scophthalmus maximus rhabdovirus

A rhabdovirus associated with a lethal hemorrhagic disease in cultured turbot Scophthalm us maximus Linnaeus was isolated. The virus induced typical cytopathogenic effects (CPE) in 9 of 15 fish cell lines examined and was then propagated and isolated from infected carp leucocyte cells (CLC). Electron microscopy observations revealed that the negatively stained virions had a typical bullet-shaped morphology with one rounded end and one flat base end. The bullet-shaped morphology was more obvious and clear in ultrathin sections of infected cells. Experimental infections also indicated that the S. maximus rhabdovirus (SMRV) was not only a viral pathogen for cultured turbot, but also had the ability to infect other fish species, such as freshwater grass carp. A partial nucleotide sequence of the SMRV polymerase gene was determined by RT-PCR using 2 pairs of degenerate primers designed according to the conserved sequences of rhabdovirus polymerase genes. Homology analysis, amino acid sequence alignment, and phylogenetic relationship analysis of the partial SMRV polymerase sequence indicated that SMRV was genetically distinct from other rhabdoviruses. Sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE) analysis of the purified SMRV revealed 5 major structural proteins, and their molecular masses were estimated to be about 250, 58, 47, 42, and 28 kDa. Significant serological reactivity differences were also observed between SMRV and its nearest neighbor, spring viremia of carp virus (SVCV). The data suggest that SMRV is likely a novel fish rhabdovirus, although it is closely related to rhabdoviruses in the genus Vesiculovirus.