Unraveling the Molecular Mechanism of Recognition of Human Interferon-Stimulated Gene Product 15 by Coronavirus Papain-Like Proteases: A Multiscale Simulation Study

Identification of inhibitors for neurodegenerative diseases targeting dual leucine zipper kinase through virtual screening and molecular dynamics simulations

Neurodegenerative diseases lead to a gradual decline in cognitive and motor functions due to the progressive loss of neurons in the central nervous system. The role of dual leucine zipper kinase (DLK) in regulating stress responses and neuronal death pathways highlights its significance as a target against neurodegenerative diseases. The non-availability of FDA-approved drugs emphasizes a need to identify novel DLK-inhibitors. We screened NPAtlas (Natural products) and MedChemExpress (FDA-approved) libraries to identify potent ATP-competitive DLK inhibitors. ADMET analyses identified four compounds (two natural products and two FDA-approved) with favourable features. Subsequently, we performed molecular dynamics simulations to examine the binding-stability and ligand-induced conformational dynamics. Molecular mechanics Poisson Boltzmann surface area (MM-PBSA) calculations demonstrated CID139591660, dithranol, and danthron having greater affinity, while CID156581477 showed lower affinity than control sunitinib. PCA and network analysis results indicated structural and network alteration post-ligand binding. Furthermore, we identified an analogue of CID156581477 using the deep learning-based web server DeLA Drug which demonstrated a higher affinity than its parent compound and the control and identified several crucial interacting residues. Overall, our study provides significant theoretical guidance for designing potent novel DLK inhibitors and compounds that could emerge as promising drug candidates against DLK following laboratory validation.

Computer-aided Affinity Enhancement of a Cross-reactive Antibody against Dengue Virus Envelope Domain III

The dengue virus (DENV), composed of four distinct but serologically related Flaviviruses, causes the most important emerging viral disease, with nearly 400 million infections yearly. Currently, there are no approved therapies. Although DENV infection induces lifelong immunity against the same serotype, the antibodies raised contribute to severe disease in heterotypic infections. Therefore, understanding the mechanism of DENV neutralization by antibodies is crucial in the design of vaccines against all serotypes. This study reports a comparative structural and energetic analysis of the monoclonal antibody (mAb) 4E11 in complex with its target domain III of the envelope protein for all four DENV serotypes. We use extensive replica molecular dynamics simulations in conjunction with the binding free energy calculations. Further single point and double mutations were designed through computational site-directed mutagenesis and observed that the re-engineered antibody exhibits high affinity to binding and broadly neutralizing activity against serotypes. Our results showed improved binding affinity by the gain of enthalpy, which could be attributed to the stabilization of salt-bridge and hydrogen bond interactions at the antigen-antibody interface. The findings provide valuable results in understanding the structural dynamics and energetic contributions that will be helpful to the design of high-affinity antibodies against dengue infections.

Activation-induced cytidine deaminase an antibody diversification enzyme interacts with chromatin modifier UBN1 in B-cells

Activation-induced cytidine deaminase (AID) is the key mediator of antibody diversification in activated B-cells by the process of somatic hypermutation (SHM) and class switch recombination (CSR). Targeting AID to the Ig genes requires transcription (initiation and elongation), enhancers, and its interaction with numerous factors. Furthermore, the HIRA chaperon complex, a regulator of chromatin architecture, is indispensable for SHM. The HIRA chaperon complex consists of UBN1, ASF1a, HIRA, and CABIN1 that deposit H3.3 onto the DNA, the SHM hallmark. We explored whether UBN1 interacts with AID using computational and in-vitro experiments. Interestingly, our in-silico studies, such as molecular docking and molecular dynamics simulation results, predict that AID interacts with UBN1. Subsequently, co-immunoprecipitation and pull-down experiments established interactions between UBN1 and AID inside B-cells. Additionally, a double immunofluorescence assay confirmed that AID and UBN1 were co-localized in the human and chicken B-cell lines. Moreover, proximity ligation assay studies validated that AID interacts with UBN1. Ours is the first report on the interaction of genome mutator enzyme AID with UBN1. Nevertheless, the fate of interaction between UBN1 and AID is yet to be explored in the context of SHM or CSR.

