Novel Lead Structures for p38 MAP Kinase via FieldScreen Virtual Screening

Practical Three-Component Regioselective Synthesis of Drug-Like 3-Aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnolines as Potential Non-Covalent Multi-Targeting Inhibitors To Combat Neurodegenerative Diseases

Neurodegenerative diseases (NDs) are one of the prominent health challenges facing contemporary society, and many efforts have been made to overcome and (or) control it. In this research paper, we described a practical one-pot two-step three-component reaction between 3,4-dihydronaphthalen-1(2H)-one (1), aryl(or heteroaryl)glyoxal monohydrates (2a-h), and hydrazine monohydrate (NH2NH2•H2O) for the regioselective preparation of some 3-aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnoline derivatives (3a-h). After synthesis and characterization of the mentioned cinnolines (3a-h), the in silico multi-targeting inhibitory properties of these heterocyclic scaffolds have been investigated upon various Homo sapiens-type enzymes, including hMAO-A, hMAO-B, hAChE, hBChE, hBACE-1, hBACE-2, hNQO-1, hNQO-2, hnNOS, hiNOS, hPARP-1, hPARP-2, hLRRK-2(G2019S), hGSK-3β, hp38α MAPK, hJNK-3, hOGA, hNMDA receptor, hnSMase-2, hIDO-1, hCOMT, hLIMK-1, hLIMK-2, hRIPK-1, hUCH-L1, hPARK-7, and hDHODH, which have confirmed their functions and roles in the neurodegenerative diseases (NDs), based on molecular docking studies, and the obtained results were compared with a wide range of approved drugs and well-known (with IC50, EC50, etc.) compounds. In addition, in silico ADMET prediction analysis was performed to examine the prospective drug properties of the synthesized heterocyclic compounds (3a-h). The obtained results from the molecular docking studies and ADMET-related data demonstrated that these series of 3-aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnolines (3a-h), especially hit ones, can really be turned into the potent core of new drugs for the treatment of neurodegenerative diseases (NDs), and/or due to the having some reactionable locations, they are able to have further organic reactions (such as cross-coupling reactions), and expansion of these compounds (for example, with using other types of aryl(or heteroaryl)glyoxal monohydrates) makes a new avenue for designing novel and efficient drugs for this purpose.

Discovery of Novel BRD4 Ligand Scaffolds by Automated Navigation of the Fragment Chemical Space

Fragment-based drug discovery (FBDD) is a very effective hit identification method. However, the evolution of fragment hits into suitable leads remains challenging and largely artisanal. Fragment evolution is often scaffold-centric, meaning that its outcome depends crucially on the chemical structure of the starting fragment. Considering that fragment screening libraries cover only a small proportion of the corresponding chemical space, hits should be seen as probes highlighting privileged areas of the chemical space rather than actual starting points. We have developed an automated computational pipeline to mine the chemical space around any specific fragment hit, rapidly finding analogues that share a common interaction motif but are structurally novel and diverse. On a prospective application on the bromodomain-containing protein 4 (BRD4), starting from a known fragment, the platform yields active molecules with nonobvious scaffold changes. The procedure is fast and inexpensive and has the potential to uncover many hidden opportunities in FBDD.

Bioisosteric Replacement as a Tool in Anti-HIV Drug Design

Bioisosteric replacement is a powerful tool for modulating the drug-like properties, toxicity, and chemical space of experimental therapeutics. In this review, we focus on selected cases where bioisosteric replacement and scaffold hopping have been used in the development of new anti-HIV-1 therapeutics. Moreover, we cover field-based, computational methodologies for bioisosteric replacement, using studies from our group as an example. It is our hope that this review will serve to highlight the utility and potential of bioisosteric replacement in the continuing search for new and improved anti-HIV drugs.

5-Keto-3-cyano-2,4-diaminothiophenes as selective maternal embryonic leucine zipper kinase inhibitors

