Identification of BRAF inhibitors through in silico screening

Computational investigation of imidazo[2,1-b]oxazole derivatives as potential mutant BRAF kinase inhibitors: 3D-QSAR, molecular docking, molecular dynamics simulation, and ADMETox studies

BRAF inhibitors are known to be an effective therapeutic target for treating melanoma and other types of cancer. Using 3D-QSAR, molecular docking, and MD simulations, this study evaluated various imidazo[2,1-b]oxazole derivatives that function as mutant BRAF kinase inhibitors. Comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) were used to create the 3D-QSAR models. CoMSIA/SEHA model has solid predictive power across several models (Q2 = 0.578; R2 = 0.828; R2pred = 0.74) and is the best model according to the numerous field models generated. The created model's predictive power was evaluated through external validation using a test set. CoMSIA/SEHA contour maps collect information that can be used to identify critical regions with solid anticancer activity. We developed four inhibitors with high predicted activity due to these observations. ADMET prediction was used to assess the toxicity of the proposed imidazo[2,1-b]oxazole compounds. The predictive molecules (T1-T4) demonstrated good ADMET properties, excluding the toxic active compounds 11r from the database. Molecular docking was also used to determine the patterns and modes of interactions between imidazo[2,1-b]oxazole ligands and receptors, which revealed that the proposed imidazo[2,1-b]oxazole scaffold was stable in the receptor's active site (PDB code: 4G9C). The suggested compounds (T1-T4) were subjected to molecular dynamics simulations lasting 100 ns to determine their binding free energies. The results showed that T2 had a more favorable binding free energy (-149.552 kJ/mol) than T1 (-112.556 kJ/mol), T3 (-115.503 kJ/mol), and T4 (-102.553 kJ/mol). The results suggest that the imidazo[2,1-b]oxazole compounds investigated in this study have potential as inhibitors of BRAF kinase and could be further developed as anticancer drugs. Highlights22 imidazo[2,1-b]oxazole compounds were subjected to research on three-dimensional quantitative conformational relationships.Using contour maps from 3D-QSAR models as a guide was used to figure out the areas and strategies for structural optimization.Combined molecular docking, molecular dynamics simulations, and binding free energy calculations to verify the inhibitor activity of the proposed 22 imidazo[2,1-b]oxazole compounds.Four potential B-RAF Kinase inhibitors were discovered, providing theoretical clues for developing a highly anticancer agent.Communicated by Ramaswamy H. Sarma.

KinomeMETA: meta-learning enhanced kinome-wide polypharmacology profiling

Kinase inhibitors are crucial in cancer treatment, but drug resistance and side effects hinder the development of effective drugs. To address these challenges, it is essential to analyze the polypharmacology of kinase inhibitor and identify compound with high selectivity profile. This study presents KinomeMETA, a framework for profiling the activity of small molecule kinase inhibitors across a panel of 661 kinases. By training a meta-learner based on a graph neural network and fine-tuning it to create kinase-specific learners, KinomeMETA outperforms benchmark multi-task models and other kinase profiling models. It provides higher accuracy for understudied kinases with limited known data and broader coverage of kinase types, including important mutant kinases. Case studies on the discovery of new scaffold inhibitors for membrane-associated tyrosine- and threonine-specific cdc2-inhibitory kinase and selective inhibitors for fibroblast growth factor receptors demonstrate the role of KinomeMETA in virtual screening and kinome-wide activity profiling. Overall, KinomeMETA has the potential to accelerate kinase drug discovery by more effectively exploring the kinase polypharmacology landscape.

A novel phosphorylation site involved in dissociating RAF kinase from the scaffolding protein 14-3-3 and disrupting RAF dimerization

