In silico hit optimization toward AKT inhibition: fragment-based approach,molecular docking and molecular dynamics study

Abstract

Protein kinase B also known as AKT is a cardinal node in different signaling pathways that regulates diverse cell processes. AKT has three isoforms that share high homology. Hyperactivation of each isoform is related with different types of cancer. This work describes the computational search for new inhibitors using a hit optimization process of the previously reported AKT pan inhibitor, a 2,4,6-trisubstituted pyridine. A database of new molecules was proposed using a variant of fragment-based docking methodology and previous reported considerations. Molecular docking followed by molecular dynamics studies were performed to select the best compounds and analyze their behavior. Protein–ligand complexes energy was calculated using molecular mechanics Poisson–Boltzmann surface area protocol. Further, proposed molecules were compared with the ChEMBL database of compounds assayed against AKT. Data analysis leads to determine the structural requirements necessary for a favorable interaction of the proposed ligands with the AKT pocket. Molecular dynamics data suggested that the pK_a of the ligands is important for the stability in the AKT pocket. Molecular similarity analysis shows that proposed ligands have not been previously reported. Thus, ligands with high docking scores and favorable behavior on molecular dynamics simulations are proposed as potential AKT inhibitors.     

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.     

Selective inhibition of human cathepsin S by 2,4,6-trisubstituted 1,3,5-triazine analogs

We report herein the synthesis and biological evaluation of a new series of 2,4,6-trisubstituted 1,3,5-triazines as reversible inhibitors of human cysteine cathepsins. The desired products bearing morpholine and N-Boc piperidine, respectively, were obtained in three to four steps from commercially available trichlorotriazine. Seventeen hitherto unknown compounds were evaluated in vitro against various cathepsins for their inhibitory properties. Among them, compound 7c (4-(morpholin-4-yl)-6-[4-(trifluoromethoxy)anilino]-1,3,5-triazine-2-carbonitrile) was identified as the most potent and selective inhibitor of cathepsin S (Ki??=??2??±??0.3?nM). Also 7c impaired the autocatalytic maturation of procathepsin S. Molecular docking studies support that 7c bound within the active site of cathepsin S, by interacting with Gly23, Cys25 and Trp26 (S1 subsite), with Asn67, Gly69 and Phe70 (S2 subsite) and with Gln19 (S1′ pocket).     

Carbonic anhydrase inhibitors. The X-ray crystal structure of human isoform II in adduct with an adamantyl analogue of acetazolamide resides in a less utilized binding pocket than most hydrophobic inhibitors

We investigated the inhibitory activity of several 1,3,4-thiadiazole-sulfonamides against all catalytically active CA (EC 4.2.1.1), CA I–XV. The tail derivatizing the 5-position in the 1,3,4-thiadiazole-2-sulfonamide scaffold was observed to be critical as an inhibitory determinant of these compounds. The high resolution X-ray crystal structure of hCA II in complex with 5-(1-adamantylcarboxamido)-1,3,4-thiadiazole-2-sulfonamide, showed the adamantyl moiety of the inhibitor residing in a less utilized binding pocket than that of most hydrophobic inhibitors, lined by the amino acid residues Ile91, Val121 and Phe131. This binding site may explain the diverse inhibition profiles of 5-carboxamide- and sufonamide-derivatized 1,3,4-thiadiazole-2-sulfonamides and offers a hot spot for designing isoform selective inhibitors, considering that residues 91 and 131 are highly variable among the 13 catalytically active isoforms.     

Kinetic mechanism of AKT/PKB enzyme family

