Identification of selective MMP-9 inhibitors through multiple e-pharmacophore,ligand-based pharmacophore,molecular docking,and density functional theory approaches

Matrix metalloproteinase-9 (MMP-9) is a significant target for the development of drugs for the treatment of arthritis, CNS disorders, and cancer metastasis. The structure-based and ligand-based methods were used for the virtual screening (VS) of database compounds to obtain potent and selective MMP-9 inhibitors. Experimentally known MMP-9 inhibitors were used to grow up ligand-based three pharmacophore models utilizing Schrodinger suite. The X-ray crystallographic structures of MMP-9 with different inhibitors were used to develop five energy-optimized structure-based (e-pharmacophore) models. All developed pharmacophores were validated and applied to screen the Zinc database. Pharmacophore matched compounds were subjected to molecular docking to retrieve hits with novel scaffolds. The molecules with diverse structures, high docking scores and low binding energies for various crystal structures of MMP-9, were selected as final hits. The Induced fit docking (IFD) analysis provided significant information about the driving of inhibitor to approve a suitable bioactive conformational position in the active site of protein. Since charge transfer reaction occurs during receptor–ligand interaction, therefore, electronic features of hits (ligands) are interesting parameters to explain the binding interactions. Density functional theory (DFT) at B3LYP/6-31G* level was utilized to explore electronic features of hits. The docking study of hits using AutoDock was helpful to establish the binding interactions. The study illustrates that the combined pharmacophore approach is advantageous to identify diverse hits which have better binding affinity to the active site of the enzyme for all possible bioactive conformations. The approach used in the study is worthy to design drugs for other targets.     

Structure based de novo design of novel glycogen synthase kinase 3 inhibitors

Glycogen synthase kinase 3 (GSK-3) is a serine/threonine kinase that has captured great attention in drug discovery projects. Structure based design has been successfully carried out to find a novel class of GSK-3 inhibitors using the Ludi de novo ligand design program. A total of 15 potential leads are suggested from the study. The structures have been validated through detailed analysis of the Ludi score values and by molecular docking experiment using FlexX. The hits have been further verified through: (1) visual examination of how well the hits dock into the GSK-3beta binding site; (2) comparative analysis of their FlexX, G_Score, PMF_Score, ChemScore, and D_scores values; (3) a comparative investigation of the docking scores of the hits with those of the reported inhibitors after calibration of the docking procedure with 17 previously reported inhibitors; (4) determination of the binding mode of the hits and comparison with that of the so far known inhibitors. Hits retaining interactions with the common amino acids of GSK-3beta binding site were taken to represent potential leads. Structurally the hits designed are mainly flat nitrogen heterocycles. These hits are expected to be important additions to the search of GSK-3 inhibitors and may provide invaluable insights to further understand the structural basis of catalysis and inhibition of this kinase.     

New molecular scaffolds for the design of Mycobacterium tuberculosis type II dehydroquinase inhibitors identified using ligand and receptor based virtual screening

Using ligand and receptor based virtual screening approaches we have identified potential virtual screening hits targeting type II dehydroquinase from Mycobacterium tuberculosis, an effective and validated anti-mycobacterial target. Initially, we applied a virtual screening workflow based on a combination of 2D structural fingerprints, 3D pharmacophore and molecular docking to identify compounds that rigidly match specific aspects of ligand bioactive conformation. Subsequently, the resulting compounds were ranked and prioritized using receptor interaction fingerprint based scoring and quantitative structure activity relationship model developed using already known actives. The virtual screening hits prioritized belong to several classes of molecular scaffolds with several available substitution positions that could allow chemical modification to enhance binding affinity. Finally, identified hits may be useful to a medicinal chemist or combinatorial chemist to pick up the new molecular starting points for medicinal chemistry optimization for the design of novel type II dehydroquinase inhibitors.     

Preparation and evaluation of soluble epoxide hydrolase inhibitors with improved physical properties and potencies for treating diabetic neuropathic pain

