In search of AKT kinase inhibitors as anticancer agents: structure-based design,docking, and molecular dynamics studies of 2,4,6-trisubstituted pyridines

The AKT isoforms are a group of key kinases that play a critical role in tumorigenesis. These enzymes are overexpressed in different types of cancers, such as breast, colon, prostate, ovarian, and lung. Because of its relevance the AKT isoforms are attractive targets for the design of anticancer molecules. However, it has been found that AKT1 and AKT3 isoforms have a main role in tumor progression and metastasis; thus, the identification of AKT isoforms specific inhibitors seems to be a challenge. Previously, we identified an ATP binding pocket pan-AKT inhibitor, this compound is a 2,4,6-trisubstituted pyridine (compound 11), which represents a new interesting scaffold for the developing of AKT inhibitors. Starting from the 2,4,6-trisubstituted pyridine scaffold, and guided by structure-based design technique, 42 new inhibitors were designed and further evaluated in the three AKT isoforms by multiple docking approach and molecular dynamics. Results showed that seven compounds presented binding selectivity for AKT1 and AKT3, better than for AKT2. The binding affinities of these seven compounds on AKT1 and AKT3 isoforms were mainly determined by hydrophobic contributions between the aromatic portion at position 4 of the pyridine ring with residues Phe236/234, Phe237/235, Phe438/435, and Phe442/439 in the ATP binding pocket. Results presented in this work provide an addition knowledge leading to promising selective AKT inhibitors.     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Designing of benzothiazole derivatives as promising EGFR tyrosine kinase inhibitors: a pharmacoinformatics study

Abstract

Benzothiazole derivatives represent an important class of therapeutic chemical agents and are widely used for interesting biological activities and therapeutic functions including anticancer, antitumor and antimicrobial. In this study, we have performed similarity/substructure-based search of eMolecule database to find out promising benzothiazole derivatives as EGFR tyrosine kinase inhibitors. Several screening criteria that included molecular docking, pharmacokinetics and synthetic accessibility were used on initially derived about 7000 molecules consisting of benzothiazole as major component. Finally, four molecules were found to be promising EGFR tyrosine kinase inhibitors. The best docked pose of each molecule was considered for binding interactions followed by molecular dynamics (MD) and binding energy calculation. Molecular docking clearly showed the final proposed derivatives potential to form a number of binding interactions. MD simulation trajectories undoubtedly indicated that the EGFR protein becomes stable when proposed derivatives bind to the receptor cavity. Strong binding affinity was found for all molecules toward the EGFR which was substantiated by the binding energy calculation using the MM-PBSA approach. Therefore, proposed benzothiazole derivatives may be promising EGFR tyrosine kinase inhibitors for potential application as cancer therapy.

Communicated by Ramaswamy H. Sarma     

Validation of potential inhibitors for SrtA against Bacillus anthracis by combined approach of ligand-based and molecular dynamics simulation

The development of SrtA inhibitors targeting the biothreat organism namely Bacillus anthracis was achieved by the combined approach of pharmacophore modeling, binding interactions, electron transferring capacity, ADME, and Molecular dynamics studies. In this study, experimentally reported Ba-SrtA inhibitors (pyridazinone and pyrazolethione derivatives) were considered for the development of enhanced pharmacophoric model. The obtained AAAHR hypothesis was a pure theoretical concept that accounts for common molecular interaction network present in experimentally active pyridazinone and pyrazolethione derivatives. Pharmacophore-based screening of AAAHR hypothesis provides several new compounds, and those compounds were treated with four phases of docking protocols with combined Glide-QPLD docking approach. In this approach, scoring and charge accuracy variations were seen to be dominated by QM/MM approach through the allocation of partial charges. Finally, we reported the best compounds from binding db, Chembridge db, and Toslab based on scoring values, energy parameters, electron transfer reaction, ADME, and cell adhesion inhibition activity. The dynamic state of interaction and binding energy assess that new compounds are more active inside the binding pocket and these compounds on experimental validations will survive as better inhibitors for targeting the cell adhesion mechanism of Ba-SrtA.     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

