Repurposing of FDA-approved drugs to target MurB and MurE enzymes in Mycobacterium tuberculosis

Abstract

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb) is one amongst the top 10 causes of death worldwide. The growing rise in antibiotic resistance compounded with slow and expensive drug discovery has further aggravated the situation. ‘Drug repurposing’ is a promising approach where known drugs are examined for a new indication. In the present study, we have attempted to identify drugs that could target MurB and MurE enzymes involved in the muramic acid synthesis pathway (Mur Pathway) in Mtb. FDA-approved drugs from two repositories i.e. Drug Bank (1932 drugs) and e-LEA3D (1852 drugs) were screened against these proteins. Several criteria were applied to study the protein-drug interactions and the consensus drugs were further studied by molecular dynamics (MD) simulation. Our study found Sulfadoxine (–7.3?kcal/mol) and Pyrimethamine (–7.8?kcal/mol) to show stable interaction with MurB while Lifitegrast (–10.5?kcal/mol) and Sildenafil (–9.1?kcal/mol) showed most reliable interaction with MurE. Furthermore, binding free energy (ΔG_bind), RMSD and RMSF data and the number of hydrogen bonds corroborated the stability of interactions and hence these drugs for repurposing should be explored further.

Communicated by Ramaswamy H. Sarma     

Design of novel quinoline-aminopiperidine derivatives as Mycobacterium tuberculosis (MTB) GyrB inhibitors: an in silico study

Tuberculosis (TB) is an infectious disease that causes a number of deaths, and the development of new, safer and more adequate TB inhibitors/drugs has become a necessity as well as a great challenge. Mycobacterial DNA gyrase B subunit has been identified to be one of the potentially underexploited drug targets in the field of anti-tubercular drug discovery. To design the novel and potent Mycobacterium tuberculosis (MTB) inhibitors, we performed molecular modeling studies that combined the 3D-QSAR, molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations. Forty eight quinoline-aminopiperidine inhibitors which act on DNA gyrase B subunit were used for constructing 3D-QSAR models. The results showed that the best CoMFA model had the high performance with q²?=?0.643, r²?=?0.947, while the best CoMSIA model yielded q²?=?0.536, r²?=?0.948. The contour map was in good agreement with the docking and MD simulations which strongly demonstrated that the molecular modeling was reliable. Based on this information, several potential compounds were designed and their inhibitory activities were also verified by the accomplished models and ADME/T predictions. We hope that our research could bring new ideas to facilitate the development of novel inhibitors with higher inhibitory activity for TB.

Communicated by Ramaswamy H. Sarma     

Design of inhibitors of thymidylate kinase from Variola virus as new selective drugs against smallpox: part II

Abstract

Acknowledging the importance of studies toward the development of measures against terrorism and bioterrorism, this study aims to contribute to the design of new prototypes of potential drugs against smallpox. Based on a former study, nine synthetic feasible prototypes of selective inhibitors for thymidylate kinase from Variola virus (VarTMPK) were designed and submitted to molecular docking, molecular dynamics simulations and binding energy calculations. The compounds are simplifications of two more complex scaffolds, with a guanine connected to an amide or alcohol through a spacer containing ether and/or amide groups, formerly suggested as promising for the design of selective inhibitors of VarTMPK. Our study showed that, despite the structural simplifications, the compounds presented effective energy values in interactions with VarTMPK and HssTMPK and that the guanine could be replaced by a simpler imidazole ring linked to a –NH₂ group, without compromising the affinity for VarTMPK. It was also observed that a positive charge in the imidazole ring is important for the selectivity toward VarTMPK and that an amide group in the spacer does not contribute to selectivity. Finally, prototype 3 was pointed as the most promising to be synthesized and experimentally evaluated.

Communicated by Ramaswamy H. Sarma     

Identification of novel scaffolds to inhibit human mitotic kinesin Eg5 targeting the second allosteric binding site using in silico methods

Human mitotic kinesins are potential anticancer drug targets because of their essential role in mitotic cell division. The kinesin Eg5 (Kinesin-5, kif11) has gained much attention in this regard and has many inhibitors in different phases of clinical trials. All drug candidates considered for Eg5 so far binds to the binding site (Site 1) formed by the loop L5, helices α2 and α3 and are uncompetitive to ATP/ADP. Recently, it has been reported that Eg5 also has a second binding site (Site 2) formed by helices α4 and α6. In the current work, we have screened the compounds in the diversity set-III from National Cancer Institute (NCI) and Zinc database to identify potential inhibitors for Eg5 that specifically binds to the site 2. The compounds were ranked based on the glide extra precision docking scores and the top ranked compounds were found to have pyridazine scaffold. The top five compounds were further evaluated for other drug like properties. Stability of protein-ligand complexes were analyzed using molecular dynamic simulations. Our studies suggest that pyridazine analogs have good MDCK, permeability properties and high binding affinity to the human Eg5.     

