Current computer-aided drug design最新文献

Background: Si-Ni-San (SNS) is the formula prescription of Traditional Chinese Medicine (TCM) with anti-depression properties, but its underlying mechanisms remain unclear.

Objective: This study provides novel approaches for the study of TCM and offers new opportunities for exploring the pharmacological properties of SNS.

Methods: The ingredients in SNS implicated in the treatment of depression were identified and studied using network pharmacology. SwissTargetPrediction and molecular docking were used to study the interaction of SNS ingredients and their targets. The protective effect of these ingredients and their cocktail in rat pheochromocytoma cells (PC12) exposed to corticosterone (Cor) were evaluated using the CCK-8 assay, Hoechst 33342 staining, 2',7'-dichlorodihydro fluorescein diacetate (H2DCFDA) staining, 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging assay, and in-cell Western analysis.

Results: The network pharmacology study showed that the HIF-1 signaling pathway was the most crucial pathway implicated in the anti-depressive property of SNS. MAPK1 (ERK2), MAPK3 (ERK1), AKT1, VEGFA, STAT3, and EGF were identified as hub target proteins in the HIF-1 signaling pathway. Quercetin, naringenin, licochalcone A, and kaempferol from SNS, which targeted the six proteins mentioned above, were used to create a cocktail. This cocktail exerted protective properties, decreased the oxidative stress in PC12 exposed to Cor, and successfully regulated the expressions of AKT1, p-AKT1, ERK1, ERK2, p-ERK1/2, STAT3, p- STAT3, and VEGFA induced by Cor exposure. The SwissTargetPrediction and molecular docking study showed that the cocktail may regulate the HIF-1 signaling pathway by directly binding with AKT1 and MAPK1.

Conclusion: The cocktail from SNS comprised of quercetin, naringenin, licochalcone A, and kaempferol exerts anti-depression potentiality by modulating the HIF-1 signaling pathway via direct interactions with AKT1 and MAPK1.

Background: Virtual screening (VS) could select possible effective candidates from a large number of organic compounds, which plays an important role in network pharmacology. Virtual screening is a very important step in network pharmacology.

Objective: The accuracy of screening compounds directly determines the subsequent network construction, target determination and pathway analysis. In order to improve the accuracy of screening the important compounds in herbs for treating diabetes, a novel methodology based on complex-valued flexible neural tree (CVFNT) model and negative sample selection algorithm is presented.

Methods: In our method, diabetes-related targets were obtained by literature search. According to diabetes-related targets, active compounds were searched from the public database. The negative sample selection algorithm based on Tanimoto index was proposed to establish inactive compound set. The CVFNT model optimized was utilized to screen effective candidate compounds.

Result: Our proposed method performs better than eight classical classifiers in terms of TPR, FPR, Precision, Specificity, F1, AUC and ROC curve. Our method could also predict 18 compounds from Liangxue Sanyu Decoction, which are involved in the treatment of diabetes.

Conclusion: Our proposed method could effectively improve the identification accuracy of diabetes- related compounds, enhance the performance of network pharmacology analysis, and provide a new idea for drug research and development for the treatment and prevention of diabetes.

Background: Thermolabile hemolysin (TLH) is a key virulent protein of Vibrio alginolyticus, known for its hemolytic and phospholipase activities, leading to shrimp vibriosis disease. It has been suggested as a potential therapeutic candidate for vibriosis therapy.

Methods: Computational studies, including molecular docking, toxicity analysis, and molecular dynamics (MD) simulations, were conducted to investigate the inhibition of the phospholipase activity of TLH by phytochemicals from Emblica officinalis.

Results: Out of the twenty-nine compounds, the top three, including Ellagic acid (CID 5281855), Quercetin (CID 5280343), and Kaempferol (CID 5280863), were sorted based on their highest molecular docking scores of -9.2, -8.9, and -8.8, respectively. Subsequently, molecular dynamics (MD) simulations of these selected leads were performed to observe the structural stability of these compounds in the binding sites of TLH protein. The MD simulation outcomes indicated that all three compounds demonstrated superior stability throughout 100 nanoseconds compared to the control compound Resveratrol. The molecular simulation results suggest stable interactions, with average root-mean-square deviation (RMSD) and root-meansquare fluctuation (RMSF) values of 1-2 Å and 0-3 Å. Pharmacokinetic and toxicity analyses were conducted to evaluate the suitability and toxicity of these selected compounds. All top three compounds passed the Lipinski rule, and toxicity criteria.