Plant derived active compounds of ayurvedic neurological formulation, Saraswatharishta as a potential dual leucine zipper kinase inhibitor: an in-silico study

Recent findings have highlighted the essential role of dual leucine zipper kinase (DLK) in neuronal degeneration. Saraswatharishta (SWRT), an ayurvedic formulation utilized in traditional Indian medicine, has demonstrated effectiveness in addressing neurodegenerative diseases. Herein, we aim to delve into the atomistic details of the mode of action of phytochemicals present in SWRT against DLK. Our screening process encompassed over 500 distinct phytochemicals derived from the main ingredients of the SWRT formulation. Through a comparative analysis of docking scores and relative poses, we successfully identified four novel compounds, which underwent further investigation via 2 × 500 ns long molecular dynamics (MD) simulations. Among the top four compounds, CID16066851 sourced from the Acorus calamus displayed the most stable complex with DLK. The molecular mechanics Poisson - Boltzmann surface area (MM-PBSA) calculations highlighted the significance of electrostatic and van der Waals interactions in the binding recognition process. Additionally, we identified key residues, namely Phe192, Leu243, Val139, and Leu141, as hotspots that predominantly govern the DLK-inhibitor interaction. Notably, the leading compounds are sourced from the Acorus calamus, Syzygium aromaticum, Zingiber officinale, and Anethum sowa plants present in the SWRT formulation. Overall, the findings of our study hold promise for future drug development endeavors combating neurodegenerative conditions.Communicated by Ramaswamy H. Sarma.

Withania somnifera extract reduces gastric cancerous properties through inhibition of gankyrin in cellular milieu produced by Helicobacter pylori and Epstein Barr virus

Helicobacter pylori and Epstein Barr virus (EBV) are group1 carcinogens and their role in Gastric cancer (GC) is well established. Previously we have shown that H. pylori and EBV appears to support aggressive gastric oncogenesis through the upregulation of oncoprotein Gankyrin. Natural plant active molecules have the potential to interrupt oncogenesis. Herein, we investigated the potential of Withania somnifera root extract (WSE) as a possible chemotherapeutic agent against host oncoprotein Gankyrin whose expression was altered by H. pylori and EBV-associated modified cellular milieu. The results show that WSE does not have any inhibitory effect on H. pylori and EBV-associated gene transcripts except for the lmps (lmp1, lmp2a, and lmp2B). Moreover, the WSE exert their anticancer activity via host cellular response and decreased the expression of cell-migratory (mmp3 and mmp7); cell-cycle regulator (pcna); antiapoptotic gene (bcl2); increased the expression of the proapoptotic gene (apaf1 and bax); and tumor suppressor (p53, prb, and pten). Knockdown of Gankyrin followed by the treatment of WSE also decreases the expression of TNF-ɑ, Akt, and elevated the expression of NFkB, PARP, Casp3, and Casp9. WSE also reduces cell migration, and genomic instability and forced the cells to commit programmed cell death. Moreover, molecular simulation studies revealed that out of eight active compounds of WSE, only four compounds such as withaferin A (WFA), withanoside IV (WA4), withanolide B (WNB), and withanolide D (WND) showed direct stable interaction with Gankyrin. This article reports for the first time that treatment of WSE decreased the cancerous properties through host cellular response modulation in gastric epithelial cells coinfected with H. pylori and EBV.Communicated by Ramaswamy H. Sarma.

The exploration of phytocompounds theoretically combats SARS-CoV-2 pandemic against virus entry, viral replication and immune evasion

BACKGROUND

