New potential inhibitors of mTOR: a computational investigation integrating molecular docking,virtual screening and molecular dynamics simulation

The mTOR (mammalian or mechanistic Target Of Rapamycin), a complex metabolic pathway that involves multiple steps and regulators, is a major human metabolic pathway responsible for cell growth control in response to multiple factors and that is dysregulated in various types of cancer. The classical inhibition of the mTOR pathway is performed by rapamycin and its analogs (rapalogs). Considering that rapamycin binds to an allosteric site and performs a crucial role in the inhibition of the mTOR complex without causing the deleterious side effects common to ATP-competitive inhibitors, we employ ligand-based drug design strategies, such as virtual screening methodology, computational determination of ADME/Tox properties of selected molecules, and molecular dynamics in order to select molecules with the potential to become non-ATP-competitive inhibitors of the mTOR enzymatic complex. Our findings suggest five novel potential mTOR inhibitors, with similar or better properties than the classic inhibitor complex, rapamycin.     

Hydroxamates of para-aminobenzoic acid as selective inhibitors of HDAC8

A series of hydroxamates (4a–4l) were prepared from p-aminobenzoic acid to inhibit HDAC8. The idea is to substitute rigid aromatic ring in place of less rigid piperazine ring of hydroxamates reported earlier by our group. It is expected to increase potency retaining the selectivity. Result obtained suggested that the modifications carried out retained the selectivity towards HDAC8 isoform and increasing the potency in very few cases. Increase in potency is also associated with variation in cap aryl region. Two compounds (4f & 4l) were found to inhibit HDAC8 at concentrations (IC₅₀) less than 20 μM.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Synthesis and molecular docking studies of some 4-phthalimidobenzenesulfonamide derivatives as acetylcholinesterase and butyrylcholinesterase inhibitors

A series of 4-phthalimidobenzenesulfonamide derivatives were designed, synthesized and evaluated for the inhibitory activities against acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). Structures of the title compounds were confirmed by spectral and elemental analyses. The cholinesterase (ChE) inhibitory activity studies were carried out using Ellman’s colorimetric method. The biological activity results revealed that all of the title compounds (except for compound 8) displayed high selectivity against AChE. Among the tested compounds, compound 7 was found to be the most potent against AChE (IC₅₀=?1.35?±?0.08?μM), while compound 3 exhibited the highest inhibition against BuChE (IC₅₀=?13.41?±?0.62?μM). Molecular docking studies of the most active compound 7 in AChE showed that this compound can interact with both the catalytic active site (CAS) and the peripheral anionic site (PAS) of AChE.     

Shape-based virtual screening,docking, and molecular dynamics simulations to identify Mtb-ASADH inhibitors

Aspartate β-semialdehyde dehydrogenase (ASADH) is a key enzyme for the biosynthesis of essential amino acids and several important metabolites in microbes. Inhibition of ASADH enzyme is a promising drug target strategy against Mycobacterium tuberculosis (Mtb). In this work, in silico approach was used to identify potent inhibitors of Mtb-ASADH. Aspartyl β-difluorophosphonate (β-AFP), a known lead compound, was used to understand the molecular recognition interactions (using molecular docking and molecular dynamics analysis). This analysis helped in validating the computational protocol and established the participation of Arg99, Glu224, Cys130, Arg249, and His256 amino acids as the key amino acids in stabilizing ligand–enzyme interactions for effective binding, an essential feature is H-bonding interactions with the two arginyl residues at the two ends of the ligand. Best binding conformation of β-AFP was selected as a template for shape-based virtual screening (ZINC and NCI databases) to identify compounds that competitively inhibit the Mtb-ASADH. The top rank hits were further subjected to ADME and toxicity filters. Final filter was based on molecular docking analysis. Each screened molecule carries the characteristics of the highly electronegative groups on both sides separated by an average distance of 6?Å. Finally, the best predicted 20 compounds exhibited minimum three H-bonding interactions with Arg99 and Arg249. These identified hits can be further used for designing the more potent inhibitors against ASADH family. MD simulations were also performed on two selected compounds (NSC4862 and ZINC02534243) for further validation. During the MD simulations, both compounds showed same H-bonding interactions and remained bound to key active residues of Mtb-ASADH.     

