Computational Analysis and Binding Site Identification of Type III Secretion System ATPase from Pseudomonas aeruginosa

In silico design of epitope-based peptide vaccine against non-typhoidal Salmonella through immunoinformatic approaches

Non-typhoidal Salmonella (NTS) is one of the leading bacterial causes of many invasive human infections with a high antibiotic resistance profile. With this concern, the current study aimed to design an effective epitope-based peptide vaccine against NTS species as a successive and substitutive protective measure of invasive NTS disease. To design rationally, the current study considered a comprehensive in silico workflow combination of both immunoinformatics and molecular modeling approaches, including molecular docking and molecular dynamics (MD) simulation. We identified the two most promising T cell epitopes KVLYGIFAI and YGIFAITAL, and three B cell epitopes AAPVQVGEAAGS, TGGGDGSNT, and TGGGDGSNTGTTTT, in the outer membrane of NTS. Using these epitopes, a multiepitope vaccine was subsequently constructed along with appropriate adjuvant and linkers, which showed a good binding affinity and stability with toll-like receptor 2 (TLR2) in both molecular docking and MD simulation. Furthermore, in silico immune simulation described a strong immune response with a high number of antibodies, interferon-γ, and activated B and T cells. This study collectively suggests that predicted vaccine constructs could be considered potential vaccine candidates against common NTS species.Communicated by Ramaswamy H. Sarma.

Structure and dynamics of membrane protein in SARS-CoV-2

SARS-CoV-2 membrane (M) protein performs a variety of critical functions in virus infection cycle. However, the expression and purification of membrane protein structure is difficult despite tremendous progress. In this study, the 3 D structure is modeled followed by intensive validation and molecular dynamics simulation. The lack of suitable homologous templates (>30% sequence identities) leads us to construct the membrane protein models using template-free modeling (de novo or ab initio) approach with Robetta and trRosetta servers. Comparing with other model structures, it is evident that trRosetta (TM-score: 0.64; TM region RMSD: 2 Å) can provide the best model than Robetta (TM-score: 0.61; TM region RMSD: 3.3 Å) and I-TASSER (TM-score: 0.45; TM region RMSD: 6.5 Å). 100 ns molecular dynamics simulations are performed on the model structures by incorporating membrane environment. Moreover, secondary structure elements and principal component analysis (PCA) have also been performed on MD simulation data. Finally, trRosetta model is utilized for interpretation and visualization of interacting residues during protein-protein interactions. The common interacting residues including Phe103, Arg107, Met109, Trp110, Arg131, and Glu135 in the C-terminal domain of M protein are identified in membrane-spike and membrane-nucleocapsid protein complexes. The active site residues are also predicted for potential drug and peptide binding. Overall, this study might be helpful to design drugs and peptides against the modeled membrane protein of SARS-CoV-2 to accelerate further investigation. Communicated by Ramaswamy H. Sarma.

Dynamic insights into the effects of nonsynonymous polymorphisms (nsSNPs) on loss of TREM2 function

Single nucleotide variations in Triggering Receptor Expressed on Myeloid Cells 2 (TREM2) are associated with many neurodegenerative diseases, including Nasu-Hakola disease (NHD), frontotemporal dementia (FTD), and late-onset Alzheimer's disease because they disrupt ligand binding to the extracellular domain of TREM2. However, the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) in TREM2 on disease progression remain unknown. In this study, we identified several high-risk nsSNPs in the TREM2 gene using various deleterious SNP predicting algorithms and analyzed their destabilizing effects on the ligand recognizing region of the TREM2 immunoglobulin (Ig) domain by molecular dynamics (MD) simulation. Cumulative prediction by all tools employed suggested the three most deleterious nsSNPs involved in loss of TREM2 function are rs549402254 (W50S), rs749358844 (R52C), and rs1409131974 (D104G). MD simulation showed that these three variants cause substantial structural alterations and conformational remodeling of the apical loops of the TREM2 Ig domain, which is responsible for ligand recognition. Detailed analysis revealed that these variants substantially increased distances between apical loops and induced conformation remodeling by changing inter-loop nonbonded contacts. Moreover, all nsSNPs changed the electrostatic potentials near the putative ligand-interacting region (PLIR), which suggested they might reduce specificity or loss of binding affinity for TREM2 ligands. Overall, this study identifies three potential high-risk nsSNPs in the TREM2 gene. We propose further studies on the molecular mechanisms responsible for loss of TREM2 function and the associations between TREM2 nsSNPs and neurodegenerative diseases.