Rapidly accelerated fibrosarcoma (ARAF, BRAF, CRAF) kinase is central to the MAPK pathway (RAS-RAF-MEK-ERK). Inactive RAF kinase is believed to be monomeric, autoinhibited, and cytosolic, while activated RAF is recruited to the membrane via RAS-GTP, leading to the relief of autoinhibition, phosphorylation of key regulatory sites, and dimerization of RAF protomers. Although it is well known that active and inactive BRAF have differential phosphorylation sites that play a crucial role in regulating BRAF, key details are still missing. In this study, we report the characterization of a novel phosphorylation site, BRAFS732 (equivalent in CRAFS624), located in proximity to the C-terminus binding motif for the 14-3-3 scaffolding protein. At the C terminus, 14-3-3 binds to BRAFpS729 (CRAFpS621) and enhances RAF dimerization. We conducted mutational analysis of BRAFS732A/E and CRAFS624A/E and revealed that the phosphomimetic S→E mutant decreases 14-3-3 association and RAF dimerization. In normal cell signaling, dimerized RAF phosphorylates MEK1/2, which is observed in the phospho-deficient S→A mutant. Our results suggest that phosphorylation and dephosphorylation of this site fine-tune the association of 14-3-3 and RAF dimerization, ultimately impacting MEK phosphorylation. We further characterized the BRAF homodimer and BRAF:CRAF heterodimer and identified a correlation between phosphorylation of this site with drug sensitivity. Our work reveals a novel negative regulatory role for phosphorylation of BRAFS732 and CRAFS624 in decreasing 14-3-3 association, dimerization, and MEK phosphorylation. These findings provide insight into the regulation of the MAPK pathway and may have implications for cancers driven by mutations in the pathway.

Molecular Modeling Unveils the Effective Interaction of B-RAF Inhibitors with Rare B-RAF Insertion Variants

The Food and Drug Administration (FDA) has approved MAPK inhibitors as a treatment for melanoma patients carrying a mutation in codon V600 of the BRAF gene exclusively. However, BRAF mutations outside the V600 codon may occur in a small percentage of melanomas. Although these rare variants may cause B-RAF activation, their predictive response to B-RAF inhibitor treatments is still poorly understood. We exploited an integrated approach for mutation detection, tumor evolution tracking, and assessment of response to treatment in a metastatic melanoma patient carrying the rare p.T599dup B-RAF mutation. He was addressed to Dabrafenib/Trametinib targeted therapy, showing an initial dramatic response. In parallel, in-silico ligand-based homology modeling was set up and performed on this and an additional B-RAF rare variant (p.A598_T599insV) to unveil and justify the success of the B-RAF inhibitory activity of Dabrafenib, showing that it could adeptly bind both these variants in a similar manner to how it binds and inhibits the V600E mutant. These findings open up the possibility of broadening the spectrum of BRAF inhibitor-sensitive mutations beyond mutations at codon V600, suggesting that B-RAF V600 WT melanomas should undergo more specific investigations before ruling out the possibility of targeted therapy.

Computer-Aided Drug Design Boosts RAS Inhibitor Discovery

The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design.

Target-Based Small Molecule Drug Discovery for Colorectal Cancer: A Review of Molecular Pathways and In Silico Studies

Colorectal cancer is one of the most prevalent cancer types. Although there have been breakthroughs in its treatments, a better understanding of the molecular mechanisms and genetic involvement in colorectal cancer will have a substantial role in producing novel and targeted treatments with better safety profiles. In this review, the main molecular pathways and driver genes that are responsible for initiating and propagating the cascade of signaling molecules reaching carcinoma and the aggressive metastatic stages of colorectal cancer were presented. Protein kinases involved in colorectal cancer, as much as other cancers, have seen much focus and committed efforts due to their crucial role in subsidizing, inhibiting, or changing the disease course. Moreover, notable improvements in colorectal cancer treatments with in silico studies and the enhanced selectivity on specific macromolecular targets were discussed. Besides, the selective multi-target agents have been made easier by employing in silico methods in molecular de novo synthesis or target identification and drug repurposing.

Comparative differential cytotoxicity of clinically used SERMs in human cancer lines of different origin and its predictive molecular docking studies of key target genes involved in cancer progression and treatment responses

