Histone deacetylase (HDAC) inhibitors are known to induce multiple epigenetic modifications affecting signaling networks and act synergistically with phosphatidylinositol 3-kinase (PI3K) inhibitors for the treatment of cancer. Herein we present a novel design approach for cancer drug development by incorporating HDAC inhibitory functionality into a PI3K inhibitor pharmacophore to construct dual-acting inhibitors. The designed compounds were synthesized and showed inhibitory activities against PI3K and HDAC. The representative dual PI3K/HDAC inhibitors, compounds 12a–j, showed potent antiproliferative activities against K562 and Hut78 in cellular assays. This work may lay the foundation for developing novel dual PI3K/HDAC inhibitors as potential anticancer therapeutics.

Covalent drugs have attracted increasing attention in recent years due to good inhibitory activity and selectivity. Targeting noncatalytic cysteines with irreversible inhibitors is a powerful approach for enhancing pharmacological potency and selectivity because cysteines can form covalent bonds with inhibitors through their nucleophilic thiol groups. However, most human kinases have multiple noncatalytic cysteines within the active site; to accurately predict which cysteine is most likely to form covalent bonds is of great importance but remains a challenge when designing irreversible inhibitors. In this work, FTMap was first applied to check its ability in predicting covalent binding site defined as the region where covalent bonds are formed between cysteines and irreversible inhibitors. Results show that it has excellent performance in detecting the hot spots within the binding pocket, and its hydrogen bond interaction frequency analysis could give us some interesting instructions for identification of covalent binding cysteines. Furthermore, we proposed a simple but useful covalent fragment probing approach and showed that it successfully predicted the covalent binding site of seven targets. By adopting a distance-based method, we observed that the closer the nucleophiles of covalent warheads are to the thiol group of a cysteine, the higher the possibility that a cysteine is prone to form a covalent bond. We believe that the combination of FTMap and our distance-based covalent fragment probing method can become a useful tool in detecting the covalent binding site of these targets.Keywords: covalent binding sites; cysteines; fragment probing; irreversible inhibitors;

The bromodomain protein module and histone deacetylase (HDAC), which recognize and remove acetylated lysine, respectively, have emerged as important epigenetic therapeutic targets in cancer treatments. Herein we presented a novel design approach for cancer drug development by combination of bromodomain and HDAC inhibitory activity in one molecule. The designed compounds were synthesized which showed inhibitory activity against bromodomain 4 and HDAC1. The representative dual bromodomain/HDAC inhibitors, compound 11 and 12, showed potent antiproliferative activities against human leukaemia cell line K562 and MV4-11 in cellular assays. This work may lay the foundation for developing dual bromodomain/HDAC inhibitors as potential anticancer therapeutics.

Intermolecular S···O interactions are very common and are important in biological systems, but until recently, the presence of these contacts in protein–ligand systems largely depended on serendipitous discovery instead of rational design. Here we provide insight into the phenomenon of intermolecular S···O contacts by focusing on three sulfur-containing aromatic rings. Quantum mechanics is employed to characterize the strength and directionality of the S···O interactions and to determine their energy dependence on their geometric parameters. Protein Data Bank mining is performed to systematically determine the occurrence and geometry of intermolecular S···O interactions, and several representative examples are discussed. Three typical cases are investigated using a combined quantum mechanics/molecular mechanics approach to demonstrate the potential of these interactions in improving binding affinities and physiochemical properties. Overall, our work elucidates the structures and energy features of intermolecular S···O interactions and addresses their use in molecular design.

•Compound 6a was designed by combining the scaffolds of UI-125 and Sorafenib.•SAR was explored from the 2-amino-3-purinylpyridine derivatives.•All compounds have B-RafV600E inhibitory activities at nanomolar ranges.•Some compounds showed potent antiproliferation against melanoma A375 cell line.•Compound 20g displayed promising antitumor efficacy in vivo.By combining the scaffolds of UI-125 and Sorafenib, a series of bis-aryl ureas and amides based on 2-amino-3-purinylpyridine moiety were designed and synthesized as novel DFG-out B-RafV600E inhibitors. Among them, 20c–e, 20g and 21h displayed potent antiproliferative activities against melanoma A375 (B-RafV600E) cell lines with IC50 values of 3.190, 2.276, 1.856, 1.632 μM and 1.839 μM, respectively, comparable with the positive control Vemurafenib (IC50 = 3.32 μM). Selected compounds were tested for the ERK inhibition in human melanoma A375 (B-RafV600E) and SK-MEL-2 (B-RafWT) cell lines by Western blot. The results revealed that our compounds inhibited the proliferation of melanoma A375 cells (B-RafV600E) through ERK pathway, without paradoxical activation of ERK in melanoma SK-MEL-2 cells (B-RafWT). Eventually, 20g and 21h were selected to confirm their inhibitory effects on tumor growth in A375 xenograft models in mice. Compound 20g exhibited equivalent antitumor efficacy in vivo (T/C = 44.37%), compared to Sorafenib (T/C = 37.35%), by 23-day repetitive administration of a single dose of 50 mg/kg without significant body weight loss.A series of 2-amino-3-purinylpyridine derivatives were synthesized and screened against B-RafV600E and A375 cell line, leading to the discovery of compound 20g, which exhibited potent A375 antiproliferative effects in vivo.