Spermine as a possible endogenous allosteric activator of carboxypeptidase A: multispectroscopic and molecular simulation studies

Abstract

Carboxypeptidase A (EC.3.4.17.1) is a zinc-containing proteolytic enzyme that removes the C-terminal amino acid from a peptide chain with the free carboxylate-terminal. In this study, the effect of spermine interaction on the structure and thermal stability of Carboxypeptidase A was investigated by ultraviolet???visible spectroscopy, fluorescence spectroscopy, circular dichroism, Kinetic measurement, molecular docking and simulation studies have also been followed at the pH of 7.5. The transition temperature of Carboxypeptidase A, as a criterion of protein thermal stability, in the presence of spermine was enhanced by increasing the concentration of spermine. The results of fluorescence intensity changes, at two temperatures of 308 and 318?K, also suggested that spermine had a great ability to quench the fluorescence of Carboxypeptidase A through the static quenching procedure. The thermodynamic parameters changes, including standard Gibbs free-energy, entropy and enthalpy, showed that the binding of spermine to Carboxypeptidase A was spontaneous and the hydrogen bonding and van der Waals interactions played a major role in stabilizing the Carboxypeptidase A–spermine complex. The changes in the content of the α-helix and the β-sheet of the Carboxypeptidase A with binding to spermine were shown by the CD spectra method. Further, kinetic studies revealed that by increasing concentration of spermine, the activity of Carboxypeptidase A was enhanced. Also, the docking study revealed that the hydrogen bonding and van der Waals interactions played a major role in stabilizing the Carboxypeptidase A–spermine complex. As a result, spermine could be considered as an activator and a stabilizer for Carboxypeptidase A.

Communicated by Ramaswamy H. Sarma     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

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.     

An in silico pharmacological approach toward the discovery of potent inhibitors to combat drug resistance HIV-1 protease variants

Protease inhibitors (PIs) are crucial drugs in highly active antiretroviral therapy for human immunodeficiency virus-1 (HIV-1) infections. However, resistance owing to mutations challenge the long-term efficacy in the medication of HIV-1-infected individuals. Lopinavir (LPV) and darunavir (DRV), two second-generation drugs are the most potent among PIs, hustling the drug resistance when mutations occur in the active and nonactive site of the protease (PR). Herein, we strive for compounds that can stifle the function of wild-type (WT) HIV-1 PR along with four major single mutants (I54M, V82T, I84V, and L90M) instigating resistance to the PIs using in silico approach. Six common compounds are retrieved from six databases using combined pharmacophore-based and structure-based virtual screening methodology. LPV and DRV are docked and the binding free energy is calculated to set the cut-off value for selecting compounds. Further, to gain insight into the stability of the complexes the molecular dynamics simulation (MDS) is carried out, which uncovers two lead molecules namely NCI-524545 and ZINC12866729. Both the lead molecules connect with WT and mutant HIV-1 PRs through strong and stable hydrogen bond interactions when compared with LPV and DRV throughout the trajectory analysis. Interestingly, NCI-524545 and ZINC12866729 exhibit direct interactions with I50/50′ by replacing the conserved water molecule as evidenced by MDS, which indicates the credible potency of these compounds. Hence, we concluded that NCI-524545 and ZINC12866729 have great puissant to restrain the role of drug resistance HIV-1 PR variants, which can also show better activity through in vivo and in vitro conditions.     

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     

Identification and validation of l-asparaginase as a potential metabolic target against Mycobacterium tuberculosis

Multidrug-resistant Mycobacterium tuberculosis (Mtb) has emerged as a major health challenge, necessitating the search for new molecular targets. A secretory amidohydrolase, l -asparaginase of Mtb (MtA), originally implicated in nitrogen assimilation and neutralization of acidic microenvironment inside human alveolar macrophages, has been proposed as a crucial metabolic enzyme. To investigate whether this enzyme could serve as a potential drug target, it was studied for structural details and active site–specific inhibitors were tested on cultured Mycobacterial strain. The structural details of MtA obtained through comparative modeling and molecular dynamics simulations provided insights about the orchestration of an alternate reaction mechanism at the active site. This was contrary to the critical Tyr flipping mechanism reported in other asparaginases. We report the novel finding of Tyr to Val replacement in catalytic triad I along with the structural reorganization of a β-hairpin loop upon substrate binding in MtA active site. Further, 5 MtA-specific, active-site–based inhibitors were obtained by following a rigorous differential screening protocol. When tested on Mycobacterium culture, 3 of these, M3 (ZINC 4740895), M26 (ZINC 33535), and doxorubicin showed promising results with inhibitory concentrations (IC ₅₀) of 431, 100, and 56 µM, respectively. Based on our findings and considering stark differences with human asparaginase, we project MtA as a promising molecular target against which the selected inhibitors may be used to counteract Mtb infection effectively.     

