BACKGROUND: Breast cancer stem cells (BCSCs) are a small subset of cells within breast tumors with characteristics similar to normal stem cells. Despite advancements in chemotherapy and targeted therapy for breast cancer, the prognosis for breast cancer patients has remained poor due to drug resistance, reoccurrence, and metastasis. Growing evidence suggests that deregulation of the self-renewal pathways, like the Wnt signaling pathway mediated by β-catenin, plays a crucial role in the survival of breast cancer stem cells. Targeting the Wnt signaling pathway in breast cancer stem cells offers a promising avenue for developing effective therapeutic strategies targeting these cells, potentially leading to improved patient outcomes and reduced tumor recurrence.
METHODS: For this purpose, we have screened a 1615 FDA-approved drug library against our target protein, β-catenin, which is involved in the Wnt signaling pathway using molecular docking analysis, molecular dynamics (MD) simulations, and molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations.
RESULTS: Molecular docking studies showed that the Lumacaftor- β-catenin complex had the lowest docking score of - 8.7 kcal/mol towards β-catenin protein than the reference inhibitor. Molecular dynamic simulations and MM/PBSA calculations were also performed for the Lumacaftor-β-catenin complex to establish the stability of the interactions involved. Considering its promising attributes and encouraging results, Lumacaftor holds significant potential as a novel therapeutic option to target BCSCs. This study opens avenues for further investigation and may pave the way for developing therapeutic potential in breast cancer treatment. Further confirmation is warranted through in vitro and clinical studies to validate the findings of this study.

Osmolytes are small organic molecules that are known to stabilize proteins and other biological macromolecules under various stressful conditions. They belong to various categories such as amino acids, methylamines, and polyols. These substances are commonly known as \'compatible solutes\' because they do not disrupt cellular processes and help regulate the osmotic balance within cells. In the case of ribonuclease A (RNase A), which is prone to aggregation, the presence of osmolytes can help to maintain its structural stability and prevent unwanted interactions leading to protein aggregation. In this study, we investigated the interaction between RNase A and several osmolytes using molecular docking and molecular dynamics (MD) simulations. We performed molecular docking to predict the binding mode and binding affinity of each osmolyte with RNase A. MD simulations were then carried out to investigate the dynamics and stability of the RNase A-osmolyte complexes. Our results show that two osmolytes, glucosylglycerol and sucrose have favorable binding affinities with RNase A. The possible role of these osmolytes in stabilizing the RNase A and prevention of aggregation is also explored. By providing computational insights into the interaction between RNase A and osmolytes, the study offers valuable information that could aid in comprehending the mechanisms by which osmolytes protect proteins and help in designing therapeutics for protein-related disorders based on osmolytes. These findings may have significant implications for the development of novel strategies aimed at preventing protein misfolding and aggregation in diverse disease conditions.Communicated by Ramaswamy H. Sarma.

Cancer is a serious threat to human life and social development and the use of scientific methods for cancer prevention and control is necessary. In this study, HQSAR, CoMFA, CoMSIA and TopomerCoMFA methods are used to establish models of 65 imidazo[4,5-b]pyridine derivatives to explore the quantitative structure-activity relationship between their anticancer activities and molecular conformations. The results show that the cross-validation coefficients q2 of HQSAR, CoMFA, CoMSIA and TopomerCoMFA are 0.892, 0.866, 0.877 and 0.905, respectively. The non-cross-validation coefficients r2 are 0.948, 0.983, 0.995 and 0.971, respectively. The externally validated complex correlation coefficients r2pred of external validation are 0.814, 0.829, 0.758 and 0.855, respectively. The PLS analysis verifies that the QSAR models have the highest prediction ability and stability. Based on these statistics, virtual screening based on R group is performed using the ZINC database by the Topomer search technology. Finally, 10 new compounds with higher activity are designed with the screened new fragments. In order to explore the binding modes and targets between ligands and protein receptors, these newly designed compounds are conjugated with macromolecular protein (PDB ID: 1MQ4) by molecular docking technology. Furthermore, to study the nature of the newly designed compound in dynamic states and the stability of the protein-ligand complex, molecular dynamics simulation is carried out for N3, N4, N5 and N7 docked with 1MQ4 protease structure for 50 ns. A free energy landscape is computed to search for the most stable conformation. These results prove the efficient and stability of the newly designed compounds. Finally, ADMET is used to predict the pharmacology and toxicity of the 10 designed drug molecules.

