Discovery, design and synthesis of 6H-anthra[1,9-cd]isoxazol-6-one scaffold as G9a inhibitor through a combination of shape-based virtual screening and structure-based molecular modification

Carbazole and tetrahydro-carboline derivatives as dopamine D3 receptor antagonists with the multiple antipsychotic-like properties

Dopamine D3 receptor (D3R) is implicated in multiple psychotic symptoms. Increasing the D3R selectivity over dopamine D2 receptor (D2R) would facilitate the antipsychotic treatments. Herein, novel carbazole and tetrahydro-carboline derivatives were reported as D3R selective ligands. Through a structure-based virtual screen, ZLG-25 (D3R Ki = 685 nmol/L; D2R Ki > 10,000 nmol/L) was identified as a novel D3R selective bitopic ligand with a carbazole scaffold. Scaffolds hopping led to the discovery of novel D3R-selective analogs with tetrahydro-β-carboline or tetrahydro-γ-carboline core. Further functional studies showed that most derivatives acted as hD3R-selective antagonists. Several lead compounds could dose-dependently inhibit the MK-801-induced hyperactivity. Additional investigation revealed that 23j and 36b could decrease the apomorphine-induced climbing without cataleptic reaction. Furthermore, 36b demonstrated unusual antidepressant-like activity in the forced swimming tests and the tail suspension tests, and alleviated the MK-801-induced disruption of novel object recognition in mice. Additionally, preliminary studies confirmed the favorable PK/PD profiles, no weight gain and limited serum prolactin levels in mice. These results revealed that 36b provided potential opportunities to new antipsychotic drugs with the multiple antipsychotic-like properties.

3D QSAR pharmacophore based lead identification of G9a lysine methyltransferase towards epigenetic therapeutics

The G9a, Lysine Methyltransferase that methylates the histone 3 lysine 9 (H3K9) of the nucleosome, is an excellent epigenetic target having no clinically passed inhibitor currently owing to adverse in vivo ADMET toxicities. In this work, we have carried out detailed computational investigations to find novel and safer lead against the target using advanced 3 D QSAR pharmacophore screening of databases containing more than 400000 entrees of natural compounds. The screening was conducted at different levels at increasing stringencies by employing pharmacophore mapping, druglikenesses and interaction profiles of the selected to identify potential hit compounds. The potential hits were further screened by advanced flexible docking, ADME and toxicity analysis to eight hit compounds. Based on the comparative analysis of the hits with the reference inhibitor, we identified one lead inhibitor against the G9a, having better binding efficacy and a safer ADMET profile than the reference inhibitor. Finally, the results were further verified using robust molecular dynamics simulation and MM-GBSA binding energy calculation. The natural compounds are generally considered benign due to their long human uses and this is the first attempt of in silico screening of a large natural compound library against G9a to our best knowledge. Therefore, the finding of this study may add value towards the development of epigenetic therapeutics against the G9a.Communicated by Ramaswamy H. Sarma.

Structure based virtual screening, in vitro and in vivo analysis revealed novel potent methyltransferase G9a inhibitors as prospective anti-Alzheimer's agents

G9a is a lysine methyltransferase able to di‐methylate lysine 9 of histone H3, promoting the repression of genes involved in learning and memory. Novel strategies based on synthesizing epigenetic drugs could regulate gene expression through histone post‐translational modifications and effectively treat neurodegenerative diseases, like Alzheimer's disease (AD). Here, potential G9a inhibitors were identified using a structure‐based virtual screening against G9a, followed by in vitro and in vivo screenings. First, screening methods with the AD transgenic Caenorhabditis elegans strain CL2006, showed that the toxicity/function range was safe and recovered age‐dependent paralysis. Likewise, we demonstrated that the best candidates direct target G9a by reducing H3 K9me2 in the CL2006 strain. Further characterization of these compounds involved the assessment of the blood‐brain barrier‐permeability and impact on amyloid‐β aggregation, showing promising results. Thus, we present a G9a inhibitor candidate, F, with a novel and potent structure, providing both leads in G9a inhibitor design and demonstrating their participation in reducing AD pathology.

