Galaxy InteractoMIX: An Integrated Computational Platform for the Study of Protein-Protein Interaction Data

Designing a multi-epitope subunit vaccine against Orf virus using molecular docking and molecular dynamics

Orf virus (ORFV) is an acute contact, epitheliotropic, zoonotic, and double-stranded DNA virus that causes significant economic losses in the livestock industry. The objective of this study is to design an immunoinformatics-based multi-epitope subunit vaccine against ORFV. Various immunodominant cytotoxic T lymphocytes (CTL), helper T lymphocytes (HTL), and B-cell epitopes from the B2L, F1L, and 080 protein of ORFV were selected and linked by short connectors to construct a multi-epitope subunit vaccine. Immunogenicity was enhanced by adding an adjuvant β-defensin to the N-terminal of the vaccine using the EAAAK linker. The vaccine exhibited a significant degree of antigenicity and solubility, without allergenicity or toxicity. The 3D formation of the vaccine was subsequently anticipated, improved, and verified. The optimized model exhibited a lower Z-score of -4.33, indicating higher quality. Molecular docking results demonstrated that the vaccine strongly binds to TLR2 and TLR4. Molecular dynamics results indicated that the docked vaccine-TLR complexes were stable. Immune simulation analyses further confirmed that the vaccine can induce a marked increase in IgG and IgM antibody titers, and elevated levels of IFN-γ and IL-2. Finally, the optimized DNA sequence of the vaccine was cloned into the vector pET28a (+) for high expression in the E.coli expression system. Overall, the designed multi-epitope subunit vaccine is highly stable and can induce robust humoral and cellular immunity, making it a promising vaccine candidate against ORFV.

Interaction Between SARS-CoV-2 Spike Protein S1 Subunit and Oyster Heat Shock Protein 70

There is growing evidence that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contaminates the marine environment and is bioaccumulated in filter-feeding shellfish. Previous study shows the Pacific oyster tissues can bioaccumulate the SARS-CoV-2, and the oyster heat shock protein 70 (oHSP70) may play as the primary attachment receptor to bind SARS-CoV-2's recombinant spike protein S1 subunit (rS1). However, detailed information about the interaction between rS1 and oHSP70 is still unknown. In this study, we confirmed that the affinity of recombinant oHSP70 (roHSP70) for rS1 (KD = 20.4 nM) is comparable to the receptor-binding affinity of rACE2 for rS1 (KD = 16.7 nM) by surface plasmon resonance (SPR)-based Biacore and further validated by enzyme-linked immunosorbent assay (ELISA). Three truncated proteins (roHSP70-N/C/M) and five mutated proteins (p.I229del, p.D457del, p.V491_K495del, p.K556I, and p.ΣroHSP70) were constructed according to the molecular docking results. All three truncated proteins have significantly lower affinity for rS1 than the full-length roHSP70, indicating that all three segments of roHSP70 are involved in binding to rS1. Further, the results of SPR and ELISA showed that all five mutant proteins had significantly lower affinity for rS1 than roHSP70, suggesting that amino acids at these sites are involved in binding to rS1. This study provides a preliminary theoretical basis for the bioaccumulation of SARS-CoV-2 in oyster tissues or using roHSP70 as the capture unit to selectively enrich virus particles for detection.

Structure-based learning to predict and model protein-DNA interactions and transcription-factor co-operativity in cis-regulatory elements

Transcription factor (TF) binding is a key component of genomic regulation. There are numerous high-throughput experimental methods to characterize TF-DNA binding specificities. Their application, however, is both laborious and expensive, which makes profiling all TFs challenging. For instance, the binding preferences of ∼25% human TFs remain unknown; they neither have been determined experimentally nor inferred computationally. We introduce a structure-based learning approach to predict the binding preferences of TFs and the automated modelling of TF regulatory complexes. We show the advantage of using our approach over the classical nearest-neighbor prediction in the limits of remote homology. Starting from a TF sequence or structure, we predict binding preferences in the form of motifs that are then used to scan a DNA sequence for occurrences. The best matches are either profiled with a binding score or collected for their subsequent modeling into a higher-order regulatory complex with DNA. Co-operativity is modelled by: (i) the co-localization of TFs and (ii) the structural modeling of protein-protein interactions between TFs and with co-factors. We have applied our approach to automatically model the interferon-β enhanceosome and the pioneering complexes of OCT4, SOX2 (or SOX11) and KLF4 with a nucleosome, which are compared with the experimentally known structures.

SBILib: a handle for protein modeling and engineering

SUMMARY