SERMS like Tamoxifene, 5-hydroxy tamoxifene, raloxifene and endoxifene has been used for the treatment of hormonal imbalances and dependent cancers owing to their action via Estrogen receptors as in the treatment of estrogen sensitive breast cancers. Due to the adverse side effects, modifications and development of the existing or newer SERMS has always been of immense interest. Ormeloxifene, a SERM molecule manufactured by HLL Lifecare Ltd, India as birth control under the trade names Saheli, Novex, and Novex-DS which is also investigated against mastalgia, fibro-adenoma and abnormal uterine bleeding. Anti-cancer effects have been reported in estrogen dependent and independent cancers which shows its wide scope to be implemented in cancer therapy. Current investigation is a comprehensive effort to find the cytotoxic potential of Ormeloxifene in comparison with clinically used four SERMS in twenty six cancer cell lines of different origin using Adriamycin as positive control. Also the computational studies pertaining to selected target/ligand with respect to tumor progression, development, treatment responses and apoptosis. The studies proved effective cytotoxicity of Ormeloxifene on cancer cell lines with lower TGI, GI50 and LC50 values which are significantly comparable. Also the in silico studies proved that the docking score of the compound suggests the interaction of the compound which could tightly regulate key target genes controlling cancer like ER, EGFR kinase, EGFR-cSRC, HDAC-2, PARP-1 and BRAF. This study brings out the superior efficacy of Ormeloxifene compared to other SERMS with proven safety profile to be repositioned as an anti-cancer drug to treat diverse cancer types.

Computational evaluation of potent 2-(1H-imidazol-2-yl) pyridine derivatives as potential V600E-BRAF inhibitors

Background

V600E-BRAF is a major protein target involved in various types of human cancers. However, the acquired resistance of the V600E-BRAF kinase to the vemurafenib and the side effects of other identified drugs initiate the search for efficient inhibitors. In the current paper, virtual docking screening combined with drug likeness and ADMET properties predictions were jointly applied to evaluate potent 2-(1H-imidazol-2-yl) pyridines as V600E-BRAF kinase inhibitors.

Results

Most of the studied compounds showed better docking scores and favorable interactions with theiV600E-BRAF target. Among the screened compounds, the two most potent (14 and 30) with good rerank scores (−124.079 and − 122.290) emerged as the most effective, and potent V600E-BRAF kinase inhibitors which performed better than vemurafenib (−116.174), an approved V600E-BRAF kinase inhibitor. Thus, the docking studies exhibited that these compounds have shown competing inhibition of V600E-BRAF kinase with vemurafenib at the active site and revealed better pharmacological properties based on Lipinski’s and Veber’s drug-likeness rules for oral bioavailability and ADMET properties.

Conclusion

The docking result, drug-likeness rules, and ADMET parameters identified compounds (14 and 30) as the best hits against V600E-BRAF kinase with better pharmacological properties. This suggests that these compounds may be developed as potent V600E-BRAF inhibitors.

In silico evaluation of some 4-(quinolin-2-yl)pyrimidin-2-amine derivatives as potent V600E-BRAF inhibitors with pharmacokinetics ADMET and drug-likeness predictions

Background

The resistance of V600E-BRAF to the vemurafenib and the side effects of the identified inhibitors trigger the research for a novel and more potent anti-melanoma agents. In this study, virtual docking screening along with pharmacokinetics ADMET and drug-likeness predictions were combined to evaluate some 4-(quinolin-2-yl)pyrimidin-2-amine derivatives as potent V600E-BRAF inhibitors.

Results

Some of the selected compounds exhibited better binding scores and favorable interaction with the V600E-BRAF enzyme. Out of the screened compounds, two most potent (5 and 9) having good Rerank scores (− 128.011 and − 126.258) emerged as effective and potent V600E-BRAF inhibitors that outperformed the FDA-approved V600E-BRAF inhibitor (vemurafenib, − 118.607). Thus, the molecular docking studies revealed that the studied compounds showed competing for inhibition of V600E-BRAF with vemurafenib at the binding site and possessed better pharmacological parameters based on the drug-likeness rules filters for the oral bioavailability, and ADMET risk parameters.

Conclusion

The docking analysis, drug-likeness rules filters, and ADMET study identified compounds (5 and 9) as the best hits against V600E-BRAF kinase with enhanced pharmacological properties. This recommends that these compounds may be developed as potent anti-melanoma agents.