Maternal embryonic leucine zipper kinase (MELK) is involved in several key cellular processes and displays increased levels of expression in numerous cancer classes (colon, breast, brain, ovary, prostate and lung). Although no selective MELK inhibitors have yet been approved, increasing evidence suggest that inhibition of MELK would constitute a promising approach for cancer therapy. A weak high-throughput screening hit (17, IC₅₀ ≈ 5 μM) with lead-like properties was optimized for MELK inhibition. The early identification of a plausible binding mode by molecular modeling offered guidance in the choice of modifications towards compound 52 which displayed a 98 nM IC₅₀. A good selectivity profile was achieved for a representative member of the series (29) in a 486 protein kinase panel. Future elaboration of 52 has the potential to deliver compounds for further development with chemotherapeutic aims.

FABP4 inhibitors 3D-QSAR model and isosteric replacement of BMS309403 datasets

The data have been obtained from FABP4 inhibitor molecules previously published. The 120 compounds were used to build a 3D-QSAR model. The development of the QSAR model has been undertaken with the use of Forge software using the PM3 optimized structure and the experimental IC₅₀ of each compound. The QSAR model was also employed to predict the activity of 3000 new isosteric derivatives of BMS309403. The isosteric replacement was also validated by the synthesis and the biological screening of three new compounds reported in the related research article “3D-QSAR assisted identification of FABP4 inhibitors: An effective scaffold hopping analysis/QSAR evaluation” (Floresta et al., 2019).

Identification of new pyrrolo[2,3-d]pyrimidines as potent VEGFR-2 tyrosine kinase inhibitors: Design, synthesis, biological evaluation and molecular modeling

Vascular endothelial growth factor receptor-2 (VEGFR-2) plays a crucial role in cancer angiogenesis. In the current study, a series of novel pyrrolo[2,3-d]pyrimidine based-compounds was designed and synthesized as VEGFR-2 inhibitors, in accordance to the structure activity relationship (SAR) studies of known type II VEGFR-2 inhibitors. The newly synthesized compounds were evaluated for their ability to inhibit VEGFR-2 kinase enzyme in vitro. All the tested compounds demonstrated highly potent dose-related VEGFR-2 inhibition with IC₅₀ values in nanomolar range. Among these compounds, pyrrolo[2,3-d]pyrimidine derivatives carrying biaryl urea moieties (12d and 15c) exhibited IC₅₀ values of 11.9 and 13.6 nM respectively. Additionally, most of the newly synthesized final compounds were tested on 60 human cancer cell lines. Docking of these compounds into the inactive conformation of VEGFR-2 was performed which showed comparable binding modes to that of the FDA approved VEGFR-2 kinase inhibitors. These newly discovered potent kinase inhibitors could be considered as potential candidates for the development of new targeted anticancer agent.

Computer-aided drug design, synthesis and identification of disulfide compounds as novel and potential allosteric PAK1 inhibitors

p21-activated kinase 1 (PAK1) is an evolutionarily conserved serine/threonine protein kinase, which has been considered as one of the key regulatory factors in signaling network of tumor cells. Therefore, inhibition of PAK1 may be a potential approach to treat many types of solid tumors. Several allosteric inhibitors of PAK1 have been identified, and the most well known one is IPA-3. But its biological activity is not satisfied, and the structure activity relationship (SAR) of PAK1 allosteric inhibitors is unclear. In this study, we designed and synthesized 13 potential allosteric inhibitors by using computer-aided drug design based on the structure of the existing PAK1 allosteric inhibitors. All the compounds were characterized by ¹H-NMR and ¹³C-NMR, among which six were not reported previously. SAR was investigated by pharmacological studies and In03 and In06 showed increased PAK1 inhibition than previously reported IPA-3. These findings could guide further structure optimization of PAK1 inhibitors.

p21-activated kinase 1 (PAK1) is an evolutionarily conserved serine/threonine protein kinase, which has been considered as one of the key regulatory factors in signaling network of tumor cells.

