ABSTRACT
The far-reaching effects of the SARS-CoV-2 pandemic have crippled the progress of the world today. With the introduction of newer and newer mutated variants of the virus, it has become necessary to have a vaccine that remains useful against all the mutated strains of SARS-CoV-2. In this regard, peptide vaccines turn out to be a cheap alternative to the traditionally designed vaccines owing to their much quicker and computationally easier, and more robust design procedures. Here, in this article, we hypothesize that there are three possible peptide vaccine regions that can be targeted to prevent the surge of SARS-CoV-2. The candidates that were selected, were surface-exposed and were not sequestered by any neighbouring amino acids. They were also found to be capable of generating both B-cell and T-cell immune responses. Most importantly, none of them contains any spike protein mutation of the currently prevailing variants of SARS-CoV-2. From these findings, we have therefore concluded that these three regions can be used in wet labs for peptide vaccine design against the upcoming strains of SARS-CoV-2.
ABSTRACT
The exponential growth in the depositories of biological sequence data has generated an urgent need to store, retrieve and analyse the data efficiently and effectively for which the standard practice of using alignment procedures are not adequate due to high demand on computing resources and time. Graphical representation of sequences has become one of the most popular alignment-free strategies to analyse the biological sequences where each basic unit of the sequences - the bases adenine, cytosine, guanine and thymine for DNA/RNA, and the 20 amino acids for proteins - are plotted on a multi-dimensional grid. The resulting curve in 2D and 3D space and the implied graph in higher dimensions provide a perception of the underlying information of the sequences through visual inspection; numerical analyses, in geometrical or matrix terms, of the plots provide a measure of comparison between sequences and thus enable study of sequence hierarchies. The new approach has also enabled studies of comparisons of DNA sequences over many thousands of bases and provided new insights into the structure of the base compositions of DNA sequences. In this article we review in brief the origins and applications of graphical representations and highlight the future perspectives in this field.
Subject(s)
DNA , RNA , DNA/genetics , Sequence Analysis, DNA/methodsABSTRACT
The design for vaccines using in silico analysis of genomic data of different viruses has taken many different paths, but lack of any precise computational approach has constrained them to alignment methods and some alignment-free techniques. In this work, a precise computational approach has been established wherein two new mathematical parameters have been suggested to identify the highly conserved and surface-exposed regions which are spread over a large region of the surface protein of the virus so that one can determine possible peptide vaccine candidates from those regions. The first parameter, w, is the sum of the normalized values of the measure of surface accessibility and the normalized measure of conservativeness, and the second parameter is the area of a triangle formed by a mathematical model named 2D Polygon Representation. This method has been, therefore, used to determine possible vaccine targets against SARS-CoV-2 by considering its surface-situated spike glycoprotein. The results of this model have been verified by a parallel analysis using the older approach of manually estimating the graphs describing the variation of conservativeness and surface-exposure across the protein sequence. Furthermore, the working of the method has been tested by applying it to find out peptide vaccine candidates for Zika and Hendra viruses respectively. A satisfactory consistency of the model results with pre-established results for both the test cases shows that this in silico alignment-free analysis proposed by the model is suitable not only to determine vaccine targets against SARS-CoV-2 but also ready to extend against other viruses.
ABSTRACT
INTRODUCTION: Coronaviruses comprise a group of enveloped, positive-sense single-stranded RNA viruses that infect humans as well as a wide range of animals. The study was performed on a set of 573 sequences belonging to SARS, MERS and SARS-CoV-2 (CoVID-19) viruses. The sequences were represented with alignment-free sequence descriptors and analyzed with different chemometric methods: Euclidean/Mahalanobis distances, principal component analysis and self-organizing maps (Kohonen networks). We report the cluster structures of the data. The sequences are well-clustered regarding the type of virus; however, some of them show the tendency to belong to more than one virus type. BACKGROUND: This is a study of 573 genome sequences belonging to SARS, MERS and SARS-- CoV-2 (CoVID-19) coronaviruses. OBJECTIVES: The aim was to compare the virus sequences, which originate from different places around the world. METHODS: The study used alignment free sequence descriptors for the representation of sequences and chemometric methods for analyzing clusters. RESULTS: Majority of genome sequences are clustered with respect to the virus type, but some of them are outliers. CONCLUSION: We indicate 71 sequences, which tend to belong to more than one cluster.
Subject(s)
COVID-19 , SARS-CoV-2 , Animals , Cluster Analysis , HumansABSTRACT
BACKGROUND: In this report, we consider a data set, which consists of 310 Zika virus genome sequences taken from different continents, Africa, Asia and South America. The sequences, which were compiled from GenBank, were derived from the host cells of different mammalian species (Simiiformes, Aedes opok, Aedes africanus, Aedes luteocephalus, Aedes dalzieli, Aedes aegypti, and Homo sapiens). METHODS: For chemometrical treatment, the sequences have been represented by sequence descriptors derived from their graphs or neighborhood matrices. The set was analyzed with three chemometrical methods: Mahalanobis distances, principal component analysis (PCA) and self organizing maps (SOM). A good separation of samples with respect to the region of origin was observed using these three methods. RESULTS: Study of 310 Zika virus genome sequences from different continents. To characterize and compare Zika virus sequences from around the world using alignment-free sequence comparison and chemometrical methods. CONCLUSION: Mahalanobis distance analysis, self organizing maps, principal components were used to carry out the chemometrical analyses of the Zika sequence data. Genome sequences are clustered with respect to the region of origin (continent, country). Africa samples are well separated from Asian and South American ones.
Subject(s)
Computer Simulation , Databases, Genetic , Sequence Analysis, RNA/methods , Zika Virus Infection/epidemiology , Zika Virus Infection/genetics , Zika Virus/genetics , Africa/epidemiology , Animals , Asia/epidemiology , Cluster Analysis , Humans , South America/epidemiologyABSTRACT
SARS-CoV-2 infection has become a worldwide pandemic and is spreading rapidly to people across the globe. To combat the situation, vaccine design is the essential solution. Mutation in the virus genome plays an important role in limiting the working life of a vaccine. In this study, we have identified several mutated clusters in the structural proteins of the virus through our novel 2D Polar plot and qR characterization descriptor. We have also studied several biochemical properties of the proteins to explore the dynamics of evolution of these mutations. This study would be helpful to understand further new mutations in the virus and would facilitate the process of designing a sustainable vaccine against the deadly virus.
Subject(s)
Computational Biology , Mutation , SARS-CoV-2/genetics , Algorithms , Genes, Viral , Quantitative Structure-Activity RelationshipABSTRACT
With the increasing frequency of viral epidemics, vaccines to augment the human immune response system have been the medium of choice to combat viral infections. The tragic consequences of the Zika virus pandemic in South and Central America a few years ago brought the issues into sharper focus. While traditional vaccine development is time-consuming and expensive, recent advances in information technology, immunoinformatics, genetics, bioinformatics, and related sciences have opened the doors to new paradigms in vaccine design and applications.Peptide vaccines are one group of the new approaches to vaccine formulation. In this chapter, we discuss the various issues involved in the design of peptide vaccines and their advantages and shortcomings, with special reference to the Zika virus for which no drugs or vaccines are as yet available. In the process, we outline our work in this field giving a detailed step-by-step description of the protocol we follow for such vaccine design so that interested researchers can easily follow them and do their own designing. Several flowcharts and figures are included to provide a background of the software to be used and results to be anticipated.
Subject(s)
Computational Biology/methods , Vaccines, Subunit/genetics , Viral Proteins/chemistry , Zika Virus/immunology , Humans , Mutation , Vaccines, Subunit/immunology , Viral Proteins/genetics , Viral Proteins/immunology , Zika Virus/geneticsABSTRACT
Human life has been at the edge of catastrophe for millennia due diseases which emerge and reemerge at random. The recent outbreak of the Zika virus (ZIKV) is one such menace that shook the global public health community abruptly. Modern technologies, including computational tools as well as experimental approaches, need to be harnessed fast and effectively in a coordinated manner in order to properly address such challenges. In this paper, based on our earlier research, we have proposed a four-pronged approach to tackle the emerging pathogens like ZIKV: (a) Epidemiological modelling of spread mechanisms of ZIKV; (b) assessment of the public health risk of newly emerging strains of the pathogens by comparing them with existing strains/pathogens using fast computational sequence comparison methods; (c) implementation of vaccine design methods in order to produce a set of probable peptide vaccine candidates for quick synthesis/production and testing in the laboratory; and (d) designing of novel therapeutic molecules and their laboratory testing as well as validation of new drugs or repurposing of drugs for use against ZIKV. For each of these stages, we provide an extensive review of the technical challenges and current state-of-the-art. Further, we outline the future areas of research and discuss how they can work together to proactively combat ZIKV or future emerging pathogens.
ABSTRACT
INTRODUCTION: Among the mosquito-borne human-infecting flavivirus species that include Zika, West Nile, yellow fever, Japanese encephalitis and Dengue viruses, the Zika virus is found to be closest to Dengue virus, sharing the same clade in the Flavivirus phylogenetic tree. We consider these five flaviviruses and on closer examination in our analyses, the nucleotide sequences of the Dengue viral genes (envelope and NS5) and genomes are seen to be quite widely different from the other four flaviviruses. We consider the extent of this distinction and determine the advantage and/or disadvantage such differences may confer upon the Dengue viral pathogenesis. METHODS: We have primarily used a 2D graphical representation technique to show the differences in base distributions in these five flaviviruses and subsequently, obtained quantitative estimates of the differences. Similarity/dissimilarity between the viruses based on the genes were also determined which showed that the differences with the Dengue genes are more pronounced. RESULTS: We found that the Dengue viruses compared to the other four flaviviruses spread rapidly worldwide and became endemic in various regions with small alterations in sequence composition relative to the host populations as revealed by codon usage biases and phylogenetic examination. CONCLUSION: We conclude that the Dengue genes are indeed more widely separated from the other aforementioned mosquito-borne human-infecting flaviviruses due to excess adenine component, a feature that is sparse in the literature. Such excesses have a bearing on drug and vaccine, especially peptide vaccine, development and should be considered appropriately.
Subject(s)
Adenine , Dengue Virus/genetics , Flavivirus/genetics , Base Sequence/genetics , Evolution, MolecularABSTRACT
We briefly review the situations arising out of epidemics that erupt rather suddenly, threatening life and livelihoods of humans. Ebola, Zika and the Nipah virus outbreaks are recent examples where the viral epidemics have led to considerably high degree of fatalities or debilitating consequences. The problems are accentuated by a lack of drugs or vaccines effective against the new and emergent viruses, and the inordinate amount of temporal and financial resources that are required to combat the novel pathogens. Progress in computational, biological and informational sciences have made it possible to consider design of synthetic vaccines that can be rapidly developed and deployed to help stem the damages. In this review, we consider the pros and cons of this new paradigm and suggest a new system where the manufacturing process can be decentralized to provide more targeted vaccines to meet the urgent needs of protection in case of a rampaging epidemic.
Subject(s)
Hemorrhagic Fever, Ebola/epidemiology , Henipavirus Infections/epidemiology , Peptides/immunology , Viral Vaccines/immunology , Zika Virus Infection/epidemiology , Animals , Hemorrhagic Fever, Ebola/immunology , Hemorrhagic Fever, Ebola/prevention & control , Hemorrhagic Fever, Ebola/virology , Henipavirus Infections/immunology , Henipavirus Infections/prevention & control , Humans , Nipah Virus/immunology , Zika Virus Infection/immunology , Zika Virus Infection/prevention & control , Zika Virus Infection/virologyABSTRACT
Graphical representation methods constitute a class of alignment-free techniques for comparative study of biomolecular sequences. In this brief commentary, we study how well some of these methods can discriminate among closely related genes.
Subject(s)
Computer-Aided Design , Influenza A Virus, H7N9 Subtype/immunology , Influenza Vaccines/immunology , Influenza, Human/prevention & control , Epitopes/immunology , Humans , Influenza A Virus, H7N9 Subtype/genetics , Influenza, Human/epidemiology , Influenza, Human/immunology , Influenza, Human/virology , Mutation , SoftwareABSTRACT
The Zika virus infections have reached epidemic proportions in the Latin American countries causing severe birth defects and neurological disorders. While several organizations have begun research into design of prophylactic vaccines and therapeutic drugs, computer assisted methods with adequate data resources can be expected to assist in these measures to reduce lead times through bioinformatics approaches. Using 60 sequences of the Zika virus envelope protein available in the GenBank database, our analysis with numerical characterization techniques and several web based bioinformatics servers identified four peptide stretches on the Zika virus envelope protein that are well conserved and surface exposed and are predicted to have reasonable epitope binding efficiency. These peptides can be expected to form the basis for a nascent peptide vaccine which, enhanced by incorporation of suitable adjuvants, can elicit immune response against the Zika virus infections.
Subject(s)
Computational Biology , Drug Design , Peptides/immunology , Viral Vaccines/immunology , Zika Virus Infection/immunology , Zika Virus/immunology , Adjuvants, Immunologic , Epitopes/chemistry , Epitopes/immunology , Peptides/chemical synthesis , Peptides/chemistry , Viral Envelope Proteins/chemistry , Viral Envelope Proteins/genetics , Viral Envelope Proteins/immunology , Viral Vaccines/chemical synthesis , Viral Vaccines/chemistry , Zika Virus/chemistry , Zika Virus/genetics , Zika Virus Infection/prevention & control , Zika Virus Infection/virologyABSTRACT
Human papillomavirus (HPV) occurs in many types, some of which cause cervical, genital, and other cancers. While vaccination is available against the major cancer-causing HPV types, many others are not covered by these preventive measures. Herein, we present a bioinformatics study for the designing of multivalent peptide vaccines against multiple HPV types as an alternative strategy to the virus-like particle vaccines being used now. Our technique of rational design of peptide vaccines is expected to ensure stability of the vaccine against many cycles of mutational changes, elicit immune response, and negate autoimmune possibilities. Using the L1 capsid protein sequences, we identified several peptides for potential vaccine design for HPV 16, 18, 33, 35, 45, and 11 types. Although there are concerns about the epitope-binding affinities for the peptides identified in this process, the technique indicates possibilities of multivalent, adjuvanted, peptide vaccines against a wider range of HPV types, and tailor-made different combinations of the peptides to address frequency variations of types over different population groups as required for prophylaxis and at lower cost than are in use at the present time.
ABSTRACT
A large number of alignment-free techniques of graphical representation and numerical characterization (GRANCH) of bio-molecular sequences have been proposed in the recent past years, but the relative efficacy of these methods in determining the degree of similarities and dissimilarities of such sequences have not been ascertained. OBJECTIVE: Our objective is to make an assessment of the relative efficacy of these methods in determining the degree of similarities and dissimilarities of bio-molecular sequences. METHOD: We have chosen 7 published/communicated methods that represent various classes of GRANCH techniques and computed the descriptors that are expected to characterize similarities and dissimilarities in several sets of gene sequences. We critically appraise the different methods and determine which of these yield non-redundant structural information that could be used to compute different properties of the sequences, and which are correlated enough to one another so that using the simplest representative of the group would suffice. We also do a principal component analysis (PCA) to determine how the variances in the calculated sequence descriptors are explained by the computed principal components (PCs). RESULTS: We found that some of the descriptors are strongly correlated implying a commonality of structural information encoded by them while others are distinctly separate. The PCA results show that the first three PC's explain >97% of the variances. CONCLUSION: We found that some mathematical DNA descriptors calculated by a few of these techniques correlate strongly with one another implying a redundancy in the structural information quantified by those descriptors; others are not strongly correlated with one another suggesting that they encode non-redundant sequence information. From this and our PCA results, our recommendation would be to use minimally correlated set of descriptors or orthogonal descriptors like PCs derived from the descriptor set for the characterization of nucleic acid structure and function.
Subject(s)
DNA/genetics , RNA/genetics , Animals , Base Sequence , DNA/chemistry , Data Display , Exons , Humans , Principal Component Analysis , RNA/chemistry , Statistics as Topic , beta-Globins/geneticsABSTRACT
The growing incidences of new viral diseases and increasingly frequent viral epidemics have strained therapeutic and preventive measures; the high mutability of viral genes puts additional strains on developmental efforts. Given the high cost and time requirements for new drugs development, vaccines remain as a viable alternative, but there too traditional techniques of live-attenuated or inactivated vaccines have the danger of allergenic reactions and others. Peptide vaccines have, over the last several years, begun to be looked on as more appropriate alternatives, which are economically affordable, require less time for development and hold the promise of multi-valent dosages. The developments in bioinformatics, proteomics, immunogenomics, structural biology and other sciences have spurred the growth of vaccinomics where computer assisted approaches serve to identify suitable peptide targets for eventual development of vaccines. In this mini-review we give a brief overview of some of the recent trends in computer assisted vaccine development with emphasis on the primary selection procedures of probable peptide candidates for vaccine development.
Subject(s)
Drug Discovery/methods , Molecular Docking Simulation/methods , Quantitative Structure-Activity Relationship , Vaccines, Subunit/chemistry , Vaccines, Subunit/immunologyABSTRACT
BACKGROUND: The influenza genome is highly variable due primarily to two mechanisms: antigenic drift and antigenic shift. A third mechanism for genetic change, known as copy choice or template switching, can arise during replication when, if two viral strains infect a cell, a part of a gene from the second viral strain can be copied into the growing progeny of a gene of the first viral strain as replacement leading to a new variety of the virus. This template switching between the same genes of the two strains is known as homologous recombination. While genetic drift and shift are well-understood, the presence or absence of intra-segment homologous recombination in influenza genomes is controversial. CONTEXT AND PURPOSE OF STUDY: We are interested to study the possibility of subunit-wise homologous recombination. The idea is that where well-defined subunits are separated by consensus sequences, it might be possible for template switching to take place at such junctions. The influenza hemagglutinin gene has basically two subunits, HA1 and HA2, with HA1 being mostly surface exposed and containing the active site for binding to cells, while HA2 secures the hemagglutinin to the viral coat. We undertook a thorough search of the major human infecting influenza hemagglutinin gene sequences, viz., the H1N1, H5N1, H3N2 and H7N9 subtypes, over the period 2010-2014 in Asia to determine if certain sequences could be identified that had HA1 from a previous strain and HA2 from another. RESULTS: Our search yielded several instances where sequence identities between segments of various strains could be interpreted as indicating possibilities of segment exchange. In some cases, on closer examination they turn out to differ by a few mutations in each segment, due perhaps to the short time span of our database. CONCLUSIONS AND POTENTIAL IMPLICATIONS: The study reported here, and in combination with our earlier observations on the neuraminidase, shows that subunit-wise recombination-like events in the influenza genes may be occurring more often than have been accounted for and merits further detailed studies.
Subject(s)
Computational Biology/statistics & numerical data , Hemagglutinin Glycoproteins, Influenza Virus/genetics , Homologous Recombination , Influenza A virus/genetics , Sequence Analysis, DNA/statistics & numerical data , Algorithms , Animals , Computational Biology/methods , Evolution, Molecular , Genes, Viral/genetics , Genome, Viral/genetics , Humans , Influenza A virus/classification , Influenza A virus/metabolism , Models, Genetic , Sequence Analysis, DNA/methodsABSTRACT
BACKGROUND: The recent epidemic of Zika virus infections in South and Latin America have raised serious concern on its ramifications for the population in the Americas and spread of the virus worldwide. The Zika virus disease is a relatively new phenomenon for which sufficient and comprehensive data and investigative reports have not been available to date. OBJECTIVE: To carry out a bioinformatics study of the available Zika virus genomic sequences to characterize the virus. METHOD: 2D graphical representation method is used for visual rendering and compute sequence parameters and descriptors of the African and Asian-American groups of the Zika viruses to characterize the sequences. We also used MEGA5.2 and other software to compute various biological properties of interest like phylogenetic relationships, transition-transversion ratios, amino acid usage, codon usage bias and hydropathy index of the Zika genomes and virions. RESULTS: The phylogenetic relationships show that the African and Asian-American Zika virus genomes are grouped in two clades. The 2D plots of typical genomes of these types also show dramatic differences indicating that the gene sequences at the 5'-end coding regions for the structural proteins are rather strongly conserved. Among other characteristics, the transition/transversion ratio matrices for the sequences in each of the two clades show that analogous to the dengue virus, the transition rates are about 10 to 15 times the transversion rates. CONCLUSION: These findings are important for computer-assisted approaches towards surveillance of emerging Zika virus strains as well as in the design of drugs and vaccines to combat the growth and spread of the Zika virus.