Identifying the novel inhibitors of lysine-specific demethylase 1 (LSD1) combining pharmacophore-based and structure-based virtual screening

Abstract

Lysine-specific demethylase 1 (LSD1) has been reported to connect with a range of solid tumors. Thus, the exploration of LSD1 inhibitors has emerged as an effective strategy for cancer treatment. In this study, we constructed a pharmacophore model based on a series of flavin adenine dinucleotide (FAD)-competing inhibitors bearing triazole???dithiocarbamate scaffold combining docking, structure–activity relationship (SAR) study, and molecular dynamic (MD) simulation. Meanwhile, another pharmacophore model was also constructed manually, relying on several speculated substrate-competing inhibitors and reported putative vital interactions with LSD1. On the basis of the two pharmacophore models, multi-step virtual screenings (VSs) were performed against substrate-binding pocket and FAD-binding pocket, respectively, combining pharmacophore-based and structure-based strategy to exploit novel LSD1 inhibitors. After bioassay evaluation, four compounds among 21 hits with diverse and novel scaffolds exhibited inhibition activity at the range of 3.63–101.43?μM. Furthermore, substructure-based enrichment was performed, and four compounds with a more potent activity were identified. After that, the time-dependent assay proved that the most potent compound with IC₅₀ 2.21?μM inhibits LSD1 activity in a manner of time-independent. In addition, the compound exhibited a cellular inhibitory effect against LSD1 in MGC-803 cells and may inhibit cell migration and invasion by reversing EMT in cultured gastric cancer cells. Considering the binding mode and SAR of the series of compounds, we could roughly deem that these compounds containing 3-methylxanthine scaffold act through occupying substrate-binding pocket competitively. This study presented a new starting point to develop novel LSD1 inhibitors.     

Discovery of potential novel microsomal triglyceride transfer protein inhibitors via virtual screening of pharmacophore modelling and molecular docking

In order to identify potential natural inhibitors against the microsomal triglyceride transfer protein (MTP), HipHop models were generated using 20 known inhibitors from the Binding Database. Using evaluation indicators, the best hypothesis model, Hypo1, was selected and utilised to screen the Traditional Chinese Medicine Database, which resulted in a hit list of 58 drug-like compounds. A homology model of MTP was built by MODELLER and was minimised by CHARMm force field. It was then validated by Ramachandran plot and Verify-3D so as to obtain a stable structure, which was further used to refine the 58 hits using molecular docking studies. Then, five compounds with higher docking scores which satisfied the docking requirements were discovered. Among them, Ginkgetin and Dauricine were most likely to be candidates that exhibition inhibiting effect on MTP. The screening strategy in this study is relatively new to the discovery of MTP inhibitors in medicinal chemistry. Moreover, it is important to note that, lomitapide, an approved MTP inhibitor, fits well with Hypo1 as well as our homology model of MTP, which confirmed the rationality of our studies. The results indicated the applicability of molecular modeling for the discovery of potential natural MTP inhibitors from traditional Chinese herbs.     

The plasma membrane of Entamoeba histolytica: structure and dynamics.

The plasma membrane components of the parasitic protozoan Entamoeba histolytica, the causative agent of human invasive amebiasis, have been biochemically and immunologically characterized during the last decade. In addition, genes coding for certain surface proteins have been cloned. In spite of these advances, a unified characterization of plasma membrane antigenic components of the parasite is still required for badly needed advancements in the design of useful diagnostic, epidemiologic, and immunoprophylactic tools. Here we review current knowledge on this issue and address the problem of the considerable variation in the electrophoretic profiles of plasma membrane proteins obtained by different groups. In addition, the differences in the degree of recognition of reported membrane antigens with human immune sera, and the diverse interpretations concerning the possible functions of the surface molecules characterized are discussed. A comparative analysis of plasma membrane proteins of E histolytica trophozoites using three different isolation methods revealed that it is possible to select for specific membrane proteins, depending on the lysis conditions. In our view, the method of Calderón and Avila preserves more proteins than other methods tested. Using sera from recent cases of invasive amebiasis studied by several laboratories in various geographical areas, a basic antigenic pattern of 11 principal proteins with molecular weights of 220, 170, 150, 125, 97, 80, 60, 45, 20 and 9 kDA was established for the pathogenic E histolytica strain HM1:IMSS, used by most research groups.     