Catalytic Dearomative Spirocyclization via Gold Carbene Species Derived from Ynamides: Efficient Synthesis of 2-Azaspiro[4.5]decan-3-ones

An intramolecular catalytic dearomatization of phenols via gold carbene species proceeded to provide 2-azaspiro[4.5]decan-3-ones. The use of NHC ligand and water as a co-solvent was critical for achieving high reactivity. This reaction did not require hazardous diazo compounds as carbene sources and proceeded even under air. The obtained spirocyclic product could be readily transformed into a gabapentin derivative by hydrogenation and deprotection.

Discovery of Compounds that Positively Modulate the High Affinity Choline Transporter

Cholinergic hypofunction is associated with decreased attention and cognitive deficits in the central nervous system in addition to compromised motor function. Consequently, stimulation of cholinergic neurotransmission is a rational therapeutic approach for the potential treatment of a variety of neurological conditions. High affinity choline uptake (HACU) into acetylcholine (ACh)-synthesizing neurons is critically mediated by the sodium- and pH-dependent high-affinity choline transporter (CHT, encoded by the SLC5A7 gene). This transporter is comparatively well-characterized but otherwise unexplored as a potential drug target. We therefore sought to identify small molecules that would enable testing of the hypothesis that positive modulation of CHT mediated transport would enhance activity-dependent cholinergic signaling. We utilized existing and novel screening techniques for their ability to reveal both positive and negative modulation of CHT using literature tools. A screening campaign was initiated with a bespoke compound library comprising both the Pfizer Chemogenomic Library (CGL) of 2,753 molecules designed specifically to help enable the elucidation of new mechanisms in phenotypic screens and 887 compounds from a virtual screening campaign to select molecules with field-based similarities to reported negative and positive allosteric modulators. We identified a number of previously unknown active and structurally distinct molecules that could be used as tools to further explore CHT biology or as a starting point for further medicinal chemistry.

Discovery of Non-ATP-Competitive Inhibitors of Polo-like Kinase 1

Polo-like kinase 1 (Plk1) is an evolutionarily conserved serine/threonine kinase, and its N-terminal kinase domain (KD) controls cell signaling through phosphorylation. Inhibitors of Plk1 are potential anticancer drugs. Most known Plk1 KD inhibitors are ATP-competitive compounds, which may suffer from low selectivity. In this study we discovered novel non-ATP-competitive Plk1 KD inhibitors by virtual screening and experimental studies. Potential binding sites in Plk1 KD were identified by using the protein binding site detection program Cavity. The identified site was subjected to molecular-docking-based virtual screening. The activities of top-ranking compounds were evaluated by in vitro enzyme assay with full-length Plk1 and direct binding assay with Plk1 KD. Several compounds showed inhibitory activity, and the most potent was found to be 3-((2-oxo-2-(thiophen-2-yl)ethyl)thio)-6-(pyridin-3-ylmethyl)-1,2,4-triazin-5(4H)-one (compound 4) with an IC50 value of 13.1 ± 1.7 μm. Our work provides new insight into the design of kinase inhibitors that target non-ATP binding sites.

Core chemotype diversification in the HIV-1 entry inhibitor class using field-based bioisosteric replacement

Demand remains for new inhibitors of HIV-1 replication and the inhibition of HIV-1 entry is an extremely attractive therapeutic approach. Using field-based bioisosteric replacements, we have further extended the chemotypes available for development in the HIV-1 entry inhibitor class. Moreover, using field-based disparity analysis of the compounds, 3D structure-activity relationships were derived that will be useful in the further development of these inhibitors towards clinical utility.

Discovery and optimization of novel small-molecule HIV-1 entry inhibitors using field-based virtual screening and bioisosteric replacement