Conclusion: Therefore, these compounds have the potential to serve as effective therapeutics for controlling Vibrio alginolyticus infection in shrimp.

Introduction: It has been reported that the extension of conjugation in chalcone scaffolds considerably enhanced the potency, selectivity, reversibility, and competitive mode of MAO-B inhibition. In this study, using the experimental results of IC50 values of fifteen halogenated conjugated dienone derivatives (MK1-MK15) against MAO-B, we developed a 3DQSAR model.

Methods: Further, we created a 3D pharmacophore model in active compounds in the series. The built model selected three variables (G2U, RDF115m, RDF155m) among the 653 AlvaDesc molecular descriptors, with a r² value of 0.87 and a Q² _cv for cross-validation equal to 0.82. The three variables were mostly associated with the direction of symmetry and the likelihood of discovering massive atoms at great distances. The evaluated molecules exhibited a good correlation between experimental and predicted data, indicating that the IC₅₀ value of the structure MK2 was related to the interatomic distances of 15.5 Å between bromine and chloro substituents. Furthermore, the molecules in the series with the highest activity were those with enhanced second component symmetry directional index from the 3D representation, which included the structures MK5 and MK6.

Results: Additionally, a pharmacophore hypothesis was developed and validated using the decoy Schrodinger dataset, with an ROC score of 0.87 and an HHRR 1 fitness score that ranged from 2.783 to 3.00. The MK series exhibited a significant blood-brain barrier (BBB) permeability, according to exploratory analyses and in silico projections, and almost all analogues were expected to have strong BBB permeability.

Conclusion: Further DFT research revealed that electrostatics were important in the interactions with MAO-B.

Objectives: The study aimed to explore the crucial genes involved in cancer-related biological processes, including EMT, autophagy, apoptosis, anoikis, and metastasis. It also sought to identify common genes among the pathways linked to these biological processes, determine the level of Bcl-2 expression in various types of cancers, and find a potent inhibitor of Bcl-2 among natural compounds.

Methods: Common genes involved in the pathways related to EMT, autophagy, apoptosis, anoikis, and metastasis were explored, and the level of the most frequently overexpressed gene that was Bcl-2, in various types of cancers was analyzed by gene expression analysis. A set of 102 natural compounds was sorted according to their docking scores using molecular docking and filtering. The top-ranked molecule was chosen for additional molecular dynamics (MD) simulation for 100 ns. Differential gene expression analysis was performed for Dioscin using GEO2R.

Results: The study identified four common genes, Bcl-2, Bax, BIRC3, and CHUK, among the pathways linked to EMT, autophagy, apoptosis, anoikis, and metastasis. Bcl-2 was highly overexpressed in many cancers, including Acute Myeloid Leukemia, Diffuse large B cell lymphoma, and Thymoma. The Dioscin structure in the Bcl-2 binding site received the highest docking score and the most relevant interactions. Dioscin's determined binding free energy by MM/GBSA was -52.21 kcal/mol, while the same calculated by MM/PBSA was -9.18 kcal/mol. A p-value of less than 0.05 was used to determine the statistical significance of the analysis performed using GEO2R. It was observed that Dioscin downregulates Bcl-2, BIRC3, and CHUK and upregulates the pro-apoptotic protein Bax.

Conclusion: The study concluded that Dioscin has the potential to act as a protein inhibitor, with a noteworthy value of binding free energy and relevant interactions with the Bcl-2 binding site. Dioscin might be a good alternative for targeting multiple cancer pathways through a single target.

Background: Sanguinarine (SAN) has been reported to have antioxidant, antiinflammatory, and antimicrobial activities with potential for the treatment of osteoporosis (OP).

Objective: This work purposed to unravel the molecular mechanisms of SAN in the treatment of OP.