Pharmacophore modeling and virtual screening to identify potential RET kinase inhibitors

Chemical features based 3D pharmacophore model for REarranged during Transfection (RET) tyrosine kinase were developed by using a training set of 26 structurally diverse known RET inhibitors. The best pharmacophore hypothesis, which identified inhibitors with an associated correlation coefficient of 0.90 between their experimental and estimated anti-RET values, contained one hydrogen-bond acceptor, one hydrogen-bond donor, one hydrophobic, and one ring aromatic features. The model was further validated by a testing set, Fischer’s randomization test, and goodness of hit (GH) test. We applied this pharmacophore model to screen NCI database for potential RET inhibitors. The hits were docked to RET with GOLD and CDOCKER after filtering by Lipinski’s rules. Ultimately, 24 molecules were selected as potential RET inhibitors for further investigation.     

Elucidation of the DNA Synthetic Cycle of Entamoeba Spp. Using Flow Cytometry and Mathematical Modeling

ABSTRACT. We developed a method to study the DNA synthetic cycles of Entamoeba histolytica and Entamoeba invadens by flow cytometry (FCM) based on a preparative procedure to reduce both high levels of natural fluorescence and non-specific adsorption of fluorochromes. We modeled G₁, S, and G₂ phases as a series of overlapping Gaussian curves. Both E. histolytica and E. invadens displayed G₁, S, and G₂ proportions that are consistent with eukaryotic cell populations in exponential or stationary growth phase. Exponential phase E. histolytica populations contained a hypodiploid subset with a mass of about 20% less than the diploid value which we estimate by FCM to be 24 × 10^-14 g DNA/cell. Exponential phase E. invadens populations contained a hypodiploid subset with a mass of about 6% less than the diploid value which we estimate by FCM to be 30 × 10^-14 g DNA/cell.     

Identification of novel ASK1 inhibitors using virtual screening

Apoptosis signal-regulating kinase 1 (ASK1, also called MAP3K5) is a mitogen-activated protein kinase kinase kinase (MAP3K) that plays important roles in stress-induced cell death and inflammation, and is expected as a new therapeutic target for cancer, cardiovascular diseases, and neurodegenerative diseases. We identified novel ASK1 inhibitors by virtual screening from the public chemical library collected by Chemical Biology Research Initiative (CBRI) at the University of Tokyo.     

Identification of novel urease inhibitors: pharmacophore modeling,virtual screening and molecular docking studies

Abstract

Pharmacophore modeling and atom-based three-dimensional quantitative structure–activity relationship (3D-QSAR) have been developed on N-acylglycino- and hippurohydroxamic acid derivatives, which are known potential inhibitors of urease. This is followed by virtual screening and ADMET (absorption, distribution, metabolism, excretion and toxicity) studies on a large library of known drugs in order to get lead molecules as Helicobacter pylori urease inhibitors. A suitable three-featured pharmacophore model comprising one H-bond acceptor and two H-bond donor features (ADD.10) has been found to be the best QSAR model. An external library of compounds (～3000 molecules), pre-filtered using Lipinski’s rule of five, has been further screened using the pharmacophore model ADD.10. By analyzing the fitness of the hits with respect to the pharmacophore model and their binding interaction inside the urease active site, four molecules have been predicted to be extremely good urease inhibitors. Two of these have significant potential and should be taken up for further drug-designing process.     

Identification of potent virtual leads to design novel PLK1 inhibitors: pharmacophore modelling,virtual screening and molecular docking studies

The aim of this study was to identify novel scaffolds and utilise them in designing potent PLK1 inhibitors. Three-dimensional pharmacophore models on the basis of chemical features were developed for PLK1 on the basis of the known inhibitors. The best pharmacophore model, Hypo 1, which has the highest correlation (0.96), the highest cost difference (75.7494), the lowest total cost and RMSD (75.7494, 0.5458), contains two hydrophobics, one ring aromatic and one hydrogen donor. Hypo 1 was validated by the test set, decoy set and the Fischer's randomisation method. Then it was used for chemical database virtual screening. The hit compounds were filtered by Lipinski's rule of five and absorption, distribution, metabolism, elimination and toxicity properties. Finally, 24 compounds with good estimated activity values were used for docking studies. These results will be used to develop new inhibitors of PLK1 as leads.     