With the emergence of drug-resistant strains and the cumulative toxicities associated with current therapies, demand remains for new inhibitors of HIV-1 replication. The inhibition of HIV-1 entry is an attractive, yet underexploited therapeutic approach with implications for salvage and preexposure prophylactic regimens, as well as topical microbicides. Using the combination of a field-derived bioactive conformation template to perform virtual screening and iterative bioisosteric replacements, coupled with in silico predictions of absorption, distribution, metabolism, and excretion, we have identified new leads for HIV-1 entry inhibitors.

Insight into the Structural Determinants of Imidazole Scaffold-Based Derivatives as TNF-α Release Inhibitors by in Silico Explorations

Presently, 151 widely-diverse pyridinylimidazole-based compounds that show inhibitory activities at the TNF-α release were investigated. By using the distance comparison technique (DISCOtech), comparative molecular field analysis (CoMFA), and comparative molecular similarity index analysis (CoMSIA) methods, the pharmacophore models and the three-dimensional quantitative structure-activity relationships (3D-QSAR) of the compounds were explored. The proposed pharmacophore model, including two hydrophobic sites, two aromatic centers, two H-bond donor atoms, two H-bond acceptor atoms, and two H-bond donor sites characterizes the necessary structural features of TNF-α release inhibitors. Both the resultant CoMFA and CoMSIA models exhibited satisfactory predictability (with Q² (cross-validated correlation coefficient) = 0.557, R²_ncv (non-cross-validated correlation coefficient) = 0.740, R²_pre (predicted correlation coefficient) = 0.749 and Q² = 0.598, R²_ncv = 0.767, R²_pre = 0.860, respectively). Good consistency was observed between the 3D-QSAR models and the pharmacophore model that the hydrophobic interaction and hydrogen bonds play crucial roles in the mechanism of actions. The corresponding contour maps generated by these models provide more diverse information about the key intermolecular interactions of inhibitors with the surrounding environment. All these models have extended the understanding of imidazole-based compounds in the structure-activity relationship, and are useful for rational design and screening of novel 2-thioimidazole-based TNF-α release inhibitors.

Accurate calculation of mutational effects on the thermodynamics of inhibitor binding to p38α MAP kinase: a combined computational and experimental study

A major current challenge for drug design efforts focused on protein kinases is the development of drug resistance caused by spontaneous mutations in the kinase catalytic domain. The ubiquity of this problem means that it would be advantageous to develop fast, effective computational methods that could be used to determine the effects of potential resistance-causing mutations before they arise in a clinical setting. With this long-term goal in mind, we have conducted a combined experimental and computational study of the thermodynamic effects of active-site mutations on a well-characterized and high-affinity interaction between a protein kinase and a small-molecule inhibitor. Specifically, we developed a fluorescence-based assay to measure the binding free energy of the small-molecule inhibitor, SB203580, to the p38α MAP kinase and used it measure the inhibitor's affinity for five different kinase mutants involving two residues (Val38 and Ala51) that contact the inhibitor in the crystal structure of the inhibitor-kinase complex. We then conducted long, explicit-solvent thermodynamic integration (TI) simulations in an attempt to reproduce the experimental relative binding affinities of the inhibitor for the five mutants; in total, a combined simulation time of 18.5 μs was obtained. Two widely used force fields - OPLS-AA/L and Amber ff99SB-ILDN - were tested in the TI simulations. Both force fields produced excellent agreement with experiment for three of the five mutants; simulations performed with the OPLS-AA/L force field, however, produced qualitatively incorrect results for the constructs that contained an A51V mutation. Interestingly, the discrepancies with the OPLS-AA/L force field could be rectified by the imposition of position restraints on the atoms of the protein backbone and the inhibitor without destroying the agreement for other mutations; the ability to reproduce experiment depended, however, upon the strength of the restraints' force constant. Imposition of position restraints in corresponding simulations that used the Amber ff99SB-ILDN force field had little effect on their ability to match experiment. Overall, the study shows that both force fields can work well for predicting the effects of active-site mutations on small molecule binding affinities and demonstrates how a direct combination of experiment and computation can be a powerful strategy for developing an understanding of protein-inhibitor interactions.

