Coevolution based immunoinformatics approach considering variability of epitopes to combat different strains: A case study using spike protein of SARS-CoV-2.

In the recent past several vaccines were developed to combat the COVID-19 disease. Unfortunately, the protective efficacy of the current vaccines has been reduced due to the high mutation rate in SARS-CoV-2. Here, we successfully implemented a coevolution based immunoinformatics approach to design an epitope-based peptide vaccine considering variability in spike protein of SARS-CoV-2. The spike glycoprotein was investigated for B- and T-cell epitope prediction. Identified T-cell epitopes were mapped on previously reported coevolving amino acids in the spike protein to introduce mutation. The non-mutated and mutated vaccine components were constructed by selecting epitopes showing overlapping with the predicted B-cell epitopes and highest antigenicity. Selected epitopes were linked with the help of a linker to construct a single vaccine component. Non-mutated and mutated vaccine component sequences were modelled and validated. The in-silico expression level of the vaccine constructs (non-mutated and mutated) in E. coli K12 shows promising results. The molecular docking analysis of vaccine components with toll-like receptor 5 (TLR5) demonstrated strong binding affinity. The time series calculations including root mean square deviation (RMSD), radius of gyration (RGYR), and energy of the system over 100 ns trajectory obtained from all atom molecular dynamics simulation showed stability of the system. The combined coevolutionary and immunoinformatics approach used in this study will certainly help to design an effective peptide vaccine that may work against different strains of SARS-CoV-2. Moreover, the strategy used in this study can be implemented on other pathogens.

A strategy to optimize the peptide-based inhibitors against different mutants of the spike protein of SARS-CoV-2.

The SARS-CoV-2 virus has caused high-priority health concerns at a global level. Vaccines have stalled the proliferation of viruses to some extent. Yet, the emergence of newer, potentially more infectious, and dangerous mutants such as Delta and Omicron are among the major challenges in finding a more permanent solution for this pandemic. The effectiveness of antivirals Molnupiravir and Paxlovid, authorized for emergency use by the FDA, are yet to be assessed on a larger population. Patients with a high risk of disease progression have received treatment with antibody-cocktail. Most of the mutations leading to the new lineage of SARS-CoV-2 are found in the spike protein of this virus that plays a key role in facilitating host entry. The current study has investigated how to modify a promising peptide-based inhibitor of spike protein, LCB3, against common mutations, N501Y and K417N in the target protein so that it retains its efficacy against the spike protein. LCB3 being a prototype for protein-based inhibitors is an ideal testing system to learn about protein-based inhibitors. This study proposes the substitutions of amino acid residues of LCB3 inhibitor using a structure-based approach that considers free energy decomposition of residues, the distance between atoms, and the interaction among amino acids. The binding free energy calculations suggest a possible improvement in the binding affinity of existing inhibitor LCB3 to the mutant forms of the S-protein using simple substitutions at specific positions of the inhibitor. This approach, being general, can be used in different inhibitors and other mutations and help in fighting against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Coevolutionary forces shaping the fitness of SARS-CoV-2 spike glycoprotein against human receptor ACE2.

The current global health problem caused by SARS-CoV-2 has challenged the scientific community in various ways. Therefore, worldwide several scientific groups are exploring SARS-CoV-2 from different aspects including its origin, spread, severe infectivity, and also to find a cure. It is now well known that spike glycoprotein helps SARS-CoV-2 to enter inside the human host through a cellular receptor ACE2. However, the role of coevolutionary forces that makes SARS-CoV-2 spike glycoprotein more fit towards its human host remains unexplored. Therefore, in present bioinformatics study we identify coevolving amino acids in spike glycoprotein. Additionally, the effects of coevolution on the stability of the spike glycoprotein as well as its binding with receptor ACE2 were predicted. The results clearly indicate that coevolutionary forces play a pivotal role in increasing the fitness of spike glycoprotein against ACE2.

Physical interaction between nuclear accumulated CC-NB-ARC-LRR protein and WRKY64 promotes EDS1 dependent Fusarium wilt resistance in chickpea.

Fusarium wilt is one of the most serious diseases affecting chickpea (Cicer arietinum L.). Here, we identified a putative Resistance Gene Analog (CaRGA) from chickpea, encoding a coiled-coil (CC) nucleotide-binding oligomerization domain (NB-ARC) containing leucine-rich repeat (LRR) protein (CC-NLR protein) that confers resistance against Fusarium oxysporum f. sp. ciceri race1 (Foc1). Over-expression and silencing of CaRGA in chickpea resulted in enhanced resistance and hyper-susceptibility, respectively against Foc1. Furthermore, defense response to Foc1 depends on CC-NLR interaction with WRKY64 transcription factor. CaRGA mediated wilt resistance largely compromised when WRKY64 was silenced. We also determined in planta intramolecular interactions and self-association of chickpea CC-NLR protein. The study shows CC domain suppressing auto-activation of the full-length CC-NLR protein in the absence of pathogen through self-inhibitory intramolecular interaction with NB-ARC domain, which is attenuated by self-interactions to LRR domain. Chickpea CC-NLR protein forms homocomplexes and then interacts with WRKY64. CC-NLR protein further phosphorylates WRKY64 thereby, ubiquitination and proteasome mediated degradation are protected. Phosphorylated WRKY64 with increased stability binds to EDS1 promoter and stimulates its transcription that induces in planta ectopic cell-death. The detailed analysis of CC-NLR and WRKY interactions provide a better understanding of the immune regulation by NLR proteins under biotic stresses.