Developing imidazole analogues as potential inhibitor for Leishmania donovani trypanothione reductase: virtual screening,molecular docking,dynamics and ADMET approach

Visceral leishmaniasis (VL) affects Indian subcontinent, African and South American continent, and it covers 70 countries worldwide. Visceral form of leishmaniasis is caused by Leishmania donovani in Indian subcontinent which is lethal if left untreated. Extensive resistance to antileishmanial drugs such as sodium stibogluconate, pentamidine and miltefosine and their decreased efficacy has been reported in the endemic region. Amphotericin B drug has shown good antileishmanial activity with significant toxicity, but its cost of treatment has limited the outreach of this treatment to affected people living in endemic zone. So, there is an urgent need to identify new antileishmanial drugs with excellent activity and minimal toxicity issues. Trypanothione reductase, a component of antioxidant system, is necessary for parasite growth and survival to raise infection. To develop potential inhibitor, we docked nine hundred and eighty-four 5-nitroimidazole analogues along with clomipramine which is a well-known inhibitor for TR. Total one hundred and forty-seven 5-nitroimidazole analogues with better docking score than clomipramine were chosen for ADMET and QikProp studies. Among these imidazole analogues, total twenty-four imidazole analogues and clomipramine were chosen on the basis of their ADMET, QikProp, and prime MM-GBSA study. Later on, two analogues with best MM-GBSA dG bind were undergone molecular dynamic simulation to ensure protein–ligand interactions. Using above approach, we confirm that ethyl 2-acetyl-5-[4-butyl-2-(3-hydroxypentyl)-5-nitro-1H-imidazol-1-yl]pent-2-enoate can be a drug candidate against L. donovani for the treatment of VL in the Indian subcontinent.     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

Molecular docking studies on DMDP derivatives as human DHFR inhibitors

Molecular docking is routinely used for understanding drug‐receptor interaction in modern drug design. Here, we describe the docking of 2, 4-diamino-5-methyl-5-deazapteridine (DMDP) derivatives as inhibitors to human dihydrofolate reductase (DHFR). We docked 78 DMDP derivates collected from literature to DHFR and studied their specific interactions with DHFR. A new shape-based method, LigandFit, was used for docking DMDP derivatives into DHFR active sites. The result indicates that the molecular docking approach is reliable and produces a good correlation coefficient (r² = 0.499) for the 73 compounds between docking score and IC₅₀ values (Inhibitory Activity). The chloro substituted naphthyl ring of compound 63 makes significant hydrophobic contact with Leu 22, Phe 31 and Pro 61 of the DHFR active site leading to enhanced inhibition of the enzyme. The docked complexes provide better insights to design more potent DHFR inhibitors prior to their synthesis.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Exploring the binding mechanism of HDAC8 selective inhibitors: Lessons from the modification of Cap group

The abnormal expression of histone deacetylase 8 (HDAC8) has been reported to associate with various cancer entities (colon, breast cancer, pancreas, etc.) as well as parasitic diseases, making HDAC8 gradually develop into an attractive and potential therapeutic target. Among the various design strategies of selective HDAC8 inhibitors (modification of Cap, Linker, or zinc binding group regions), the optimization of Cap region has aroused great interest among the researchers. However, the detailed information underlying how the modification of Cap region influences the inhibitory activities is still unclear, and in this study, compounds 2c, 3g, and 3n were selected to explore the differences in binding mechanisms brought by Cap modifications via various computational approaches at the atomic level. Five residues (Y293, H167, D254, D165, and M261) have a large difference in energy contributions to the constructed systems, and the subpocket formed by Y293 and M261 could interact with Cap groups, triggering the differences in the energy contributions of the residues (H167, D254, and D165) located in metal-catalytic center. In summary, the compounds 2c, 3g, and 3n were selected as molecular probes to explore the binding mechanism, and the residues (Y293 and M261) forming the subpocket should be paid special attention in the design and synthesis of novel selective HDAC8 inhibitors.     

Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.     

Synthesis of quinoline derivatives as diabetic II inhibitors and molecular docking studies

