QSAR modeling,pharmacophore-based virtual screening,and ensemble docking insights into predicting potential epigallocatechin gallate (EGCG) analogs against epidermal growth factor receptor

Epigallocatechin gallate (EGCG) is a major polyphenols of green tea may have the possibility to inhibit epidermal growth factor receptor (EGFR) activity and lead to reduce non-small cell lung cancer (NSCLC) progression. However, EGCG has some toxic features; moreover, there is a lack of explorations into the molecular interaction mechanisms of EGCG and the EGFR. In this examination, integration of quantitative structure–activity relationship (QSAR) modeling, pharmacophore-based virtual screening, and ensemble docking approaches were used to predict potential novel EGCG analogs as effective EGFR inhibitors. QSAR modeling of logP and logS predictions and toxicity endpoint investigation for a set of 82 compounds were shown good predictive ability and robustness from the applicability domain and confusion matrix elucidations. Virtual screening and docking studies revealed that seven high potential EGCG analogs as strong EGFR binders. Molecular interactions interpretations indicated some insights into the structural features of ligands that efficiently interfere with mutation possible residues (Gly⁷¹⁹ and Thr⁷⁹⁰) of the EGFR. The hydrogen bonds, hydrophobic interactions, atomic π-cation interactions and salt bridges of ligands are contributing additional stability to receptor structure, which can lead to blocking the intracellular protein-tyrosine kinase activity, including EGFR associated pathways activation in NSCLC. Therefore, this can characterize as a block-cluster mechanism between EGCG analogs and EGFR complexes. In silico anti-EGFR and anticancer activity predictions suggested that, ligands could act as promising pharmacological, anticancer, and drug-like templates of EGFR towards moderating the NSCLC progressions. These results and provided pinpoints could be beneficial to recognize probable therapeutic targets for NSCLC therapy.     

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.     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Screening from TCM Database@Taiwan and QSAR model for identifying HER2 inhibitors

Human epidermal growth factor receptor type 2 (HER2) overexpression, which has been reported to increase the malignancy of human ovarian cancer cells and the metastatic potential of human breast cancer cells, is an important factor in tumour formation and tumour growth. However, currently available HER2 inhibitors, such as Lapatinib, for cancer therapy cause adverse side effects including diarrhoea, rash and possible liver toxicity. We hoped to find novel agents that cause less adverse side effects by performing virtual screening process on the world's largest traditional Chinese medicine compound database. The results thus obtained were then validated using 3D quantitative structure–activity relationship model. Top three candidates were selected from the docking results. The top three candidates and the control both formed a hydrogen bond with the key residue, Lys724. This showed that the candidates and the control have similar binding effects to HER2. These candidates were investigated using comparative molecular field analysis and comparative molecular similarity indices analysis models. The results from these models showed high correlation coefficients (r ²) of 0.9547 and 0.9226, respectively. All top three candidates had high docking scores, favourable pharmacophores and functional groups forming stable hydrogen bonds with HER2. These properties suggested stable binding affinities and favourable interaction with HER2. We concluded that these candidates may be further investigated as potential HER2 inhibitors.     

In silico docking and molecular dynamics simulation of 3-dehydroquinate synthase (DHQS) from <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis>

The shikimate pathway is as an attractive target because it is present in bacteria, algae, fungi, and plants but does not occur in mammals. In Mycobacterium tuberculosis (MTB), the shikimate pathway is integral to the biosynthesis of naphthoquinones, menaquinones, and mycobactin. In these study, novel inhibitors of 3-dehydroquinate synthase (DHQS), an enzyme that catalyzes the second step of the shikimate pathway in MTB, were determined. 12,165 compounds were selected from two public databases through virtual screening and molecular docking analysis using PyRx 8.0 and Autodock 4.2, respectively. A total of 18 compounds with the best binding energies (?13.23 to ?8.22 kcal/mol) were then selected and screened for absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis, and nine of those compounds were found to satisfy all of the ADME and toxicity criteria. Among those nine, the three compounds—ZINC633887?(binding energy =??10.29 kcal/mol), ZINC08983432?(?9.34 kcal/mol), and PubChem73393?(?8.61 kcal/mol)—with the best binding energies were further selected for molecular dynamics (MD) simulation analysis. The results of the 50-ns MD simulations showed that the two compounds ZINC633887 and PubChem73393 formed stable complexes with DHQS and that the structures of those two ligands remained largely unchanged at the ligand-binding site during the simulations. These two compounds identified through docking and MD simulation are potential candidates for the treatment of TB, and should undergo validation in vivo and in vitro.     

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with ?259.41, ?110.09, ?188.25, ?163.30, ?202.10, and ?158.79 kJ mol^?1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.     