Methods: OP-related genes and SAN-related targets were predicted from public databases. Differential expression analysis and VennDiagram were adopted to detect SAN-related targets against OP. Protein-protein interaction (PPI) network was served for core target identification. Molecular docking and DeepPurpose algorithm were further adopted to investigate the binding ability between core targets and SAN. Gene pathway scoring of these targets was calculated utilizing gene set variation analysis (GSVA). Finally, we explored the effect of SAN on the expressions of core targets in preosteoblastic MC3T3-E1 cells.

Results: A total of 21 candidate targets of SAN against OP were acquired. Furthermore, six core targets were identified, among which CASP3, CTNNB1, and ERBB2 were remarkably differentially expressed in OP and healthy individuals. The binding energies of SAN with CASP3, CTNNB1, and ERBB2 were -6, -6.731, and -7.162 kcal/mol, respectively. Moreover, the GSVA scores of the Wnt/calcium signaling pathway were significantly lower in OP cases than in healthy individuals. In addition, the expression of CASP3 was positively associated with Wnt/calcium signaling pathway. CASP3 and ERBB2 were significantly lower expressed in SAN group than in DMSO group, whereas the expression of CTNNB1 was in contrast.

Conclusion: CASP3, CTNNB1, and ERBB2 emerge as potential targets of SAN in OP prevention and treatment.

Objectives: Cancer poses a great threat to human health, and effective drugs to treat it are always needed. Several compounds containing a 2-aminopyrazine framework have been identified as antitumor agents with SHP2 inhibition activities. This current work aimed to search for more potent novel compounds possessing a 2-aminopyrazine moiety with antitumor activities.

Methods: A series of 12 novel 2-aminopyrazine derivatives was synthesized, and their structures were confirmed by spectroscopic techniques. The inhibitory activities of all the synthesized compounds against MDA-MB-231 and H1975 cancer cell lines were evaluated by an MTT assay. The most potent compound 3e was analyzed by flow cytometry. Subsequently, computational studies were performed to investigate the possible antitumor mechanisms of compound 3e.

Results: The results indicated that compound 3e exhibited potent antitumor activities with IC50 values of 11.84 ± 0.83 μM against H1975 cells and 5.66 ± 2.39 μM against MDA-MB-231 cells, which were more potent than the SHP2 inhibitor GS493 (IC₅₀ = 19.08 ± 1.01 μM against H1975 cells and IC₅₀ = 25.02 ± 1.47 μM against MDA-MB-231 cells). Further analysis by flow cytometry demonstrated that compound 3e induced cell apoptosis in H1975 cells. The results of the molecular docking and MD simulations, including RMSD, RMSF, PCA, DCCM and binding energy and decomposition analyses, revealed that compound 3e probably selectively inhibited SHP2.

Conclusion: A new compound having a 2-aminopyrazine substructure with potent inhibitory activities against the H1975 and MDA-MB-231 cancer cells was obtained, meriting further investigation as an antitumor drug.

Background: Blocking Kv 1.2 and Kv 1.3 potassium channels using scorpion venom- derived toxins holds potential therapeutic value. These channels are implicated in autoimmune diseases such as neurodegenerative diseases, multiple sclerosis, rheumatoid arthritis, and type 1 diabetes.

Objectives: The present work aims at the discovery and in silico activity analysis of potassium channel blockers (KTxs) from the cDNA library derived from the venom gland of Iranian scorpion Hemiscorpius lepturus (H. lepturus).

Methods: The sequence regarding potassium channel blockers were extracted based on Gene Ontology for H. lepturus venom gland. Homology analyses, superfamily, family, and evolutionary signatures of H. lepturus KTxs (H.L KTxs) were determined by using BLASTP, COBALT, PROSITE, and InterPro servers. The predicted 3D structures of H.L KTxs were superimposed against their homologs to predict structure activity relationship. Molecular docking analysis was also performed to predict the binding affinity of H.L KTxs to Kv 1.2 and Kv 1.3 channels. Finally, the toxicity was predicted.

