Pharmacoinformatics-based identification of chemically active molecules against Ebola virus

Abstract

Ebola is a dangerous virus transmitted by animals and humans and to date there is no curable agent for such a deadly infectious disease. In this study, pharmacoinformatics-based methods were adopted to find effective novel chemical entities against Ebola virus. A well predictive and statistical robust pharmacophore model was developed from known Ebola virus inhibitors collected from the literature. The model explained the significance of each of hydrogen bond acceptor and donor, and two hydrophobic regions for activity. The National Cancer Institute and Asinex (Antiviral library) databases were screened using the final validated pharmacophore model. Initial hits were further screened with a set of criteria and finally eight molecules from both databases were proposed as promising anti Ebola agents. Further molecular docking and molecular dynamics studies were carried out and it was found that the proposed molecules possessed capability to interact with amino residues of Ebola protein as well as retaining equilibrium of protein-ligand systems. Finally, the binding energies were calculated using molecular mechanics Poisson–Boltzmann surface area approach and all proposed molecules showed strong binding affinity towards the Ebola protein receptor.

Communicated by Ramaswamy H. Sarma     

Identification potential inhibitors against the Streptococcus quorum-sensing signal pathway

Abstract

Streptococcal infections are common in human and antibiotics are frequently prescribed in clinical practice. However, infections caused by drug-resistant strains are particularly difficult to treat using common antibiotics. Hence, there is an urgent need for new antibiotics. Quorum sensing is a regulatory mechanism involving cell communication that is thought to play an important role in various bacterial infections, including those caused by Streptococcus. The ATP-binding cassette transporter ComA of Streptococcus is essential for quorum-sensing signal production. The inhibition of the ComA peptidase domain (ComA PEP) suppresses the quorum-sensing pathway and resulting changes in phenotype and/or behavior. Using virtual screening and molecular dynamics simulations, two promising candidate compounds, ZINC32918029 and ZINC6751571, were found. These compounds had similar binding modes and interactions to the experimentally determined reference inhibitor 6CH. However, a significantly stronger negative binding energy was achieved (?113.501?±?15.312?KJ/mol and ?103.153?±?11.912?KJ/mol for ZINC32918029 and ZINC6751571, respectively). Molecular dynamics simulations also revealed that ZINC32918029 and ZINC6751571 had a strong affinity for ComA PEP. These results indicate that ZINC32918029 and ZINC6751571 are promising candidate inhibitors of the Streptococcus quorum-sensing pathway.

Communicated by Ramaswamy H. Sarma     

Homology modeling,docking and structure-based virtual screening for new inhibitor identification of Klebsiella pneumoniae heptosyltransferase-III

Abstract

Klebsiella pneumoniae (K. pneumoniae) is a Gram-negative opportunistic pathogen commonly associated with hospital-acquired infections that are often resistant even to antibiotics. Heptosyltransferase (HEP) belongs to the family of glycosyltransferase-B (GT-B) and plays an important in the synthesis of lipopolysaccharides (LPS) essential for the formation of bacterial cell membrane. HEP-III participates in the transfer of heptose sugar to the outer surface of bacteria to synthesize LPS. LPS truncation increases the bacterial sensitivity to hydrophobic antibiotics and detergents, making the HEP as a novel drug target. In the present study, we report the 3D homology model of K. pneumoniae HEP-III and its structure validation. Active site was identified based on similarities with known structures using Dali server, and structure-based pharmacophore model was developed for the active site substrate ADP. The generated pharmacophore model was used as a 3D search query for virtual screening of the ASINEX database. The hit compounds were further filtered based on fit value, molecular docking, docking scores, molecular dynamics (MD) simulations of HEP-III complexed with hit molecules, followed by binding free energy calculations using Molecular Mechanics-Poisson–Boltzmann Surface Area (MM-PBSA). The insights obtained in this work provide the rationale for design of novel inhibitors targeting K. pneumoniae HEP-III and the mechanistic aspects of their binding.

Communicated by Ramaswamy H. Sarma     

Docking and molecular dynamics studies of new potential inhibitors of the human epidermal receptor 2

Compounds similar to lapatinib and gefitinib have been investigated as potential inhibitors of the intracellular receptor tyrosine kinase (RTK) domain of the human epidermal receptor 2 (HER2), which is a promising molecular target to the drug design of new chemotherapies for breast, lung, ovarian and colorectal cancers. In this study, we have searched potential HER2 inhibitors used for treatment of other illnesses such as hepatitis, bacterial infections and sexual impotence screened in the DrugBank. The compounds selected were subjected to virtual screening docking in order to evaluate the main interactions between them and the RTK domain of HER2. The selected compounds were investigated by flexible docking, molecular dynamics studies and ΔG _bind calculations. The results suggest that antrafenine, saprisartan, reserpine, irinotecan and udenafil are potential candidates to inhibit the RTK domain of HER2.     