Cathepsin K is a type of cysteine proteinase that is primarily expressed in osteoclasts and has a key role in the breakdown of bone matrix protein during bone resorption. Many studies suggest that the deficiency of cathepsin K is concomitant with a suppression of osteoclast functioning, therefore rendering the resorptive properties of cathepsin K the most prominent target for osteoporosis. This innovative work has identified a novel anti-osteoporotic agent against Cathepsin K by using a comparison of machine learning and deep learning-based virtual screening followed by their biological evaluation. Out of ten shortlisted compounds, five of the compounds (JFD02945, JFD02944, RJC01981, KM08968 and SB01934) exhibit more than 50% inhibition of the Cathepsin K activity at 0.1 μM concentration and are considered to have a promising inhibitory effect against Cathepsin K. The comprehensive docking, MD simulation, and MM/PBSA investigations affirm the stable and effective interaction of these compounds with Cathepsin K to inhibit its function. Furthermore, the compounds RJC01981, KM08968 and SB01934 are represented to have promising anti-osteoporotic properties for the management of osteoporosis owing to their significantly well predicted ADMET properties.

DPP4 inhibitors can control glucose homeostasis by increasing the level of GLP-1 incretins hormone due to dipeptidase mimicking. Despite the potent effects of DPP4 inhibitors, these compounds cause unwanted toxicity attributable to their effect on other enzymes. As a result, it seems essential to find novel and DPP4 selective compounds. In this study, we introduce a potent and selective DPP4 inhibitor via structure-based virtual screening, molecular docking, molecular dynamics simulation, MM/PBSA calculations, DFT analysis, and ADMET profile. The screened compounds based on similarity with FDA-approved DPP4 inhibitors were docked towards the DPP4 enzyme. The compound with the highest docking score, ZINC000003015356, was selected. For further considerations, molecular docking studies were performed on selected ligands and FDA-approved drugs for DPP8 and DPP9 enzymes. Molecular dynamics simulation was run during 200 ns and the analysis of RMSD, RMSF, Rg, PCA, and hydrogen bonding were performed. The MD outputs showed stability of the ligand-protein complex compared to available drugs in the market. The total free binding energy obtained for the proposed DPP4 inhibitor was more negative than its co-crystal ligand (N7F). ZINC000003015356 confirmed the role of the five Lipinski rule and also, have low toxicity parameter according to properties. Finally, DFT calculations indicated that this compound is sufficiently soft.

Lupus Nephritis (LN) is an autoimmune disease affecting the kidneys, and conventional drug studies have limitations due to its imprecise and complex pathogenesis. Therefore, the aim of this study was to design a novel Lupus Nephritis-targeted drug with good clinical due potential, high potency and selectivity by computer-assisted approach.NIK belongs to the serine/threonine protein kinase, which is gaining attention as a drug target for Lupus Nephritis. we used bioinformatics, homology modelling and sequence comparison analysis, small molecule ab initio design, ADMET analysis, molecular docking, molecular dynamics simulation, and MM/PBSA analysis to design and explore the selectivity and efficiency of a novel Lupus Nephritis-targeting drug, ClImYnib, and a classical NIK inhibitor, NIK SMI1. We used bioinformatics techniques to determine the correlation between lupus nephritis and the NF-κB signaling pathway. De novo drugs design was used to create a NIK-targeted inhibitor, ClImYnib, with lower toxicity, after which we used molecular dynamics to simulate NIK SMI1 against ClImYnib, and the simulation results showed that ClImYnib had better selectivity and efficiency. Our research delves into the molecular mechanism of protein ligands, and we have designed and validated an excellent NIK inhibitor using multiple computational simulation methods. More importantly, it provides an idea of target designing small molecules.Communicated by Ramaswamy H. Sarma.