The long unstructured region of Bcl-xl modulates its structural dynamics.

Bcl-xl protein has a long unstructured loop attached to its structured region which joins two helices. The necessity to have this unstructured segment in Bcl-xl is not yet well understood. To what extent the unstructured segment can influence the dynamics of the structured region of protein, with potential to influence the function, has been investigated in this work. Molecular dynamics simulation and principal component analysis show how the loop affects the internal motions of the protein, particularly its ligand binding pocket. Generally an unstructured region in the structure would promote flexibility resulting entropic stability but in contrary, here it narrows down the conformational space of the structured region of protein that could be hypothesized to impact the functional precision. Effects of the loop propagate to the binding pocket through structural rearrangements of polar side chains. The immediate suspicion of possible impact of phosphorylation to modulate the function of the protein is proven to be a fact, as the phosphorylated S49 and S62 located on the large unstructured region are seen to perturb the electrostatic network of the structure; an observation that validates and clarifies the role of loop as a modulator through biophysical and biochemical mechanisms. Proteins 2017; 85:1567-1579. © 2017 Wiley Periodicals, Inc.

Interactions between Bcl-xl and its inhibitors: Insights into ligand design from molecular dynamics simulation.

The Bcl-xl protein is a potential drug target for cancer, and it has a relatively flat and flexible binding pocket. ABT263 is one of the most promising molecules that inhibit Bcl-xl, and it was developed from its precursor ABT737 with suitable substitutions. However, the structural and mechanistic implications of those changes have not yet been reported. Molecular dynamics simulation has revealed that the conformational microstates of the complex of Bcl-xl and ABT263 shows heterogeneity at the binding interface with Bcl-xl in contrast to the precise interactions witnessed in case of ABT737. This occurs because not all the functional groups of ABT263 are able to anchor into the binding pocket simultaneously at the time of complexation; leaving at least one group weakly associated every time. The insight into the mechanism shows that, in spite of such mutual exclusivity, the resultant effect becomes beneficial, i.e., becomes more effective than ABT737. Going against the traditional belief, the calculations also confirm that there is no benefit of reshaping the highly flexible binding pocket to allow the ligand to be comfortably accommodated and avoid conflicting orientations of the functional groups, as the destabilization becomes active from other sources. These structural clues and in-silico tests suggest possible avenues for improving the binding affinity of ABT263 through further in-vitro and in-vivo tests.

Adaptability in protein structures: structural dynamics and implications in ligand design.

The basic framework of understanding the mechanisms of protein functions is achieved from the knowledge of their structures which can model the molecular recognition. Recent advancement in the structural biology has revealed that in spite of the availability of the structural data, it is nontrivial to predict the mechanism of the molecular recognition which progresses via situation-dependent structural adaptation. The mutual selectivity of protein-protein and protein-ligand interactions often depends on the modulations of conformations empowered by their inherent flexibility, which in turn regulates the function. The mechanism of a protein's function, which used to be explained by the ideas of 'lock and key' has evolved today as the concept of 'induced fit' as well as the 'population shift' models. It is felt that the 'dynamics' is an essential feature to take into account for understanding the mechanism of protein's function. The design principles of therapeutic molecules suffer from the problems of plasticity of the receptors whose binding conformations are accurately not predictable from the prior knowledge of a template structure. On the other hand, flexibility of the receptors provides the opportunity to improve the binding affinity of a ligand by suitable substitution that will maximize the binding by modulating the receptors surface. In this paper, we discuss with example how the protein's flexibility is correlated with its functions in various systems, revealing the importance of its understanding and for making applications. We also highlight the methodological challenges to investigate it computationally and to account for the flexible nature of the molecules in drug design.

New insight to structure-function relationship of GalNAc mediated primary interaction between insecticidal Cry1Ac toxin and HaALP receptor of Helicoverpa armigera.

Over the last few decades Cry1Ac toxin has been widely used in controlling the insect attack due to its high specificity towards target insects. The pore-forming toxin undergoes a complex mechanism in the insect midgut involving sequential interaction with specific glycosylated receptors in which terminal GalNAc molecule plays a vital role. Recent studies on Cry toxins interactions with specific receptors revealed the importance of several amino acid residues in domain III of Cry1Ac, namely Q509, N510, R511, Y513 and W545, serve as potential binding sites that surround the putative GalNAc binding pocket and mediate the toxin-receptor interaction. In the present study, alanine substitution mutations were generated in the Cry1Ac domain III region and functional significance of those key residues was monitored by insect bioassay on Helicoverpa armigera larvae. In addition, ligand blot analysis and SPR binding assay was performed to monitor the binding characteristics of Cry1Ac wild type and mutant toxins towards HaALP receptor isolated from Helicoverpa armigera. Mutagenesis data revealed that, alanine substitutions in R511, Y513 and W545 substantially impacted the relative affinity towards HaALP receptor and toxicity toward target insect. Furthermore, in silico study of GalNAc-mediated interaction also confirmed the important roles of these residues. This structural analysis will provide a detail insight for evaluating and engineering new generation Cry toxins to address the problem of change in insect behavioral patterns.