Extra precision glide docking, free energy calculation and molecular dynamics studies of 1,2-diarylethane derivatives as potent urease inhibitors

Leveraging shape screening and molecular dynamics simulations to optimize PARP1-Specific chemo/radio-potentiators for antitumor drug design

PARP1 plays a pivotal role in DNA repair within the base excision pathway, making it a promising therapeutic target for cancers involving BRCA mutations. Current study is focused on the discovery of PARP inhibitors with enhanced selectivity for PARP1. Concurrent inhibition of PARP1 with PARP2 and PARP3 affects cellular functions, potentially causing DNA damage accumulation and disrupting immune responses. In step 1, a virtual library of 593 million compounds has been screened using a shape-based screening approach to narrow down the promising scaffolds. In step 2, hierarchical docking approach embedded in Schrödinger suite was employed to select compounds with good dock score, drug-likeness and MMGBSA score. Analysis supplemented with decomposition energy, molecular dynamics (MD) simulations and hydrogen bond frequency analysis, pinpointed that active site residues; H862, G863, R878, M890, Y896 and F897 are crucial for specific binding of ZINC001258189808 and ZINC000092332196 with PARP1 as compared to PARP2 and PARP3. The binding of ZINC000656130962, ZINC000762230673, ZINC001332491123, and ZINC000579446675 also revealed interaction involving two additional active site residues of PARP1, namely N767 and E988. Weaker or no interaction was observed for these residues with PARP2 and PARP3. This approach advances our understanding of PARP-1 specific inhibitors and their mechanisms of action, facilitating the development of targeted therapeutics.

Molecular modelling of SdiA protein by selected flavonoid and terpenes compounds to attenuate virulence in Klebsiella pneumoniae

Klebsiella pneumoniae is one of the perturbing multidrug resistant (MDR) and ESKAPE pathogens contributing to the mounting morbidity, mortality and extended rate of hospitalization. Its virulence, often regulated by quorum sensing (QS) reinforces the need to explore alternative and prospective antivirulence agents, relatively from plants secondary metabolites. Computer aided drug discovery using molecular modelling techniques offers advantage to investigate prospective drugs to combat MDR pathogens. Thus, this study employed virtual screening of selected terpenes and flavonoids from medicinal plants to interrupt the QS associated SdiA protein in K. pneumoniae to attenuate its virulence. 4LFU was used as a template to model the structure of SdiA. ProCheck, Verify3D, Ramachandran plot scores, and ProSA-Web all attested to the model's good quality. Since SdiA protein in K. pneumoniae leads to the expression of virulence, 31 prospective bioactive compounds were docked for antagonistic potential. The stability of the protein-ligand complex, atomic motions and inter-atomic interactions were further investigated through molecular dynamics simulations (MDS) at 100 ns production runs. The binding free energy was estimated using the molecular mechanics/poisson-boltzmann surface area (MM/PB-SA). Furthermore, the drug-likeness properties of the studied compounds were validated. Docking studies showed phytol possesses the highest binding affinity (-9.205 kcal/mol) while glycitein had -9.752 kcal/mol highest docking score. The MDS of the protein in complex with the best-docked compounds revealed phytol with the highest binding energy of -44.2625 kcal/mol, a low root-mean-square deviation (RMSD) value of 1.54 Å and root-mean-square fluctuation (RMSF) score of 1.78 Å. Analysis of the drug-likeness properties prediction and bioavailability of these compounds revealed their conformed activity to lipinski's rules with bioavailability scores of 0.55 F. The studied terpenes and flavonoids compounds effectively thwart SdiA protein, therefore regulate inter- or intra cellular communication and associated in virulence Enterobacteriaceae, serving as prospective antivirulence drugs.Communicated by Ramaswamy H. Sarma.

Atomic level and structural understanding of natural ligands inhibiting Helicobacter pylori peptide deformylase through ligand and receptor based screening, SIFT, molecular dynamics and DFT - a structural computational approach

Helicobacter pylori is a Gram-negative microaerophilic gastric pathogen, responsible for the cause of peptic ulcer around half of the global population. Although several antibiotics and combination therapies have been employed for H. pylori-related gastric ulcer and cancer regiments, identifying potent inhibitors for specific targets of this bacterium will help assessing better treatment periodicity and methods to eradicate H. pylori. Herein, 1,000,000 natural compounds were virtually screened against Helicobacter pylori Peptide deformylase (HpPDF). Pharmacophore hypotheses were created using ligand and receptor-based pharmacophore modeling of GLIDE. Stringent HTVS and IFD docking protocol of GLIDE predicted leads with stable intermolecular bonds and scores. Molecular dynamics simulation of HpPDF was carried out for 100 ns using GROMACS. Hits ZINC00225109 and ZINC44896875 came up with a glide score of -9.967 kcal/mol and -12.114 kcal/mol whereas; reference compound actinonin produced a glide score of -9.730 kcal/mol. Binding energy values of these hits revealed the involvement of significant Van der Waals and Coulomb forces and the deduction of lipophilic forces that portray the deep hydrophobic residues in the S1pocket of H. pylori. The DFT analysis established the electron density-based features of the molecules and observed that the results correlate with intermolecular docking interactions. Analysis of the MD trajectories revealed the crucial residues involved in HpPDF - ligand binding and the conformational changes in the receptor. We have identified and deciphered the crucial features necessary for the potent ligand binding at catalytic site of HpPDF. The resulting ZINC natural compound hits from the study could be further employed for potent drug development.Communicated by Ramaswamy H. Sarma.