EHMT2/G9a as an Epigenetic Target in Pediatric and Adult Brain Tumors

Epigenetic mechanisms, including post-translational modifications of DNA and histones that influence chromatin structure, regulate gene expression during normal development and are also involved in carcinogenesis and cancer progression. The histone methyltransferase G9a (euchromatic histone lysine methyltransferase 2, EHMT2), which mostly mediates mono- and dimethylation by histone H3 lysine 9 (H3K9), influences gene expression involved in embryonic development and tissue differentiation. Overexpression of G9a has been observed in several cancer types, and different classes of G9a inhibitors have been developed as potential anticancer agents. Here, we review the emerging evidence suggesting the involvement of changes in G9a activity in brain tumors, namely glioblastoma (GBM), the main type of primary malignant brain cancer in adults, and medulloblastoma (MB), the most common type of malignant brain cancer in children. We also discuss the role of G9a in neuroblastoma (NB) and the drug development of G9a inhibitors.

Recent progress in histone methyltransferase (G9a) inhibitors as anticancer agents

Epigenetics is the study of heritable changes in gene expression without changing the DNA sequence - a change in phenotype without a change in genotype. Epigenetic abnormalities can lead to serious diseases such as cancer in organisms. Histone methylation is one of the several manifestations of epigenetics, and requires specific enzymes to catalyze, for example, G9a, which is a histone methyl transferase. G9a catalyzes the methylation of histone 3 lysine 9 (H3K9) and histone 3 lysine 27 (H3K27). In addition, G9a also plays an essential role in DNA replication, damage and repair, and gene expression by regulating DNA methylation. Moreover, G9a has been found to be overexpressed in many tumor cells and is associated with the occurrence and development of tumors. Because of its unique characteristics, G9a has become a very promising target for anti-cancer agents. Over the last decade, dozens of G9a inhibitors have been discovered as potential anticancer therapeutic agents. In this review, we summarize and classify current G9a inhibitors, the challenges and future direction are also discussed in detail.

Insights for the design of protein lysine methyltransferase G9a inhibitors

The epigenetic control of gene expression could be affected by addition and/or removal of post-translational modifications such as phosphorylation, acetylation and methylation of histone proteins, as well as methylation of DNA (5-methylation on cytosines). Misregulation of these modifications is associated with altered gene expression, resulting in various disease conditions. G9a belongs to the protein lysine methyltransferases that specifically methylates the K9 residue of histone H3, leading to suppression of several tumor suppressor genes. In this review, G9a functions, role in various diseases, structural biology aspects for inhibitor design, structure-activity relationship among the reported inhibitors are discussed which could aid in the design and development of potent G9a inhibitors for cancer treatment in the future.

Discovery of a Novel Chemotype of Histone Lysine Methyltransferase EHMT1/2 (GLP/G9a) Inhibitors: Rational Design, Synthesis, Biological Evaluation, and Co-crystal Structure

Since the discovery of compound BIX01294 over 10 years ago, only a very limited number of nonquinazoline inhibitors of H3K9-specific methyltransferases G9a and G9a-like protein (GLP) have been reported. Herein, we report the identification of a novel chemotype for G9a/GLP inhibitors, based on the underinvestigated 2-alkyl-5-amino- and 2-aryl-5-amino-substituted 3 H-benzo[ e][1,4]diazepine scaffold. Our research efforts resulted in the identification 12a (EML741), which not only maintained the high in vitro and cellular potency of its quinazoline counterpart, but also displayed improved inhibitory potency against DNA methyltransferase 1, improved selectivity against other methyltransferases, low cell toxicity, and improved apparent permeability values in both parallel artificial membrane permeability assay (PAMPA) and blood-brain barrier-specific PAMPA, and therefore might potentially be a better candidate for animal studies. Finally, the co-crystal structure of GLP in complex with 12a provides the basis for the further development of benzodiazepine-based G9a/GLP inhibitors.

The Distribution of Standard Deviations Applied to High Throughput Screening