In silico insights on tankyrase protein: A potential target for colorectal cancer

Abstract

The Wnt/β-catenin pathway plays an important regulatory role in cancer signaling and cell regenerative mechanisms. Its suppression has long been considered as an important challenge of anticancer treatment and management. The poly(ADP-ribose) polymerase (PARP) family represented as a new class of therapeutic targets with diverse potential disease indications. Tankyrase (TNKS) is considered to be a potential target for the intervention of various cancers. The main objective of the work is to explore the molecular and quantum mechanics of the drug-like compounds and to identify the potential inhibitors for TNKS protein using the structure and ligand-based virtual screening from several databases and to explore the binding pocket and interactions of active residues. The screened compounds were further filtered using binding-free energy calculation and molecular dynamics simulation studies. The results have provided a strong molecular knowledge of TNKS and offered top hit potent inhibitors. The identified lead compounds LC_40781, LC_40777, LC_39767, LC_8346, NCI_682438, and NCI_721141 were observed to have potent activity against TNKS protein. The hydrogen bonding of compounds with Asp1198, His1201, Tyr1203 in TNKS1 and Gly1032, Ser1068 in TNKS2 are the key interactions plays a major role in binding energy. Therefore, the outcome of the study would help for further validation and provides valuable information to guide the future TNKS-specific inhibitor designing.

Communicated by Ramaswamy H. Sarma     

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.     

Molecular modeling,dynamics, and an insight into the structural inhibition of cofactor independent phosphoglycerate mutase isoform 1 from Wuchereria bancrofti using cheminformatics and mutational studies

Phosphoglycerate mutase catalyzes the interconversion between 2-phosphoglycerate and 3-phosphoglycerate in the glycolytic and gluconeogenic pathways. They exist in two unrelated forms, that is either cofactor (2,3-diphosphoglycerate) dependent or cofactor-independent. These two enzymes have no similarity in amino acid sequence, tertiary structure, and in catalytic mechanism. Wuchereria bancrofti (WB) contains the cofactor-independent form, whereas other organisms can possess the dependent form or both. Since, independent phosphoglycerate mutase (iPGM) is an essential gene for the survival of nematodes, and it has no sequence or structural similarity to the cofactor-dependent phosphoglycerate mutase found in mammals, it represents an attractive drug target for the filarial nematodes. In this current study, a putative cofactor-iPGM gene was identified in the protein sequence of the WB. In the absence of crystal structure, a three-dimensional structure was determined using the homology modeling approximation, and the most stable protein conformation was identified through the molecular dynamics simulation studies, using GROMACS 4.5. Further, the functional or characteristic residues were identified through the sequence analysis, potential inhibitors were short-listed and validated, and potential inhibitors were ranked using the cheminformatics and molecular dynamics simulations studies, Prime MM-GBSA approach, respectively.     

Virtual screening,molecular dynamics,and binding free energy calculations on human carbonic anhydrase IX catalytic domain for deciphering potential leads

Carbonic anhydrase IX is a tumor-associated membrane-bound metallo-enzyme which catalyzes the reversible hydration of carbon dioxide (CO₂) to bicarbonate (HCO₃^?) and proton (H⁺) ions. It is a hypoxia-inducible enzyme and plays a critical role in tumor pH homeostasis favoring tumor cell invasiveness and drug resistance. Over expression of CAIX is documented in cancers of breast, lung, kidney, colon/rectum, etc. Chemical inhibition of CAIX activity has proven to be an effective therapeutic modality towards targeting cancer. Hence, in this study, we intend to identify potential molecules from NCI (National Cancer Institute) and Maybridge databases implementing high-throughput virtual screening. CAIX co-crystallized with acetazolamide (a known inhibitor of CAIX) (PDB ID: 3IAI) was used for reference-guided docking protocol. The potential inhibitors among the coupled data sets were finalized based on Glide docking score, Prime/MMGBSA scoring, significant intermolecular interactions, ADMET (absorption, distribution, metabolism and excretion, toxicity) prediction and stability of complex formation, molecular dynamics simulation, and comparative analysis. By this study, we propose NSC_93618, NSC_170253, NSC_93618, JFD03677, SEW06488, and BTB09372 to be highly significant, as all these compounds were found to qualify as potential leads surpassing all the stringent filtering process. However, NSC_93618 was found to be the most potential, as it featured with higher complex stability with strong bonded interactions, binding affinity synonymous to acetazolamide. Hence, these proposed compounds shall prove to be effective in targeting CAIX towards modulating carcinogenesis.     