The bromodomain is a key protein–protein interaction module that specifically reads the acetylation marks of histones in epigenetic regulation. Currently, lots of inhibitors targeting the bromodomain have been reported as therapeutic agents. To better understand the interaction mechanism of bromodomain inhibitors, 20 diverse bromodomain inhibitors were studied using a combination of computational methods, including molecular docking, interaction fingerprinting, molecular dynamics simulation and binding free energy calculation. As a result, interactions important for the activity were critically analyzed, and the energy contribution in terms of individual residues was explored. These integrated results provided insights into two hot spots in the active site of the bromodomain, where the hydrophobic hot spot formed by Trp81, Val87, Leu92 and Ile146 played a central role in the interaction, and the hydrogen-bond hot spot mediated by Asn140 exhibited a moderate contribution to the binding affinity of the bromodomain inhibitors. This interaction mechanism study may facilitate the rational design of novel small-molecule bromodomain inhibitors.

The discovery of novel scaffolds against a specific target has long been one of the most significant but challengeable goals in discovering lead compounds. A scaffold that binds in important regions of the active pocket is more favorable as a starting point because scaffolds generally possess greater optimization possibilities. However, due to the lack of sufficient chemical space diversity of the databases and the ineffectiveness of the screening methods, it still remains a great challenge to discover novel active scaffolds. Since the strengths and weaknesses of both fragment-based drug design and traditional virtual screening (VS), we proposed a fragment VS concept based on Bayesian categorization for the discovery of novel scaffolds. This work investigated the proposal through an application on VEGFR-2 target. Firstly, scaffold and structural diversity of chemical space for 10 compound databases were explicitly evaluated. Simultaneously, a robust Bayesian classification model was constructed for screening not only compound databases but also their corresponding fragment databases. Although analysis of the scaffold diversity demonstrated a very unevenly distribution of scaffolds over molecules, results showed that our Bayesian model behaved better in screening fragments than molecules. Through a literature retrospective research, several generated fragments with relatively high Bayesian scores indeed exhibit VEGFR-2 biological activity, which strongly proved the effectiveness of fragment VS based on Bayesian categorization models. This investigation of Bayesian-based fragment VS can further emphasize the necessity for enrichment of compound databases employed in lead discovery by amplifying the diversity of databases with novel structures.

c-Met has been considered as an attractive target for developing antitumor agents. The highly selective c-Met inhibitors provide invaluable opportunities for the combination with other therapies safely to achieve the optimal efficacy. In this work, a series of triazolopyrazine c-Met inhibitors with exquisitely selectivity were investigated using a combination of molecular docking, three-dimensional quantitative structure–activity relationship (3D-QSAR), and molecular dynamics simulation. Comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) models were developed to reveal the structural determinants for c-Met inhibition. Both models were validated to have high reliability and predictability, and contour map analysis suggested feature requirements for different substituents on the scaffold. It is worth noting that an important hydrogen bond rich region was identified in the unique narrow channel, which is distinct from other kinases. Molecular dynamics simulations and binding free energy calculations provided further support that suitable groups in this hydrogen bond rich region made great contributions to the binding of ligands. Moreover, hydrogen bonds with residues of the narrow channel were also indicated to be essential to improve the activity and selectivity. This study will facilitate the discovery and optimization of novel c-Met inhibitors with higher activity and selectivity.

A series of novel indole derivatives were designed and synthesized and their inhibitory activity against B-Raf and HepG2 cell were also described. Among them, compounds 7a and 7b exhibited excellent potency, which showed the potential for further research as lead compounds.A series of novel indole derivatives were designed and synthesized and their inhibitory activity against B-Raf and HepG2 cells were also described. Among them, compounds 7a and 7b exhibited excellent potency, which showed the potential for further research as lead compounds.

In recent years, various virtual screening (VS) tools have been developed, and many successful screening campaigns have been showcased. However, whether by conventional molecular docking or pharmacophore screening, the selection of virtual hits is based on the ranking of compounds by scoring functions or fit values, which remains the bottleneck of VS due to insufficient accuracy. As the limitations of individual methods persist, a comprehensive comparison and integration of different methods may provide insights into selecting suitable methods for VS. Here, we evaluated the performance of molecular docking, fingerprint-based 2D similarity and multicomplex pharmacophore in an individual and a combined manner, through a retrospective VS study on VEGFR-2 inhibitors. An integrated two-layer workflow was developed and validated through VS of VEGFR-2 inhibitors against the DUD-E database, which demonstrated improved VS performance through a ligand-based method ECFP_4, followed by molecular docking, and then a strict multicomplex pharmacophore. Through a retrospective comparison with six published papers, this integrated approach outperformed 43 out of 45 methods, indicating a great effectiveness. This kind of integrated VS approach can be extended to other targets for the screening and discovery of inhibitors.