Discovery of VEGFR inhibitors through virtual screening and energy assessment

Vascular endothelial growth factor receptor-2 (VEGFR-2) is crucial in promoting tumor angiogenesis and cancer metastasis. Thus, inhibition of VEGFR-2 has appeared as a good tactic for cancer treatment. To find out novel VEGFR-2 inhibitors, first, the PDB structure of VEGFR-2, 6GQO, was selected based on atomic nonlocal environment assessment (ANOLEA) and PROCHECK assessment. 6GQO was then further used for structure-based virtual screening (SBVS) of different molecular databases, including US-FDA approved drugs, US-FDA withdrawn drugs, may bridge, MDPI, and Specs databases using Glide. Based on SBVS, receptor fit, drug-like filters, and absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis of 427877 compounds, the best 22 hits were selected. From the 22 hits, hit 5 complex with 6GQO was put through molecular mechanics/generalized born surface area (MM/GBSA) study and hERG binding. The MM/GBSA study revealed that hit 5 possesses lesser binding free energy with more inferior stability in the receptor pocket than the reference compound. The VEGFR-2 inhibition assay of hit 5 disclosed an IC₅₀ of 165.23 nM against VEGFR-2, which can be possibly enhanced through structural modifications.     

Exploring dual inhibitory role of febrifugine analogues against Plasmodium utilizing structure-based virtual screening and molecular dynamic simulation

Malaria is an endemic disease caused by the protozoan parasite Plasomodium falciparum. Febrifugine analogues are natural compound obtained from the traditional Chinese herbs have shown significant antimalarial and anticancerous efficacy in experimental model. Development of resistance against the existing antimalarial drug has alarmed the scientific innovators to find a potential antimalarial molecule which can be further used by endemic countries for the elimination of this disease. In this study, structure-based virtual screening and molecular dynamics (MD) base approaches were used to generate potential antimalarial compound against plasmepsin II and prolyl-tRNA synthetase of Plasmodium. Here, we have docked series of febrifugine analogues (n = 11,395) against plasmepsin II in three different docking modes and then it was compared with previously reported target prolyl-tRNA synthetase. Extra precision docking resulted into 235 ligands having better docking score were subject for QikProp analysis. Better ligands (n = 39) obtained from QikProp analysis were subject for ADMET prediction and docking protocol validation through the estimation of receiver operator characteristics. In the later stage, 24 ligands obtained from ADMET study were subject for the estimation of binding energy through MM-GBSA and same were also docked against prolyl-tRNA synthetase to get compounds with dual inhibitor role. Finally, MD simulation and 2D fingerprint MACCS study of two best ligands have shown significant interaction with plasmepsin II and homology against known active ligand with noteworthy MACCS index, respectively. This study concludes that FA12 could be potential drug candidate to fight against Plasmodium falciparum parasites.     

Discovery and molecular docking of quinolyl-thienyl chalcones as anti-angiogenic agents targeting VEGFR-2 tyrosine kinase

Vascular endothelial growth factor Receptor-2 (VEGFR-2) kinase inhibition is one of the well established strategies to promptly tackle tumor growth by suppression of angiogenesis. In the current study, structure-based virtual screening methodology of a series of quinolyl-thienyl chalcones indicated their strong potential as VEGFR-2 kinase inhibitors. In vitro VEGFR-2 kinase inhibitory activity was found to be significant (compound 19, IC(50): 73.41nM). All compounds showed significant inhibition of human umbilical vein endothelial cells (HUVEC) proliferation (compound 19, IC(50): 21.78nM). Molecular interactions of the compounds were studied using molecular docking studies.     

Similarity search combined with docking and molecular dynamics for novel hAChE inhibitor scaffolds

The main purpose of this study was to address the performance of virtual screening methods based on ligands and the protein structure of acetylcholinesterase (AChE) in order to retrieve novel human AChE (hAChE) inhibitors. In addition, a protocol was developed to identify novel hit compounds and propose new promising AChE inhibitors from the ZINC database with 10 million commercially available compounds. In this sense, 3D similarity searches using rapid overlay of chemical structures and similarity analysis through comparison of electrostatic overlay of docked hits were used to retrieve AChE inhibitors from collected databases. Molecular dynamics simulation of 100 ns was carried out to study the best docked compounds from similarity searches. Some key residues were identified as crucial for the dual binding mode of inhibitor with the interaction site. All results indicated the relevant use of EON and docking strategy for identifying novel hit compounds as promising potential anticholinesterase candidates, and seven new structures were selected as potential hAChE inhibitors.

