BLASTP Algorithm Research Articles

Are zinc fingers of α-proteobacteria an important molecular mechanism?

The zinc-finger proteins in α-proteobacteria are considered to be the precursors of zinc fingers in higher Eukaryotes, due to their structural similarities and the crucial role played by α-proteobacteria during eukaryogenesis. Since they interact with macromolecules, their emergence in Bacteria, and later in the Eukarya, provided those organisms with major functional upgrades and evolutionary advantages. Hence, the aim of this research is to prove the important role of α-proteobacterial zinc fingers in this taxon, and the significance of their relics, present in higher Eukaryotes. A part of this investigation involved identifying similar zinc-finger proteins between α-proteobacteria and humans using the blastp algorithm, finding their shared domains in the CDD database, and modelling their structures in the ChimeraX programme. The remaining research objectives were reached by analysing data from pre-existing studies, in order to come up with an additional set of conclusions, relevant to this investigation. It was established that the Ros/MucR zinc-finger proteins regulate many pathways, crucial for the survival of α-proteobacteria and their interactions with Eukaryotes. Additionally, it was found that many zinc-finger proteins supply α-proteobacteria with eukaryotic mechanisms, which differentiate them from other bacterial taxa. It was also concluded that α-proteobacterial zinc fingers may be responsible for the resistance of α-proteobacteria to certain heavy metals. This investigation also proposes a new evolutionary hypothesis for the emergence of zinc fingers in Proteobacteria, and presents further arguments in favour of the theory that the Eukaryota acquired zinc fingers from α-proteobacteria during eukaryogenesis.

Genome and Transcriptome Sequence Resources and Effector Repertoire of Pythium myriotylum Drechsler.

Genome and Transcriptome Sequence Resources and Effector Repertoire of Pythium myriotylum Drechsler.

Genome-wide exploration of sugar transporter (sweet) family proteins in Fabaceae for Sustainable protein and carbon source.

Sugar transporter proteins (STPs) are membrane proteins required for sugar transport throughout cellular membranes. They plays an imperative role in sugar transmission across the plant and determinants of crop yield. However, the analysis of these important STPs Sugars Will Eventually be Exported Transporters (SWEET) family in legumes is still not well-documented and remains unclear. Therefore, the in-silico analysis of STPs has been performed to unravel their cellular, molecular, and structural composition in legume species. This study conducted a systematic search for STPs in Cajanus cajan using the Blastp algorithm to understand its molecular basis. Here, we performed a comprehensive analysis of 155 identified SWEET proteins across 12 legumes species, namely (Cajanus cajan, Glycine max, Vigna radiate, Vigna angularis, Medicago truncatula, Lupinus angustifolius, Glycine soja, Spatholobus suberectus, Cicer arietinum, Arachis ipaensis, Arachis hypogaea, Arachis duranensis). The amino acid composition and motif analysis revealed that SWEET proteins are rich in essential amino acids such as leucine, valine, isoleucine, phenylalanine, and serine while less profuse in glutamine, tryptophan, cysteine, and histidine. A total of four main conserved motifs of SWEET proteins are also highly abundant in these amino acids. The present study deciphered the details on primary physicochemical properties, secondary, tertiary structure, and phylogenetic analysis of SWEETs protein. Majorities of SWEET proteins (72.26%) are in stable form with an average instability index of 36.5%, and it comprises a higher fraction of positively charged amino acid Arg + Lys residues. Secondary structure analysis shown that these proteins are richer in alpha-helix (40%) than extended strand (30%) and random coil (25%), respectively. Furthermore, to infer their mechanism at a structural and functional level which play an essential roles in growth, development, and stress responses. This study will be useful to examine photosynthetic productivity, embryo sugar content, seed quality, and yield enhancement in Fabaceae for a sustainable source of essential amino acids and carbon source.

Evaluation of proteome complexes normalizing osmoregulation in salt stressed Luffa acutangula (L.) Roxb.