Soluble epoxide hydrolase (sEH), a novel therapeutic target for neuropathic pain, is a largely cytosolic enzyme that degrades epoxy-fatty acids (EpFAs), an important class of lipid signaling molecules. Many inhibitors of sEH have been reported, and to date, the 1,3-disubstituted urea has the highest affinity reported for the sEH among the central pharmacophores evaluated. An earlier somewhat water soluble sEH inhibitor taken to the clinic for blood pressure control had mediocre potency (both affinity and kinetics) and a short in vivo half-life. We undertook a study to overcome these difficulties, but the sEH inhibitors carrying a 1,3-disubstituted urea often suffer poor physical properties that hinder their formulation. In this report, we described new strategies to improve the physical properties of sEH inhibitors with a 1,3-disubstituted urea while maintaining their potency and drug-target residence time (a complementary in vitro parameter) against sEH. To our surprise, we identified two structural modifications that substantially improve the potency and physical properties of sEH inhibitors carrying a 1,3-disubstituted urea pharmacophore. Such improvements will greatly facilitate the movement of sEH inhibitors to the clinic.     

Discovery of novel class 1 phosphatidylinositide 3-kinases (PI3K) fragment inhibitors through structure-based virtual screening

The discovery of ligand efficient and lipophilicity efficient fragment inhibitors of class 1 phosphatidylinositide 3-kinases (PI3K) is reported. A fragment version of the AstraZeneca compound bank was docked to a homology model of the PI3K p110β isoform. Interaction-based scoring of the predicted binding poses served to further prioritise the virtual fragment hits. Experimental screening confirmed potency for a total of 18 fragment inhibitors, belonging to five different structural classes.     

Design of Benzene-1,2-diamines as selective inducible nitric oxide synthase inhibitors: a combined de novo design and docking analysis

Selective inhibition of inducible nitric oxide synthases (iNOS) has been a challenging problem for researchers pursuing work in finding methods to treat inflammatory disorders, shock, etc. Though many inhibitors have been studied to date, all are associated with selectivity or potency problems. Additionally, most of the reported compounds have several similarities and fewer number of novel structures are being tried. There is an increasing need to design novel molecules for this target. In this work, de novo design using LUDI, combined with docking analysis using FlexX has been employed in an attempt to identify novel scaffolds. Benzene-1,2-diamines were identified which could mimic the interactions of the substrate analogs and other inhibitors. Comparative docking scores in each of the isoforms of nitric oxide synthase were employed to recognize hits for iNOS selectivity. Figure Figure shows the docked poses of the ligand M226 along with that of the reference GW274150. (FlexX analysis)     

High-throughput identification of inhibitors of human mitochondrial peptide deformylase

The human mitochondrial peptide deformylase (HsPDF) provides a potential new target for broadly acting antiproliferative agents. To identify novel nonpeptidomimetic and nonhydroxamic acid-based inhibitors of HsPDF, the authors have developed a high-throughput screening (HTS) strategy using a fluorescence polarization (FP)-based binding assay as the primary assay for screening chemical libraries, followed by an enzymatic-based assay to confirm hits, prior to characterization of their antiproliferative activity against established tumor cell lines. The authors present the results and performance of the established strategy tested in a pilot screen of 2880 compounds and the identification of the 1st inhibitors. Two common scaffolds were identified within the hits. Furthermore, cytotoxicity studies revealed that most of the confirmed hits have antiproliferative activity. These findings demonstrate that the designed strategy can identify novel functional inhibitors and provide a powerful alternative to the use of functional assays in HTS and support the hypothesis that HsPDF inhibitors may constitute a new class of antiproliferative agent.     

Human soluble epoxide hydrolase: structural basis of inhibition by 4-(3-cyclohexylureido)-carboxylic acids