The opportunistic bacterium Acinetobacter baumannii, which belongs to ESKAPE group of pathogenic bacteria, is leading cause of infections associated with gram-negative bacteria. Acinetobacter baumannii causes severe diseases, such as VAP (ventilator-associated pneumonia), meningitis, and UTI (urinary tract infections) among the nosocomial infections contracted in hospitals. The high infection rate and growing resistance to the vast array of antibiotics makes it paramount to look for new therapeutic strategies against this pathogen. The most promising therapeutic targets are the proteins involved in the synthesis of peptidoglycan which is chief component of bacterial cell wall, MurE is one of those enzymes and is responsible for the addition of one unit of meso-diaminopimelic acid (meso-A2pm) to the nucleotide precursor, UDPMurNAc-L-Ala-D-Glu, and aids in the formation of crosslinker pentapeptide chain. The three-dimensional structure of MurE was modelled using homology modelling technique and then vHTS was performed using this model against Approved Drug Library on DrugRep server using AutoDock Vina. Out of 500 drug molecules, two were selected based on estimated binding affinity, interaction pattern, interacting residues, etc. The selected drug molecules are DB12887 (Tazemetostat) and DB13879 (Glecaprevir). Then, MD simulations were performed on native MurE and its complexes with ligands to examine their dynamical behaviour, stability, integrity, compactness, and folding properties. The protein-ligand complexes were then subjected to binding free energy calculations using the MM/PBSA-based binding free energy analysis and the values are -109.788 ± 8.03 and -152.753 ± 11.98 kcal for MurE-DB12887 and MurE-DB13879 complex, respectively. All the analysis performed on MD trajectories for the complexes of these ligands with protein provided plenty of dependable evidences to consider these molecules for inhibition of MurE enzyme from A. baumannii. Communicated by Ramaswamy H. Sarma.

Alzheimer\'s is characterized as a progressive neurodegenerative disease due to beta-amyloid accumulation in the brain. Some previous studies reported that RXR agonists could be effective in the treatment of Alzheimer\'s disease. There are currently numerous attempts being made to discover a natural RXR agonist that is more potent than 9-cis-retinoic acid (9CR). One of the most efficient resources for finding high-potential compounds is natural databases. In this study, 81215 compounds from the IB screen library as natural databases were docked against the RXR-alpha binding site. The best compounds discovered interact with the RXR-alpha binding site with a lower binding energy (-11 to -13 kcal/mol) than the binding energy of -10.94 kcal/mol for 9-cis, which means that these compounds could interact stronger with RXR-alpha than 9CR. All selected compounds could pass the blood-brain barrier. Physiochemical properties assessment indicated that all compounds passed Lipinski\'s rule and had the potential to be oral drug candidates. The stability of protein-ligand complexes during a timescale of 100 ns by Molecular Dynamics simulation demonstrated that all compounds could effectively interact with the RXR binding site. The molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) represented that all selected hit compounds had a better binding affinity to the alpha RXR binding site compared to 9CR, which means these hit compounds had potential drug candidates for the treatment of Alzheimer\'s disease. However, experimental assessment is needed to validate this result.Communicated by Ramaswamy H. Sarma.

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop\'s flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI\'s binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Acinetobacter baumannii is one of the emerging causes of hospital acquired infections and this bacterium, due to multi-drug resistant and Extensive Drug resistant has been able to develop resistance against the antimicrobial agents that are being used to eliminate it. A.baumannii has been the cause of death in immune compromised patients in hospitals. Hence it is the urgent need of time to find potential inhibitors for this bacterium to cease its virulence and affect its survival inside host organisms. The Dihydrofolate reductase enzyme, which is an important biocatalyst in the conversion of Dihydrofolate to Tetrahydrofolate, is an important drug target protein. In the present study high throughput screening is used to identify the inhibitors of this enzyme. The prioritized ligand molecular candidates identified through virtual screening for the substrate binding site of the predicted model are Z1447621107, Z2604448220 and Z1830442365. The Molecular Dynamics Simulation study suggests that potential inhibitor of the Dihydrofolate reductase enzyme would prevent bacteria from completing its life cycle, affecting its survival. Finally the complexes were analysed for binding free energy of the Dihydrofolate reductase enzyme complexes with the ligands.