Modern-day agriculture is facing the challenge of sustaining global food security. However, the rapid increase in salinity stress among arable areas poses a major threat to crop health and yield. Salinity stress is one of the most common and rapidly spreading stress that has a detrimental effect on the productivity of edible plant family i.e. Cucurbitaceae. The present study endeavors to evaluate the Osmoregulators (anti-oxidants and proteins), that supports the growth of two varieties of Luffa acutangula (L.) Roxb. under salt stress. The 2-3 weeks old saplings were exposed to salt stress (up to 200 mM NaCl) for one week. Post-treatment the osmoregulatory metabolites like Trehalose, Proline & enzymic anti-oxidants like peroxidase (POD), Superoxide dismutase (SOD) and proteins using LC-MS/MS were analyzed. In both the varieties, Trehalose increased with increasing salt concentration, while the level of Proline increased in Variety 1 and decreased in Variety 2. With increasing salt concentrations, the POD activity decreased in both varieties whereas that of SOD levels increased in Variety 2 and decreased in Variety 1. The protein identified by LC-MS/MS and functional annotation analysis employing Uniport database & BlastP algorithm, aided in the detection of differentially expressed proteins in response to salt stress. This was followed by metabolic interaction annotation enrichment analysis by FunRich 3.0 tool, enabling characterization of proteins to be involved in the Calvin cycle, amino acids biosynthesis, carbohydrate and energy metabolism, ROS defence, hormonal biosynthesis and signal transduction. The augmentation of the metabolic activities of the Calvin cycle, biosynthesis of amino acids, carotenoids and peroxisomes, glycolytic pathway and the tricarboxylic acid cycle will conceivably influence the photosynthetic capacity in L. acutangula varieties under salt stress. The upsurge of key enzymes involved in these above described biological processes possibly appears to play an important role in the enhancement of salt tolerance.

The mitochondrial aspartate/glutamate carrier (AGC or Aralar1) isoforms in D. melanogaster: biochemical characterization, gene structure, and evolutionary analysis

BackgroundIn man two mitochondrial aspartate/glutamate carrier (AGC) isoforms, known as aralar and citrin, are required to accomplish several metabolic pathways. In order to fill the existing gap of knowledge in Drosophila melanogaster, we have studied aralar1 gene, orthologue of human AGC-encoding genes in this organism. MethodsThe blastp algorithm and the “reciprocal best hit” approach have been used to identify the human orthologue of AGCs in Drosophilidae and non-Drosophilidae. Aralar1 proteins have been overexpressed in Escherichia coli and functionally reconstituted into liposomes for transport assays. ResultsThe transcriptional organization of aralar1 comprises six isoforms, three constitutively expressed (aralar1-RA, RD and RF), and the remaining three distributed during the development or in different tissues (aralar1-RB, RC and RE). Aralar1-PA and Aralar1-PE, representative of all isoforms, have been biochemically characterized. Recombinant Aralar1-PA and Aralar1-PE proteins share similar efficiency to exchange glutamate against aspartate, and same substrate affinities than the human isoforms. Interestingly, although Aralar1-PA and Aralar1-PE diverge only in their EF-hand 8, they greatly differ in their specific activities and substrate specificity. ConclusionsThe tight regulation of aralar1 transcripts expression and the high request of aspartate and glutamate during early embryogenesis suggest a crucial role of Aralar1 in this Drosophila developmental stage. Furthermore, biochemical characterization and calcium sensitivity have identified Aralar1-PA and Aralar1-PE as the human aralar and citrin counterparts, respectively. General significanceThe functional characterization of the fruit fly mitochondrial AGC transporter represents a crucial step toward a complete understanding of the metabolic events acting during early embryogenesis.

Characteristic of Glutamate Cysteine Ligase Gene and its response to the salinity and temperature stress in Chlamydomonas sp. ICE-L from Antarctica

Chlamydomonas sp. ICE-L, an Antarctic ice alga, has high tolerance ability to freezing and salinity. Glutathione (GSH) is an important small antioxidative molecule in the growth and stress responses of plants including algae. We cloned a full-length cDNA encoding glutamate cysteine ligase (ICE-LGCL), the key enzyme of GSH synthesis, from Chlamydomonas sp. ICE-L by RT-PCR and rapid amplification of cDNA ends technique (RACE). The cDNA has 2199 bp nucleotides with an open reading frame (ORF) of 1452 bp encoding a polypeptide of 453 amino acids. BLASTP algorithm results showed that ICE-LGCL shared 51%-70% amino acid sequence identity with the reported GCLs and shared the highest identity with Chlamydomonas reinhardtii. We also successfully made ICE-LGCL protein express in E. coli BL21. The optimum expression conditions are the induced reagent IPTG of 0.2 mmol/L, temperature of 37 °C, and the induced time of 4 h. The expression patterns of ICE-LGCL in mRNA by real-time PCR analysis showed that ICE-LGCL expressed under the different temperature and salinity challenges. Low temperature and low salinities stimulated the accumulation of ICE-LGCL mRNA in ICE-L cells. These results indicate that GCL might play an important role in Antarctic Chlamydomonas sp. ICE-L.