Synthesis and biological evaluation of novel xanthine derivatives as potential apoptotic antitumor agents

A series of novel xanthine/NO donor hybrids containing 1,3,8-trisubstituted or 1,8-disubstituted xanthine derivatives were designed and synthesized. The synthesized compounds were tested in a cell viability assay using human mammary gland epithelial cell line (MCF-10A) where all the compounds exhibited no cytotoxic effects and more than 90% cell viability at a concentration of 50 μM. The oxime containing compounds 7a-b and 17-24 were more active as antiproliferative agents than their non-oxime congeners 6a-b and 9-16. Hydroxyimino-phenethyl scaffold compounds 17-24 were more active than the hydroxyimino-ethyl phenyl acetamide 7a-b derivatives. Compounds 18-20 and 22-24 exhibited inhibition of EGFR with IC₅₀ ranging from 0.32 to 2.88 μM. Compounds 18-20 and 22-24 increased the level of active caspase 3 by 4-8 folds, compared to the control cells in Panc-1 cell lines compared to doxorubicin as a reference drug. Compounds 18, 22 and 23 were the most caspase-3 inducers. Compounds 22 and 23 increased the levels of caspase-8 and 9 indicating activation of both intrinsic and extrinsic pathways and showed potent induction of Bax, down-regulation of Bcl-2 protein levels and over-expression of cytochrome c levels in Panc-1 human pancreas cancer cells. Compound 23 exhibited mainly cell cycle arrest at the Pre-G1 and G2/M phases in the cell cycle analysis of Panc-1 cell line. The drug likeness profiles of compounds 18-20 and 22-24 were predicted to have good to excellent drug likeness profiles specially compounds 18-20 and 23. Finally molecular docking study was performed at the EGFR active site to suggest thier possible binding mode. The hydroxyimino-phenethyl scaffold compounds 17-24 represent an interesting starting point to optimize their pharmacokinetics and pharmacodynamics profiles.

Mechanism of BRAF Activation through Biochemical Characterization of the Recombinant Full-Length Protein

BRAF kinase plays an important role in mitogen-activated protein kinase (MAPK) signaling and harbors activating mutations in about half of melanomas and in a smaller percentage in many other cancers. Despite its importance, few in vitro studies have been performed to characterize the biochemical properties of full-length BRAF. Herein, a strategy to generate an active, intact form of BRAF protein suitable for in vitro enzyme kinetics is described. It is shown that purified, intact BRAF protein autophosphorylates the kinase activation loop and this can be enhanced by binding the MEK protein substrate through an allosteric mechanism. These studies provide in vitro evidence that BRAF selectively binds to active RAS and that the BRAF/CRAF heterodimer is the most active form, relative to their respective homodimers. Full-length BRAF analysis with small-molecule BRAF inhibitors shows that two drugs, dabrafenib and vemurafenib, can modestly enhance kinase activity of BRAF at low concentration. Taken together, this characterization of intact BRAF contributes to a framework for understanding its role in cell signaling.

Identification of Specific Inhibitors of Trypanosoma cruzi Malic Enzyme Isoforms by Target-Based HTS

Trypanosoma cruzi is the causative agent of Chagas disease. The lack of an efficient and safe treatment supports the research into novel metabolic targets, with the malic enzyme (ME) representing one such potential candidate. T. cruzi expresses a cytosolic (TcMEc) and a mitochondrial (TcMEm) ME isoform, with these activities functioning to generate NADPH, a key source of reducing equivalents that drives a range of anabolic and protective processes. To identify specific inhibitors that target TcMEs, two independent high-throughput screening strategies using a diversity library containing 30,000 compounds were employed. IC₅₀ values of 262 molecules were determined for both TcMEs, as well as for three human ME isoforms, with the inhibitors clustered into six groups according to their chemical similarity. The most potent hits belonged to a sulfonamide group that specifically target TcMEc. Moreover, several selected inhibitors of both TcMEs showed a trypanocidal effect against the replicative forms of T. cruzi. The chemical diversity observed among those compounds that inhibit TcMEs activity emphasizes the druggability of these enzymes, with a sulfonamide-based subset of compounds readily able to block TcMEc function at a low nanomolar range.