Discovery of novel small-molecule Src kinase inhibitors via a kinase-focused druglikeness rule and structure-based virtual screening

Overexpression of the non-receptor tyrosine kinase Src is implicated in the development and progression of various human cancers. Blocking signalling pathways mediated by Src is believed to be a promising anticancer strategy. We report herein the discovery of novel small-molecule Src inhibitors by crystal structure-based virtual screening. A kinase-focused druglikeness rule was proposed and used in the design of compound library. Combination of large-scale docking with DOCK and rescoring with GOLD resulted in 6 hits with moderate to potent inhibitory activity against Src. Among them, compound 1 with an IC₅₀ of 1.2 μM shows the most potent inhibitory activity. By using molecular docking, binding models of the top 3 hits (ranked by potency and ligand efficiency) with Src were constructed to provide a rational strategy that simultaneously exploits hydrogen bonding interaction in the hinge region and hydrophobic stacking in the back pocket. This approach is instructive and meaningful to further structure-based drug design of Src inhibitors.     

Pharmacophore modeling and structure-based virtual screening to identify potent inhibitors targeting LuxP of Vibrio harveyi

The main aim of the study is to identify molecules that can disrupt quorum sensing (QS) system of Vibrio harveyi and therefore perhaps the production of toxins. Recently, a novel class of dioxazaborocane derivatives has been found to block AI-2 QS by targeting LuxPQ, but the mechanism of protein inhibition is still unclear. In order to investigate the possible binding modes, all the derivatives were docked into the binding site of LuxP using induced fit docking (IFD). The computed binding affinity is in good agreement with the experimental data. Resultant protein–ligand complexes were simulated using Desmond module and the result revealed better binding of ligands in the binding site of LuxP. Both pharmacophore- and structure-based virtual screening was performed to identify novel hits against LuxP. A filtering protocol, including lipinski filters, number of rotatable bonds and three levels of docking precisions were used for the selection of hits with specific properties. The virtual screening results were then combined and analyzed, which retrieved six hits with significant Glide score, binding affinity toward LuxP. The pharmacokinetic properties of the retrieved hits are in the acceptable range. Enrichment calculation was performed to validate the final hits, to discriminate the active compounds from the inactive compounds. The identified hits could serve as a base for further drug development against LuxP of Vibrio harveyi.     

Discovery of the first SecA inhibitors using structure-based virtual screening

Bacterial protein secretion is a critical and complex process. The Sec machinery provides a major pathway for protein translocation across and integration into the cellular membrane in bacteria. Small molecule probes that perturb the functions of individual member proteins within the Sec machinery will be very important research tools as well as leads for future antimicrobial agent development. Herein we describe the discovery of inhibitors, through virtual screening, that specifically act on SecA ATPase, which is a critical member of the Sec system. These are the very first inhibitors reported for intrinsic SecA ATPase.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.     

Toward the virtual screening of Cdc25A phosphatase inhibitors with the homology modeled protein structure

Cdc25 phosphatases have been considered as attractive drug targets for anticancer therapy due to the correlation of their overexpression with a wide variety of cancers. As a method for the discovery of novel inhibitors of Cdc25 phosphatases, we have evaluated the computer-aided drug design protocol involving the homology modeling of Cdc25A and virtual screening with the two docking tools: FlexX and the modified AutoDock program implementing the effects of ligand solvation in the scoring function. The homology modeling with the X-ray crystal structure of Cdc25B as a template provides a high-quality structure of Cdc25A that enables the structure-based inhibitor design. Of the two docking programs under consideration, AutoDock is found to be more accurate than FlexX in terms of scoring putative ligands. A detailed binding mode analysis of the known inhibitors shows that they can be stabilized in the active site of Cdc25A through the simultaneous establishment of the multiple hydrogen bonds and the hydrophobic interactions. The present study demonstrates the usefulness of the modified AutoDock program as a docking tool for virtual screening of new Cdc25 phosphatase inhibitors as well as for binding mode analysis to elucidate the activities of known inhibitors. Figure Structures and available IC50 values (in μM) of the twenty Cdc25 phosphatase inhibitors seeded in docking library     