Open image in new window

Graphical abstract Compound N01 in the 4M0E hAChE crystallography structure from docking results. Yellow dashed lines Hydrogen bonds, blue dashed lines π-stacking interactions, green dashed lines cation-π interactions.

     

IQGAP1-Dependent Signaling Pathway Regulates Endothelial Cell Proliferation and Angiogenesis

Background

Vascular endothelial growth factor receptor-2 (VEGFR-2) signaling is an obligate requirement for normal development and pathological angiogenesis such as cancer and age-related macular degeneration. Although autophosphorylation of tyrosine 1173 (Y1173) of VEGFR-2 is considered a focal point for its angiogenic signal relay, however, the mechanism of phosphorylation of Y1173, signaling proteins that are recruited to this residue and their role in angiogenesis is not fully understood.

Methodology/Principal Findings

In this study we demonstrate that c-Src kinase directly through its Src homology 2 (SH2) domain and indirectly via c-Cbl binds to phospho-Y1057 of VEGFR-2. Activation of c-Src kinase by a positive feedback mechanism phosphorylates VEGFR-2 at multi-docking site, Y1173. c-Src also catalyzes tyrosine phosphorylation of IQGAP1 and acts as an adaptor to bridge IQGAP1 to VEGFR-2. In turn, IQGAP1 activates b-Raf and mediates proliferation of endothelial cells. Silencing expression of IQGAP1 and b-Raf revealed that their activity is essential for VEGF to stimulate angiogenesis in an in vivo angiogenesis model of chicken chorioallantoic membrane (CAM).

Conclusions/Significance

Angiogenesis contributes to the pathology of numerous human diseases ranging from cancer to age-related macular degeneration. Determining molecular mechanism of tyrosine phosphorylation of VEGFR-2 and identification of molecules that are relaying its angiogenic signaling may identify novel targets for therapeutic intervention against angiogenesis-associated diseases. Our study shows that recruitment and activation of c-Src by VEGFR-2 plays a pivotal role in relaying angiogenic signaling of VEGFR-2; it phosphorylates VEGFR-2 at Y1173, facilitates association and activation of IQGAP1 and other signaling proteins to VEGFR-2. IQGAP1-dependent signaling, in part, is critically required for endothelial cell proliferation, a key step in angiogenesis. Thus, Y1057 of VEGFR-2 serves to regulate VEGFR-2 function in a combinatorial manner by supporting both diversity of recruitment of angiogenic signaling proteins to VEGFR-2, and its ability to promote angiogenesis.     

Design,synthesis, docking,ADMET studies,and anticancer evaluation of new 3-methylquinoxaline derivatives as VEGFR-2 inhibitors and apoptosis inducers

Vascular endothelial growth factor receptor-2 (VEGFR-2) plays a critical role in cancer angiogenesis. Inhibition of VEGFR-2 activity proved effective suppression of tumour propagation. Accordingly, two series of new 3-methylquinoxaline derivatives have been designed and synthesised as VEGFR-2 inhibitors. The synthesised derivatives were evaluated in vitro for their cytotoxic activities against MCF-7and HepG2 cell lines. In addition, the VEGFR-2 inhibitory activities of the target compounds were estimated to indicate the potential mechanism of their cytotoxicity. To a great extent, the results of VEGFR-2 inhibition were highly correlated with that of cytotoxicity. Compound 27a was the most potent VEGFR-2 inhibitor with IC₅₀ of 3.2 nM very close to positive control sorafenib (IC₅₀ = 3.12 nM). Such compound exhibited a strong cytotoxic effect against MCF-7 and HepG2, respectively with IC₅₀ of 7.7 and 4.5 µM in comparison to sorafenib (IC₅₀ = 3.51 and 2.17 µM). In addition, compounds 28, 30f, 30i, and 31b exhibited excellent VEGFR-2 inhibition activities (IC₅₀ range from 4.2 to 6.1 nM) with promising cytotoxic activity. Cell cycle progression and apoptosis induction were investigated for the most active member 27a. Also, the effect of 27a on the level of caspase-3, caspase-9, and BAX/Bcl-2 ratio was determined. Molecular docking studies were implemented to interpret the binding mode of the target compounds with the VEGFR-2 pocket. Furthermore, toxicity and ADMET calculations were performed for the synthesised compounds to study their pharmacokinetic profiles     