Results: Seven H.L KTxs, designated as Leptukalin, were extracted from the cDNA library of H. lepturus venom gland. Homology analyses proved that they can act as potassium channel blockers and they belong to the superfamily and family of Scorpion Toxin-like and Short-chain scorpion toxins, respectively. Structural alignment results confirmed the activity of H.L KTxs. Binding affinity of all H.L KTxs to Kv 1.2 and Kv 1.3 channels ranged from -4.4 to -5.5 and -4 to -5.7 Kcal/mol, respectively. In silico toxicity assay showed that Leptukalin 3, Leptukalin 5, and Leptukalin 7 were non-toxic.

Conclusion: Three non-toxic KTxs, Leptukalin 3, 5, and 7, were successfully discovered from the cDNA library of H. lepturus venom gland. Gathering all data together, the discovered peptides are promising potassium channel blockers. Accordingly, Leptukalin 3, 5, and 7 could be suggested for complementary in vitro studies and mouse model of autoimmune diseases.

Background: Acinetobacter baumannii is one of the main causes of nosocomial infections. No vaccine has yet been licensed for use in humans, and efforts are still ongoing.

Objective: In the present study, we have predicted the B-cell epitopes of A. baumannii's outer membrane protein K (OMPK) by using epitope prediction algorithms as possible vaccine candidates for future studies.

Methods: The linear B-cell epitopes were predicted by seven different prediction tools. The 3D structure of OMPK was modeled and used for discontinuous epitope prediction by ElliPro and DiscoTope 2.0 tools. The final linear epitopes and the discontinuous epitope segments were checked for potential allergenicity, toxicity, human similarity, and experimental records. The structure and physicochemical features of the final epitopic peptide were assessed by numerous bioinformatics tools.

Results: Many B-cell epitopes were detected that could be assessed for possible antigenicity and immunogenicity. Also, an epitopic 22-mer region (peptide) of OMPK was found that contained both linear and discontinuous B-cell epitopes. This epitopic peptide has been found to possess appropriate physicochemical and structural properties to be an A. baumannii vaccine candidate.

Conclusion: Altogether, here, the high immunogenic B-cell epitopes of OMPK have been identified, and a high immunogenic 22-mer peptide as an A. baumannii vaccine candidate has been introduced. The in vitro/in vivo studies of this peptide are recommended to decide its real efficacy and efficiency.

Aims and objectives: This study aimed to evaluate the pharmacological mechanism of Hederagenin (HD) combined with oxaliplatin (L-OHP) in treating gastric cancer (GC) through network pharmacology combined with experimental verification.

Material and methods: Network pharmacology methods were used to screen potential targets for HD, L-OHP, and GC-related targets from public databases, and the intersection of the three gene sets was taken. Cross genes were analyzed through protein-protein interaction (PPI) networks to predict core targets, and related pathways were predicted through GO and KEGG enrichment analysis. The experimental results were verified by the in vitro experiments. HD was applied on AGS/L-OHP cells, and then cellular chemosensitivity and the expressions of P-gp, Survivin, Bcl-2, p-Akt, and p-PI3K genes were detected. Wound assay and Transwell Chamber assay were employed to detect the effect of HD on AGS/L-OHP cells. Nude mice xenograft models transfected using AGS/L-OHP cells were also treated with HD in order to verify the results. The size and weight of the tumor, as well as the expressions of P-gp, Survivin, Bcl-2, p- Akt and p-PI3K genes, were also measured.

Results: KEGG analysis showed that the anti-gastric cancer effect of HD was mediated mainly by PI3K-Akt signaling pathways. The PI3K-Akt signaling pathway containing more enriched genes may play a greater role in anti-gastric cancer. It was observed that for AGS/L-OHP cells jointly treated with HD and L-OHP, their activity, migration and invasion were significantly lower than those treated only using HD or L-OHP group. Moreover, expressions of p-Akt, p- PI3K, Bcl-2, P-gp, and Survivin for the HD+L-OHP group decreased significantly. Results of the in vivo experiments showed that the sizes and weights of tumors in the HD+L-OHP group were the lowest compared to the HD group and L-OHP group.

Conclusion: Our findings suggest that HD may reduce the resistance of AGS/L-OHP cells to LOHP by regulating the PI3K/Akt signaling pathway.