Design, synthesis and biological evaluation of novel benzo-α-pyrone containing piperazine derivatives as potential BRAFV600E inhibitors

The increasingly acquired resistance to vemurafenib and side effects of known inhibitors motivate the search for new and more effective anti-melanoma drugs. In this Letter, virtual screening and scaffold growth were combined together to achieve new molecules as BRAF^V600E inhibitors. Along with docking simulation, a primary screen in vitro was performed to filter the modifications for the lead compound, which was then substituted, synthesized and evaluated for their inhibitory activity against BRAF^V600E and several melanoma cell lines. Out of the obtained compounds, derivative 3l was identified as a potent BRAF^V600E inhibitor and exerted an anticancer effect through BRAF^V600E inhibition. The following biological evaluation assays confirmed that 3l could induce cell apoptosis and marked DNA fragmentation. Furthermore, 3l could arrest the cell cycle at the G0/G1 phase in melanoma cells. The docking simulation displayed that 3l could tightly bind with the crystal structure of BRAF^V600E at the active site. Overall, the biological profile of 3l suggests that this compound may be developed as a potential anticancer agent.

Insight into molecular dynamics simulation of BRAF(V600E) and potent novel inhibitors for malignant melanoma

BRAF inhibitors have changed the standard therapeutic protocol for advanced or metastatic melanoma which harbored notorious BRAF(V600E) single mutation. However, drug resistance to BRAF inhibitors happens just like other cancer treatment. In this study, we constructed the ideal BRAF(V600E)-modeled structure through homology modeling and introduced the method of structure-based docking or virtual screening from the large compound database. Through certain methods of molecular dynamics simulation, we realized that BRAF(V600E) had quite prominent difference of molecular character or structural variation from the wild-type BRAF protein. It might confer the metamorphic character of advanced melanoma for the patients who harbored BRAF(V600E) mutation. By the methods of ligand-based quantitative structure-activity relationship and molecular dynamics simulation, we further recommend that aknadicine and 16beta-hydroxy-19s-vindolinine N-oxide from the traditional Chinese medicine are potent novel inhibitors for the management of malignant melanoma in the future.

Structural optimization and biological screening of a steroidal scaffold possessing cucurbitacin-like functionalities as B-Raf inhibitors

Inhibition of the mitogen-activated protein kinase (MAPK) pathway by targeting the commonly occurring mutated B-Raf in melanoma has become a practical method for the development of drugs and drug candidates. In order to expand upon the currently reported structural scaffolds used to target the MAPK pathway, molecular docking studies led to the installation an α,β-unsaturated ketone side chain, related to the cucurbitacin class of natural products, on to an estrone core via an aldol condensation reaction, along with installation of the Δ(9,11) olefin to assemble what has been defined as a pseudo-cis configuration at the B/C ring juncture. Combination of these cucurbitacin-like features resulted in a compound with an enhanced biological profile against the A-375 mutant B-Raf cell line, in regards to their cytotoxicity and inhibitory activity toward phosphorylated extracellular-signal-regulated kinase (ERK).

Insight into the structural features of pyrazolopyrimidine- and pyrazolopyridine-based B-Raf(V600E) kinase inhibitors by computational explorations

Presently, both ligand-based and receptor-based 3D-QSAR modelings were performed on 107 pyrazolopyrimidine- and pyrazolopyridine-based inhibitors of B-Raf(V600E) kinase. The optimal model is successful to predict the inhibitors' activity with Q(2) of 0.504, R(2) ncv of 0.960, and R(2) pred of 0.872. Besides, the 3D contour maps explain well the structural requirements of the interaction between the ligand and the receptor. Furthermore, molecular docking and MD were also carried out to study the binding mode. Our findings are the following: (i) Bulky substituents at position 3, 10 and ring D improve the inhibitory activity, but impair the activity at position 5, 11, and 19. (ii) Electropositive groups at position 10, 13 and 20 and electronegative groups at position 2 increase the biological activity. (iii) Hydrophobic substituents at ring C are beneficial to improve the biological activity, while hydrophilic substituents at position 11 and ring D are good for the activity. (4) This scaffold of inhibitors may bind to the B-Raf kinase with an 'L' conformation and belong to type III binding mode, which is fixed by hydrophobic interaction and hydrogen bonds with residues from hinge region and DFG motif. These results may be a guidance to develop new B-Raf(V600E) kinase inhibitors.