Fast batch searching for protein homology based on compression and clustering

BackgroundIn bioinformatics community, many tasks associate with matching a set of protein query sequences in large sequence datasets. To conduct multiple queries in the database, a common used method is to run BLAST on each original querey or on the concatenated queries. It is inefficient since it doesn’t exploit the common subsequences shared by queries.ResultsWe propose a compression and cluster based BLASTP (C2-BLASTP) algorithm to further exploit the joint information among the query sequences and the database. Firstly, the queries and database are compressed in turn by procedures of redundancy analysis, redundancy removal and distinction record. Secondly, the database is clustered according to Hamming distance among the subsequences. To improve the sensitivity and selectivity of sequence alignments, ten groups of reduced amino acid alphabets are used. Following this, the hits finding operator is implemented on the clustered database. Furthermore, an execution database is constructed based on the found potential hits, with the objective of mitigating the effect of increasing scale of the sequence database. Finally, the homology search is performed in the execution database. Experiments on NCBI NR database demonstrate the effectiveness of the proposed C2-BLASTP for batch searching of homology in sequence database. The results are evaluated in terms of homology accuracy, search speed and memory usage.ConclusionsIt can be seen that the C2-BLASTP achieves competitive results as compared with some state-of-the-art methods.

MuBLASTP: database-indexed protein sequence search on multicore CPUs.

BackgroundThe Basic Local Alignment Search Tool (BLAST) is a fundamental program in the life sciences that searches databases for sequences that are most similar to a query sequence. Currently, the BLAST algorithm utilizes a query-indexed approach. Although many approaches suggest that sequence search with a database index can achieve much higher throughput (e.g., BLAT, SSAHA, and CAFE), they cannot deliver the same level of sensitivity as the query-indexed BLAST, i.e., NCBI BLAST, or they can only support nucleotide sequence search, e.g., MegaBLAST. Due to different challenges and characteristics between query indexing and database indexing, the existing techniques for query-indexed search cannot be used into database indexed search.ResultsmuBLASTP, a novel database-indexed BLAST for protein sequence search, delivers identical hits returned to NCBI BLAST. On Intel Haswell multicore CPUs, for a single query, the single-threaded muBLASTP achieves up to a 4.41-fold speedup for alignment stages, and up to a 1.75-fold end-to-end speedup over single-threaded NCBI BLAST. For a batch of queries, the multithreaded muBLASTP achieves up to a 5.7-fold speedups for alignment stages, and up to a 4.56-fold end-to-end speedup over multithreaded NCBI BLAST.ConclusionsWith a newly designed index structure for protein database and associated optimizations in BLASTP algorithm, we re-factored BLASTP algorithm for modern multicore processors that achieves much higher throughput with acceptable memory footprint for the database index.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1302-4) contains supplementary material, which is available to authorized users.

Exploring the Molecular Basis for Binding of Inhibitors by Threonyl-tRNA Synthetase from Brucella abortus: A Virtual Screening Study.

Targeting threonyl-tRNA synthetase (ThrRS) of Brucella abortus is a promising approach to developing small-molecule drugs against bovine brucellosis. Using the BLASTp algorithm, we identified ThrRS from Escherichia coli (EThrRS, PDB ID 1QF6), which is 51% identical to ThrRS from Brucella abortus (BaThrRS) at the amino acid sequence level. EThrRS was used as the template to construct a BaThrRS homology model which was optimized using molecular dynamics simulations. To determine the residues important for substrate ATP binding, we identified the ATP-binding regions of BaThrRS, docked ATP to the protein, and identified the residues whose side chains surrounded bound ATP. We then used the binding site of ATP to virtually screen for BaThrRS inhibitors and got seven leads. We further characterized the BaThrRS-binding site of the compound with the highest predicted inhibitory activity. Our results should facilitate future experimental effects to find novel drugs for use against bovine brucellosis.

