Design and synthesis of novel macrocyclic 2-amino-6-arylpyrimidine Hsp90 inhibitors

Drug discovery and optimization based on the co-crystal structure of natural product with target

Due to their structural diversities and prevalent biological activities, natural products (NPs) are momentous resources for drug discovery. Although NPs have a wide range of biological activities, many exhibit structural complexity that leads to synthetic difficulties, which combines with inefficient biological activity, toxicity, and unfavorable pharmacokinetic characteristics and ultimately imparts poor safety and efficacy outcomes. Progress in crystallization and computational techniques allow crystallography to have a seasonable influences on drug discovery. By co-crystallizing with proteins, therapeutic targets of NPs in specific diseases can be identified. By analyzing the co-crystal information, the structure-activity relationships (SARs) of NPs targeting specific proteins can be grasped. Under the guidance of co-crystal information, directional structural modification and simplification are powerful strategies for overcoming limitations of NPs, improving the success rate of NP-based drug discovery, and obtaining NP-based drugs with high selectivity, low toxicity and favorable pharmacokinetic characteristics. Here, we review the co-crystal information of a selection of NPs, focusing on the SARs of NPs reflected by co-crystal information and the modification and simplification strategies of NPs, and discuss how to apply co-crystal information in the optimization of NP-based lead compound.

Heterocyclic Compounds as Hsp90 Inhibitors: A Perspective on Anticancer Applications

Heat shock proteins (Hsps) have garnered special attention in cancer therapy as molecular chaperones with regulatory/mediatory effects on folding, maintenance/stability, maturation, and conformation of proteins as well as their effects on prevention of protein aggregation. Hsp90 ensures the stability of various client proteins needed for the growth of cells or the survival of tumor cells; therefore, they are overexpressed in tumor cells and play key roles in carcinogenesis. Accordingly, Hsp90 inhibitors are recognized as attractive therapeutic agents for investigations pertaining to tumor suppression. Natural Hsp90 inhibitors comprising geldanamycin (GM), reclaimed analogs of GM including 17-AAG and DMAG, and radicicol, a natural macrocyclic antifungal, are among the first potent Hsp90 inhibitors. Herein, recently synthesized heterocyclic compounds recognized as potent Hsp90 inhibitors are reviewed along with the anticancer effects of heterocyclic compounds, comprising purine, pyrazole, triazine, quinolines, coumarin, and isoxazoles molecules.

An Insight into the Medicinal Chemistry Perspective of Macrocyclic Derivatives with Antitumor Activity: A Systematic Review

The profound pharmacological properties of macrocyclic compounds have led to their development as drugs. In conformationally pre-organized ring structures, the multiple functions and stereochemical complexity provided by the macrocycle result in high affinity and selectivity of protein targets while maintaining sufficient bioavailability to reach intracellular locations. Therefore, the construction of macrocycles is an ideal choice to solve the problem of “undruggable” targets. Inspection of 68 macrocyclic drugs on the market showed that 10 of them were used to treat cancer, but this structural class still has been poorly explored within drug discovery. This perspective considers the macrocyclic compounds used for anti-tumor with different targets, their advantages and disadvantages, and the various synthetic methods of them.

In vivo-in vitro correlation of antitumor activity of heat shock protein 90 (HSP90) inhibitors with a pharmacokinetics/pharmacodynamics analysis using NCI-N87 xenograft mice

The in vitro antitumor activity (e.g. IC₅₀) of anticancer drugs is important for selecting candidate compounds for in vivo drug efficacy study in the early stage of drug discovery. In this study, we investigated the relationship between in vitro IC₅₀ and in vivo EC₅₀ using six heat shock protein 90 (HSP90) inhibitors.IC₅₀ of each compound was calculated from in vitro cell proliferation assays using the NCI-N87 cancer cell line. Each compound was administered to NCI-N87 xenograft mice, and EC₅₀ and the maximum tumour-killing rate constant were calculated from pharmacokinetics/pharmacodynamics analyses using plasma concentrations and tumour volumes.IC₅₀ obtained in vitro was poorly correlated with EC₅₀ obtained in vivo, while a good correlation (r = 0.856) was observed between them when corrected with the unbound fraction ratio.The results of this study using of HSP90 inhibitors as model compounds suggest importance of the consideration of an unbound fraction to evaluate the relationship between IC₅₀ and EC₅₀. These results will contribute to improvement in the prediction accuracy of in vivo drug efficacy from in vitro activity and the efficiency of drug discovery research.

Scaffold Hopping Transformations Using Auxiliary Restraints for Calculating Accurate Relative Binding Free Energies