The novel coronavirus disease-2019 (COVID-19) that emerged in China, is an extremely contagious and pathogenic viral infection caused by the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) that has sparked a global pandemic. The few and limited availability of approved therapeutic agents or vaccines is of great concern. Urgently, Remdesivir, Nirmatrelvir, Molnupiravir, and some phytochemicals including polyphenol, flavonoid, alkaloid, and triterpenoid are applied to develop as repurposing drugs against the SARS-CoV-2 invasion.

METHODS

This study was conducted to perform molecular docking and absorption, distribution, metabolism, excretion and toxicity (ADMET) analysis of the potential phytocompounds and repurposing drugs against three targets of SARS-CoV-2 proteins (RNA dependent RNA polymerase, RdRp, Endoribonclease, S-protein of ACE2-RBD).

RESULTS

The docking data illustrated Arachidonic acid, Rutin, Quercetin, and Curcumin were highly bound with coronavirus polyprotein replicase and Ebolavirus envelope protein. Furthermore, anti- Ebolavirus molecule Remedesivir, anti-HIV molecule Chloroquine, and Darunavir were repurposed with coronavirus polyprotein replicase as well as Ebolavirus envelope protein. The strongest binding interaction of each targets are Rutin with RdRp, Endoribonclease with Amentoflavone, and ACE2-RBD with Epigallocatechin gallate.

CONCLUSIONS

Taken altogether, these results shed a light on that phytocompounds have a therapeutic potential for the treatment of anti-SARS-CoV-2 may base on multi-target effects or cocktail formulation for blocking viral infection through invasion/activation, transcription/reproduction, and posttranslational cleavage to battle COVID-19 pandemic.

Finding inhibitors and deciphering inhibitor-induced conformational plasticity in the Janus kinase via multiscale simulations

The Janus kinase (JAK) is a master regulator of the JAK/STAT pathway. Dysregulation of this signalling cascade causes neuroinflammation and autoimmune disorders. Therefore, JAKs have been characterized as an attractive target for developing anti-inflammatory drugs. Nowadays, designing efficient, effective, and specific targeted therapeutics without being cytotoxic has gained interest. We performed the virtual screening of natural products in combination with pharmacological analyses. Subsequently, we performed molecular dynamics simulations to study the stability of the ligand-bound complexes and ligand-induced inactive conformations. Notably, inactive kinases display remarkable conformational plasticity; however, ligand-induced molecular mechanisms of these conformations are still poorly understood. Herein, we performed a free energy landscape analysis to explore the conformational plasticity of the JAK1 kinase. Leonurine, STOCK1N-68642, STOCK1N-82656, and STOCK1N-85809 bound JAK1 exhibited a smooth transition from an active (αC-in) to a completely inactive conformation (αC-out). Ligand binding induces disorders in the αC-helix. Molecular mechanics Poisson Boltzmann surface area (MM/PBSA) calculation suggested three phytochemicals, namely STOCK1N-68642, Epicatechin, and STOCK1N-98615, have higher binding affinity compared to other ligand molecules. The ligand-induced conformational plasticity revealed by our simulations differs significantly from the available crystal structures, which might help in designing allosteric drugs.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Plant-derived active compounds as a potential nucleocapsid protein inhibitor of SARS-CoV-2: an in-silico study

The coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2. This virus has a high mismatch repair proofreading ability due to its unique exonuclease activity, making it knotty to treat. The nucleocapsid protein can serve as a potential antiviral drug target, as this protein is responsible for multiple captious functions during the viral life cycle. Herein, we have investigated the potential to repurpose active antiviral compounds of plant origins for treating the SARS-CoV-2 infection. In the present study, we followed the molecular docking methodology to screen druggable natural plants' active compounds against the nucleocapsid protein of SARS-CoV-2. The virtual screening of all 68 compounds revealed that the top seven active compounds, such as withanolide D, hypericin, silymarin, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein, exhibit good binding affinity with druggable ADME properties, toxicity, and Pass prediction. The stability of the docked complexes was studied by conducting molecular simulations of 100 ns. MM-GBSA calculated the binding free energy uncovered that withanolide D, hypericin, and silymarin result in highly stable binding conformations in three different sites of the nucleocapsid protein. However, further investigation is needed in order to validate the candidacy of these inhibitors for clinical trials. HighlightsNatural plants' active compounds may aid in the inhibition of SARS-CoV-2 replication and COVID-19 therapeutics.Hypericin, silymarin, withanolide D, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein are effective against SARS-CoV-2 N protein.Studied natural plants' active compounds could be useful against COVID-19 and its associated organs comorbidities.ADMET properties of selected compounds favor these compounds as druggable candidates.Communicated by Ramaswamy H. Sarma.