Design, synthesis and biological evaluation of (1,3-diphenyl-1H-pyrazol-4-yl) methyl benzoate derivatives as potential BRAFV600E inhibitors

A series of (1,3-diphenyl-1H-pyrazol-4-yl) methyl benzoate derivatives (6a–10d) were designed and synthesized and evaluated as BRAF^V600 inhibitors. Among them, compound 10a showed the most potent inhibitory activity against A375, WM266.4 and BRAF^V600Ein vitro with IC₅₀ values of 1.36 μM, 0.94 μM and 0.11 μM, respectively.

Identification and synthesis of N-(thiophen-2-yl) benzamide derivatives as BRAF(V600E) inhibitors

The V600E BRAF kinase mutation, which activates the downstream MAPK signaling pathway, commonly occurs in about 8% of all human malignancies and about 50% of all melanomas. In this study, we employed virtual screening and chemical synthesis to identify a series of N-(thiophen-2-yl) benzamide derivatives as potent BRAF(V600E) inhibitors. Structure-activity relationship studies of these derivatives revealed that compounds b40 and b47 are the two most potent BRAF(V600E) inhibitors in this series.

Development of a novel class of B-Raf(V600E)-selective inhibitors through virtual screening and hierarchical hit optimization

Oncogenic mutations in critical nodes of cellular signaling pathways have been associated with tumorigenesis and progression. The B-Raf protein kinase, a key hub in the canonical MAPK signaling cascade, is mutated in a broad range of human cancers and especially in malignant melanoma. The most prevalent B-Raf(V600E) mutant exhibits elevated kinase activity and results in constitutive activation of the MAPK pathway, thus making it a promising drug target for cancer therapy. Herein, we describe the development of novel B-Raf(V600E) selective inhibitors via multi-step virtual screening and hierarchical hit optimization. Nine hit compounds with low micromolar IC(50) values were identified as B-Raf(V600E) inhibitors through virtual screening. Subsequent scaffold-based analogue searching and medicinal chemistry efforts significantly improved both the inhibitor potency and oncogene selectivity. In particular, compounds 22f and 22q possess nanomolar IC(50) values with selectivity for B-Raf(V600E)in vitro and exclusive cytotoxicity against B-Raf(V600E) harboring cancer cells.

Pharmacophore-based virtual screening and biological evaluation of small molecule inhibitors for protein arginine methylation

Protein arginine methyltransferases (PRMTs) are proved to play vital roles in chromatin remodeling, RNA metabolism, and signal transduction. Aberrant regulation of PRMT activity is associated with various pathological states such as cancer and cardiovascular disorders. Development and application of small molecule PRMT inhibitors will provide new avenues for therapeutic discovery. The combination of pharmacophore-based virtual screening methods with radioactive methylation assays provided six hits identified as inhibitors against the predominant arginine methyltransferase PRMT1 within micromolar potency. Two potent compounds, A9 and A36, exhibited the inhibitory effect by directly targeting substrate H4 other than PRMT1 and displayed even higher inhibition activity than the well-known PRMT inhibitors AMI-1. A9 significantly inhibits proliferation of castrate-resistant prostate cancer cells. Together, A9 may be a potential inhibitor against advanced hormone-independent cancers, and the work will provide clues for the future development of specific compounds that block the interaction of PRMTs with their targets.

Identification of a novel family of BRAF(V600E) inhibitors