Results of Intracecal Inoculation of Germfree and Conventional Guinea Pigs and Germfree Rats with Axenically Cultivated Entamoeba histolytica

SYNOPSIS. Germfree and conventional guinea pigs and germfree rats were inoculated with large numbers of axenically cultivated Entamoeba histolytica. Amebic lesions were not found in any of the animals, and there was no indication that the ameba had established lumen infections in any instance. The axenic amebae appeared to have lost their pathogenicity. This loss was believed to depend, in part at least, upon the fact that they did not encyst and thus did not complete their life cycle in axenic culture.     

Homology modeling of human histone deacetylase 10 and design of potential selective inhibitors

Abstract

Histone deacetylases (HDACs) are implicated in the pathology of various cancers, and their pharmacological blockade has proven to be promising in reversing the malignant phenotypes. However, lack of crystal structures of some of the human HDAC isoforms (e.g., HDAC10) hinders the design of the isoform-selective inhibitor. Here, the recently solved X-ray crystal structure of Danio rerio (zebrafish) HDAC10 (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) was retrieved from the PDB and used as a template structure to model the three-dimensional structure of human HDAC10. The overall quality of the best model (M0017) was assessed by computing its z-score—a measure of the deviation of the total energy of the structure with respect to an energy distribution derived from random conformations and by docking of known HDAC10 inhibitors to its catalytic cavity. Furthermore, to identify potential HDAC10-selective inhibitor ligand-based virtual screening was carried out against the ZINC database. The free modeled structure of HDAC10 and its complexes with quisinostat and the highest-ranked compound ZINC19749069 were submitted to molecular dynamics simulation. The comparative analysis of root-mean-squared deviation, root-mean-squared fluctuation, radius of gyration (Rg), and potential energy of these systems showed that HDAC10-ZINC19749069 complex remained the most stable over time. Thus, M0017 could be potentially used for structure-based inhibitor against HDAC10, and ZINC19749069 may provide a scaffold for further optimization.

Communicated by Ramaswamy H. Sarma     

Attachment of Entamoeba histolytica to Glass in a Defined Maintenance Medium: Specific Requirement for Cysteine and Ascorbic Acid*

SYNOPSIS. Cysteine and ascorbic acid were previously shown to be required by Entamoeba histolytica trophozoites for attachment to glass, elongation, and ameboid movement as well as for short-term (12–24 h) survival in a balanced salt solution containing bovine serum albumin and a vitamin solution (Maintenance Medium 1). If the only function of cysteine and ascorbate was to decrease the redox potential, other reducing agents should be effective. However, the requirement for cysteine in the presence of ascorbic acid was highly specific. Equally effective were D- and L-cysteine; however, of many other compounds tested, only thioglycolic acid, ascorbic acid, or L-cystine (in decreasing order) were somewhat active. Under N₂ atmosphere, cysteine and ascorbic acid were still required, although their concentrations could be halved. The ability to attach in the maintenance medium was irreversibly lost after only 5 min of cysteine-ascorbic acid deprivation; however, there was no decrease in viability when the amebae were transferred to growth medium within 30 min. Cysteine thiol groups in the medium were oxidized rapidly regardless of the concentration of ascorbic acid or the presence of amebae; however, ascorbic acid prolonged attachment of amebae.     

Discovery of novel HCV polymerase inhibitors using pharmacophore-based virtual screening

We report the use of pharmacophore-based virtual screening as an efficient tool for the discovery of novel HCV polymerase inhibitors. A three-dimensional pharmacophore model for the HCV-796 binding site, NNI site IV inhibitor, to the enzyme was built by means of the structure-based focusing module in Cerius2 program. Using these models as a query for virtual screening, we produced a successful example of using pharmacophore-based virtual screening to identify novel compounds with HCV replicon assay through inhibition of HCV polymerization. Among the hit compounds, compounds 1 and 2 showed 56% and 48% inhibition of NS5B polymerization activity at 20 μM, respectively. In addition, compound 1 also exhibited replicon activity with EC₅₀ value of 2.16 μM. Following up the initial hit, we obtained derivatives of compound 1 and evaluated polymerization inhibition activity and HCV replicon assay. These results provide information necessary for the development of more potent NS5B inhibitors.