Molecular docking, 3D-QSAR,molecular dynamics,synthesis and anticancer activity of tyrosine kinase 2 (TYK 2) inhibitors

Abstract

A therapeutic rationale is proposed by selectively targeting tyrosine kinase 2 (TYK 2) to obtain potent TYK 2 inhibitors by molecular modeling studies. In the present study, we have taken tyrosine kinase (TYK 2) inhibitors and carried out molecular docking, 3?D quantitative structure–activity relationship (3D-QSAR) analysis and molecular dynamics (MD). Based on the 3D-QSAR results thirteen new compounds (R-1 to R-13) were designed and synthesized in good yields. The synthesized molecules were evaluated for their in vitro anticancer activity against LnCap and A549 cell lines. The molecules R-1, R-3, R-5, R-7, and R-10 exhibited considerable anti cancer activity.     

Identification of natural compound inhibitors against PfDXR: A hybrid structure-based molecular modeling approach and molecular dynamics simulation studies

In the present contribution, multicomplex-based pharmacophore studies were carried out on the structural proteome of Plasmodium falciparum 1-deoxy-D -xylulose-5-phosphate reductoisomerase. Among the constructed models, a representative model with complementary features, accountable for the inhibition was used as a primary filter for the screening of database molecules. Auxiliary evaluations of the screened molecules were performed via drug-likeness and molecular docking studies. Subsequently, the stability of the docked inhibitors was envisioned by molecular dynamics simulations, principle component analysis, and molecular mechanics-Poisson-Boltzmann surface area-based free binding energy calculations. The stability assessment of the hits was done by comparing with the reference (beta-substituted fosmidomycin analog, LC5) to prioritize more potent candidates. All the complexes showed stable dynamic behavior while three of them displayed higher binding free energy compared with the reference. The work resulted in the identification of the compounds with diverse scaffolds, which could be used as initial leads for the design of novel PfDXR inhibitors.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．     

Insight into ligand selectivity in HCV NS5B polymerase: molecular dynamics simulations, free energy decomposition and docking

Modeling studies were performed on HCV NS5B polymerase in an effort to design new inhibitors. The binding models of five different scaffold inhibitors were investigated and compared by using molecular dynamics simulations, free energy calculation and decomposition. Our results show Tyr448 plays the most critical role in the binding of most inhibitors. In addition, favorable contributions of residues Pro197, Arg200, Cys366, Met414 and Tyr448 in a deep hydrophobic pocket prove to be important for the selectivity of inhibitors. Furthermore, an optimized docking protocol was presented based on cross-docking the five inhibitors in the palm binding site of this enzyme using the Autodock program. This protocol was used later to virtually screen NCI and Maybridge diversity set libraries. The binding site was profiled via the statistics and analysis of the hydrogen bond networks formed between the receptor and the top-ranked diversity set compounds. Based on our detailed binding site analysis two useful rules were proposed to guide the selection of promising hits.     

Virtual screening against metalloenzymes for inhibitors and substrates

Molecular docking uses the three-dimensional structure of a receptor to screen databases of small molecules for potential ligands, often based on energetic complementarity. For many docking scoring functions, which calculate nonbonded interactions, metalloenzymes are challenging because of the partial covalent nature of metal-ligand interactions. To investigate how well molecular docking can identify potential ligands of metalloenzymes using a "standard" scoring function, we have docked the MDL Drug Data Report (MDDR), a functionally annotated database of 95,000 small molecules, against the X-ray crystal structures of five metalloenzymes. These enzymes included three zinc proteases, the nickel analogue of an iron enzyme, and a molybdenum metalloenzyme. The ability of the docking program to retrospectively enrich the annotated ligands as high-scoring hits for each enzyme and to calculate proper geometries was evaluated. In all five systems, the annotated ligands within the MDDR were enriched at least 20 times over random. To test the approach prospectively, a sixth target, the zinc beta-lactamase from Bacteroides fragilis, was screened against the fragment-like subset of the ZINC database. We purchased and tested 15 compounds from among the top 50 top-ranked ligands from docking, and found 5 inhibitors with apparent K(i) values less than 120 microM, the best of which was 2 microM. A more ambitious test still was predicting actual substrates for a seventh target, a Zn-dependent phosphotriesterase from Pseudomonas diminuta. Screening the Available Chemicals Directory (ACD) identified 25 thiophosphate esters as potential substrates within the top 100 ranked compounds. Eight of these, all previously uncharacterized for this enzyme, were acquired and tested, and all were confirmed experimentally as substrates. These results suggest that a simple, noncovalent scoring function may be used to identify inhibitors of at least some metalloenzymes.     