Commentary: Plant Auxin Biosynthesis Did Not Originate in Charophytes.

Commentary: Plant Auxin Biosynthesis Did Not Originate in Charophytes.

A Chemical Proteomics Approach to Profiling the ATP-binding Proteome of Mycobacterium tuberculosis

Tuberculosis, caused by Mycobacterium tuberculosis, remains one of the leading causes of death worldwide despite extensive research, directly observed therapy using multidrug regimens, and the widespread use of a vaccine. The majority of patients harbor the bacterium in a state of metabolic dormancy. New drugs with novel modes of action are needed to target essential metabolic pathways in M. tuberculosis; ATP-competitive enzyme inhibitors are one such class. Previous screening efforts for ATP-competitive enzyme inhibitors identified several classes of lead compounds that demonstrated potent anti-mycobacterial efficacy as well as tolerable levels of toxicity in cell culture. In this report, a probe-based chemoproteomic approach was used to selectively profile the M. tuberculosis ATP-binding proteome in normally growing and hypoxic M. tuberculosis. From these studies, 122 ATP-binding proteins were identified in either metabolic state, and roughly 60% of these are reported to be essential for survival in vitro. These data are available through ProteomeXchange with identifier PXD000141. Protein families vital to the survival of the tubercle bacillus during hypoxia emerged from our studies. Specifically, along with members of the DosR regulon, several proteins involved in energy metabolism (Icl/Rv0468 and Mdh/Rv1240) and lipid biosynthesis (UmaA/Rv0469, DesA1/Rv0824c, and DesA2/Rv1094) were found to be differentially abundant in hypoxic versus normal growing cultures. These pathways represent a subset of proteins that may be relevant therapeutic targets for development of novel ATP-competitive antibiotics.

Chaperone-usher fimbriae of Escherichia coli.

Chaperone-usher (CU) fimbriae are adhesive surface organelles common to many Gram-negative bacteria. Escherichia coli genomes contain a large variety of characterised and putative CU fimbrial operons, however, the classification and annotation of individual loci remains problematic. Here we describe a classification model based on usher phylogeny and genomic locus position to categorise the CU fimbrial types of E. coli. Using the BLASTp algorithm, an iterative usher protein search was performed to identify CU fimbrial operons from 35 E. coli (and one Escherichia fergusonnii) genomes representing different pathogenic and phylogenic lineages, as well as 132 Escherichia spp. plasmids. A total of 458 CU fimbrial operons were identified, which represent 38 distinct fimbrial types based on genomic locus position and usher phylogeny. The majority of fimbrial operon types occupied a specific locus position on the E. coli chromosome; exceptions were associated with mobile genetic elements. A group of core-associated E. coli CU fimbriae were defined and include the Type 1, Yad, Yeh, Yfc, Mat, F9 and Ybg fimbriae. These genes were present as intact or disrupted operons at the same genetic locus in almost all genomes examined. Evaluation of the distribution and prevalence of CU fimbrial types among different pathogenic and phylogenic groups provides an overview of group specific fimbrial profiles and insight into the ancestry and evolution of CU fimbriae in E. coli.

CUDA-BLASTP: Accelerating BLASTP on CUDA-Enabled Graphics Hardware

Scanning protein sequence database is an often repeated task in computational biology and bioinformatics. However, scanning large protein databases, such as GenBank, with popular tools such as BLASTP requires long runtimes on sequential architectures. Due to the continuing rapid growth of sequence databases, there is a high demand to accelerate this task. In this paper, we demonstrate how GPUs, powered by the Compute Unified Device Architecture (CUDA), can be used as an efficient computational platform to accelerate the BLASTP algorithm. In order to exploit the GPU’s capabilities for accelerating BLASTP, we have used a compressed deterministic finite state automaton for hit detection as well as a hybrid parallelization scheme. Our implementation achieves speedups up to 10.0 on an NVIDIA GeForce GTX 295 GPU compared to the sequential NCBI BLASTP 2.2.22. CUDA-BLASTP source code which is available at https://sites.google.com/site/liuweiguohome/software.

BlastR—fast and accurate database searches for non-coding RNAs

