VIRALpro: a tool to identify viral capsid and tail sequences.
Ontology highlight
ABSTRACT: Motivation
Not only sequence data continue to outpace annotation information, but also the problem is further exacerbated when organisms are underrepresented in the annotation databases. This is the case with non-human-pathogenic viruses which occur frequently in metagenomic projects. Thus, there is a need for tools capable of detecting and classifying viral sequences.Results
We describe VIRALpro a new effective tool for identifying capsid and tail protein sequences, which are the cornerstones toward viral sequence annotation and viral genome classification.Availability and implementation
The data, software and corresponding web server are available from http://scratch.proteomics.ics.uci.edu as part of the SCRATCH suite.Contact
clovis.galiez@inria.fr or pfbaldi@uci.eduSupplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Galiez C
PROVIDER: S-EPMC5860506 | biostudies-literature | 2016 May
REPOSITORIES: biostudies-literature
Publications
VIRALpro: a tool to identify viral capsid and tail sequences.
Galiez Clovis C Magnan Christophe N CN Coste Francois F Baldi Pierre P
Bioinformatics (Oxford, England) 20160105 9
<h4>Motivation</h4>Not only sequence data continue to outpace annotation information, but also the problem is further exacerbated when organisms are underrepresented in the annotation databases. This is the case with non-human-pathogenic viruses which occur frequently in metagenomic projects. Thus, there is a need for tools capable of detecting and classifying viral sequences.<h4>Results</h4>We describe VIRALpro a new effective tool for identifying capsid and tail protein sequences, which are th ...[more]
Similar Datasets
Project description:We now have access to the sequences of tens of millions of proteins. These protein sequences are essential for modern molecular biology and computational biology. The vast majority of protein sequences are derived from gene prediction tools and have no experimental supporting evidence for their translation. Despite the increasing accuracy of gene prediction tools there likely exists a large number of spurious protein predictions in the sequence databases. We have developed the Spurio tool to help identify spurious protein predictions in prokaryotes. Spurio searches the query protein sequence against a prokaryotic nucleotide database using tblastn and identifies homologous sequences. The tblastn matches are used to score the query sequence's likelihood of being a spurious protein prediction using a Gaussian process model. The most informative feature is the appearance of stop codons within the presumed translation of homologous DNA sequences. Benchmarking shows that the Spurio tool is able to distinguish spurious from true proteins. However, transposon proteins are prone to be predicted as spurious because of the frequency of degraded homologs found in the DNA sequence databases. Our initial experiments suggest that less than 1% of the proteins in the UniProtKB sequence database are likely to be spurious and that Spurio is able to identify over 60 times more spurious proteins than the AntiFam resource. The Spurio software and source code is available under an MIT license at the following URL: https://bitbucket.org/bateman-group/spurio.
| S-EPMC5897793 | biostudies-literature
Project description:Capturing group-specific sequences between two groups of genomic/metagenomic sequences is critical for the follow-up identifications of singular nucleotide variants (SNVs), gene families, microbial species or other elements associated with each group. A sequence that is present, or rich, in one group, but absent, or scarce, in another group is considered a "group-specific" sequence in our study. We developed a user-friendly tool, KmerGO, to identify group-specific sequences between two groups of genomic/metagenomic long sequences or high-throughput sequencing datasets. Compared with other tools, KmerGO captures group-specific k-mers (k up to 40 bps) with much lower requirements for computing resources in much shorter running time. For a 1.05 TB dataset (.fasta), it takes KmerGO about 21.5 h (including k-mer counting) to return assembled group-specific sequences on a regular stand-alone workstation with no more than 1 GB memory. Furthermore, KmerGO can also be applied to capture trait-associated sequences for continuous-trait. Through multi-process parallel computing, KmerGO is implemented with both graphic user interface and command line on Linux and Windows free from any pre-installed supporting environments, packages, and complex configurations. The output group-specific k-mers or sequences from KmerGO could be the inputs of other tools for the downstream discovery of biomarkers, such as genetic variants, species, or genes. KmerGO is available at https://github.com/ChnMasterOG/KmerGO.
| S-EPMC7477287 | biostudies-literature
Project description:I present a review of the theoretical and computational methodologies that have been used to model the assembly of viral capsids. I discuss the capabilities and limitations of approaches ranging from equilibrium continuum theories to molecular dynamics simulations, and I give an overview of some of the important conclusions about virus assembly that have resulted from these modeling efforts. Topics include the assembly of empty viral shells, assembly around single-stranded nucleic acids to form viral particles, and assembly around synthetic polymers or charged nanoparticles for nanotechnology or biomedical applications. I present some examples in which modeling efforts have promoted experimental breakthroughs, as well as directions in which the connection between modeling and experiment can be strengthened.
| S-EPMC4318123 | biostudies-literature
Project description:Bacteriophages isolated on Mycobacterium smegmatis mc(2)155 represent many distinct genomes sharing little or no DNA sequence similarity. The genomes are architecturally mosaic and are replete with genes of unknown function. A new group of genomes sharing substantial nucleotide sequences constitute Cluster J. The six mycobacteriophages forming Cluster J are morphologically members of the Siphoviridae, but have unusually long genomes ranging from 106.3 to 117 kbp. Reconstruction of the capsid by cryo-electron microscopy of mycobacteriophage BAKA reveals an icosahedral structure with a triangulation number of 13. All six phages are temperate and homoimmune, and prophage establishment involves integration into a tRNA-Leu gene not previously identified as a mycobacterial attB site for phage integration. The Cluster J genomes provide two examples of intron splicing within the virion structural genes, one in a major capsid subunit gene, and one in a tail gene. These genomes also contain numerous free-standing HNH homing endonuclease, and comparative analysis reveals how these could contribute to genome mosaicism. The unusual Cluster J genomes provide new insights into phage genome architecture, gene function, capsid structure, gene mobility, intron splicing, and evolution.
| S-EPMC3706429 | biostudies-literature
Project description:Identifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.We have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder's performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.VirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder's potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.This innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.
| S-EPMC5501583 | biostudies-literature
Project description:The HIV-1 capsid, a conical shell encasing viral nucleoprotein complexes, is involved in multiple post-entry processes during viral replication. Many host factors can directly bind to the HIV-1 capsid protein (CA) and either promote or prevent HIV-1 infection. The viral capsid is currently being explored as a novel target for therapeutic interventions. In the past few decades, significant progress has been made in our understanding of the capsid-host interactions and mechanisms of action of capsid-targeting antivirals. At the same time, a large number of different viral capsids, which derive from many HIV-1 mutants, naturally occurring variants, or diverse lentiviruses, have been characterized for their interactions with capsid-binding molecules in great detail utilizing various experimental techniques. This review provides an overview of how sequence variation in CA influences phenotypic properties of HIV-1. We will focus on sequence differences that alter capsid-host interactions and give a brief account of drug resistant mutations in CA and their mutational effects on viral phenotypes. Increased knowledge of the sequence-function relationship of CA helps us deepen our understanding of the adaptive potential of the viral capsid.
| S-EPMC8549334 | biostudies-literature
Project description:Most of the existing research in assembly pathway prediction/analysis of viral capsids makes the simplifying assumption that the configuration of the intermediate states can be extracted directly from the final configuration of the entire capsid. This assumption does not take into account the conformational changes of the constituent proteins as well as minor changes to the binding interfaces that continue throughout the assembly process until stabilization. This article presents a statistical-ensemble-based approach that samples the configurational space for each monomer with the relative local orientation between monomers, to capture the uncertainties in binding and conformations. Further, instead of using larger capsomers (trimers, pentamers) as building blocks, we allow all possible subassemblies to bind in all possible combinations. We represent the resulting assembly graph in two different ways: First, we use the Wilcoxon signed-rank measure to compare the distributions of binding free energy computed on the sampled conformations to predict likely pathways. Second, we represent chemical equilibrium aspects of the transitions as a Bayesian Factor graph where both associations and dissociations are modeled based on concentrations and the binding free energies. We applied these protocols on the feline panleukopenia virus and the Nudaurelia capensis virus. Results from these experiments showed a significant departure from those that one would obtain if only the static configurations of the proteins were considered. Hence, we establish the importance of an uncertainty-aware protocol for pathway analysis, and we provide a statistical framework as an important first step toward assembly pathway prediction with high statistical confidence.
| S-EPMC5757074 | biostudies-literature
Project description:Lambda-like double-stranded (ds) DNA bacteriophage undergo massive conformational changes in their capsid shell during the packaging of their viral genomes. Capsid shells are complex organizations of hundreds of protein subunits that assemble into intricate quaternary complexes that ultimately are able to withstand over 50 atm of pressure during genome packaging. The extensive integration between subunits in capsids requires the formation of an intermediate complex, termed a procapsid, from which individual subunits can undergo the necessary refolding and structural rearrangements needed to transition to the more stable capsid. Although various mature capsids have been characterized at atomic resolution, no such procapsid structure is available for a dsDNA virus or bacteriophage. Here we present a procapsid X-ray structure at 3.65 A resolution, termed prohead II, of the lambda-like bacteriophage HK97, the mature capsid structure of which was previously solved to 3.44 A (ref. 2). A comparison of the two largely different capsid forms has unveiled an unprecedented expansion mechanism that describes the transition. Crystallographic and hydrogen/deuterium exchange data presented here demonstrate that the subunit tertiary structures are significantly different between the two states, with twisting and bending motions occurring in both helical and beta-sheet regions. We also identified subunit interactions at each three-fold axis of the capsid that are maintained throughout maturation. The interactions sustain capsid integrity during subunit refolding and provide a fixed hinge from which subunits undergo rotational and translational motions during maturation. Previously published calorimetric data of a closely related bacteriophage, P22, showed that capsid maturation was an exothermic process that resulted in a release of 90 kJ mol(-1) of energy. We propose that the major tertiary changes presented in this study reveal a structural basis for an exothermic maturation process probably present in many dsDNA bacteriophage and possibly viruses such as herpesvirus, which share the HK97 subunit fold.
| S-EPMC2765791 | biostudies-literature
Project description:Stable structures of icosahedral symmetry can serve numerous functional roles, including chemical microencapsulation and delivery of drugs and biomolecules, epitope presentation to allow for an efficient immunization process, synthesis of nanoparticles of uniform size, observation of encapsulated reactive intermediates, formation of structural elements for supramolecular constructs, and molecular computing. By examining physical models of spherical virus assembly we have arrived at a general synthetic strategy for producing chemical capsids at size scales between fullerenes and spherical viruses. Such capsids can be formed by self-assembly from a class of molecules developed from a symmetric pentagonal core. By designing chemical complementarity into the five interface edges of the molecule, we can produce self-assembling stable structures of icosahedral symmetry. We considered three different binding mechanisms: hydrogen bonding, metal binding, and formation of disulfide bonds. These structures can be designed to assemble and disassemble under controlled environmental conditions. We have conducted molecular dynamics simulation on a class of corannulene-based molecules to demonstrate the characteristics of self-assembly and to aid in the design of the molecular subunits. The edge complementarities can be of diverse structure, and they need not reflect the fivefold symmetry of the molecular core. Thus, self-assembling capsids formed from coded subunits can serve as addressable nanocontainers or custom-made structural elements.
| S-EPMC2410070 | biostudies-literature
Project description:R-loops are RNA-DNA hybrid sequences that are emerging players in various biological processes, occurring in both prokaryotic and eukaryotic cells. In viruses, R-loop investigation is limited and functional importance is poorly understood. Here, we performed a computational approach to investigate prevalence, distribution, and location of R-loop forming sequences (RLFS) across more than 6000 viral genomes. A total of 14637 RLFS loci were identified in 1586 viral genomes. Over 70% of RLFS-positive genomes are dsDNA viruses. In the order Herpesvirales, RLFS were presented in all members whereas no RLFS was predicted in the order Ligamenvirales. Analysis of RLFS density in all RLFS-positive genomes revealed unusually high RLFS densities in herpesvirus genomes, with RLFS densities particularly enriched within repeat regions such as the terminal repeats (TRs). RLFS in TRs are positionally conserved between herpesviruses. Validating the computationally-identified RLFS, R-loop formation was experimentally confirmed in the TR and viral Bcl-2 promoter of Kaposi sarcoma-associated herpesvirus (KSHV). These predictions and validations support future analysis of RLFS in regulating the replication, transcription, and genome maintenance of herpesviruses.
| S-EPMC7156643 | biostudies-literature