Refined homology model of cytochrome bcc complex B subunit for virtual screening of potential anti-tuberculosis agents

Abstract

Cytochrome bcc complex is important for ATP synthesis and cellular activity, as a crucial step in the terminal reduction of oxygen in aerobic electron transport chains. The b subunit of cytochrome bcc complex (QcrB) has been reported as a promising anti-tuberculosis target, with many novel anti-tuberculosis scaffolds reported. However, the 3D structure of mycobacterium tuberculosis (M. tuberculosis) QcrB has not been released, making it hard to understand the interactions between QcrB and its inhibitors as well as to develop novel anti-tuberculosis scaffolds. Herein we built the optimal homology model of M. tuberculosis QcrB using the M. smegmatis QcrB structure as template, which was refined through all-atom molecular dynamics simulation. Then, the binding modes of known inhibitors were predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the accuracy of docking results and stability of the protein-inhibitor complexes. The informative key residues within QcrB site enabled us to perform structure-based virtual library screening to obtain potential M. tuberculosis QcrB inhibitors, which were validated through molecular dynamics simulation and MM-GBSA calculation and analyzed through pharmacokinetic properties prediction. Our research would provide a deeper insight into the interactions between M. tuberculosis QcrB and its inhibitors, which boosts to develop novel therapy against tuberculosis.

Communicated by Ramaswamy H. Sarma     

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.     

Identification of novel PPARα/γ dual agonists by virtual screening,ADMET prediction and molecular dynamics simulations

PPARα and PPARγ have been the most widely studied Peroxisome proliferator-activated receptor (PPAR) subtypes due to their important roles in regulating glucose, lipids, and cholesterol metabolism. By combining the lowering serum triglyceride levels benefit of PPARα agonists (such as fibrates) with the glycemic advantages of the PPARγ agonists (such as TZD), the dual PPAR agonists approach can both improve the metabolic effects and minimize the side effects caused by either agent alone, and hence, has become a promising strategy for designing effective drugs against type-2 diabetes. In this study, by means of virtual screening, ADMET prediction and molecular dynamics (MD) simulations techniques, one compound-ASN15761007 with high binding score, low toxicity were gained. It was observed by MD simulations that ASN15761007 not only possessed the same function as AZ242 did in activating PPARα and BRL did in activating PPARγ, but also had more favorable conformation for binding to the two receptors. Our results provided an approach to rapidly produce novel PPARα/γ dual agonists which might be a potential lead compound to develop against insulin resistance and hyperlipidemia.     

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.     

3D QSAR pharmacophore-based virtual screening for the identification of potential inhibitors of tyrosinase

Abstract

Tyrosinase plays an important role in melanin biosynthesis and protects skin against ultraviolet radiations. Functional deficiency of tyrosinase results in serious dermatological diseases. Tyrosinase also participates in neuromelanin formation in the human brain, which leads to neurodegeneration resulting in Parkinson’s disease. In fruits and vegetables, tyrosinase plays a critical role in senescence, causing undesired browning that results in faster deterioration and shorter shelf lines. The only commercially available tyrosinase is mushroom tyrosinase and it shows the highest homology to the mammalian tyrosinase. Although kojic acid is currently used as a tyrosinase inhibitor, they have serious side effects such as dermatitis, carcinogenesis and hepatotoxicity. Therefore, in order to develop a more active and safer tyrosinase inhibitor, 3D QSAR pharmacophore models were generated based on experimentally known inhibitors. The pharmacophore model, Hypo1, was developed with a large cost difference, high correlation coefficient and low RMS deviation. Hypo1 showed a good spatial arrangement; consisting of five-point features including two hydrogen bond acceptor, one hydrogen bond donor and two hydrophobic features. Hypo1 was further validated by cost analysis, test set and Fisher’s randomisation method. Hypo1 was used as a 3D query for screening the in-house drug-like databases, and the hits were further selected by applying ADMET, Lipinski’s rule of five and fit value criteria. To identify binding conformations, the obtained hits were subjected to molecular docking. Finally, molecular dynamics simulations revealed the appropriate binding modes of hit compounds. To conclude, we propose the final three hit compounds with new structural scaffolds as a virtual candidate as tyrosinase inhibitors.

Communicated by Ramaswamy H. Sarma