In silico screening of drug-target interactions is a key part of the drug discovery process. Changes in the drug scaffold via contraction or expansion of rings, the breaking of rings, and the introduction of cyclic structures from acyclic structures are commonly applied by medicinal chemists to improve binding affinity and enhance favorable properties of candidate compounds. These processes, commonly referred to as scaffold hopping, are challenging to model computationally. Although relative binding free energy (RBFE) calculations have shown success in predicting binding affinity changes caused by perturbing R-groups attached to a common scaffold, applications of RBFE calculations to modeling scaffold hopping are relatively limited. Scaffold hopping inevitably involves breaking and forming bond interactions of quadratic functional forms, which is highly challenging. A novel method for handling ring opening/closure/contraction/expansion and linker contraction/expansion is presented here. To the best of our knowledge, RBFE calculations on linker contraction/expansion have not been previously reported. The method uses auxiliary restraints to hold the atoms at the ends of a bond in place during the breaking and forming of the bonds. The broad applicability of the method was demonstrated by examining perturbations involving small-molecule macrocycles and mutations of proline in proteins. High accuracy was obtained using the method for most of the perturbations studied. The rigor of the method was isolated from the force field by validating the method using relative and absolute hydration free energy calculations compared to standard simulation results. Unlike other methods that rely on λ-dependent functional forms for bond interactions, the method presented here can be employed using modern molecular dynamics software without modification of codes or force field functions.

Synthetic Opportunities and Challenges for Macrocyclic Kinase Inhibitors

Macrocycles are typically cyclic variants of inhibitors derived from uncyclized canonical molecules or from natural products. For medicinal chemistry, drug-like macrocycles have received increasing interest over the past few years, since it has been demonstrated that macrocyclization can favorably alter the biological and physiochemical properties as well as selectivity in comparison to the acyclic analogue. Recent drug approvals such as Lorlatinib, glecaprevir, or voxilaprevir underline the clinical relevance of drug-like macrocycles. However, the synthesis of drug-like macrocycles can be challenging, since the ring-closing reaction is generally challenging with yields depending on the size and geometry of the bridging linker. Nevertheless, macrocycles are one opportunity to expand the synthetic toolbox for medicinal chemistry to provide bioactive molecules. Therefore, we reviewed the past literature of drug-like macrocycles highlighting reactions that have been successfully used for the macrocyclization. We classified the cyclization reactions by their type, ring-size, yield, and macrocyclization efficiency index.

1,3,5-Triazine Analogs: A Potent Anticancer Scaffold

Background:

This review presents the exhaustive exploration of 1,3,5-triazine scaffold for development of analogs of anticancer drugs, over the last century. In the recent years, striazine moiety has been one of the most studied moiety, showing broad-spectrum pharmacological activities such as antibacterial, antifungal, analgesic, anti-HIV, antileishmanial, antitrypanosomal, antimalarial and antiviral. Nowadays, many boffins are have become interested in novel synthesis of s-triazine derivatives because of low cost and ease of availability.

Methods:

This scaffold has been extensively investigated mainly in the past decade. Many products have been synthesized from different starting materials and these synthetic products possess anticancer potential against various cell lines.

Results:

Many 1,3,5-triazine analogs exhibited significant anticancer activity in various models and cell lines exhibiting different mechanisms. Some analogs have also shown good pharmacokinetic parameters with less IC50 values.

Conclusion:

Various 1,3,5-triazine analogs have shown potent activities and may be regarded as clinical candidates for future anticancer formulations. This review may be helpful to those researchers seeking required information with regard to the drug design and medicinal properties of 1,3,5-triazine derivatives for selected targets. This review may also offer help to find and improve clinically viable anticancer molecules.

Design, synthesis, and molecular docking studies of novel pyrazolyl 2-aminopyrimidine derivatives as HSP90 inhibitors

A series of novel pyrazolyl 2-aminopyrimidine derivatives (7a-t) were designed based on scaffold hopping techniques, synthesized and biologically evaluated for their HSP90 inhibition and anticancer activity. Several compounds exhibited potent HSP90 inhibition with IC₅₀ values less than that of the reference standard 17-AAG (1.25 µM). The most potent compound 7t displayed excellent HSP90 inhibition with an IC₅₀ of 20 nM and in vitro antiproliferative potential against three cancer cell lines (IC₅₀  < 5 µM). 7t also induced dose dependent degradation of client proteins (pHER2 and pERK1/2) in Western blot analysis. Several structural features of 7p-t oriented the molecules to retain all the essential binding interactions with HSP90, as observed by rationalized docking studies. Therefore, the para-nitrophenyl ring on the central pyrazole ring along with the 2-amino group on the pyrimidine ring are the crucial features in the development of novel HSP90 inhibitors based on this scaffold for targeted anticancer therapy.