Comparative Structural Dynamics of Isoforms of Helicobacter pylori Adhesin BabA Bound to Lewis b Hexasaccharide via Multiple Replica Molecular Dynamics Simulations

BabA of Helicobacter pylori is the ABO blood group antigen-binding adhesin. Despite considerable diversity in the BabA sequence, it shows an extraordinary adaptation in attachment to mucosal layers. In the current study, multiple replica molecular dynamics simulations were conducted in a neutral aqueous solution to elucidate the conformational landscape of isoforms of BabA bound to Lewis b (Le^b) hexasaccharide. In addition, we also investigated the underlying molecular mechanism of the BabA-glycan complexation using the MM/GBSA scheme. The conformational dynamics of Le^b in the free and protein-bound states were also studied. The carbohydrate-binding site across the four isoforms was examined, and the conformational variability of several vital loops was observed. The cysteine–cysteine loops and the two diversity loops (DL1 and DL2) were identified to play an essential role in recognizing the glycan molecule. The flexible crown region of BabA was stabilized after association with Le^b. The outward movement of the DL2 loop vanished upon ligand binding for the Spanish specialist strain (S381). Our study revealed that the S831 strain shows a stronger affinity to Le^b than other strains due to an increased favorable intermolecular electrostatic contribution. Furthermore, we showed that the α1-2-linked fucose contributed most to the binding by forming several hydrogen bonds with key amino acids. Finally, we studied the effect of the acidic environment on the BabA-glycan complexation via constant pH MD simulations, which showed a reduction in the binding free energy in the acidic environment. Overall, our study provides a detailed understanding of the molecular mechanism of Le^b recognition by four isoforms of H. pylori that may help the development of therapeutics targeted at inhibiting H. pylori adherence to the gastric mucosa.

Phosphorylation-Induced Conformational Dynamics and Inhibition of Janus Kinase 1 by Suppressors of Cytokine Signaling 1

The dysfunction of the JAK/STAT (Janus kinase/signal transducers and activators of transcription) pathway results in several pathophysiological conditions, including autoimmune disorders. The negative feedback regulators of the JAK/STAT signaling pathway, suppressors of cytokine signaling (SOCS), act as a natural inhibitor of JAK and inhibit aberrant activity. SOCS1 is the most potent member of the SOCS family, whose kinase inhibitory region targets the substrate-binding groove of JAK with high affinity and blocks the phosphorylation of JAK kinases. Overall, we performed an aggregate of 13 μs molecular dynamics simulations on the activation loop's three different phosphorylation (double and single) states. Results from our simulations show that the single Tyr¹⁰³⁴ phosphorylation could stabilize the JAK1/SOCS1 complex as well as the flexible activation segment. The phosphate-binding loop (P-loop) shows conformational variability at dual and single phosphorylated states. Principal component analysis and protein structure network (PSN) analysis reveal that the different phosphorylation states and SOCS1 binding induce intermediate inactive conformations of JAK1, which could be a better target for future JAK1 selective drug design. PSN analysis suggests that the com-pY1034 system is stabilized due to higher values of network hubs than the other two complex systems. Moreover, the binding free energy calculations suggest that pTyr¹⁰³⁴ states show a higher affinity toward SOCS1 than the dual and pTyr¹⁰³⁵ states. We believe that the mechanistic understanding of JAK1/SOCS1 complexation will aid future studies related to peptide inhibitors based on SOCS1.