The BRAF oncoprotein is mutated in about half of malignant melanomas and other cancers, and a kinase activating single valine to glutamate substitution at residue 600 (BRAF(V600E)) accounts for over 90% of BRAF-mediated cancers. Several BRAF(V600E) inhibitors have been developed, although they harbor some liabilities, thus motivating the development of other BRAF(V600E) inhibitor options. We report here the use of an ELISA based high-throughput screen to identify a family of related quinolol/naphthol compounds that preferentially inhibit BRAF(V600E) over BRAF(WT) and other kinases. We also report the X-ray crystal structure of a BRAF/quinolol complex revealing the mode of inhibition, employ structure-based medicinal chemistry efforts to prepare naphthol analogues that inhibit BRAF(V600E) in vitro with IC(50) values in the 80-200 nM range under saturating ATP concentrations, and demonstrate that these compounds inhibit MAPK signaling in melanoma cells. Prospects for improving the potency and selectivity of these inhibitors are discussed.

Chemistry of polyhalogenated nitrobutadienes, 10: Synthesis of highly functionalized heterocycles with a rigid 6-amino-3-azabicyclo[3.1.0]hexane moiety

The nitropolychlorobutadienes 3, 4 are valuable building blocks for various amination and successive heterocyclization products. Nucleophilic substitution reactions of the partially protected, bioactive amines 1, 2 with either vinyl, imidoyl or carbonyl chlorides result in the formation of the enamines 11, 12, 13, 16, 25, the amidine 6, and the amides 20, 21, respectively. In the following, cyclization to the highly functionalized pyrazoles 27, 28, pyrimidine 26 and pyridopyrimidine 24 succeeded. Deprotection of 21, 12 and 28 proved to be only partially feasible.

Design, synthesis and biological evaluation of novel (E)-α-benzylsulfonyl chalcone derivatives as potential BRAF inhibitors

Activating mutations in the BRAF serine/threonine kinase are found in more than 70% of human melanomas, >90% of which are BRAF(V600E). It provides new therapeutic opportunities in malignant melanoma. In silico and in vitro screening of our compound collection has identified Hit 2 as BRAF(V600E) inhibitor. Based on its structure, a series of novel (E)-α-benzylsulfonyl chalcone derivatives (13-40) were designed and synthesized. Compound 38 exhibited the most potent inhibitory activity with an IC(50) value of 0.17 μM for BRAF(V600E) and GI(50) value of 0.52 μM for mutant BRAF-dependent cells. The results of cell based pERK activity and cellular selectivity suggested that those compounds could selectively inhibit proliferation of mutant BRAF-dependent melanoma cell line through inhibition of oncogenic BRAF.

Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds

N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.

Molecular dynamics simulation, free energy calculation and structure-based 3D-QSAR studies of B-RAF kinase inhibitors

(V600E)B-RAF kinase is the most frequent onco-genic protein kinase mutation in melanoma and is a promising target to treat malignant melanoma. In this work, a molecular modeling study combining QM-polarized ligand docking, molecular dynamics, free energy calculation, and three-dimensional quantitative structure-activity relationships (3D-QSAR) was performed on a series of pyridoimidazolone compounds as the inhibitors of (V600E)B-RAF kinase to understand the binding mode between the inhibitors and (V600E)B-RAF kinase and the structural requirement for the inhibiting activity. 3D-QSAR models, including CoMFA and CoMSIA, were developed from the conformations obtained by QM-polarized ligand docking strategy. The obtained models have a good predictive ability in both internal and external validation. Furthermore, molecular dynamics simulation and free energy calculations were employed to determine the detailed binding process and to compare the binding mode of the inhibitors with different activities. The binding free energies calculated by MM/PBSA gave a good correlation with the experimental biological activity. The decomposition of free energies by MM/GBSA indicates the van der Waals interaction is the major driving force for the interaction between the inhibitors and (V600E)B-RAF kinase. The hydrogen bond interactions between the inhibitors with Glu501 and Asp594 of the (V600E)B-RAF kinase help to stabilize the DFG-out conformation. The results from this study can provide some insights into the development of novel potent (V600E)B-RAF kinase inhibitors.

Development of drugs for Epstein-Barr virus using high-throughput in silico virtual screening

IMPORTANCE OF THE FIELD