In searchof the potenttherapeutic agent as an α-glucosidase inhibitor, we have synthesized twenty-five analogs (1–25) of quinoline-based Schiff bases as an inhibitoragainst α-glucosidase enzyme under positive control acarbose (IC₅₀ = 38.45 ± 0.80 µM). From the activity profile it was foundthat analogs 1, 2, 3, 4, 11, 12 and 20with IC₅₀values 12.40 ± 0.40, 9.40 ± 0.30, 14.10 ± 0.40, 6.20 ± 0.30, 14.40 ± 0.40, 7.40 ± 0.20 and 13.20 ± 0.40 µMrespectively showed most potent inhibition among the series even than standard drug acarbose (IC₅₀ = 38.45 ± 0.80 µM). Here in the present study analog 4 (IC₅₀ = 6.20 ± 0.30 µM) was found with many folds better α-glucosidase inhibitory activity than the reference drug. Eight analogs like 5, 7, 8, 16, 17, 22, 24 and 25 among the whole series displayed less than 50% inhibition. The substituents effects on phenyl ring thereby superficially established through SAR study. Binding interactions of analogs and the active site of ligands proteins were confirmed through molecular docking study. Spectroscopic techniques like ¹H NMR, ¹³C NMR and ESIMS were used for characterization.     

Synthesis,molecular docking,binding free energy calculation and molecular dynamics simulation studies of benzothiazol-2-ylcarbamodithioates as Staphylococcus aureus MurD inhibitors

Abstract

A new series of benzothiazol-2-ylcarbamodithioate functional compounds 5a-f has been designed, synthesized and characterized by spectral data. These compounds were screened for their in vitro antibacterial activity against strains of Staphylococcus aureus (NCIM 5021, NCIM 5022 and methicillin-resistant isolate 43300), Bacillus subtilis (NCIM 2545), Escherichia coli (NCIM 2567), Klebsiella pneumoniae (NCIM 2706) and Psudomonas aeruginosa (NCIM 2036). Compounds 5a and 5d exhibited significant activity against all the tested bacterial strains. Specifically, compounds 5a and 5d showed potent activity against K. pneumoniae (NCIM 2706), while compound 5a also displayed potent activity against S. aureus (NCIM 5021). Compound 5d showed minimum IC₅₀ value of 13.37?μM against S. aureus MurD enzyme. Further, the binding interactions of compounds 5a-f in the catalytic pocket have been investigated using the extra-precision molecular docking and binding free energy calculation by MM-GBSA approach. A 30?ns molecular dynamics simulation of 5d/modeled S. aureus MurD enzyme was performed to determine the stability of the predicted binding conformation.     

Pharmacophore modeling,molecular docking and molecular dynamics studies on natural products database to discover novel skeleton as non-purine xanthine oxidase inhibitors

Gout is a common inflammatory arthritis caused by the deposition of urate crystals within joints. It is increasingly in prevalence during the past few decades as shown by the epidemiological survey results. Xanthine oxidase (XO) is a key enzyme to transfer hypoxanthine and xanthine to uric acid, whose overproduction leads to gout. Therefore, inhibiting the activity of xanthine oxidase is an important way to reduce the production of urate. In the study, in order to identify the potential natural products targeting XO, pharmacophore modeling was employed to filter databases. Here, two methods, pharmacophore based on ligand and pharmacophore based on receptor-ligand, were constructed by Discovery Studio. Then GOLD was used to refine the potential compounds with higher fitness scores. Finally, molecular docking and dynamics simulations were employed to analyze the interactions between compounds and protein. The best hypothesis was set as a 3D query to screen database, returning 785 and 297 compounds respectively. A merged set of the above 1082 molecules was subjected to molecular docking, which returned 144 hits with high-fitness scores. These molecules were clustered in four main kinds depending on different backbones. What is more, molecular docking showed that the representative compounds established key interactions with the amino acid residues in the protein, and the RMSD and RMSF of molecular dynamics results showed that these compounds can stabilize the protein. The information represented in the study confirmed previous reports. And it may assist to discover and design new backbones as potential XO inhibitors based on natural products.     

Design,docking studies and molecular dynamics of new potential selective inhibitors of Plasmodium falciparum serine hydroxymethyltransferase