AKT/PKB is a phosphoinositide-dependent serine/threonine protein kinase that plays a critical role in the signal transduction of receptors. It also serves as an oncogene in the tumorigenesis of cancer cells when aberrantly activated by genetic lesions of the PTEN tumor suppressor, phosphatidylinositol 3-kinase, and receptor tyrosine kinase overexpression. Here we have characterized and compared kinetic mechanisms of the three AKT isoforms. Initial velocity studies revealed that all AKT isozymes follow the sequential kinetic mechanism by which an enzyme-substrate ternary complex forms before the product release. The empirically derived kinetic parameters are apparently different among the isoforms. AKT2 showed the highest Km value for ATP, and AKT3 showed the highest kcat value. The patterns of product inhibition of AKT1, AKT2, and AKT3 by ADP were all consistent with an ordered substrate addition mechanism with ATP binding to the enzymes prior to the peptide substrate. Further analysis of steady state kinetics of AKT1 in the presence of dead-end inhibitors supported the finding and suggested that the AKT family of kinases catalyzes reactions via an Ordered Bi Bi sequential mechanism with ATP binding to the enzyme prior to peptide substrate and ADP being released after the phosphopeptide product. These results suggest that ATP is an initiating factor for the catalysis of AKT enzymes and may play a role in the regulation AKT enzyme activity in cells.     

Synthesis and biological evaluation of novel 2,4,6-triazine derivatives as antimicrobial agents

A series of 2,4,6-trisubstituted [1,3,5]triazines were synthesized and evaluated for their antimicrobial activity against two representative Gram-positive, Gram-negative bacteria and two fungi. Biological data revealed that among all the compounds screened, compounds 3f, 3g, 3h, 3i, 3m, 3o and 3p found to have promising antimicrobial activity against all the selected pathogenic bacteria and fungi. Out of the synthesized compounds seven analogues have shown MIC in the range of 6.25-12.5 μg/mL. These compounds were generally nontoxic and may prove useful as antimicrobial agents.     

2,3,5-Trisubstituted pyridines as selective AKT inhibitors. Part II: Improved drug-like properties and kinase selectivity from azaindazoles

A novel series of AKT inhibitors containing 2,3,5-trisubstituted pyridines with novel azaindazoles as hinge binding elements are described. Among these, the 4,7-diazaindazole compound 2c has improved drug-like properties and kinase selectivity than those of indazole 1, and displays greater than 80% inhibition of GSK3β phosphorylation in a BT474 tumor xenograft model in mice.     

Identification through structure-based methods of a bacterial NAD+-dependent DNA ligase inhibitor that avoids known resistance mutations

In an attempt to identify novel inhibitors of NAD⁺-dependent DNA ligase (LigA) that are not affected by a known resistance mutation in the adenosine binding pocket, a detailed analysis of the binding sites of a variety of bacterial ligases was performed. This analysis revealed several similarities to the adenine binding region of kinases, which enabled a virtual screen of known kinase inhibitors. From this screen, a thienopyridine scaffold was identified that was shown to inhibit bacterial ligase. Further characterization through structure and enzymology revealed the compound was not affected by a previously disclosed resistance mutation in Streptococcus pneumoniae LigA, Leu75Phe. A subsequent medicinal chemistry program identified substitutions that resulted in an inhibitor with moderate activity across various Gram-positive bacterial LigA enzymes.     

Discovery of novel pyruvate dehydrogenase kinases inhibitors by screening of an in-house small molecule library for anti-lung cancer therapeutics

Pyruvate dehydrogenase kinases (PDKs) are widely over-expressed in various human solid cancers, making them attractive therapeutic targets for cancer treatment. Herein, we report the identification of structurally novel PDKs inhibitors by screening of an in-house small molecule library. Biochemical assay indicated that the identified compounds 1–4 inhibited PDK1 activity with EC₅₀ values of 0.50, 1.99, 4.64, and 0.42?µM, respectively. The ITC analysis suggested that the identified compounds 1–4 were pan-isoform PDK inhibitors, which bound to and inhibited the four PDK isoforms. Moreover, 1–4 dose-dependently reduced pyruvate dehydrogenase complex phosphorylation in NCI-H1975 cell. Molecular docking suggested that the most potent compound 4 docked well in the ATP binding pocket of the four PDK isoforms, forming direct hydrogen bond interactions with the conserved amino acids Thr and Asp in ATP binding pocket of PDKs. The cell viability assay demonstrated that 4 potently blocked NCI-H1975 cell proliferation (IC₅₀?=?3.32?µM), but had little effect on human normal lung cell MRC-5 even with the tested concentration up to 40?µM. All the data demonstrated that 4 was a promising lead for the development of structurally novel PDKs inhibitor for the cancer treatment.     