Epstein-Barr virus (EBV) is a ubiquitous human herpesvirus that is causally associated with endemic forms of Burkitt's lymphoma, nasopharyngeal carcinoma and lymphoproliferative disease in immunosuppressed individuals. On a global scale, EBV infects > 90% of the adult population and is responsible for ∼ 1% of all human cancers. To date, there is no efficacious drug or therapy for the treatment of EBV infection and EBV-related diseases.

AREAS COVERED IN THIS REVIEW

In this review, we discuss the existing anti-EBV inhibitors and those under development. We discuss the value of different molecular targets, including EBV lytic DNA replication enzymes as well as proteins that are expressed exclusively during latent infection, such as EBV nuclear antigen 1 (EBNA-1) and latent membrane protein 1. As the atomic structure of the EBNA-1 DNA binding domain has been described, it is an attractive target for in silico methods of drug design and small molecule screening. We discuss the use of computational methods that can greatly facilitate the development of novel inhibitors and how in silico screening methods can be applied to target proteins with known structures, such as EBNA-1, to treat EBV infection and disease.

WHAT THE READER WILL GAIN

The reader is familiarized with the problems in targeting of EBV for inhibition by small molecules and how computational methods can greatly facilitate this process.

TAKE HOME MESSAGE

Despite the impressive efficacy of nucleoside analogs for the treatment of herpesvirus lytic infection, there remain few effective treatments for latent infections. As EBV latent infection persists within and contributes to the formation of EBV-associated cancers, targeting EBV latent proteins is an unmet medical need. High-throughput in silico screening can accelerate the process of drug discovery for novel and selective agents that inhibit EBV latent infection and associated disease.

Models for prediction of (V)600(E)BRAF and melanoma cells growth inhibitory activities of pyridoimidazolones

Targeted inhibition of activated BRAF mutation has emerged as a most promising and putative therapeutic approach for the anticancer drug development. In the present study, an in-silico approach using decision tree and moving average analysis has been applied to a data set comprising of 43 analogues of pyridoimidazolones for development of models for prediction of both (V)600(E)BRAF and melanoma cells (BRAF WM266.4) growth inhibitory activities. A decision tree was mainly employed for determining the importance of molecular descriptors (n=46). The value of majority of these descriptors for each analogue in the dataset was computed using E-Dragon software (version 1.0). The decision tree learned the information from the input data with an accuracy of 98% and correctly predicted the cross-validated (10-fold) data with accuracy up to 79%. A total of three non-correlating descriptors, identified best by the decision tree analysis, were subsequently utilized for development of suitable models using moving average analysis. These proposed models resulted in the prediction of (V)600(E)BRAF inhibitory activity (IC50) and melanoma cells growth (SRB GI50) inhibitory activity with an overall accuracy of ≥90%. The statistical significance of models/descriptors was assessed through intercorrelation analysis, sensitivity, specificity and Matthew's correlation coefficient.

Computational study on the interaction of N1 substituted pyrazole derivatives with B-raf kinase: an unusual water wire hydrogen-bond network and novel interactions at the entrance of the active site

Docking and molecular dynamics (MD) simulations of N1 substituted pyrazole derivatives complexed with B-Raf kinase were performed to gain insight into the structural and energetic preferences of these inhibitors. First, a comparative study of fully automated docking programs AutoDock, ICM, GLIDE, and Surflex-Dock in closely approximating the X-ray crystal structure of the inhibitor (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime was performed. Afterward, the dynamics of the above-mentioned compound and the less active analogous compounds with 1-methyl-4-piperidinyl and tetrahydro-2H-pyran-4-yl groups at position N1 of pyrazole ring inside the B-Raf active site were analyzed by MD simulations. We found that the most active compound has stable interactions with residues Ile463 and His539 at the entrance of the B-Raf active site. Those interactions were in very good agreement with more reliable quantum mechanics/molecular mechanics calculations performed on the torsional angle phi between the pyrazole ring and the substituents at position N1. In addition, we identified a water wire connecting N2 of the pyrazole ring, Cys532, and Ser536, which is composed of three water molecules for the most active compound. We found some differences in the water wire hydrogen-bond network formed by less active compounds. We suggest that the differences between these structural features are responsible for the differences in activity among the studied compounds.