X-ray crystal structures of human soluble epoxide hydrolase (sEH) complexed with four different dialkylurea inhibitors bearing pendant carboxylate "tails" of varying length have been determined at 2.3-3.0 A resolution. Similarities among inhibitor binding modes reinforce the proposed roles of Y381 and/or Y465 as general acids that protonate the epoxide ring of the substrate in concert with nucleophilic attack of D333 at the electrophilic epoxide carbon. Additionally, the binding of these inhibitors allows us to model the binding mode of the endogenous substrate 14,15-epoxyeicosatrienoic acid. Contrasts among inhibitor binding modes include opposite orientations of inhibitor binding in the active-site hydrophobic tunnel. Alternative binding orientations observed for this series of inhibitors to human sEH, as well as the binding of certain dialkylurea inhibitors to human sEH and murine sEH, complicate the structure-based design of human sEH inhibitors with potential pharmaceutical applications in the treatment of hypertension. Thus, with regard to the optimization of inhibitor designs targeting human sEH, it is critical that human sEH and not murine sEH be utilized for inhibitor screening, and it is critical that structures of human sEH-inhibitor complexes be determined to verify inhibitor binding orientations that correlate with measured affinities.     

Analysis of non-peptidic compounds as potential malarial inhibitors against Plasmodial cysteine proteases via integrated virtual screening workflow

Falcipain-2 (FP-2) and falcipain-3 (FP-3), haemoglobin-degrading enzymes in Plasmodium falciparum, are validated drug targets for the development of effective inhibitors against malaria. However, no commercial drug-targeting falcipains has been developed despite their central role in the life cycle of the parasites. In this work, in silico approaches are used to identify key structural elements that control the binding and selectivity of a diverse set of non-peptidic compounds onto FP-2, FP-3 and homologues from other Plasmodium species as well as human cathepsins. Hotspot residues and the underlying non-covalent interactions, important for the binding of ligands, are identified by interaction fingerprint analysis between the proteases and 2-cyanopyridine derivatives (best hits). It is observed that the size and chemical type of substituent groups within 2-cyanopyridine derivatives determine the strength of protein–ligand interactions. This research presents novel results that can further be exploited in the structure-based molecular-guided design of more potent antimalarial drugs.     

Fragment-based screening for inhibitors of PDE4A using enthalpy arrays and X-ray crystallography

Fragment-based screening has typically relied on X-ray or nuclear magnetic resonance methods to identify low-affinity ligands that bind to therapeutic targets. These techniques are expensive in terms of material and time, so it useful to have a higher throughput method to reliably prescreen a fragment library to identify a subset of compounds for structural analysis. Calorimetry provides a label-free method to assay binding and enzymatic activity that is unaffected by the spectroscopic properties of the sample. Conventional microcalorimetry is hampered by requiring large quantities of reagents and long measurement times. Nanocalorimeters can overcome these limitations of conventional isothermal titration calorimetry. Here we have used enthalpy arrays, which are arrays of nanocalorimeters, to perform an enzyme activity-based fragment screen for competitive inhibitors of phosphodiesterase 4A (PDE4A). Several inhibitors with K ( I ) <2 mM were identified and moved to X-ray crystallization trials. Although the co-crystals did not yield high-resolution data, evidence of binding was observed, and the chemical structures of the hits were consistent with motifs of known PDE4 inhibitors. This study shows how array calorimetry can be used as a prescreening method for fragment-based lead discovery with enzyme targets and provides a list of candidate fragments for inhibition of PDE4A.     

Novel human mPGES-1 inhibitors identified through structure-based virtual screening