Fragment-based drug design has emerged as an important methodology for lead discovery and drug design. Different with other studies focused on fragment library design and active fragment identification, a fragment-based strategy was developed in combination with three-dimensional quantitative structure–activity relationship (3D-QSAR) for structural optimization in this study. Based on a validated scaffold or fragment hit, a series of structural optimization was conducted to convert it to lead compounds, including 3D-QSAR modelling, active site analysis, fragment-based structural optimization and evaluation of new molecules. 3D-QSAR models and active site analysis provided sufficient information for confirming the SAR and pharmacophoric features for fragments. This strategy was evaluated through the structural optimization on a c-Met inhibitor scaffold 5H-benzo[4,5]cyclohepta[1,2-b]pyridin-5-one, which resulted in an c-Met inhibitor with high inhibitory activity. Our study suggested the effectiveness of this fragment-based strategy and the druggability of our newly explored active region. The reliability of this strategy indicated it could also be applied to facilitate lead optimization of other targets.

β-Carboline family of compounds is a large group of alkaloids widely distributed in nature and exhibits broad-spectrum anti-tumor activities. We designed and synthesized two series of novel 1-carboxamide- and 6-sulfonamide-substituted β-carboline derivatives 7a–p and 12a–b, and their wild type B-Raf kinase inhibitory activities were described. Most compounds showed moderate to excellent inhibitory activities. Among them, 1-carboxamide-6-(N-(3-(dimethylamino)propyl)-sulfamoyl)-β-carboline, 7e exhibited potent activity (IC50 = 1.62 μM), showing the potential for further investigation as a lead compound.Two series of novel 1-carboxamide and 6-sulfonamide substituted β-carboline derivatives 7a–o, 12a–b were designed and synthesized and their wild type B-Raf kinase inhibitory activities were described. Compound 1-carboxamide-6-(N-(3-(dimethylamino)propyl)sulfamoyl)-β-carboline, 7e exhibited potent activity (IC50 = 1.62 μM), showing the potential for further investigation as a lead compound.

Vascular endothelial growth factor (VEGF) and its receptor tyrosine kinase VEGFR-2 or kinase insert domain receptor (KDR) have been identified as promising targets for novel anticancer agents. To achieve new potent inhibitors of KDR, we conducted molecular fragment replacement (MFR) studies for the understanding of 3D-QSAR modeling and the docking investigation of arylphthalazines and 2-((1H-Azol-1-yl)methyl)-N-arylbenzamides-based KDR inhibitors. Two favorable 3D-QSAR models (CoMFA with q2, 0.671; r2, 0.969; CoMSIA with q2, 0.608; r2, 0.936) have been developed to predict the biological activity of new compounds. The new molecular database generated by MFR was virtually screened using Glide (docking) and further evaluated with CoMFA prediction, protein–ligand interaction fingerprint (PLIF) and ADMET analysis. 44 N-(pyridin-4-ylmethyl)aniline derivatives as novel potential KDR inhibitors were finally obtained. In this paper, the work flow developed could be applied to de novo drug design and virtual screening potential KDR inhibitors, and use hit compounds to further optimize and design new potential KDR inhibitors.

A three dimensional (3D) chemical feature based pharmacophore model was developed for selective histone deacetylase 1 (HDAC1) inhibitors, which provides an efficient way to discuss the isoform selectivity of HDAC inhibitors. In contrast to the classical pan-HDAC pharmacophore, two hydrophobic features (HY and HYAr2) were found in the chemical feature based pharmacophore model, which might be responsible for the selectivity of HDAC1 inhibitions. Molecular docking also highlighted the two hydrophobic features, which are located in the internal cavity adjacent to the active site. The results contribute to our understanding of the molecular mechanism underlying the selectivity of HDAC1 inhibitors and suggest a possible target region to design novel selective HDAC1 inhibitors.A selective pharmacophore model was developed based on a series of selective HDAC1 inhibitors. Two hydrophobic features (HY and HYAr2) were responsible for the selectivity of HDAC1 inhibitions.

A three-dimensional pharmacophore model was developed based on 25 currently available Raf-1 kinase inhibitors. The best pharmacophore hypothesis (Hypo1), consisting of four chemical features (one hydrogen-bond acceptor, one hydrogen-bond donor, and two hydrophobic groups), has a correlation coefficient of 0.972. The results of our study provide a valuable tool in designing new leads with desired biological activity by virtual screening.The best pharmacophore hypothesis (Hypo1), consisting of four chemical features, and the high active compound (compound 6) can be well mapped onto the Hypo1 model.