Modeling of Arylamide Helix Mimetics in the p53 Peptide Binding Site of hDM2 Suggests Parallel and Anti-Parallel Conformations Are Both Stable

The design of novel α-helix mimetic inhibitors of protein-protein interactions is of interest to pharmaceuticals and chemical genetics researchers as these inhibitors provide a chemical scaffold presenting side chains in the same geometry as an α-helix. This conformational arrangement allows the design of high affinity inhibitors mimicking known peptide sequences binding specific protein substrates. We show that GAFF and AutoDock potentials do not properly capture the conformational preferences of α-helix mimetics based on arylamide oligomers and identify alternate parameters matching solution NMR data and suitable for molecular dynamics simulation of arylamide compounds. Results from both docking and molecular dynamics simulations are consistent with the arylamides binding in the p53 peptide binding pocket. Simulations of arylamides in the p53 binding pocket of hDM2 are consistent with binding, exhibiting similar structural dynamics in the pocket as simulations of known hDM2 binders Nutlin-2 and a benzodiazepinedione compound. Arylamide conformations converge towards the same region of the binding pocket on the 20 ns time scale, and most, though not all dihedrals in the binding pocket are well sampled on this timescale. We show that there are two putative classes of binding modes for arylamide compounds supported equally by the modeling evidence. In the first, the arylamide compound lies parallel to the observed p53 helix. In the second class, not previously identified or proposed, the arylamide compound lies anti-parallel to the p53 helix.     

Structural determinants of phosphoinositide 3-kinase inhibition by wortmannin, LY294002, quercetin, myricetin, and staurosporine

The specific phosphoinositide 3-kinase (PI3K) inhibitors wortmannin and LY294002 have been invaluable tools for elucidating the roles of these enzymes in signal transduction pathways. The X-ray crystallographic structures of PI3Kgamma bound to these lipid kinase inhibitors and to the broad-spectrum protein kinase inhibitors quercetin, myricetin, and staurosporine reveal how these compounds fit into the ATP binding pocket. With a nanomolar IC50, wortmannin most closely fits and fills the active site and induces a conformational change in the catalytic domain. Surprisingly, LY294002 and the lead compound on which it was designed, quercetin, as well as the closely related flavonoid myricetin bind PI3K in remarkably different orientations that are related to each other by 180 degrees rotations. Staurosporine/PI3K interactions are reminiscent of low-affinity protein kinase/staurosporine complexes. These results provide a rich basis for development of isoform-specific PI3K inhibitors with therapeutic potential.     

Mono-/dihydroxybenzoic acid esters and phenol pyridinium derivatives as inhibitors of the mammalian carbonic anhydrase isoforms I,II, VII,IX, XII and XIV

Using hydroxy-/dihydroxybenzoic acids as leads, a series of methyl, ethyl and iso-propyl esters of 4-hydroxy-benzoic acid, 2,4-, 2,5-, 2,6-, 3,4-, and 3,5-dihydroxybenzoic acids and of coumaric acid, were obtained and investigated for the inhibition of six mammalian carbonic anhydrase (CA, EC 4.2.1.1) isoforms, that is, the cytosolic CA I, II and VII, and the transmembrane CA IX, XII and XIV, many of which are established drug targets. Other compounds incorporating phenol/catechol moieties were obtained from dopamine by reaction with fluorescein isothiocyanate or with 2,4,6-trisubstituted pyrylium salts. Some aminophenols were also derivatized in a similar manner, by using pyrylium salts. Many of these compounds showed increased inhibitory action compared to the lead compounds from which they were obtained, with efficacy in the submicromolar range against most investigated CA isoforms. As phenols are a class of less investigated CA inhibitors (CAIs) compared to the sulfonamides, and their mechanism of inhibition is less well understood, compounds of the type designed here may be helpful in gaining more insights into these phenomena.     