DFT study of the Lewis acid mediated synthesis of 3-acyltetramic acids

The recent emergence and re-emergence of alphaviruses, in particular the chikungunya virus (CHIKV), in numerous countries has invoked a worldwide threat to human health, while simultaneously generating an economic burden on affected countries. There are currently no vaccines or effective drugs available for the treatment of the CHIKV, and with few lead compounds reported, the vital medicinal chemistry is significantly more challenging. This study reports on the discovery of potential inhibitors for the nsP3 macro domain of CHIKV using molecular docking, virtual screening, and molecular dynamics simulations, as well as work done to evaluate and confirm the active site of nsP3. Virtual screening was carried out based on blind docking as well as focused docking, using the database of 1541 compounds from NCI Diversity Set II, to identify hit compounds for nsP3. The top hit compounds were further subjected to molecular dynamic simulations, yielding a greater understanding of the dynamic behavior of nsP3 and its complexes with various ligands, concurrently confirming the outcomes of docking, and establishing in silico lead compounds which target the CHIKV nsP3 enzyme.

Prediction of binding for a kind of non-peptic HCV NS3 serine protease inhibitors from plants by molecular docking and MM-PBSA method

In the study, molecular dynamics simulations combined with MM-PBSA (Molecular Mechanics and Poisson-Boltzmann Surface Area) technique were applied to predict the binding mode of the polyphenol inhibitor in the binding pocket of the HCV NS3 serine protease for which the ligand-protein crystal structure is not available. The most favorable geometry of three candidates from molecular docking had a binding free energy about 3 and 6kcal/mol more favorable than the other two candidates, respectively, and was identified as the correct binding mode. In the mode, the correlation of the calculated and experimental binding affinities of all five polyphenol compounds is satisfactory indicated by r(2)=0.92. The most favorable binding mode suggests that two galloyl residues at 3 and 4 positions of the glucopyranose ring of the inhibitors interact with SER139, GLY137, ALA157, and ASP81 by hydrogen bond interaction and with ALA156 and HIE57 by hydrophobic interaction and are essential for the activities of the studied inhibitors.     

Molecular epidemiology of antibiotic-resistant genes and potent inhibitors against TEM,CTX-M-14, CTX-M-15, and SHV-1 proteins of Escherichia coli in district Peshawar,Pakistan

The Extended Spectrum Beta-Lactamases (ESBLs) producing bacteria is an issue of concern for clinicians resulting in minimize the treatment options. To overcome resistance mechanisms, novel inhibitors with good Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties must inhibit the ESBLs resistant genes. The current study aimed to identify the antibiotic resistance genes of ESBLs producing E. coli and a single inhibitor was designed to inhibit all the resistant proteins. The results showed that 42.9% ESBL producers had CTX-M (69.9%), TEM (63.4%), SHV (34.5%) and CTX-M-14 (17.5%) genes. The ESBLs producing isolates were resistant to cephalosporins, quinolones, and sulfonamide with Minimum Inhibitory Concentration (MICs) ranging from 64 to >256 μg/ml. To design multi inhibitory ligands, RECAP synthesis was used for the de-novo discovery of 1000 inhibitors database. Protein crystal structures were retrieved from Protein Data Base (PDB). Lipinski’s rules of five were applied to the novel inhibitors database to improve the ADMET properties. The novel inhibitors database was selected for docking simulations. Placement of the ligand was used by the London dG algorithm implemented in Molecular Operating Environment (MOE), while GBVI/WSA dG algorithm was used for final refinement. Based on docking score, visual inspection of ligands interaction with key residues, binding affinity, and binding energy of ligands with proteins, ten compounds were selected for ESBLs proteins with best ADMET properties, binding energy, and binding affinity the reported ones. These hits compounds have unique scaffolds and are predicted to be a starting point for developing potent inhibitors against antibiotic-resistant proteins.     