Virtual screening using a conformationally flexible target protein: models for ligand binding to p38α MAPK

We have used virtual screening to develop models for the binding of aryl substituted heterocycles to p38α MAPK. Virtual screening was conducted on a number of p38α MAPK crystal structures using a library of 46 known p38α MAPK inhibitors containing a heterocyclic core substituted by pyridine and fluorophenyl rings (structurally related to SB203580) and a set of decoy compounds. Multiple protonation states and tautomers of active and decoy compounds were considered. Each docking model was evaluated using receiver operating characteristic (ROC) curves and enrichment factors. The two best performing single crystal structures were found to be 1BL7 and 2EWA, with enrichment factors of 14.1 and 13.0 at 2% of the virtual screen respectively. Ensembles of up to four receptors of similar conformations were generated, generally giving good or very good performances with high ROC AUCs and good enrichment. The 1BL7-2EWA ensemble was able to outperform each of its constituent receptors and gave high enrichment factors of 17.3, 12.0, 8.0 at 2, 5 and 10% respectively, of the virtual screen. A ROC AUC of 0.94 was obtained for this ensemble. This method may be applied to other proteins where there are a large number of inhibitor classes with different binding site conformations.

Molecular dynamics simulation and free energy calculation studies of the binding mechanism of allosteric inhibitors with p38α MAP kinase

p38 MAP kinase is a promising target for anti-inflammatory treatment. The classical kinase inhibitors imatinib and sorafenib as well as BI-1 and BIRB-796 were reported to bind in the DFG-out form of human p38α, known as type II or allosteric kinase inhibitors. Although DFG-out conformation has attracted great interest in the design of type II kinase inhibitors, the structural requirements for binding and mechanism of stabilization of DFG-out conformation remain unclear. As allosteric inhibition is important to the selectivity of kinase inhibitor, herein the binding modes of imatinib, sorafenib, BI-1 and BIRB-796 to p38α were investigated by molecular dynamics simulation. Binding free energies were calculated by molecular mechanics/Poisson-Boltzmann surface area method. The predicted binding affinities can give a good explanation of the activity difference of the studied inhibitors. Furthermore, binding free energies decomposition analysis and further structural analysis indicate that the dominating effect of van der Waals interaction drives the binding process, and key residues, such as Lys53, Gly71, Leu75, Ile84, Thr106, Met109, Leu167, Asp168, and Phe169, play important roles by forming hydrogen bond, salt bridge, and hydrophobic interactions with the DFG-out conformation of p38α. Finally, we also conducted a detailed analysis of BI-1, imatinib, and sorafenib binding to p38α in comparison with BIRB-796 exploited for gaining potency as well as selectivity of p38 inhibitors. These results are expected to be useful for future rational design of novel type II p38 inhibitors.

Microwave assisted synthesis of novel pyrazolone derivatives attached to a pyrimidine moiety and evaluation of their anti-inflammatory, analgesic and antipyretic activities

In the present research work, the motto was to develop new chemical entities as potential anti-inflammatory, analgesic and antipyretic agents. Various 4-(2-amino-6-(substituted)pyrimidin-4-yl)-3-methyl-1-(substituted)-1H-pyrazol-5(4H)-one derivatives (5a-5j) and their Schiff bases (6a-6j) were synthesized. The newly synthesized compounds were characterized by TLC and spectral data. The compounds containing pyrazolone and amino pyrimidine as basic moieties (5a-5j), were screened for their anti-inflammatory, analgesic and antipyretic activities, compounds 5a, 5c-5f, 5h exhibited activities nearly similar to the standard. The pharmacological studies reveal that the presence of 4-hydroxy, 4-methoxy, 4-(N,N-dimethylamino) or 2-hydroxy groups on phenyl ring at C6 of amino pyrimidine exhibits anti-inflammatory, analgesic and antipyretic activities nearly similar to the standard and substitutions like 4-chloro, 2-nitro, 3-nitro or 4-nitro on same phenyl ring lead to a decrease in activities.