Structural Consequence of Non-Synonymous Single-Nucleotide Variants in the N-Terminal Domain of LIS1

Disruptive neuronal migration during early brain development causes severe brain malformation. Characterized by mislocalization of cortical neurons, this condition is a result of the loss of function of migration regulating genes. One known neuronal migration disorder is lissencephaly (LIS), which is caused by deletions or mutations of the LIS1 (PAFAH1B1) gene that has been implicated in regulating the microtubule motor protein cytoplasmic dynein. Although this class of diseases has recently received considerable attention, the roles of non-synonymous polymorphisms (nsSNPs) in LIS1 on lissencephaly progression remain elusive. Therefore, the present study employed combined bioinformatics and molecular modeling approach to identify potential damaging nsSNPs in the LIS1 gene and provide atomic insight into their roles in LIS1 loss of function. Using this approach, we identified three high-risk nsSNPs, including rs121434486 (F31S), rs587784254 (W55R), and rs757993270 (W55L) in the LIS1 gene, which are located on the N-terminal domain of LIS1. Molecular dynamics simulation highlighted that all variants decreased helical conformation, increased the intermonomeric distance, and thus disrupted intermonomeric contacts in the LIS1 dimer. Furthermore, the presence of variants also caused a loss of positive electrostatic potential and reduced dimer binding potential. Since self-dimerization is an essential aspect of LIS1 to recruit interacting partners, thus these variants are associated with the loss of LIS1 functions. As a corollary, these findings may further provide critical insights on the roles of LIS1 variants in brain malformation.

Molecular docking and dynamics of Nickel-Schiff base complexes for inhibiting β-lactamase of Mycobacterium tuberculosis

In recent years, multidrug-resistance has become a primary concern in the treatment and management of tuberculosis, an infectious disease caused by Mycobacterium tuberculosis. In this context, searching new anti-tuberculosis agents particularly targeting the β-lactamase (BlaC) is reported to be promising as this enzyme is one of the key player in the development of multidrug resistance. This study reports the design of some Nickel (Ni) based tetradentate N₂O₂ Schiff bases, employing density functional theory. All analogs are optimized at B3LYP/SDD level of theory. Dipole moment, electronic energy, enthalpy, Gibbs free energy, HOMO-LUMO gap, and softness of these modified drugs are also investigated. Molecular interactions between designed ligands and BlaC have been analyzed by molecular docking approach, followed by molecular dynamics (MD) simulation. All designed compounds show low HOMO-LUMO gap, while addition of halogen increases the dipole moment of the compounds. Docking and MD simulation investigations reveal that the designed compounds are more potent than standard inhibitor, where Ile117, Pro290, Arg236 and Thr253 residues of BlaC are found to play important role in the ligand binding. Through MD simulation study, the best binding compound is also observed to form stable complex by increasing the protein rigidness. The ADME/T analysis suggests that modified drugs are less toxic and shows an improved pharmacokinetic properties than that of the standard drug. These results further confirm the ability of Ni-directed Schiff bases to bind simultaneously to the active site of BlaC and support them as potential candidates for the future treatment of tuberculosis disease.

In silico-based vaccine design against Ebola virus glycoprotein

Ebola virus (EBOV) is one of the lethal viruses, causing more than 24 epidemic outbreaks to date. Despite having available molecular knowledge of this virus, no definite vaccine or other remedial agents have been developed yet for the management and avoidance of EBOV infections in humans. Disclosing this, the present study described an epitope-based peptide vaccine against EBOV, using a combination of B-cell and T-cell epitope predictions, followed by molecular docking and molecular dynamics simulation approach. Here, protein sequences of all glycoproteins of EBOV were collected and examined via in silico methods to determine the most immunogenic protein. From the identified antigenic protein, the peptide region ranging from 186 to 220 and the sequence HKEGAFFLY from the positions of 154-162 were considered the most potential B-cell and T-cell epitopes, correspondingly. Moreover, this peptide (HKEGAFFLY) interacted with HLA-A*32:15 with the highest binding energy and stability, and also a good conservancy of 83.85% with maximum population coverage. The results imply that the designed epitopes could manifest vigorous enduring defensive immunity against EBOV.