Chaperome heterogeneity and its implications for cancer study and treatment

The chaperome is the collection of proteins in the cell that carry out molecular chaperoning functions. Changes in the interaction strength between chaperome proteins lead to an assembly that is functionally and structurally distinct from each constituent member. In this review, we discuss the epichaperome, the cellular network that forms when the chaperome components of distinct chaperome machineries come together as stable, functionally integrated, multimeric complexes. In tumors, maintenance of the epichaperome network is vital for tumor survival, rendering them vulnerable to therapeutic interventions that target critical epichaperome network components. We discuss how the epichaperome empowers an approach for precision medicine cancer trials where a new target, biomarker, and relevant drug candidates can be correlated and integrated. We introduce chemical biology methods to investigate the heterogeneity of the chaperome in a given cellular context. Lastly, we discuss how ligand-protein binding kinetics are more appropriate than equilibrium binding parameters to characterize and unravel chaperome targeting in cancer and to gauge the selectivity of ligands for specific tumor-associated chaperome pools.

Accurate and Reliable Prediction of the Binding Affinities of Macrocycles to Their Protein Targets

Macrocycles have been emerging as a very important drug class in the past few decades largely due to their expanded chemical diversity benefiting from advances in synthetic methods. Macrocyclization has been recognized as an effective way to restrict the conformational space of acyclic small molecule inhibitors with the hope of improving potency, selectivity, and metabolic stability. Because of their relatively larger size as compared to typical small molecule drugs and the complexity of the structures, efficient sampling of the accessible macrocycle conformational space and accurate prediction of their binding affinities to their target protein receptors poses a great challenge of central importance in computational macrocycle drug design. In this article, we present a novel method for relative binding free energy calculations between macrocycles with different ring sizes and between the macrocycles and their corresponding acyclic counterparts. We have applied the method to seven pharmaceutically interesting data sets taken from recent drug discovery projects including 33 macrocyclic ligands covering a diverse chemical space. The predicted binding free energies are in good agreement with experimental data with an overall root-mean-square error (RMSE) of 0.94 kcal/mol. This is to our knowledge the first time where the free energy of the macrocyclization of linear molecules has been directly calculated with rigorous physics-based free energy calculation methods, and we anticipate the outstanding accuracy demonstrated here across a broad range of target classes may have significant implications for macrocycle drug discovery.

Computational Macrocyclization: From de novo Macrocycle Generation to Binding Affinity Estimation

Macrocycles play an increasing role in drug discovery, but their synthesis is often demanding. Computational tools that suggest macrocyclization based on a known binding mode and that estimate the binding affinity of these macrocycles could have a substantial impact on the medicinal chemistry design process. For both tasks, we established a workflow with high practical value. For five diverse pharmaceutical targets we show that the effect of macrocyclization on binding can be calculated robustly and accurately. Applying this method to macrocycles designed by LigMac, a search tool for de novo macrocyclization, our results suggest that we have a robust protocol in hand to design macrocycles and prioritize them prior to synthesis.

1,3,5-Triazines: A promising scaffold for anticancer drugs development

This review covering literature reports from the beginning of this century to 2016 describes the synthetic pathways, the antitumor activity, the structure-activity relationship and, whenever reported, the possible mechanism of action of 1,3,5-triazine derivatives as well as of their hetero-fused compounds. Many 1,3,5-triazine derivatives, both uncondensed and hetero-fused, have shown remarkable antitumor activities and some of them reached clinical development.

ForceGen 3D structure and conformer generation: from small lead-like molecules to macrocyclic drugs

We introduce the ForceGen method for 3D structure generation and conformer elaboration of drug-like small molecules. ForceGen is novel, avoiding use of distance geometry, molecular templates, or simulation-oriented stochastic sampling. The method is primarily driven by the molecular force field, implemented using an extension of MMFF94s and a partial charge estimator based on electronegativity-equalization. The force field is coupled to algorithms for direct sampling of realistic physical movements made by small molecules. Results are presented on a standard benchmark from the Cambridge Crystallographic Database of 480 drug-like small molecules, including full structure generation from SMILES strings. Reproduction of protein-bound crystallographic ligand poses is demonstrated on four carefully curated data sets: the ConfGen Set (667 ligands), the PINC cross-docking benchmark (1062 ligands), a large set of macrocyclic ligands (182 total with typical ring sizes of 12-23 atoms), and a commonly used benchmark for evaluating macrocycle conformer generation (30 ligands total). Results compare favorably to alternative methods, and performance on macrocyclic compounds approaches that observed on non-macrocycles while yielding a roughly 100-fold speed improvement over alternative MD-based methods with comparable performance.