High throughput screening (HTS) assesses compound libraries for “activity” using target assays. A subset of HTS data contains a large number of sample measurements replicated a small number of times providing an opportunity to introduce the distribution of standard deviations (DSD). Applying the DSD to some HTS data sets revealed signs of bias in some of the data and discovered a sub-population of compounds exhibiting high variability which may be difficult to screen. In the data examined, 21% of 1189 such compounds were pan-assay interference compounds. This proportion reached 57% for the most closely related compounds within the sub-population. Using the DSD, large HTS data sets can be modelled in many cases as two distributions: a large group of nearly normally distributed “inactive” compounds and a residual distribution of “active” compounds. The latter were not normally distributed, overlapped inactive distributions – on both sides –, and were larger than typically assumed. As such, a large number of compounds are being misclassified as “inactive” or are invisible to current methods which could become the next generation of drugs. Although applied here to HTS, it is applicable to data sets with a large number of samples measured a small number of times.

Shape similarity guided pose prediction: lessons from D3R Grand Challenge 3

To extend the utility of ligand 3D shape similarity into pose prediction and virtual screening, we have previously developed CDVS and PoPSS methods. Both of them utilize ligand 3D shape similarity with the crystallographic ligands to improve pose prediction. While CDVS utilizes shape similarity to select suitable receptor structures for molecular docking, PoPSS places a ligand conformation of the highest shape similarity with crystal ligands into the target protein binding pocket which is then refined by side-chain repacking and Monte Carlo energy minimization. Analyses of PoPSS revealed some drawbacks in ligand conformation generation and the scoring scheme used. Moreover, as PoPSS does not sample the ligand conformation after placing it in the binding pocket, it relies solely on conformation generation methods to produce native like conformations. To address these limitations of PoPSS method, we report here a modified approach named as PoPSS-Lite, where side-chain repacking was replaced by a simple grid-based energy minimization. This modification also allowed the sampling of terminal functional groups while keeping the core scaffold fixed. Furthermore, shape similarity calculations were improved by increasing the number of ligand conformations and using a different similarity metric. The performance of PoPSS-Lite was prospectively evaluated in D3R GC3. Comparison of PoPSS-Lite demonstrated superior performance over PoPSS and CDVS with lower mean and median RMSDs. Furthermore, comparison with other D3R GC3 pose prediction submissions revealed top performance for PoPSS-Lite. Our D3R GC3 result extends our perspective that ligand 3D shape similarity with known crystallographic information can be successfully used to predict the binding pose of ligands with unknown binding modes. Our D3R GC3 results further highlight the necessity for improvement in conformer generation methods in order to improve shape similarity guided pose prediction.

Advances in the Development of Shape Similarity Methods and Their Application in Drug Discovery

Molecular similarity is a key concept in drug discovery. It is based on the assumption that structurally similar molecules frequently have similar properties. Assessment of similarity between small molecules has been highly effective in the discovery and development of various drugs. Especially, two-dimensional (2D) similarity approaches have been quite popular due to their simplicity, accuracy and efficiency. Recently, the focus has been shifted toward the development of methods involving the representation and comparison of three-dimensional (3D) conformation of small molecules. Among the 3D similarity methods, evaluation of shape similarity is now gaining attention for its application not only in virtual screening but also in molecular target prediction, drug repurposing and scaffold hopping. A wide range of methods have been developed to describe molecular shape and to determine the shape similarity between small molecules. The most widely used methods include atom distance-based methods, surface-based approaches such as spherical harmonics and 3D Zernike descriptors, atom-centered Gaussian overlay based representations. Several of these methods demonstrated excellent virtual screening performance not only retrospectively but also prospectively. In addition to methods assessing the similarity between small molecules, shape similarity approaches have been developed to compare shapes of protein structures and binding pockets. Additionally, shape comparisons between atomic models and 3D density maps allowed the fitting of atomic models into cryo-electron microscopy maps. This review aims to summarize the methodological advances in shape similarity assessment highlighting advantages, disadvantages and their application in drug discovery.

Discovery of Reversible DNA Methyltransferase and Lysine Methyltransferase G9a Inhibitors with Antitumoral in Vivo Efficacy