In Silico and In Vitro Screening of Antipathogenic Properties of Melianthus comosus (Vahl) against Pseudomonas aeruginosa

Bacterial quorum sensing (QS) system regulates pathogenesis, virulence, and biofilm formation, and together they contribute to nosocomial infections. Opportunistic pathogens, such as Pseudomonas aeruginosa, rely on QS for regulating virulence factors. Therefore, blocking the QS system may aid management of various infectious diseases caused by human pathogens. Plant secondary metabolites can thwart bacterial colonization and virulence. As such, this study was undertaken to evaluate three extracts from the medicinal plant, Melianthus comosus, from which phytochemical compounds were identified with potential to inhibit QS-dependent virulence factors in P. aeruginosa. Chemical profiling of the three extracts identified 1,2-benzene dicarboxylic acid, diethyl ester, neophytadiene and hexadecanoic acid as the common compounds. Validation of antibacterial activity confirmed the same MIC values of 0.78 mg/mL for aqueous, methanol and dichloromethane extracts while selected guanosine showed MIC 0.031 mg/mL. Molecular docking analysis showed anti-quorum sensing (AQS) potential of guanosine binding to CviR’ and 2UV0 proteins with varying docking scores of −5.969 and −8.376 kcal/mol, respectively. Guanosine inhibited biofilm cell attachment and biofilm development at 78.88% and 34.85%, respectively. Significant swimming and swarming motility restriction of P. aeruginosa were observed at the highest concentration of plant extracts and guanosine. Overall, guanosine revealed the best swarming motility restrictions. M. comosus extracts and guanosine have shown clear antibacterial effects and subsequent reduction of QS-dependent virulence activities against P.aeruginosa. Therefore, they could be ideal candidates in the search for antipathogenic drugs to combat P.aeruginosa infections.

Drug screening and identification of key candidate genes and pathways of rheumatoid arthritis

Rheumatoid arthritis (RA), which normally manifests as a multi-joint inflammatory reaction, is a common immunological disease in clinical practice. However, the pathogenesis of RA has not yet been fully elucidated. Rituximab (RTX) is an effective drug in the treatment of RA, however its therapeutic efficacy and mechanism of action require further investigation. Thus, the present study aimed to screen the candidate key regulatory genes and explain the potential mechanisms of RA. Gene chips of RA and normal joint tissues were analyzed and, gene chips of RTX before and after treatment were investigated. In the present study, strong evidence supporting the pathogenesis of RA and mechanism of action of RTX were also revealed. Differentially expressed genes (DEGs) were analyzed using the limma package of RStudio software. A total of 1,150 DEGs were detected in RA compared with normal joint tissues. The upregulated genes were enriched in ‘interleukin-12 production’, ‘I-κB kinase/NF-κB signaling’, ‘regulation of cytokine production involved in immune response’ and ‘cytokine metabolic process’. Functional enrichment analysis showed that RTX was primarily involved in the inhibition of ‘adaptive immune response’, ‘B cell activation involved in immune response’ and ‘immune effector process’. Subsequently, leukocyte immunoglobulin-like receptor subfamily B member 1 (LILRB1), a hub gene with high connectivity degree, was selected, and traditional Chinese medicine libraries were molecularly screened according to the structure of the LILRB1 protein. The results indicated that kaempferol 3-O-β-D-glucosyl-(1→2)-β-D-glucoside exhibited the highest docking score. In the present study, the DEGs and their biological functions in RA and the pharmacological mechanism of RTX action were determined. Taken together, the results suggested that LILRB1 may be used as a molecular target for RA treatment, and kaempferol 3-O-β-D-glucosyl-(1→2)-β-D-glucoside may inhibit the pathological process of RA.

Identification of Key Candidate Genes and Pathways of Candida albicans-Infected Human Umbilical Vein Endothelial Cells and Drug Screening

Candida albicans is a common opportunistic pathogen that can cause serious infection by blood transmission. C. albicans enters the blood circulation and adheres to the endothelial cells of the vascular wall. However, the detailed mechanism of the effect of C. albicans on the endothelial cells remains unclear. In this study, the microarray expression profile of human umbilical vein endothelial cells exposed to C. albicans was analyzed. The 191 up-regulated genes were enriched in TNF, T cell receptor, and NF-kappa B signaling pathways. The 71 down-regulated genes were enriched in pyruvate metabolic, purine nucleotide metabolic, purine nucleotide biosynthetic, and humoral immune response processes. Gene set enrichment analysis showed that apoptosis, oxidative phosphorylation, IL6/JAK/STAT3 signaling pathways were enriched. Moreover, two hub genes with a high degree of connectivity, namely, MYC and IL6, were selected. Molecular screening of traditional Chinese medicine libraries was performed on the basis of the structure of MYC protein. The okanin had the highest docking score. MYC might be used as molecular targets for treatment. In addition, okanin may inhibit the infection of C. albicans. Thus, MYC can be subjected to further research.