Selectivity of kinase inhibitor fragments

A kinase-focused screening set of fragments has been assembled and has proved successful for the discovery of ligand-efficient hits against many targets. Here we present some of our general conclusions from this exercise. Notably, we present the first profiling results for literature fragments that have previously been used as starting points for optimization against individual kinases. We consider the importance of screening format and the extent to which selectivity is helpful in selecting fragments for progression. Results are also outlined for fragments targeting the DFG-out conformation and for atypical kinases such as PIM1 and lipid kinases.

The Interweaving of Cheminformatics and HTS

The aim of this chapter is to describe the stages of early drug discovery that can be assisted by techniques commonly used in the field of cheminformatics. In fact, cheminformatics tools can be applied all the way from the design of compound libraries and the analysis of HTS results, to the discovery of functional relationships between compounds and their targets.

Modeling approaches for ligand-based 3D similarity

3D ligand-based similarity approaches are widely used in the early phases of drug discovery for tasks such as hit finding by virtual screening or compound design with quantitative structure–activity relationships. Here in we review widely used software for performing such tasks. Some techniques are based on relatively mature technology, shape-based similarity for instance. Typically, these methods remained in the realm of the expert user, the experienced modeler. However, advances in implementation and speed have improved usability and allow these methods to be applied to databases comprising millions of compounds. There are now many reports of such methods impacting drug-discovery projects. As such, the medicinal chemistry community has become the intended market for some of these new tools, yet they may consider the wide array and choice of approaches somewhat disconcerting. Each method has subtle differences and is better suited to certain tasks than others. In this article we review some of the widely used computational methods via application, provide straightforward background on the underlying theory and provide examples for the interested reader to pursue in more detail. In the new era of preclinical drug discovery there will be ever more pressure to move faster and more efficiently, and computational approaches based on 3D ligand similarity will play an increasing role in in this process.

Inhibitors of Leishmania GDP-mannose pyrophosphorylase identified by high-throughput screening of small-molecule chemical library

The current treatment for leishmaniasis is based on chemotherapy, which relies on a handful of drugs with serious limitations, such as high cost, toxicity, and a lack of efficacy in regions of endemicity. Therefore, the development of new, effective, and affordable antileishmanial drugs is a global health priority. Leishmania synthesizes a range of mannose-rich glycoconjugates that are essential for parasite virulence and survival. A prerequisite for glycoconjugate biosynthesis is the conversion of monosaccharides to the activated mannose donor, GDP-mannose, the product of a reaction catalyzed by GDP-mannose pyrophosphorylase (GDP-MP). The deletion of the gene encoding GDP-MP in Leishmania led to a total loss of virulence, indicating that the enzyme is an ideal drug target. We developed a phosphate sensor-based high-throughput screening assay to quantify the activity of GDP-MP and screened a library containing approximately 80,000 lead-like compounds for GDP-MP inhibitors. On the basis of their GDP-MP inhibitory properties and chemical structures, the activities of 20 compounds which were not toxic to mammalian cells were tested against ex vivo amastigotes and in macrophage amastigote assays. The most potent compound identified in the primary screen (compound 3), a quinoline derivative, demonstrated dose-dependent activity in both assays (50% inhibitory concentration = 21.9 microM in the macrophage assay) and was shown to be nontoxic to human fibroblasts. In order to elucidate signs of an early structure-activity relationship (SAR) for this class of compounds, we obtained and tested analogues of compound 3 and undertook limited medicinal chemistry optimization, which included the use of a number of SAR probes of the piperazinyl aryl substituent of compound 3. We have identified novel candidate compounds for the design and synthesis of antileishmanial therapeutics.