Microsomal prostaglandin E synthase-1 (mPGES-1) is an inducible prostaglandin E synthase after exposure to pro-inflammatory stimuli and, therefore, represents a novel target for therapeutic treatment of acute and chronic inflammatory disorders. It is essential to identify mPGES-1 inhibitors with novel scaffolds as new leads or hits for the purpose of drug design and discovery that aim to develop the next-generation anti-inflammatory drugs. Herein we report novel mPGES-1 inhibitors identified through a combination of large-scale structure-based virtual screening, flexible docking, molecular dynamics simulations, binding free energy calculations, and in vitro assays on the actual inhibitory activity of the computationally selected compounds. The computational studies are based on our recently developed three-dimensional (3D) structural model of mPGES-1 in its open state. The combined computational and experimental studies have led to identification of new mPGES-1 inhibitors with new scaffolds. In particular, (Z)-5-benzylidene-2-iminothiazolidin-4-one is a promising novel scaffold for the further rational design and discovery of new mPGES-1 inhibitors. To our best knowledge, this is the first time a 3D structural model of the open state mPGES-1 is used in structure-based virtual screening of a large library of available compounds for the mPGES-1 inhibitor identification. The positive experimental results suggest that our recently modeled trimeric structure of mPGES-1 in its open state is ready for the structure-based drug design and discovery.     

Discovery of novel Cyclophilin D inhibitors starting from three dimensional fragments with millimolar potencies

Fragment-based screening by SPR enabled the discovery of chemical diverse fragment hits with millimolar binding affinities to the peptidyl-prolyl isomerase Cyclophilin D (CypD). The CypD protein crystal structures of 6 fragment hits provided the basis for subsequent medicinal chemistry optimization by fragment merging and linking yielding three different chemical series with either urea, oxalyl or amide linkers connecting millimolar fragments in the S1′ and S2 pockets. We successfully improved the in vitro CypD potencies in the biochemical FP and PPIase assays and in the biophysical SPR binding assay from millimolar towards the low micromolar and submicromolar range by >1000-fold for some fragment derivatives. The initial SAR together with the protein crystal structures of our novel CypD inhibitors provide a suitable basis for further hit-to-lead optimization.     

Computer-aided discovery of novel non-peptide inhibitors against amyloid-beta (Aβ) peptide aggregation for treating Alzheimer's disease

A non-peptide inhibitor that is metabolically stable, orally active and capable of crossing the blood–brain barrier has been a popular option for treating Alzheimer's disease (AD). To identify novel non-peptide inhibitors for AD drug development, a structure-based pharmacophore model (SBPM) was developed using the representative docked conformation of the recently discovered peptide inhibitor PGKLVYA in the potential binding site on the Aβ(17–42) protofibril. The best SBPM, consisting of two hydrophobic, one hydrogen bond donor, and one positive ionisable feature, was further validated using ligand pharmacophore mapping studies. The well-validated SBPM was then used as the 3D query in virtual screening to identify potential hits from the National Cancer Institute database. These hits were subsequently filtered by toxicity prediction and molecular docking, and their binding stabilities and affinities were validated by 20-ns molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area analysis, respectively. Finally, two Hits (NSC35984 and NSC102747) were identified as potential leads, which exhibited higher binding stability and affinity towards Aβ compared with PGKVYA. Our results also suggest that these two Hits have the ability to prevent Aβ adopting toxic β-sheet structures, and can be easily synthesised and have structural novelty, indicating that they are promising candidates for treating AD.     

Synthesis and evaluation of heteroaryl substituted diazaspirocycles as scaffolds to probe the ATP-binding site of protein kinases

With the success of protein kinase inhibitors as drugs to target cancer, there is a continued need for new kinase inhibitor scaffolds. We have investigated the synthesis and kinase inhibition of new heteroaryl-substituted diazaspirocyclic compounds that mimic ATP. Versatile syntheses of substituted diazaspirocycles through ring-closing metathesis were demonstrated. Diazaspirocycles directly linked to heteroaromatic hinge binder groups provided ligand efficient inhibitors of multiple kinases, suitable as starting points for further optimization. The binding modes of representative diazaspirocyclic motifs were confirmed by protein crystallography. Selectivity profiles were influenced by the hinge binder group and the interactions of basic nitrogen atoms in the scaffold with acidic side-chains of residues in the ATP pocket. The introduction of more complex substitution to the diazaspirocycles increased potency and varied the selectivity profiles of these initial hits through engagement of the P-loop and changes to the spirocycle conformation, demonstrating the potential of these core scaffolds for future application to kinase inhibitor discovery.     

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies

Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.     