Searching for anthranilic acid-based thumb pocket 2 HCV NS5B polymerase inhibitors through a combination of molecular docking, 3D-QSAR and virtual screening

A combination of the following computational methods: (i) molecular docking, (ii) 3-D Quantitative Structure Activity Relationship Comparative Molecular Field Analysis (3D-QSAR CoMFA), (iii) similarity search and (iv) virtual screening using PubChem database was applied to identify new anthranilic acid-based inhibitors of hepatitis C virus (HCV) replication. A number of known inhibitors were initially docked into the “Thumb Pocket 2” allosteric site of the crystal structure of the enzyme HCV RNA-dependent RNA polymerase (NS5B GT1b). Then, the CoMFA fields were generated through a receptor-based alignment of docking poses to build a validated and stable 3D-QSAR CoMFA model. The proposed model can be first utilized to get insight into the molecular features that promote bioactivity, and then within a virtual screening procedure, it can be used to estimate the activity of novel potential bioactive compounds prior to their synthesis and biological tests.     

The comparison of VEGFR-1-binding domain of VEGF-A with modelled VEGF-C sheds light on receptor specificity

Receptor specificity determines the role of vascular endothelial growth factors (VEGFs), which either induce angiogenesis via VEGFR-1 and VEGFR-2 receptors or lymphangiogenesis via the VEGFR-3 receptor. Among the VEGFs, VEGF-A and VEGF-B predominantly induce angiogenesis while VEGF-C and VEGF-D induce lymphangiogenesis. The answer for the question of why VEGF-C and VEGF-D are not able to bind VEGFR-1 and behave as angiogenic growth factors may hide behind the details of the tertiary structures of these proteins. In the present study, the tertiary structure of human VEGF-C protein was modelled and the model was compared with the known human VEGF-A tertiary structure. In overall, the modelled structure highly resembled the structure of VEGF-A. The respective key residues that are involved in cysteine-knot motif formation in VEGF-A are similarly located and identically oriented in VEGF-C, indicating the presence of a VEGF-A-like homodimer. However, a VEGF-C homodimer created via monomer docking did not superimpose well with the VEGF-A homodimer. Rigid docking models of VEGF-C with the VEGFR-1 receptor revealed that in the VEGF-C–VEGFR-1 complex, the receptor–protein-interacting residues were not correctly oriented to induce angiogenesis via VEGFR-1. Mapping the electrostatic surface potentials to the protein surfaces revealed noteworthy number of dissimilarities between VEGF-A and VEGF-C, indicating that overall both proteins differ in their folding properties and stability.     

Quinoxalinone (Part II). Discovery of (Z)-3-(2-(pyridin-4-yl)vinyl)quinoxalinone derivates as potent VEGFR-2 kinase inhibitors

Inhibition of VEGFR-2 kinase has been highlighted as one of the well-defined strategies to suppress tumor growth via blockade of angiogenesis. Guided by the principles of bioisosteric replacement and pharmacophoric fragment migration, a series of novel quinoxalinone derivates were designed, synthesized and evaluated for their VEGFR-2 inhibitory potencies. Among them, compounds 7c, 8b, 8c, 8e and 10b displayed antiangiogenic abilities via the in vitro tube formation assay (cellular level) and ex vivo rat aortic ring assay (tissue level) at a low concentration (0.1 μM). By means of in vivo zebrafish embryo model, two (Z)-3-(2-(pyridin-4-yl)vinyl)quinoxalinone derivates 8c and 8e showed significant antiangiogenesis effects, suggesting they have potentials to be developed into antiangiogenesis agents via further structural optimization. Moreover, these two compounds also demonstrated potent inhibition toward VEGFR-2 and B-raf kinases in a low concentration (1 μM). A possible interpretation of our evaluation result has been presented by a molecular docking study by docking representative compound 8c with VEGFR-2.     

Dietary Compound Isoliquiritigenin Inhibits Breast Cancer Neoangiogenesis via VEGF/VEGFR-2 Signaling Pathway

