Discovery of potential competitive inhibitors against With-No-Lysine kinase 1 for treating hypertension by virtual screening, inverse pharmacophore-based lead optimization, and molecular dynamics simulations

Multilayered screening for multi-targeted anti-Alzheimer's and anti-Parkinson's agents through structure-based pharmacophore modelling, MCDM, docking, molecular dynamics and DFT: a case study of HDAC4 inhibitors

Abstract

Alzheimer's disease (AD) and Parkinson's disease (PD) are neurological conditions that primarily impact the elderly having distinctive traits and some similarities in terms of symptoms and progression. The multifactorial nature of AD and PD encourages exploring potentiality of multi-target therapy for addressing these conditions to conventional, the "one drug one target" strategy. This study highlights the searching of potential HDAC4 inhibitors through multiple screening approaches. In this context, structure-based pharmacophore model, ligand profiler mapping and MCDM approaches were performed for target prioritization. Similarly, ligand profiler, MCDM and Docking studies were performed to prioritize multi-targeted HDAC4 inhibitors. These comprehensive approaches unveiled 5 common targets and 5 multi-targeted prioritized compounds consensually. MD simulations, DFT and binding free energy calculations corroborated the stability and robustness of propitious compound 774 across 5 prioritized targets. In conclusion, the screened compound 774 (ChEMBL 4063938) could be a promising multi-targeted therapy for managing AD and PD further rendering experimental validation.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00302-4.

Elucidating the conformational dynamics of histo-blood group antigens and their interactions with the rotavirus spike protein through computational lens

In the present study, we investigated the conformational dynamics of histo-blood group antigens (HBGAs) and their interactions with the VP8* domain of four rotavirus genotypes (P[4], P[6], P[19], and P[11]) utilizing all-atom molecular dynamics simulations in explicit water. Our study revealed distinct changes in the dynamic behavior of the same glycan due to linkage variations. We observed that LNFPI HBGA having a terminal β linkage shows two dominant conformations after complexation, whereas only one was obtained for LNFPI with a terminal α linkage. Interestingly, both variants displayed a single dominant structure in the free state. Similarly, LNT and LNnT show a shift in their dihedral linkage profile between their two terminal monosaccharides because of a change in the linkage from β(1-3) to β(1-4). The molecular mechanics generalized Born surface area (MM/GBSA) calculations yielded the highest binding affinity for LNFPI(β)/P[6] (-13.93 kcal/mol) due to the formation of numerous hydrogen bonds between VP8* and HBGAs. LNnT binds more strongly to P[11] (-12.88 kcal/mol) than LNT (-4.41 kcal/mol), suggesting a single change in the glycan linkage might impact its binding profile significantly. We have also identified critical amino acids and monosaccharides (Gal and GlcNAc) that contributed significantly to the protein-ligand binding through the per-residue decomposition of binding free energy. Moreover, we found that the interaction between the same glycan and different protein receptors within the same rotavirus genogroup influenced the micro-level dynamics of the glycan. Overall, our study helps a deeper understanding of the H-type HBGA and rotavirus spike protein interaction.Communicated by Ramaswamy H. Sarma.

Exploring molecular interactions of potential inhibitors against the spleen tyrosine kinase implicated in autoimmune disorders via virtual screening and molecular dynamics simulations

The spleen tyrosine kinase (Syk) plays a pivotal role in immune cells' signal transduction mechanism. While fostamatinib, an FDA-approved Syk inhibitor, is currently used to treat immune thrombocytopenia, the search for improved Syk-targeted medications to treat autoimmune diseases is still underway. Herein, we screened 38,493 compounds against Syk and selected eight leads based on the docking score and ADMET properties, and performed 3× 200 ns long molecular dynamics simulations of the apo and Syk-ligand complexes. We considered R406, the active component of fostamatinib, as a control. The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations demonstrated the lead1 (Δ G b i n d  = -30.35 kcal/mol) exhibited a similar binding free energy as the control (Δ G b i n d = -29.82 kcal/mol). The Syk stabilizing effect of lead1 was also indicated in its network features, sampling space, and residual correlation motion analysis. We further generated 100 structural analogues of lead1 using deep learning, and one of the analogues displayed a better binding free energy (Δ G b i n d = -47.58 kcal/mol) compared to the control or lead1, facilitated by more favourable van der Waals interactions and lesser binding-opposing net polar forces. This analogue may be further exploited to develop effective therapeutics against Syk-associated diseases after validation in vitro and in vivo.

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.

An update regarding the role of WNK kinases in cancer

Mammalian WNK kinases (WNKs) are serine/threonine kinases that contain four members, WNK1-4. They function to maintain ion homeostasis and regulate blood pressure in mammals. Recent studies have revealed that the dysregulation of WNKs contributes to tumor growth, metastasis, and angiogenesis through complex mechanisms, especially through phosphorylating kinase substrates SPS1-related proline/alanine-rich kinase (SPAK) and oxidative stress-responsive kinase 1 (OSR1). Here, we review and discuss the relationships between WNKs and several key factors/biological processes in cancer, including ion channels, cation chloride cotransporters, sodium bicarbonate cotransporters, signaling pathways, angiogenesis, autophagy, and non-coding RNAs. In addition, the potential drugs for targeting WNK-SPAK/OSR1 signaling have also been discussed. This review summarizes and discusses knowledge of the roles of WNKs in cancer, which provides a comprehensive reference for future studies.