Conformational selection in silico: loop latching motions and ligand binding in enzymes

Ligand binding frequently induces significant conformational changes in a protein receptor. Understanding and predicting such conformational changes represent an important challenge for computational biology, including applications to structure-based drug design. We describe an approach to this problem based on the assumption that the holo state is at least transiently populated in the absence of a ligand; this hypothesis has been referred to as "conformational selection." Here, we apply a method that tests this hypothesis on a challenging class of ligand-induced conformational changes, which we refer to as loop latching: the closing of a loop around an active site that sequesters the ligand from solvent. The method uses a combination of replica exchange molecular dynamics and a loop prediction algorithm to generate low-energy loop structures, and docking to select the conformation appropriate for binding a particular ligand. On a test set of six proteins, it yields loop structures including hololike conformations, generally below 2 A RMSD from the liganded structure, for loops that span up to 15 residues. Docking serves as a stringent test of the predictions. In five of the six cases, the predicted loop conformations improve the ranks of cognate ligands relative to using the apo structure, although the results remain, in most cases, significantly worse than using a holo structure. The poses of the cognate ligands are correct in four of the six test cases, while they are correct for five of the six using a holo structure.     

Biochemical changes induced by grape seed extract and low level laser therapy administration during intraoral wound healing in rat liver: an experimental and in silico study

In the present study, the changes that occur in rat liver tissue as a result of the use of grape seed extract (GSE) and low level laser therapy (LLLT) in intraoral wound (IW) healing are analyzed using biochemical parameters. Diode laser application groups received 8 J/cm² dose LLLT once a day for 4 days (810 nm wavelength, continuous mode, 0.25 W, 9 s). As a result of the biological parameter analysis, it was determined that the oxidative damage caused by the IWs and recovery period on 7th and 14th days could be substantially removed with GSE applications that have antioxidant capacity especially in rat liver tissue. In addition, the active compound of grape seed, catechin is studied in the active site of glycogen synthase kinase 3 (GSK3) target using molecular modeling approaches. Post-processing molecular dynamics (MD) results for catechin is compared with a standard GSK3 inhibitor. MD simulations assisted for better understanding of inhibition mechanism and the crucial amino acids contributing in the ligand binding. These results along with a through free energy analysis of ligands using sophisticated simulations methods are quite striking and it suggests a greater future role for simulation in deciphering complex patterns of molecular mechanism in combination with methods for understanding drug-receptor interactions.     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Kinetic and in silico studies of hydroxy-based inhibitors of carbonic anhydrase isoforms I and II

A series of hydroxy and phenolic compounds have been assayed for the inhibition of two physiologically relevant carbonic anhydrase (CA, EC 4.2.1.1) isozymes, the cytosolic human isozymes I and II. The investigated molecules showed inhibition constants in the range of 1.07–4003 and 0.09–31.5?μM at the hCA I and hCA II enzymes, respectively. In order to investigate the binding mechanisms of these inhibitors, in silico studies were also applied. Molecular docking scores of the studied compounds are compared using three different scoring algorithms, namely Glide/SP, Glide/XP and Glide/IFD. In addition, different ADME (absorption, distribution, metabolism and excretion) analysis was performed. All the examined compounds were found within the acceptable range of pharmacokinetic profiles.     

In silico studies on marine actinomycetes as potential inhibitors for Glioblastoma multiforme

Glioblastoma multiforme (GBM) is considered to be the most common and often deadly disorder which affects the brain. It is caused by the over expression of proteins such as ephrin type-A receptor 2 (EphA2), epidermal growth factor receptor (EGFR) and EGFRvIII. These 3 proteins are considered to be the potential therapeutic targets for GBM. Among these, EphA2 is reported to be over-expressed in ˜90% of GBM. Herein we selected 35 compounds from marine actinomycetes, 5 in vitro and in vivo studied drug candidates and 4 commercially available drugs for GBM which were identified from literature and analysed by using comparative docking studies. Based on the glide scores and other in silico parameters available in Schrödinger, two selected marine actinomycetes compounds which include Tetracenomycin D and Chartreusin exhibited better binding energy among all the compounds studied in comparative docking. In this study we have demonstrated the inhibition of the 3 selected targets by the two bioactive compounds from marine actinomycetes through in-silico docking studies. Furthermore molecular dynamics simulation were also been performed to check the stability and the amino acids interacted with the 3 molecular targets (EphA2 receptor, EGFR, EGFRvIII) for GBM. Our results suggest that Tetracinomycin D and Chartreusin are the novel and potential inhibitor for the treatment of GBM.