Angiogenesis is crucial for cancer initiation, development and metastasis. Identifying natural botanicals targeting angiogenesis has been paid much attention for drug discovery in recent years, with the advantage of increased safety. Isoliquiritigenin (ISL) is a dietary chalcone-type flavonoid with various anti-cancer activities. However, little is known about the anti-angiogenic activity of isoliquiritigenin and its underlying mechanisms. Herein, we found that ISL significantly inhibited the VEGF-induced proliferation of human umbilical vein endothelial cells (HUVECs) at non-toxic concentration. A series of angiogenesis processes including tube formation, invasion and migration abilities of HUVECs were also interrupted by ISL in vitro. Furthermore, ISL suppressed sprout formation from VEGF-treated aortic rings in an ex-vivo model. Molecular mechanisms study demonstrated that ISL could significantly inhibit VEGF expression in breast cancer cells via promoting HIF-1α (Hypoxia inducible factor-1α) proteasome degradation and directly interacted with VEGFR-2 to block its kinase activity. In vivo studies further showed that ISL administration could inhibit breast cancer growth and neoangiogenesis accompanying with suppressed VEGF/VEGFR-2 signaling, elevated apoptosis ratio and little toxicity effects. Molecular docking simulation indicated that ISL could stably form hydrogen bonds and aromatic interactions within the ATP-binding region of VEGFR-2. Taken together, our study shed light on the potential application of ISL as a novel natural inhibitor for cancer angiogenesis via the VEGF/VEGFR-2 pathway. Future studies of ISL for chemoprevention or chemosensitization against breast cancer are thus warranted.     

Identification of novel serotonin 2C receptor ligands by sequential virtual screening

Pharmacophore modelling, docking and virtual screening have become important tool in drug discovery process. Serotonin 2C (5-HT_2C) receptor ligands have got major attention for their therapeutic uses as antidepressant and anorectic agents. Two step pharmacophore and docking based virtual screening was done using 5-HT_2C agonists. Two common feature pharmacophore directed virtual hits had submicromolar activity. Refined pharmacophore with excluded volumes was constructed and combined with homology model based docking. Best hit from this virtual screening showed IC₅₀ of 20.1 nM. Similarity search of this hit compound resulted more active ligand with 7.8 nM activity.     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.     

In silico screening of anti-inflammatory compounds from Lichen by targeting cyclooxygenase-2

Abstract

Non-steroidal anti-inflammatory drugs (NSAID) targeting cyclooxygenase-2 are clinically effective. However, they lack anti-thrombotic activity resulting in incidences of adverse effects like myocardial infarction, gastrointestinal and abdominal discomfort which necessitate for discovering new drug candidates with improved therapeutic effects and tolerability. Various recent researches have suggested that many lichens offer a vast reservoir for anti-inflammatory drug candidates which are natural as well as safe for human consumption. Drug discovery is a very complex and time-consuming process; however, in silico techniques can make this process simple and economic. Hence to find out natural anti-inflammatory compounds, we have carried out the virtual screening of 412 lichen compounds by molecular docking with human Cox-2 enzyme and validated the docking score by X-Score followed by ADMET and Drug-likeness analysis. The resulting 6 top-scored compounds were subjected to Molecular dynamics simulation (MDS) to analyze the stability of docked protein-ligand complex, to assess the fluctuation and conformational changes during protein-ligand interaction. The values of RMSD, Rg, and interaction energy after 30?ns of MDS revealed the good stability of these Lichen compounds in the active site pocket of Cox-2 in compare to reference, JMS. Additionally, we have done the pharmacophore analysis which found many common pharmacophore features between Lichen compounds and well known anti-inflammatory compounds. Our result shows that these lichen compounds are potential anti-inflammatory candidates and could be further modified and evaluated to develop more effective anti-inflammatory drugs with fewer side effects for the treatment of inflammatory diseases.

Communicated by Ramaswamy H. Sarma     

Uroporphyrinogen Decarboxylase as a Potential Target for Specific Components of Traditional Chinese Medicine: A Virtual Screening and Molecular Dynamics Study

Uroporphyrinogen decarboxylase (UROD) has been suggested as a protectant against radiation for head and neck cancer (HNC). In this study, we employed traditional Chinese medicine (TCM) compounds from TCM nawiaT@esabataD (http://tcm.cmu.edu.tw/) to screen for drug-like candidates with potential UROD inhibition characteristics using virtual screening techniques. Isopraeroside IV, scopolin, and nodakenin exhibited the highest Dock Scores, and were predicted to have good Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties. Two common moieties, 2H-chromen-2-one and glucoside, were observed among the top TCM candidates. Cross comparison of the docking poses indicated that candidates formed stable interactions with key binding and catalytic residues of UROD through these two moieties. The 2H-chromen-2-one moiety enabled pi-cation interactions with Arg37 and H-bonds with Tyr164. The glucoside moiety was involved in forming H-bonds with Arg37 and Asp86. From our computational results, we propose isopraeroside IV, scopolin, and nodakenin as ligands that might exhibit drug-like inhibitory effects on UROD. The glucoside and 2H-chromen-2-one moieties may potentially be used for designing inhibitors of UROD.