A Combination of Receptor-Based Pharmacophore Modeling & QM Techniques for Identification of Human Chymase Inhibitors

Inhibition of chymase is likely to divulge therapeutic ways for the treatment of cardiovascular diseases, and fibrotic disorders. To find novel and potent chymase inhibitors and to provide a new idea for drug design, we used both ligand-based and structure-based methods to perform the virtual screening(VS) of commercially available databases. Different pharmacophore models generated from various crystal structures of enzyme may depict diverse inhibitor binding modes. Therefore, multiple pharmacophore-based approach is applied in this study. X-ray crystallographic data of chymase in complex with different inhibitors were used to generate four structure–based pharmacophore models. One ligand–based pharmacophore model was also developed from experimentally known inhibitors. After successful validation, all pharmacophore models were employed in database screening to retrieve hits with novel chemical scaffolds. Drug-like hit compounds were subjected to molecular docking using GOLD and AutoDock. Finally four structurally diverse compounds with high GOLD score and binding affinity for several crystal structures of chymase were selected as final hits. Identification of final hits by three different pharmacophore models necessitates the use of multiple pharmacophore-based approach in VS process. Quantum mechanical calculation is also conducted for analysis of electrostatic characteristics of compounds which illustrates their significant role in driving the inhibitor to adopt a suitable bioactive conformation oriented in the active site of enzyme. In general, this study is used as example to illustrate how multiple pharmacophore approach can be useful in identifying structurally diverse hits which may bind to all possible bioactive conformations available in the active site of enzyme. The strategy used in the current study could be appropriate to design drugs for other enzymes as well.     

Discovery of novel inhibitors of the ZipA/FtsZ complex by NMR fragment screening coupled with structure-based design

ZipA is a membrane anchored protein in Escherichia coli that interacts with FtsZ, a homolog of eukaryotic tubulins, forming a septal ring structure that mediates bacterial cell division. Thus, the ZipA/FtsZ protein-protein interaction is a potential target for an antibacterial agent. We report here an NMR-based fragment screening approach which identified several hits that bind to the C-terminal region of ZipA. The screen was performed by 1H-15N HSQC experiments on a library of 825 fragments that are small, lead-like, and highly soluble. Seven hits were identified, and the binding mode of the best one was revealed in the X-ray crystal structure. Similar to the ZipA/FtsZ contacts, the driving force in the binding of the small molecule ligands to ZipA is achieved through hydrophobic interactions. Analogs of this hit were also evaluated by NMR and X-ray crystal structures of these analogs with ZipA were obtained, providing structural information to help guide the medicinal chemistry efforts.     

Identification of B. anthracis N5-carboxyaminoimidazole ribonucleotide mutase (PurE) active site binding compounds via fragment library screening

The de novo purine biosynthesis pathway is an attractive target for antibacterial drug design, and PurE from this pathway has been identified to be crucial for Bacillus anthracis survival in serum. In this study we adopted a fragment-based hit discovery approach, using three screening methods—saturation transfer difference nucleus magnetic resonance (STD-NMR), water-ligand observed via gradient spectroscopy (WaterLOGSY) NMR, and surface plasmon resonance (SPR), against B. anthracis PurE (BaPurE) to identify active site binding fragments by initially testing 352 compounds in a Zenobia fragment library. Competition STD NMR with the BaPurE product effectively eliminated non-active site binding hits from the primary hits, selecting active site binders only. Binding affinities (dissociation constant, K_D) of these compounds varied between 234 and 301 μM. Based on test results from the Zenobia compounds, we subsequently developed and applied a streamlined fragment screening strategy to screen a much larger library consisting of 3000 computationally pre-selected fragments. Thirteen final fragment hits were confirmed to exhibit binding affinities varying from 14 μM to 700 μM, which were categorized into five different basic scaffolds. All thirteen fragment hits have ligand efficiencies higher than 0.30. We demonstrated that at least two fragments from two different scaffolds exhibit inhibitory activity against the BaPurE enzyme.