The SBILib Python library provides an integrated platform for the analysis of macromolecular structures and interactions. It combines simple 3D file parsing and workup methods with more advanced analytical tools. SBILib includes modules for macromolecular interactions, loops, super-secondary structures, and biological sequences, as well as wrappers for external tools with which to integrate their results and facilitate the comparative analysis of protein structures and their complexes. The library can handle macromolecular complexes formed by proteins and/or nucleic acid molecules (i.e. DNA and RNA). It is uniquely capable of parsing and calculating protein super-secondary structure and loop geometry. We have compiled a list of example scenarios which SBILib may be applied to and provided access to these within the library.

AVAILABILITY AND IMPLEMENTATION

SBILib is made available on Github at https://github.com/structuralbioinformatics/SBILib.

CM2D3: Furnishing the human interactome with structural models of protein complexes derived by comparative modeling and docking

The human interactome is composed of around half a million interactions according to recent estimations and it is only for a small fraction of those that three-dimensional structural information is available. Indeed, the structural coverage of the human interactome is very low and given the complexity and time-consuming requirements of solving protein structures this problem will remain for the foreseeable future. Structural models, or predictions, of protein complexes can provide valuable information when the experimentally determined 3D structures are not available. Here we present CM2D3, a relational database containing structural models of the whole human interactome derived both from comparative modeling and data-driven docking. Starting from a consensus interactome derived from integrating several interactomics databases, a strategy was devised to derive structural models by computational means. Currently, CM2D3 includes 33338 structural models of which 5121 derived from comparative modeling and the remaining from docking. Of the latter, the structures of 14554 complexes were derived from monomers modeled by M4T while the rest were modeled with structures as predicted by AlphaFold2. Lastly, CM2D3 complements existing resources by focusing on models derived from both free-docking, as opposed to template-based docking, and hence expanding the available structural information on protein complexes to the scientific community. Database URL:http://www.bioinsilico.org/CM2D3.

Reverse vaccinology approach to design a vaccine targeting membrane lipoproteins of Salmonella typhi

Typhoid fever caused by Salmonella is one of the major health issues worldwide, resulting in millions of cases and has very high rates of morbidities. The therapeutic approaches need to be updated for the effective elimination of the bacterial pathogen. The designing of the multiepitope vaccine against Salmonella using comparative proteomics and reverse vaccinology has covered up all the epitopes that induce sufficient immune responses in the host body. Out of the 4293 proteins, 15 outer membrane proteins have been selected based on their antigenicity, low transmembrane helix (<1), and virulence-associated factors. With the help of the reverse vaccinology approach, the epitopes of MHC Class I, Class II, and B-cell with antigenic, low toxicity, and that have the potential to generate immunogenic response have been identified. Based on the comparative analysis of all the epitopes, a multiepitope-based construct has been designed. Based on physicochemical properties and docking scores for HLA and TLR4, the VC5 construct has been selected, and the molecular dynamic simulation studies have confirmed their interaction. The dissociation constant of the VC5 and TLR4 was found to be 3.1 x 10^-9. Different immune cell activation has been analyzed, representing the potentiality of the VC5 construct as an effective vaccine target. In silico cloning of VC5 in pET28a has also been performed, which requires experimental validation. Therefore, the present study designs a multi-epitope vaccine VC5 targeted to the membrane lipoproteins of Salmonella typhi.Communicated by Ramaswamy H. Sarma.

Theoretical 3D Modeling of NLRP3 Inflammasome Complex

The NOD-like receptor pyrin domain containing 3 (NLRP3) is a multidomain protein that plays a key role in innate immune response. Structures of NLRP3 in different conformational states and bound to cognate partners are available. In this chapter we present an approach to model the oligomeric structure of NLRP3 by homology modeling using multiple templates, symmetry, and refinement. The overall process presented here represents advanced exercise in structural modeling that provides unique insights into the biological role and activation of NLRP3 oligomer. Finally, the same approach can be easily adapted to the rest of the members of the NLRP family.

Construction and contextualization approaches for protein-protein interaction networks

Protein-protein interaction network (PPIN) analysis is a widely used method to study the contextual role of proteins of interest, to predict novel disease genes, disease or functional modules, and to identify novel drug targets. PPIN-based analysis uses both generic and context-specific networks. Multiple contextualization methodologies have been described, such as shortest-path algorithms, neighborhood-based methods, and diffusion/propagation algorithms. This review discusses these methods, provides intuitive representations of PPIN contextualization, and also examines how the quality of such context-specific networks could be improved by considering additional sources of evidence. As a heuristic, we observe that tasks such as identifying disease genes, drug targets, and protein complexes should consider local neighborhoods, while uncovering disease mechanisms and discovering disease-pathways would gain from diffusion-based construction.