Detecting species-site dependencies in large multiple sequence alignments

Hotspots for mutations in the SARS-CoV-2 spike glycoprotein: a correspondence analysis

Spike glycoprotein (Sgp) is liable for binding of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to the host receptors. Since Sgp is the main target for vaccine and drug designing, elucidating its mutation pattern could help in this regard. This study is aimed at investigating the correspondence of specific residues to the SgpSARS-CoV-2 functionality by explorative interpretation of sequence alignments. Centrality analysis of the Sgp dissects the importance of these residues in the interaction network of the RBD-ACE2 (receptor-binding domain) complex and furin cleavage site. Correspondence of RBD to threonine500 and asparagine501 and furin cleavage site to glutamine675, glutamine677, threonine678, and alanine684 was observed; all residues are exactly located at the interaction interfaces. The harmonious location of residues dictates the RBD binding property and the flexibility, hydrophobicity, and accessibility of the furin cleavage site. These species-specific residues can be assumed as real targets of evolution, while other substitutions tend to support them. Moreover, all these residues are parts of experimentally identified epitopes. Therefore, their substitution may affect vaccine efficacy. Higher rate of RBD maintenance than furin cleavage site was predicted. The accumulation of substitutions reinforces the probability of the multi-host circulation of the virus and emphasizes the enduring evolutionary events.

Warfarin and vitamin K epoxide reductase: a molecular accounting for observed inhibition

Vitamin K epoxide reductase (VKOR), an endoplasmic reticulum membrane protein, is the key enzyme for vitamin K-dependent carboxylation, a posttranslational modification that is essential for the biological functions of coagulation factors. VKOR is the target of the most widely prescribed oral anticoagulant, warfarin. However, the topological structure of VKOR and the mechanism of warfarin's inhibition of VKOR remain elusive. Additionally, it is not clear why warfarin-resistant VKOR mutations identified in patients significantly decrease warfarin's binding affinity, but have only a minor effect on vitamin K binding. Here, we used immunofluorescence confocal imaging of VKOR in live mammalian cells and PEGylation of VKOR's endogenous cytoplasmic-accessible cysteines in intact microsomes to probe the membrane topology of human VKOR. Our results show that the disputed loop sequence between the first and second transmembrane (TM) domain of VKOR is located in the cytoplasm, supporting a 3-TM topological structure of human VKOR. Using molecular dynamics (MD) simulations, a T-shaped stacking interaction between warfarin and tyrosine residue 139, within the proposed TY₁₃₉A warfarin-binding motif, was observed. Furthermore, a reversible dynamic warfarin-binding pocket opening and conformational changes were observed when warfarin binds to VKOR. Several residues (Y25, A26, and Y139) were found essential for warfarin binding to VKOR by MD simulations, and these were confirmed by the functional study of VKOR and its mutants in their native milieu using a cell-based assay. Our findings provide new insights into the dynamics of the binding of warfarin to VKOR, as well as into warfarin's mechanism of anticoagulation.

ALVIS: interactive non-aggregative visualization and explorative analysis of multiple sequence alignments