We present and validate BlastR, a method for efficiently and accurately searching non-coding RNAs. Our approach relies on the comparison of di-nucleotides using BlosumR, a new log-odd substitution matrix. In order to use BlosumR for comparison, we recoded RNA sequences into protein-like sequences. We then showed that BlosumR can be used along with the BlastP algorithm in order to search non-coding RNA sequences. Using Rfam as a gold standard, we benchmarked this approach and show BlastR to be more sensitive than BlastN. We also show that BlastR is both faster and more sensitive than BlastP used with a single nucleotide log-odd substitution matrix. BlastR, when used in combination with WU-BlastP, is about 5% more accurate than WU-BlastN and about 50 times slower. The approach shown here is equally effective when combined with the NCBI-Blast package. The software is an open source freeware available from www.tcoffee.org/blastr.html.

Massively parallel FPGA-based implementation of BLASTp with the two-hit method

Abstract Protein database search requests are generally being performed using the BLASTp algorithm, introduced by NCBI [1] . Since it is computationally intensive, it becomes more and more ine_ective with today's growth of sequence database sizes. The needs for an e_cient parallel implementation arise. In this paper, we focus on a massive parallelization using the FPGA-based hardware architecture RIVYERA [2] . The aim is to reach speedups in orders of magnitude with a flexible implementation while saving energy costs compared to PC-based database search. We keep our implementation close to the structure published by Kasap et al. [3] , [4] and include ideas from Sotiriades et al. [5] such that all parts of the algorithm are organized in components of a long pipeline. We also use the idea of the two-hit method [6] to keep the computational e_ort small. Besides the related work, we perform the very last step of the algorithm to produce a gapped alignment with the Needleman-Wunsch algorithm in software, only with the option of hardware processing after reconfiguration. This saves FPGA-resources and allows an even higher degree of parallelism.

Pathogen Proteins Eliciting Antibodies Do Not Share Epitopes with Host Proteins: A Bioinformatics Approach

The best way to prevent diseases caused by pathogens is by the use of vaccines. The advent of genomics enables genome-wide searches of new vaccine candidates, called reverse vaccinology. The most common strategy to apply reverse vaccinology is by designing subunit recombinant vaccines, which usually generate an humoral immune response due to B-cell epitopes in proteins. A major problem for this strategy is the identification of protective immunogenic proteins from the surfome of the pathogen. Epitope mimicry may lead to auto-immune phenomena related to several human diseases. A sequence-based computational analysis has been carried out applying the BLASTP algorithm. Therefore, two huge databases have been created, one with the most complete and current linear B-cell epitopes, and the other one with the surface-protein sequences of the main human respiratory bacterial pathogens. We found that none of the 7353 linear B-cell epitopes analysed shares any sequence identity region with human proteins capable of generating antibodies, and that only 1% of the 2175 exposed proteins analysed contain a stretch of shared sequence with the human proteome. These findings suggest the existence of a mechanism to avoid autoimmunity. We also propose a strategy for corroborating or warning about the viability of a protein linear B-cell epitope as a putative vaccine candidate in a reverse vaccinology study; so, epitopes without any sequence identity with human proteins should be very good vaccine candidates, and the other way around.

Identification and Functional Characterization of a Monofunctional Peroxisomal Enoyl-CoA Hydratase 2 That Participates in the Degradation of Even cis-Unsaturated Fatty Acids in Arabidopsis thaliana

A gene, named AtECH2, has been identified in Arabidopsis thaliana to encode a monofunctional peroxisomal enoyl-CoA hydratase 2. Homologues of AtECH2 are present in several angiosperms belonging to the Monocotyledon and Dicotyledon classes, as well as in a gymnosperm. In vitro enzyme assays demonstrated that AtECH2 catalyzed the reversible conversion of 2E-enoyl-CoA to 3R-hydroxyacyl-CoA. AtECH2 was also demonstrated to have enoyl-CoA hydratase 2 activity in an in vivo assay relying on the synthesis of polyhydroxyalkanoate from the polymerization of 3R-hydroxyacyl-CoA in the peroxisomes of Saccharomyces cerevisiae. AtECH2 contained a peroxisome targeting signal at the C-terminal end, was addressed to the peroxisome in S. cerevisiae, and a fusion protein between AtECH2 and a fluorescent protein was targeted to peroxisomes in onion cells. AtECH2 gene expression was strongest in tissues with high beta-oxidation activity, such as germinating seedlings and senescing leaves. The contribution of AtECH2 to the degradation of unsaturated fatty acids was assessed by analyzing the carbon flux through the beta-oxidation cycle in plants that synthesize peroxisomal polyhydroxyalkanoate and that were over- or underexpressing the AtECH2 gene. These studies revealed that AtECH2 participates in vivo to the conversion of the intermediate 3R-hydroxyacyl-CoA, generated by the metabolism of fatty acids with a cis (Z)-unsaturated bond on an even-numbered carbon, to the 2E-enoyl-CoA for further degradation through the core beta-oxidation cycle.

PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database

BackgroundEqually important and challenging as genome annotation, is the subsequent classification of predicted genes into their respective pathways. The Kyoto Encyclopedia of Genes and Genomes (KEGG) represents a database consisting of known genes and their respective biochemical functionalities. Although accessible online, analyses of multiple genes are time consuming and are not suitable for analyzing data sets that are proprietary.ResultsPresented here is a new software solution that utilizes the KEGG online database for pathway mapping of partial and whole prokaryotic genomes. PathwayVoyager retrieves user-defined subsets of the KEGG database and stores the data as local, blast-formatted databases. Previously selected datasets can be re-used, reducing run-time significantly. Whole or partial genomes can be automatically analyzed using NCBI's BlastP algorithm and ORFs with similarities below the user-defined threshold will be marked on pathway maps. Multiple gene hits are sorted by similarity. Since no sequence information is transmitted over the Internet, PathwayVoyager is an ideal solution for pathway mapping and reconstruction of confidential DNA sequence data.ConclusionPathwayVoyager represents an alternative approach to many already existing, more complex pathway reconstructions software solutions. This software does not require any dedicated hardware or software and is flexible and straightforward to use. It is ideally suited for environments where analyses on variable datasets are desired.

Elicitors of Plant Defense Responses from Biocontrol Strains of Trichoderma viren

ABSTRACT Effective biocontrol strains of Trichoderma virens can induce the production of defense-related compounds in the roots of cotton. Ineffective strains do not induce these compounds to significant levels. This elicittation was found to be heat stable, insoluble in chloroform, passed through a 5K molecular weight cut-off (MWCO) filter, but not a 3K MWCO filter, and was sensitive to treatment by proteinase K. When the active material was subjected to sodium dodecyl sulfate-polyacrylamide gel electrophoresis, several bands were present in the material from biocontrol-active strains that were lacking in inactive strains. When eluted and tested for elicitation activity, with or without renaturation, four bands stimulated cotton terpenoid production. One band showed cross-reaction with an antibody to the ethylene-inducing xylanase from T. viride. Another band of approximately 18 kDa, gave significant stimulation of cotton terpenoid production and increased peroxidase activity in cotton radicles in all tests, with or without renaturation. The 18-kDa protein was subjected to amino-terminal sequence analysis, and the first 19 amino acids at the amino terminus were determined to be DTVSYDTGYDNGSRSLNDV. A database homology search using the BLASTp algorithm showed the highest similarity to a serine proteinase from Fusarium sporotrichioides.

Transcriptional Down-regulation of Peroxisome Numbers Affects Selective Peroxisome Degradation in Hansenula polymorpha

We have isolated and characterized a novel transcription factor of Hansenula polymorpha that is involved in the regulation of peroxisomal protein levels. This protein, designated Mpp1p, belongs to the family of Zn(II)2Cys6 proteins. In cells deleted for the function of Mpp1p the levels of various proteins involved in peroxisome biogenesis (peroxins) and function (enzymes) are reduced compared with wild type or, in the case of the matrix protein dihydroxyacetone synthase, fully absent. Also, upon induction of mpp1 cells on methanol, the number of peroxisomes was strongly reduced relative to wild type cells and generally amounted to one organelle per cell. Remarkably, this single organelle was not susceptible to selective peroxisome degradation (pexophagy) and remained unaffected during exposure of methanol-induced cells to excess glucose conditions. We show that this mechanism is a general phenomenon in H. polymorpha in the case of cells that contain only a single peroxisome.