N-(thiazol-2-yl)-2-thiophene carboxamide derivatives as Abl inhibitors identified by a pharmacophore-based database screening of commercially available compounds

Suggestions derived from a previous ligand-based ligand design approach and docking calculations aimed at finding compound with affinity toward Abl and molecular scaffolds previously untested as Abl inhibitors, led to the identification of commercially available N-(thiazol-2-yl)-2-thiophene carboxamide derivatives with affinity in a cell-free assay up to low nanomolar concentrations, significantly enhanced with respect to that of their parent compounds previously reported. In particular, among compounds of the Asinex database, molecular docking simulations guided the choice of high-affinity ligands, predicting their binding mode and their interaction pattern with the Abl catalytic binding site. Moreover, affinity of the new compounds was also rationalized in terms of their interactions with the enzyme.     

Screening of hepatitis C NS5B polymerase inhibitors containing benzothiadiazine core: a steered molecular dynamics

Hepatic C virus (HCV) is a global health problem, resulting in liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and fatality. The NS5B polymerase of HCV plays an important role in viral RNA replication process, making it an attractive therapeutic target for design and development of anti-HCV drugs. To search new potent compounds against the HCV NS5B polymerase, the molecular docking and the steered molecular dynamics (SMD) simulation techniques were performed. The potential potent inhibitors of the NS5B polymerase were screened out from the ZINC database using structural similarity search and molecular docking technique. Five top-hit compounds (the ZINC compounds 49888724, 49054741, 49777239, 49793673, and 49780355) were then studied by the SMD simulations based on the hypothesis that a high rupture force relates to a high binding efficiency. The results demonstrated that the ZINC compound 49888724 had a greater maximum rupture force, reflecting a good binding strength and inhibitory potency than known inhibitors and the rest four ZINC compounds. Therefore, our finding indicated that the ZINC compound 49888724 is a potential candidate to be a novel NS5B inhibitor for further design. Besides, the van der Waals interaction could be considered as the main contribution for stabilizing the NS5B-ligand complex.     

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.     

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca2? dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 ?, 5.5 ?, 6 ?, 6.5 ?, 7.0 ? and 7.5 ? were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.     

Identification of molecular fingerprints of natural products for the inhibition of breast cancer resistance protein (BCRP)

BackgroundExtensive research over the past several decades, explored that the natural compounds contain different plant secondary metabolites and have the potential to inhibit breast cancer resistance protein (BCRP).PurposeTo identify crucial molecular fingerprints of some natural products for the inhibition of breast cancer resistance protein and also to screen out some potent natural BCRP inhibitors.Study designMultiple modelling strategies were applied with three main mottos: (a) Generation of robust classification models to identify the linear and non-linear relationships among the natural compounds and the inhibition of BCRP, (b) Identification of important structural fingerprints that modulate BCRP inhibition and screening of natural database to find the probable hit molecules, (c) Comprehensive ligand-receptor interactions analysis of those against the putative breast cancer resistant protein through molecular docking analysis.MethodsMonte Carlo optimization and SPCI analysis was used to identify important structural fingerprints. QSARCo. and swissADME analysis were used for screening and prediction of hits. Finally, docking analysis was performed for interaction study.ResultsIn this study, some important structural fingerprints of BCRP inhibitors were identified. Additionally, eleven natural anti-cancer compounds were predicted to be active against the BCRP and also satisfy the different drug-likeliness properties. Among them, apigenin was found to have better binding affinities against the putative target as obtained from molecular docking analysis.ConclusionThis study is an attempt to understand about the molecular fingerprints of natural compounds for the inhibition of BCRP and also to dig out some novel natural inhibitors against BCRP.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.