Binding-affinity predictions of HSP90 in the D3R Grand Challenge 2015 with docking, MM/GBSA, QM/MM, and free-energy simulations

We have estimated the binding affinity of three sets of ligands of the heat-shock protein 90 in the D3R grand challenge blind test competition. We have employed four different methods, based on five different crystal structures: first, we docked the ligands to the proteins with induced-fit docking with the Glide software and calculated binding affinities with three energy functions. Second, the docked structures were minimised in a continuum solvent and binding affinities were calculated with the MM/GBSA method (molecular mechanics combined with generalised Born and solvent-accessible surface area solvation). Third, the docked structures were re-optimised by combined quantum mechanics and molecular mechanics (QM/MM) calculations. Then, interaction energies were calculated with quantum mechanical calculations employing 970-1160 atoms in a continuum solvent, combined with energy corrections for dispersion, zero-point energy and entropy, ligand distortion, ligand solvation, and an increase of the basis set to quadruple-zeta quality. Fourth, relative binding affinities were estimated by free-energy simulations, using the multi-state Bennett acceptance-ratio approach. Unfortunately, the results were varying and rather poor, with only one calculation giving a correlation to the experimental affinities larger than 0.7, and with no consistent difference in the quality of the predictions from the various methods. For one set of ligands, the results could be strongly improved (after experimental data were revealed) if it was recognised that one of the ligands displaced one or two water molecules. For the other two sets, the problem is probably that the ligands bind in different modes than in the crystal structures employed or that the conformation of the ligand-binding site or the whole protein changes.

Enhanced antitumor activity of erlotinib in combination with the Hsp90 inhibitor CH5164840 against non-small-cell lung cancer

Inhibition of heat shock protein 90 (Hsp90) can lead to degradation of multiple client proteins, which are involved in tumor progression. Epidermal growth factor receptor (EGFR) is one of the most potent oncogenic client proteins of Hsp90. Targeted inhibition of EGFR has shown clinical efficacy in the treatment of patients with non-small-cell lung cancer (NSCLC). However, primary and acquired resistance to the existing EGFR inhibitors is a major clinical problem. In the present study, we investigated the effect of the novel Hsp90 inhibitor CH5164840 on the antitumor activity of erlotinib. The NSCLC cell lines and xenograft models were treated with CH5164840 and erlotinib to examine their mechanisms of action and cell growth inhibition. We found that CH5164840 showed remarkable antitumor activity against NSCLC cell lines and xenograft models. The addition of CH5164840 enhanced the antitumor activity of erlotinib against NCI-H292 EGFR-overexpressing xenograft models. Phosphorylation of Stat3 increased with erlotinib treatment in NCI-H292 cells, which was abrogated by Hsp90 inhibition. Furthermore, in a NCI-H1975 T790M mutation erlotinib-resistant model, CH5164840 enhanced the antitumor activity of erlotinib despite the low efficacy of erlotinib treatment alone. In addition, ERK signaling was effectively suppressed by combination treatment with erlotinib and CH5164840 in a NCI-H1975 erlotinib-resistant model. Taken together, these data indicate that CH5164840 has potent antitumor activity and is highly effective in combination with erlotinib against NSCLC tumors with EGFR overexpression and mutations. Our results support the therapeutic potential of CH5164840 as a Hsp90 inhibitor for combination therapy with EGFR-targeting agents against EGFR-addicted NSCLC.

Targeting protein-protein interactions as an anticancer strategy

The emergence and convergence of cancer genomics, targeted therapies, and network oncology have significantly expanded the landscape of protein-protein interaction (PPI) networks in cancer for therapeutic discovery. Extensive biological and clinical investigations have led to the identification of protein interaction hubs and nodes that are critical for the acquisition and maintenance of characteristics of cancer essential for cell transformation. Such cancer-enabling PPIs have become promising therapeutic targets. With technological advances in PPI modulator discovery and validation of PPI-targeting agents in clinical settings, targeting of PPI interfaces as an anticancer strategy has become a reality. Future research directed at genomics-based PPI target discovery, PPI interface characterization, PPI-focused chemical library design, and patient-genomic subpopulation-driven clinical studies is expected to accelerate the development of the next generation of PPI-based anticancer agents for personalized precision medicine. Here we briefly review prominent PPIs that mediate cancer-acquired properties, highlight recognized challenges and promising clinical results in targeting PPIs, and outline emerging opportunities.