Using knowledge- and structure-based approaches, we designed and synthesized reversible chemical probes that simultaneously inhibit the activity of two epigenetic targets, histone 3 lysine 9 methyltransferase (G9a) and DNA methyltransferases (DNMT), at nanomolar ranges. Enzymatic competition assays confirmed our design strategy: substrate competitive inhibitors. Next, an initial exploration around our hit 11 was pursued to identify an adequate tool compound for in vivo testing. In vitro treatment of different hematological neoplasia cell lines led to the identification of molecules with clear antiproliferative efficacies (GI₅₀ values in the nanomolar range). On the basis of epigenetic functional cellular responses (levels of lysine 9 methylation and 5-methylcytosine), an acceptable therapeutic window (around 1 log unit) and a suitable pharmacokinetic profile, 12 was selected for in vivo proof-of-concept ( Nat. Commun. 2017 , 8 , 15424 ). Herein, 12 achieved a significant in vivo efficacy: 70% overall tumor growth inhibition of a human acute myeloid leukemia (AML) xenograft in a mouse model.

Computer-Aided Drug Design in Epigenetics

Epigenetic dysfunction has been widely implicated in several diseases especially cancers thus highlights the therapeutic potential for chemical interventions in this field. With rapid development of computational methodologies and high-performance computational resources, computer-aided drug design has emerged as a promising strategy to speed up epigenetic drug discovery. Herein, we make a brief overview of major computational methods reported in the literature including druggability prediction, virtual screening, homology modeling, scaffold hopping, pharmacophore modeling, molecular dynamics simulations, quantum chemistry calculation, and 3D quantitative structure activity relationship that have been successfully applied in the design and discovery of epi-drugs and epi-probes. Finally, we discuss about major limitations of current virtual drug design strategies in epigenetics drug discovery and future directions in this field.

Structure-activity relationship studies of G9a-like protein (GLP) inhibitors

Given the high homology between the protein lysine methyltransferases G9a-like protein (GLP) and G9a, it has been challenging to develop potent and selective inhibitors for either enzyme. Recently, we reported two quinazoline compounds, MS0124 and MS012, as GLP selective inhibitors. To further investigate the structure-activity relationships (SAR) of the quinazoline scaffold, we designed and synthesized a range of analogs bearing different 2-amino substitutions and evaluated their inhibition potencies against both GLP and G9a. These studies led to the identification of two new GLP selective inhibitors, 13 (MS3748) and 17 (MS3745), with 59- and 65-fold higher potency for GLP over G9a, which were confirmed by isothermal titration calorimetry (ITC). Crystal structures of GLP and G9a in complex with 13 and 17 provide insight into the interactions of the inhibitors with both proteins. In addition, we generated GLP selective inhibitors bearing a quinoline core instead of the quinazoline core.

Discovery of Potent and Selective Inhibitors for G9a-Like Protein (GLP) Lysine Methyltransferase

G9a-like protein (GLP) and G9a are highly homologous protein lysine methyltransferases (PKMTs) sharing approximately 80% sequence identity in their catalytic domains. GLP and G9a form a heterodimer complex and catalyze mono- and dimethylation of histone H3 lysine 9 and nonhistone substrates. Although they are closely related, GLP and G9a possess distinct physiological and pathophysiological functions. Thus, GLP or G9a selective small-molecule inhibitors are useful tools to dissect their distinct biological functions. We previously reported potent and selective G9a/GLP dual inhibitors including UNC0638 and UNC0642. Here we report the discovery of potent and selective GLP inhibitors including 4 (MS0124) and 18 (MS012), which are >30-fold and 140-fold selective for GLP over G9a and other methyltransferases, respectively. The cocrystal structures of GLP and G9a in the complex with either 4 or 18 displayed virtually identical binding modes and interactions, highlighting the challenges in structure-based design of selective inhibitors for either enzyme.

Use of AlCl3 in Friedel Crafts arylation type reactions and beyond: an overview on the development of unique methodologies leading to N-heteroarenes

An overview on the development of unique methodologies that highlight the use of AlCl₃ in reactions leading to new N-heteroarenes of biological significance is presented.