In this paper, former studies on the interactions of the natural substrate and potential inhibitors of Plasmodium falciparum serine hydroxymethyltransferase (PfSHMT) were used to design five new potential selective inhibitors to this enzyme. Results of the docking energies calculations of these structures inside the active sites of PfSHMT and human SHMT were used to select a more suitable structure as a potential selective inhibitor to PfSHMT. Further molecular dynamics studies of this molecule and 5-formyl-6-hydrofolic acid (natural substrate) docked inside these enzymes' active sites revealed important features for additional refinements of this structure and also additional residues in the PfSHMT active site to be considered further for designing selective inhibitors.     

Quercetin derivatives as non-nucleoside inhibitors for dengue polymerase: molecular docking,molecular dynamics simulation,and binding free energy calculation

Dengue is an important public health problem in tropical and subtropical regions of the world. Neither vaccine nor an antiviral medication is available to treat dengue. This insists the need of drug discovery for dengue. In order to find a potent lead molecule, RNA-dependent RNA polymerase which is essential for dengue viral replication is chosen as a drug target. As Quercetin showed antiviral activity against several viruses, quercetin derivatives developed by combinatorial library synthesis and mined from PubChem databases were screened for a potent anti-dengue viral agent. Our study predicted Quercetin 3-(6″-(E)-p-coumaroylsophoroside)-7-rhamnoside as a dengue polymerase inhibitor. The results were validated by molecular dynamics simulation studies which reveal water bridges and hydrogen bonds as major contributors for the stability of the polymerase-lead complex. Interactions formed by this compound with residues Trp795, Arg792 and Glu351 are found to be essential for the stability of the polymerase-lead complex. Our study demonstrates Quercetin 3-(6″-(E)-p-coumaroylsophoroside)-7-rhamnoside as a potent non-nucleoside inhibitor for dengue polymerase.     

Novel pyrrolopyrimidine derivatives: design,synthesis, molecular docking,molecular simulations and biological evaluations as antioxidant and anti-inflammatory agents

Current medical approaches to control the Covid-19 pandemic are either to directly target the SARS-CoV-2 via innovate a defined drug and a safe vaccine or indirectly target the medical complications of the virus. One of the indirect strategies for fighting this virus has been mainly dependent on using anti‐inflammatory drugs to control cytokines storm responsible for severe health complications. We revealed the discovery of novel fused pyrrolopyrimidine derivatives as promising antioxidant and anti-inflammatory agents. The newly synthesised compounds were evaluated for their in vitro anti-inflammatory activity using RAW264.7 cells after stimulation with lipopolysaccharides (LPS). The results revealed that 3a, 4b, and 8e were the most potent analogues. Molecular docking and simulations of these compounds against COX-2, TLR-2 and TLR-4 respectively was performed. The former results were in line with the biological data and proved that 3a, 4b and 8e have potential antioxidant and anti-inflammatory effects.     

Structure-based virtual screening,molecular docking,ADMET and molecular simulations to develop benzoxaborole analogs as potential inhibitor against Leishmania donovani trypanothione reductase

Visceral leishmaniasis (VL) is the most fatal form of leishmaniasis and it affects 70 countries worldwide. Increasing drug resistant for antileishmanial drugs such as miltefosine, sodium stibogluconate and pentamidine has been reported in the VL endemic region. Amphotericin B has shown potential antileishmanial activity in different formulations but its cost of treatment and associated nephrotoxicity have limited its use by affected people living in the endemic zone. To control the VL infection in the affected countries, it is necessary to develop new antileishmanial compounds with high efficacy and negligible toxicity. Computer aided programs such as binding free energy estimation; ADMET prediction and molecular dynamics simulation can be used to investigate novel antileishmanial molecules in shorter duration. To develop antileishmanial lead molecule, we performed standard precision (SP) docking for 1160 benzoxaborole analogs along with reference inhibitors against trypanothione reductase of Leishmania parasite. Furthermore, extra precision (XP) docking, ADMET prediction, prime MM-GBSA was conducted over 115 ligands, showing better docking score than reference inhibitors to get potential antileishmanial compounds. Simultaneously, area under the curve (AUC) was estimated using ROC plot to validate the SP and XP docking protocol. Later on, two benzoxaborole analogs with best MM-GBSA ΔG-bind were subjected to molecular simulation and docking confirmation to ensure the ligand interaction with TR. The presented drug discovery based on computational study confirms that BOB27 can be used as a potential drug candidate and warrants further experimental investigation to fight against VL in endemic areas.