Sequence Logos and its variants are the most commonly used method for visualization of multiple sequence alignments (MSAs) and sequence motifs. They provide consensus-based summaries of the sequences in the alignment. Consequently, individual sequences cannot be identified in the visualization and covariant sites are not easily discernible. We recently proposed Sequence Bundles, a motif visualization technique that maintains a one-to-one relationship between sequences and their graphical representation and visualizes covariant sites. We here present Alvis, an open-source platform for the joint explorative analysis of MSAs and phylogenetic trees, employing Sequence Bundles as its main visualization method. Alvis combines the power of the visualization method with an interactive toolkit allowing detection of covariant sites, annotation of trees with synapomorphies and homoplasies, and motif detection. It also offers numerical analysis functionality, such as dimension reduction and classification. Alvis is user-friendly, highly customizable and can export results in publication-quality figures. It is available as a full-featured standalone version (http://www.bitbucket.org/rfs/alvis) and its Sequence Bundles visualization module is further available as a web application (http://science-practice.com/projects/sequence-bundles).

Phylogeny of the Vitamin K 2,3-Epoxide Reductase (VKOR) Family and Evolutionary Relationship to the Disulfide Bond Formation Protein B (DsbB) Family

In humans and other vertebrate animals, vitamin K 2,3-epoxide reductase (VKOR) family enzymes are the gatekeepers between nutritionally acquired K vitamins and the vitamin K cycle responsible for posttranslational modifications that confer biological activity upon vitamin K-dependent proteins with crucial roles in hemostasis, bone development and homeostasis, hormonal carbohydrate regulation and fertility. We report a phylogenetic analysis of the VKOR family that identifies five major clades. Combined phylogenetic and site-specific conservation analyses point to clade-specific similarities and differences in structure and function. We discovered a single-site determinant uniquely identifying VKOR homologs belonging to human pathogenic, obligate intracellular prokaryotes and protists. Building on previous work by Sevier et al. (Protein Science 14:1630), we analyzed structural data from both VKOR and prokaryotic disulfide bond formation protein B (DsbB) families and hypothesize an ancient evolutionary relationship between the two families where one family arose from the other through a gene duplication/deletion event. This has resulted in circular permutation of primary sequence threading through the four-helical bundle protein folds of both families. This is the first report of circular permutation relating distant α-helical membrane protein sequences and folds. In conclusion, we suggest a chronology for the evolution of the five extant VKOR clades.

Sequence Bundles: a novel method for visualising, discovering and exploring sequence motifs

Background

We introduce Sequence Bundles--a novel data visualisation method for representing multiple sequence alignments (MSAs). We identify and address key limitations of the existing bioinformatics data visualisation methods (i.e. the Sequence Logo) by enabling Sequence Bundles to give salient visual expression to sequence motifs and other data features, which would otherwise remain hidden.

Methods

For the development of Sequence Bundles we employed research-led information design methodologies. Sequences are encoded as uninterrupted, semi-opaque lines plotted on a 2-dimensional reconfigurable grid. Each line represents a single sequence. The thickness and opacity of the stack at each residue in each position indicates the level of conservation and the lines' curved paths expose patterns in correlation and functionality. Several MSAs can be visualised in a composite image. The Sequence Bundles method is designed to favour a tangible, continuous and intuitive display of information.

Results

We have developed a software demonstration application for generating a Sequence Bundles visualisation of MSAs provided for the BioVis 2013 redesign contest. A subsequent exploration of the visualised line patterns allowed for the discovery of a number of interesting features in the dataset. Reported features include the extreme conservation of sequences displaying a specific residue and bifurcations of the consensus sequence.

Conclusions

Sequence Bundles is a novel method for visualisation of MSAs and the discovery of sequence motifs. It can aid in generating new insight and hypothesis making. Sequence Bundles is well disposed for future implementation as an interactive visual analytics software, which can complement existing visualisation tools.

Role of factor H binding protein in Neisseria meningitidis virulence and its potential as a vaccine candidate to broadly protect against meningococcal disease

Neisseria meningitidis is a Gram-negative microorganism that exists exclusively in humans and can cause devastating invasive disease. Although capsular polysaccharide-based vaccines against serogroups A, C, Y, and W135 are widely available, the pathway to a broadly protective vaccine against serogroup B has been more complex. The last 11 years has seen the discovery and development of the N. meningitidis serogroup B (MnB) outer membrane protein factor H binding protein (fHBP) as a vaccine component. Since the initial discovery of fHBP, a tremendous amount of work has accumulated on the diversity, structure, and regulation of this important protein. fHBP has proved to be a virulence factor for N. meningitidis and a target for functional bactericidal antibodies. fHBP is critical for survival of meningococci in the human host, as it is responsible for the primary interaction with human factor H (fH). Binding of hfH by the meningococcus serves to downregulate the host alternative complement pathway and helps the organism evade host innate immunity. Preclinical studies have shown that an fHBP-based vaccine can elicit serum bactericidal antibodies capable of killing MnB, and the vaccine has shown very encouraging results in human clinical trials. This report reviews our current knowledge of fHBP. In particular, we discuss the recent advances in our understanding of fHBP, its importance to N. meningitidis, and its potential role as a vaccine for preventing MnB disease.

A multi-country evaluation of Neisseria meningitidis serogroup B factor H-binding proteins and implications for vaccine coverage in different age groups

BACKGROUND

Recombinant vaccines containing factor H-binding protein (fHBP) have been developed for the purpose of protection from invasive meningococcal serogroup B disease. Neisseria meningitidis fHBP sequences can be divided into 2 genetically and immunologically distinct subfamilies (A and B); thus, cross protection is conferred within but not between subfamilies. A comprehensive understanding of fHBP epidemiology is required to accurately assess the potential vaccine impact when considering different vaccination implementation strategies.

METHODS

Systematically collected invasive meningococcal serogroup B isolates from England, Wales, Northern Ireland, the United States, Norway, France and the Czech Republic were previously characterized for fHBP sequence. This study expanded the evaluation with additional meningococcal serogroup B disease isolates from Spain (n = 346) and Germany (n = 205). This expanded set (n = 1841), collected over a 6-year period (2001 to 2006), was evaluated for fHBP sequence and fHBP subfamily relative to patient age.

RESULTS

All 1841 isolates contained fhbp. fHBP sequences from Spain and Germany fell within the previously described subfamilies, with 69% of isolates belonging to subfamily B and 31% to subfamily A; prevalent sequence variants were also similar. Stratification of data by age indicated that disease in infants <1 year of age was caused by a significantly higher proportion of isolates with fHBP subfamily A variants than that seen in adolescents and young adults 11-25 years (47.7% versus 19.5%, P < 0.0001, respectively).

CONCLUSIONS

These observations highlight a difference in epidemiology of fHBP subfamilies in different age groups, with fHBP subfamily A strains causing more disease in vulnerable populations, such as infants, than in adolescents.

A hetero-dimer model for concerted action of vitamin K carboxylase and vitamin K reductase in vitamin K cycle

Vitamin K carboxylase (VKC) is believed to convert vitamin K, in the vitamin K cycle, to an alkoxide-epoxide form which then reacts with CO(2) and glutamate to generate γ-carboxyglutamic acid (Gla). Subsequently, vitamin K epoxide reductase (VKOR) is thought to convert the alkoxide-epoxide to a hydroquinone form. By recycling vitamin K, the two integral-membrane proteins, VKC and VKOR, maintain vitamin K levels and sustain the blood coagulation cascade. Unfortunately, NMR or X-ray crystal structures of the two proteins have not been characterized. Thus, our understanding of the vitamin K cycle is only partial at the molecular level. In this study, based on prior biochemical experiments on VKC and VKOR, we propose a hetero-dimeric form of VKC and VKOR that may explain the efficient oxidation and reduction of vitamin K during the vitamin K cycle.

Concept and application of a computational vaccinology workflow

BACKGROUND

The last years have seen a renaissance of the vaccine area, driven by clinical needs in infectious diseases but also chronic diseases such as cancer and autoimmune disorders. Equally important are technological improvements involving nano-scale delivery platforms as well as third generation adjuvants. In parallel immunoinformatics routines have reached essential maturity for supporting central aspects in vaccinology going beyond prediction of antigenic determinants. On this basis computational vaccinology has emerged as a discipline aimed at ab-initio rational vaccine design.Here we present a computational workflow for implementing computational vaccinology covering aspects from vaccine target identification to functional characterization and epitope selection supported by a Systems Biology assessment of central aspects in host-pathogen interaction. We exemplify the procedures for Epstein Barr Virus (EBV), a clinically relevant pathogen causing chronic infection and suspected of triggering malignancies and autoimmune disorders.

RESULTS

We introduce pBone/pView as a computational workflow supporting design and execution of immunoinformatics workflow modules, additionally involving aspects of results visualization, knowledge sharing and re-use. Specific elements of the workflow involve identification of vaccine targets in the realm of a Systems Biology assessment of host-pathogen interaction for identifying functionally relevant targets, as well as various methodologies for delineating B- and T-cell epitopes with particular emphasis on broad coverage of viral isolates as well as MHC alleles.Applying the workflow on EBV specifically proposes sequences from the viral proteins LMP2, EBNA2 and BALF4 as vaccine targets holding specific B- and T-cell epitopes promising broad strain and allele coverage.

CONCLUSION

Based on advancements in the experimental assessment of genomes, transcriptomes and proteomes for both, pathogen and (human) host, the fundaments for rational design of vaccines have been laid out. In parallel, immunoinformatics modules have been designed and successfully applied for supporting specific aspects in vaccine design. Joining these advancements, further complemented by novel vaccine formulation and delivery aspects, have paved the way for implementing computational vaccinology for rational vaccine design tackling presently unmet vaccine challenges.

Visualization of multiple alignments, phylogenies and gene family evolution

Software for visualizing sequence alignments and trees are essential tools for life scientists. In this review, we describe the major features and capabilities of a selection of stand-alone and web-based applications useful when investigating the function and evolution of a gene family. These range from simple viewers, to systems that provide sophisticated editing and analysis functions. We conclude with a discussion of the challenges that these tools now face due to the flood of next generation sequence data and the increasingly complex network of bioinformatics information sources.