Discovery of 2,3,5-trisubstituted pyridine derivatives as potent Akt1 and Akt2 dual inhibitors

This letter describes the discovery of a novel series of dual Akt1/Akt2 kinase inhibitors, based on a 2,3,5-trisubstituted pyridine scaffold. Compounds from this series, which contain a 5-tetrazolyl moiety, exhibit more potent inhibition of Akt2 than Akt1.     

p38α MAP Kinase C-Terminal Domain Binding Pocket Characterized by Crystallographic and Computational Analyses

The mitogen-activated protein (MAP) kinase protein family has a critical role in cellular signaling events, with MAP kinase p38α acting in inflammatory processes and being an important drug discovery target. MAP kinase drug design efforts have focused on small-molecule inhibitors of the ATP catalytic site, which exhibit dose-limiting adverse effects. Therefore, characterizing other potential sites that bind substrates, inhibitors, or allosteric effectors is of great interest. Here, we present the crystal structure of human p38α MAP kinase, which has a lead compound bound both in the active site and in the lipid-binding site of the C-terminal cap. This C-terminal cap is formed from an extension to the kinase fold, unique to the MAP kinase and cyclin-dependent kinase families and glycogen synthase kinase 3. Binding of this lead, 4-[3-(4-fluorophenyl)-1H-pyrazol-4-yl]pyridine, to wild-type p38α induces movement of the C-terminal cap region, creating a hydrophobic pocket centered around residue Trp197. Computational analysis of this C-terminal domain pocket indicates notable flexibility for potentially binding different-shaped compounds, including lipids, oxidized arachidonic acid species such as leukotrienes, and small-molecule effectors. Furthermore, our structural results defining the open p38α C-lobe pocket provide a detailed framework for the design of novel small molecules with affinities comparable to active-site binders: to bind and potentially modulate the shape and activity of p38α in predetermined ways. Moreover, these results and analyses of p38α suggest strategies for designing specific binding compounds applicable to other MAP kinases, as well as the cyclin-dependent kinase family and glycogen synthase kinase 3β that also utilize the C-terminal insert in their interactions.     

Benzenesulfonamides with pyrimidine moiety as inhibitors of human carbonic anhydrases I,II, VI,VII, XII,and XIII

Two groups of benzenesulfonamide derivatives, bearing pyrimidine moieties, were designed and synthesized as inhibitors of carbonic anhydrases (CA). Their binding affinities to six recombinant human CA isoforms I, II, VI, VII, XII, and XIII were determined by the thermal shift assay (TSA). The binding of several inhibitors was measured by isothermal titration calorimetry (ITC). Direct demonstration of compound inhibition was achieved by determining the inhibition constant by stopped-flow CO₂ hydration assay. The most potent compounds demonstrated selectivity towards isoform I and affinities of 0.5 nM. The crystal structures of selected compounds in complex with CA II, XII, and XIII were determined to atomic resolution. Compounds described here were compared with previously published pyrimidinebenzenesulfonamides.¹ Systematic structure–activity analysis of 40 compound interactions with six isoforms yields clues for the design of compounds with greater affinities and selectivities towards target CA isoforms.     

X-ray structure of a novel matrix metalloproteinase inhibitor complexed to stromelysin

A new class of matrix metalloproteinase (MMP) inhibitors has been identified by screening a collection of compounds against stromelysin. The inhibitors, 2,4,6-pyrimidine triones, have proven to be potent inhibitors of gelatinases A and B. An X-ray crystal structure of one representative compound bound to the catalytic domain of stromelysin shows that the compounds bind at the active site and ligand the active-site zinc. The pyrimidine triones mimic substrates in forming hydrogen bonds to key residues in the active site, and provide opportunities for placing appropriately chosen groups into the S1' specificity pocket of MMPS: A number of compounds have been synthesized and assayed against stromelysin, and the variations in potency are explained in terms of the binding mode revealed in the X-ray crystal structure.     

Discovery of PDK1 Kinase Inhibitors with a Novel Mechanism of Action by Ultrahigh Throughput Screening

The phosphoinositide 3-kinase/AKT signaling pathway plays a key role in cancer cell growth, survival, and angiogenesis. Phosphoinositide-dependent protein kinase-1 (PDK1) acts at a focal point in this pathway immediately downstream of phosphoinositide 3-kinase and PTEN, where it phosphorylates numerous AGC kinases. The PDK1 kinase domain has at least three ligand-binding sites: the ATP-binding pocket, the peptide substrate-binding site, and a groove in the N-terminal lobe that binds the C-terminal hydrophobic motif of its kinase substrates. Based on the unique PDK1 substrate recognition system, ultrahigh throughput TR-FRET and Alphascreen® screening assays were developed using a biotinylated version of the PDK1-tide substrate containing the activation loop of AKT fused to a pseudo-activated hydrophobic motif peptide. Using full-length PDK1, K_m values were determined as 5.6 μm for ATP and 40 nm for the fusion peptide, revealing 50-fold higher affinity compared with the classical AKT(Thr-308)-tide. Kinetic and biophysical studies confirmed the PDK1 catalytic mechanism as a rapid equilibrium random bireactant reaction. Following an ultrahigh throughput screen of a large library, 2,000 compounds were selected from the reconfirmed hits by computational analysis with a focus on novel scaffolds. ATP-competitive hits were deconvoluted by dose-response studies at 1× and 10× K_m concentrations of ATP, and specificity of binding was assessed in thermal shift assay. Inhibition studies using fusion PDK1-tide1 substrate versus AKT(Thr-308)-tide and kinase selectivity profiling revealed a novel selective alkaloid scaffold that evidently binds to the PDK1-interacting fragment pocket. Molecular modeling suggests a structural paradigm for the design of inhibitory versus activating allosteric ligands of PDK1.     

Antimetabolites: Design,synthesis, and cytotoxic evaluation of novel dihydropyridine thioglycosides and pyridine thioglycosides

A convenient synthesis of a novel series of dihydropyridine and pyridine thioglycosides was developed. The evaluation of anti-proliferative activity against HepG-2 cell lines (liver carcinoma cell lines) shows that most of the compounds have antitumor activity, especially 5b, 5f, 5j, 5n, 7b, 7f, 7j, 7n, 8b, 8f, and 8j. The results of molecular docking reveal that these compounds have high binding affinity by hydrogen bond formation with the binding pocket of thymidylate synthase dihydrofolate reductase (TS-DHFR).     

Azole-based inhibitors of AKT/PKB for the treatment of cancer

Through a combination of screening and structure-based rational design, we have discovered a series of N¹-(5-(heterocyclyl)-thiazol-2-yl)-3-(4-trifluoromethylphenyl)-1,2-propanediamines that were developed into potent ATP competitive inhibitors of AKT. Studies of linker strand-binding adenine isosteres identified SAR trends in potency and selectivity that were consistent with binding interactions observed in structures of the inhibitors bound to AKT1 and to the counter-screening target PKA. One compound was shown to have acceptable pharmacokinetic properties and to be a potent inhibitor of AKT signaling and of in vivo xenograft tumor growth in a preclinical model of glioblastoma.     

Structure-based drug design of 1,3,6-trisubstituted 1,4-diazepan-7-ones as selective human kallikrein 7 inhibitors

A novel series of 1,3,6-trisubstituted 1,4-diazepan-7-ones were investigated as human kallikrein 7 (KLK7, stratum corneum chymotryptic enzyme) inhibitors. Based on the X-ray co-crystal structure of compound 1 bound to human KLK7, the derivatives of this scaffold were designed, synthesized, and evaluated. Through structure-activity relationship studies focused on the side chain located in the prime site region of the enzyme, representative compounds 15, 33a, and 35a were identified as highly potent and selective inhibitors of human KLK7.