A biophysical framework for double-drugging kinases.

Selective orthosteric inhibition of kinases has been challenging due to the conserved active site architecture of kinases and emergence of resistance mutants. Simultaneous inhibition of distant orthosteric and allosteric sites, which we refer to as "double-drugging", has recently been shown to be effective in overcoming drug resistance. However, detailed biophysical characterization of the cooperative nature between orthosteric and allosteric modulators has not been undertaken. Here, we provide a quantitative framework for double-drugging of kinases employing isothermal titration calorimetry, Förster resonance energy transfer, coupled-enzyme assays, and X-ray crystallography. We discern positive and negative cooperativity for Aurora A kinase (AurA) and Abelson kinase (Abl) with different combinations of orthosteric and allosteric modulators. We find that a conformational equilibrium shift is the main principle governing cooperativity. Notably, for both kinases, we find a synergistic decrease of the required orthosteric and allosteric drug dosages when used in combination to inhibit kinase activities to clinically relevant inhibition levels. X-ray crystal structures of the double-drugged kinase complexes reveal the molecular principles underlying the cooperative nature of double-drugging AurA and Abl with orthosteric and allosteric inhibitors. Finally, we observe a fully closed conformation of Abl when bound to a pair of positively cooperative orthosteric and allosteric modulators, shedding light on the puzzling abnormality of previously solved closed Abl structures. Collectively, our data provide mechanistic and structural insights into rational design and evaluation of double-drugging strategies.

Understanding the Differences of Danusertib's Residence Time in Aurora Kinases A/B: Dissociation Paths and Key Residues Identified using Conventional and Enhanced Molecular Dynamics Simulations.

The accurate experimental estimation of protein-ligand systems' residence time (τ) has become very relevant in drug design projects due to its importance in the last stages of refinement of the drug's pharmacodynamics and pharmacokinetics. It is now well-known that it is not sufficient to estimate the affinity of a protein-drug complex in the thermodynamic equilibrium process in in vitro experiments (closed systems), where the concentrations of the drug and protein remain constant. On the contrary, it is mandatory to consider the conformational dynamics of the system in terms of the binding and unbinding processes between protein and drugs in in vivo experiments (open systems), where their concentrations are in constant flux. This last model has been proven to dictate much of several drugs' pharmacological activities in vivo. At the atomistic level, molecular dynamics simulations can explain why some drugs are more effective than others or unveil the molecular aspects that make some drugs work better in one molecular target. Here, the protein kinases Aurora A/B, complexed with its inhibitor Danusertib, were studied using conventional and enhanced molecular dynamics (MD) simulations to estimate the dissociation paths and, therefore, the computational τ values and their comparison with experimental ones. Using classical molecular dynamics (cMD), three differential residues within the Aurora A/B active site, which seems to play an essential role in the observed experimental Danusertib's residence time against these kinases, were characterized. Then, using WT-MetaD, the relative Danusertib's residence times against Aurora A/B kinases were measured in a nanosecond time scale and were compared to those τ values observed experimentally. In addition, the potential dissociation paths of Danusertib in Aurora A and B were characterized, and differences that might be explained by the differential residues in the enzyme's active sites were found. In perspective, it is expected that this computational protocol can be applied to other protein-ligand complexes to understand, at the molecular level, the differences in residence times and amino acids that may contribute to it.

Allosteric modulation of a human protein kinase with monobodies.

Despite being the subject of intense effort and scrutiny, kinases have proven to be consistently challenging targets in inhibitor drug design. A key obstacle has been promiscuity and consequent adverse effects of drugs targeting the ATP binding site. Here we introduce an approach to controlling kinase activity by using monobodies that bind to the highly specific regulatory allosteric pocket of the oncoprotein Aurora A (AurA) kinase, thereby offering the potential for more specific kinase modulators. Strikingly, we identify a series of highly specific monobodies acting either as strong kinase inhibitors or activators via differential recognition of structural motifs in the allosteric pocket. X-ray crystal structures comparing AurA bound to activating vs inhibiting monobodies reveal the atomistic mechanism underlying allosteric modulation. The results reveal 3 major advantages of targeting allosteric vs orthosteric sites: extreme selectivity, ability to inhibit as well as activate, and avoidance of competing with ATP that is present at high concentrations in the cells. We envision that exploiting allosteric networks for inhibition or activation will provide a general, powerful pathway toward rational drug design.

Novel Aurora A Kinase Inhibitor Fangchinoline Enhances Cisplatin-DNA Adducts and Cisplatin Therapeutic Efficacy in OVCAR-3 Ovarian Cancer Cells-Derived Xenograft Model.

Aurora A kinase (Aurora A) is a serine/threonine kinase regulating control of multiple events during cell-cycle progression. Playing roles in promoting proliferation and inhibiting cell death in cancer cells leads Aurora A to become a target for cancer therapy. It is overexpressed and associated with a poor prognosis in ovarian cancer. Improving cisplatin therapy outcomes remains an important issue for advanced-stage ovarian cancer treatment, and Aurora A inhibitors may improve it. In the present study, we identified natural compounds with higher docking scores than the known Aurora A ligand through structure-based virtual screening, including the natural compound fangchinoline, which has been associated with anticancer activities but not yet investigated in ovarian cancer. The binding and inhibition of Aurora A by fangchinoline were verified using cellular thermal shift and enzyme activity assays. Fangchinoline reduced viability and proliferation in ovarian cancer cell lines. Combination fangchinoline and cisplatin treatment enhanced cisplatin-DNA adduct levels, and the combination index revealed synergistic effects on cell viability. An in vivo study showed that fangchinoline significantly enhanced cisplatin therapeutic effects in OVCAR-3 ovarian cancer-bearing mice. Fangchinoline may inhibit tumor growth and enhance cisplatin therapy in ovarian cancer. This study reveals a novel Aurora A inhibitor, fangchinoline, as a potentially viable adjuvant for ovarian cancer therapy.

Role of AURKA in the hypothalamus-pituitary-testicular axis in Tibetan sheep from Tianzhu.

The mitosis-associated protein aurora kinase A (AURKA) regulates the maturation of germ cells. We have previously reported using transcriptome analysis that AURKA is expressed in yak testes. Although Tibetan sheep possess an immense economic value, their reproductive rate is low. Herein, the expression and functions of AURKA in the hypothalamus-pituitary-testicular (HPT) axis in Tibetan sheep from Tianzhu were investigated. The cDNA sequence of sheep AURKA was cloned and bioinformatics techniques were used to predict its structure. Tissue expression of AURKA was determined by qPCR, immunoblotting, immunostaining, and immunohistochemistry. The AURKA coding sequence was found to be 1218 bp in length, encoding a 405-amino acid polypeptide chain. Furthermore, the highest sequence similarity of AURKA with the corresponding sequence in other species was seen in goat and cattle; the least degree of similarity was seen in the domestic cat. In addition, AURKA expression was elevated in the testes compared to that in the hypothalamus and pituitary (p < 0.01). Moreover, AURKA was mainly localized in the hypothalamic paraventricular nucleus (magnocellular), chromophobe cells of the pituitary, and spermatogenic cells of the testis. These results indicated that AURKA might participate in sheep reproductive regulation, thus providing a reference for the study of AURKA function in the reproductive process of Tibetan sheep from Tianzhu.

Protein kinase CK2 phosphorylation of SAPS3 subunit increases PP6 phosphatase activity with Aurora A kinase.

Protein Ser/Thr phosphatase-6 (PP6) regulates pathways for activation of NF-kB, YAP1 and Aurora A kinase (AURKA). PP6 is a heterotrimer comprised of a catalytic subunit, one of three different SAPS subunits and one of three different ankyrin-repeat ANKRD subunits. Here, we show FLAG-PP6C expressed in cells preferentially binds endogenous SAPS3, and the complex is active with the chemical substrate DiFMUP. SAPS3 has multiple acidic sequence motifs recognized by protein kinase CK2 (CK2) and SAPS3 is phosphorylated by purified CK2, without affecting its associated PP6 phosphatase activity. However, HA3-SAPS3-PP6 phosphatase activity using pT288 AURKA as substrate is significantly increased by phosphorylation with CK2. The substitution of Ala in nine putative phosphorylation sites in SAPS3 was required to prevent CK2 activation of the phosphatase. Different CK2 chemical inhibitors equally increased phosphorylation of endogenous AURKA in living cells, consistent with reduction in PP6 activity. CRISPR/Cas9 deletion or siRNA knockdown of SAPS3 resulted in highly activated endogenous AURKA, and a high proportion of cells with abnormal nuclei. Activation of PP6 by CK2 can form a feedback loop with bistable changes in substrates.

PHA-680626 Is an Effective Inhibitor of the Interaction between Aurora-A and N-Myc.

Neuroblastoma is a severe childhood disease, accounting for ~10% of all infant cancers. The amplification of the MYCN gene, coding for the N-Myc transcription factor, is an essential marker correlated with tumor progression and poor prognosis. In neuroblastoma cells, the mitotic kinase Aurora-A (AURKA), also frequently overexpressed in cancer, prevents N-Myc degradation by directly binding to a highly conserved N-Myc region. As a result, elevated levels of N-Myc are observed. During recent years, it has been demonstrated that some ATP competitive inhibitors of AURKA also cause essential conformational changes in the structure of the activation loop of the kinase that prevents N-Myc binding, thus impairing the formation of the AURKA/N-Myc complex. In this study, starting from a screening of crystal structures of AURKA in complexes with known inhibitors, we identified additional compounds affecting the conformation of the kinase activation loop. We assessed the ability of such compounds to disrupt the interaction between AURKA and N-Myc in vitro, using Surface Plasmon Resonance competition assays, and in tumor cell lines overexpressing MYCN, by performing Proximity Ligation Assays. Finally, their effects on N-Myc cellular levels and cell viability were investigated. Our results identify PHA-680626 as an amphosteric inhibitor both in vitro and in MYCN overexpressing cell lines, thus expanding the repertoire of known conformational disrupting inhibitors of the AURKA/N-Myc complex and confirming that altering the conformation of the activation loop of AURKA with a small molecule is an effective strategy to destabilize the AURKA/N-Myc interaction in neuroblastoma cancer cells.

Quantitative Lysine Reactivity Profiling Reveals Conformational Inhibition Dynamics and Potency of Aurora A Kinase Inhibitors.

Small-molecule inhibitors of protein kinases attract widespread interest in the field of disease therapy because of their high specificity and ease of administration. However, dissecting the conformational inhibition dynamics of kinase inhibitors is still challenging. Here, simultaneously monitoring the conformational inhibition details and potency of Aurora A kinase inhibitors has been achieved by active isotope dimethyl labeling coupled with mass spectrometry-based quantitative lysine reactivity profiling. The conformational effects of inhibitors on lysine reactivity can be globally quantified to feasibly reveal the regions involved in the kinase dynamic inhibition. The half-maximum disturbance concentrations (DC50 values) of the conformation-specific lysine residues could directly represent the conformational selectivity and potency of kinase inhibitors. Further, K309 is discovered as a novel hotspot contributing to the inhibition of Aurora A kinase via the specific rotation of kinase activation loop. This quantitative lysine reactivity profiling strategy might greatly promote the development of targeted drugs.

Synthesis and identification of 2,4-bisanilinopyrimidines bearing 2,2,6,6-tetramethylpiperidine-N-oxyl as potential Aurora A inhibitors.

The Aurora kinases are a family of serine/threonine kinases that interact with components of the mitotic apparatus and serve as potential therapeutic targets in oncology. Herein, we reported a series of 2,4-bisanilinopyrimidines bearing 2,2,6,6-tetramethylpiperidine-N-oxyl with selective inhibition of Aurora A in either enzymatic assays or cellular phenotypic assays, and displaying more potent anti-proliferation compared with that of VX-680. The most potent compound 10a forms better interaction with Aurora A than Aurora B in molecular docking. Mechanistic studies revealed that 10a disrupt the spindle formation, block the cell cycle progression in the G2/M phase and induce apoptosis in HeLa cell. These results suggested that the produced series of compounds are potential anticancer agents for further development as selective Aurora A inhibitors.

Hydrophobic Core Variations Provide a Structural Framework for Tyrosine Kinase Evolution and Functional Specialization.

Protein tyrosine kinases (PTKs) are a group of closely related enzymes that have evolutionarily diverged from serine/threonine kinases (STKs) to regulate pathways associated with multi-cellularity. Evolutionary divergence of PTKs from STKs has occurred through accumulation of mutations in the active site as well as in the commonly conserved hydrophobic core. While the functional significance of active site variations is well understood, relatively little is known about how hydrophobic core variations contribute to PTK evolutionary divergence. Here, using a combination of statistical sequence comparisons, molecular dynamics simulations, mutational analysis and in vitro thermostability and kinase assays, we investigate the structural and functional significance of key PTK-specific variations in the kinase core. We find that the nature of residues and interactions in the hydrophobic core of PTKs is strikingly different from other protein kinases, and PTK-specific variations in the core contribute to functional divergence by altering the stability and dynamics of the kinase domain. In particular, a functionally critical STK-conserved histidine that stabilizes the regulatory spine in STKs is selectively mutated to an alanine, serine or glutamate in PTKs, and this loss-of-function mutation is accommodated, in part, through compensatory PTK-specific interactions in the core. In particular, a PTK-conserved phenylalanine in the I-helix appears to structurally and functionally compensate for the loss of STK-histidine by interacting with the regulatory spine, which has far-reaching effects on enzyme activity, inhibitor sensing, and stability. We propose that hydrophobic core variations provide a selective advantage during PTK evolution by increasing the conformational flexibility, and therefore the allosteric potential of the kinase domain. Our studies also suggest that Tyrosine Kinase Like kinases such as RAF are intermediates in PTK evolutionary divergence inasmuch as they share features of both PTKs and STKs in the core. Finally, our studies provide an evolutionary framework for identifying and characterizing disease and drug resistance mutations in the kinase core.

Structural basis of N-Myc binding by Aurora-A and its destabilization by kinase inhibitors.

Myc family proteins promote cancer by inducing widespread changes in gene expression. Their rapid turnover by the ubiquitin-proteasome pathway is regulated through phosphorylation of Myc Box I and ubiquitination by the E3 ubiquitin ligase SCFFbxW7 However, N-Myc protein (the product of the MYCN oncogene) is stabilized in neuroblastoma by the protein kinase Aurora-A in a manner that is sensitive to certain Aurora-A-selective inhibitors. Here we identify a direct interaction between the catalytic domain of Aurora-A and a site flanking Myc Box I that also binds SCFFbxW7 We determined the crystal structure of the complex between Aurora-A and this region of N-Myc to 1.72-Å resolution. The structure indicates that the conformation of Aurora-A induced by compounds such as alisertib and CD532 is not compatible with the binding of N-Myc, explaining the activity of these compounds in neuroblastoma cells and providing a rational basis for the design of cancer therapeutics optimized for destabilization of the complex. We also propose a model for the stabilization mechanism in which binding to Aurora-A alters how N-Myc interacts with SCFFbxW7 to disfavor the generation of Lys48-linked polyubiquitin chains.

Systematic Computational Design and Identification of Low Picomolar Inhibitors of Aurora Kinase A.

Aurora kinase A (AKA) has served as an effective molecular target for the development of cancer therapeutics. A series of potent AKA inhibitors with the (4-methoxy-pyrimidin-2-yl)-phenyl-amine (MPPA) scaffold are identified using a systematic computer-aided drug design protocol involving structure-based virtual screening, de novo design, and free energy perturbation (FEP) simulations. To enhance the accuracy of the virtual screening to find a proper molecular core and de novo design to optimize biochemical potency, we preliminarily improved the scoring function by implementing a reliable hydration energy term. The overall design strategy proves successful to the extent that some inhibitors reveal exceptionally high potency at low picomolar levels; this was achieved by substituting phenyl, chlorine, and tetrazole moieties on the MPPA scaffold. The establishment of bidentate hydrogen bonds with backbone groups in the hinge region appears to be necessary for the high biochemical potency, consistent with the literature X-ray crystallographic data. The picomolar inhibitory activity also stems from the simultaneous formation of additional hydrogen bonds with the side chains of the hinge region and P-loop residues. The FEP simulation results show that the inhibitory activity surges to the low picomolar level because the interactions in the ATP-binding site of AKA become strong by structural modifications enough to overbalance the increase in dehydration cost. Because of the exceptionally high biochemical potency, the AKA inhibitors reported in this study are anticipated to serve as a new starting point for the discovery of anticancer medicine.

A moving target: structure and disorder in pursuit of Myc inhibitors.

The Myc proteins comprise a family of ubiquitous regulators of gene expression implicated in over half of all human cancers. They interact with a large number of other proteins, such as transcription factors, chromatin-modifying enzymes and kinases. Remarkably, few of these interactions have been characterized structurally. This is at least in part due to the intrinsically disordered nature of Myc proteins, which adopt a defined conformation only in the presence of binding partners. Owing to this behaviour, crystallographic studies on Myc proteins have been limited to short fragments in complex with other proteins. Most recently, we determined the crystal structure of Aurora-A kinase domain bound to a 28-amino acid fragment of the N-Myc transactivation domain. The structure reveals an α-helical segment within N-Myc capped by two tryptophan residues that recognize the surface of Aurora-A. The kinase domain acts as a molecular scaffold, independently of its catalytic activity, upon which this region of N-Myc becomes ordered. The binding site for N-Myc on Aurora-A is disrupted by certain ATP-competitive inhibitors, such as MLN8237 (alisertib) and CD532, and explains how these kinase inhibitors are able to disrupt the protein-protein interaction to affect Myc destabilization. Structural studies on this and other Myc complexes will lead to the design of protein-protein interaction inhibitors as chemical tools to dissect the complex pathways of Myc regulation and function, which may be developed into Myc inhibitors for the treatment of cancer.

SERS and MD simulation studies of a kinase inhibitor demonstrate the emergence of a potential drug discovery tool.

We demonstrate the use of surface-enhanced Raman spectroscopy (SERS) as an excellent tool for identifying the binding site of small molecules on a therapeutically important protein. As an example, we show the specific binding of the common antihypertension drug felodipine to the oncogenic Aurora A kinase protein via hydrogen bonding interactions with Tyr-212 residue to specifically inhibit its activity. Based on SERS studies, molecular docking, molecular dynamics simulation, biochemical assays, and point mutation-based validation, we demonstrate the surface-binding mode of this molecule in two similar hydrophobic pockets in the Aurora A kinase. These binding pockets comprise the same unique hydrophobic patches that may aid in distinguishing human Aurora A versus human Aurora B kinase in vivo. The application of SERS to identify the specific interactions between small molecules and therapeutically important proteins by differentiating competitive and noncompetitive inhibition demonstrates its ability as a complementary technique. We also present felodipine as a specific inhibitor for oncogenic Aurora A kinase. Felodipine retards the rate of tumor progression in a xenografted nude mice model. This study reveals a potential surface pocket that may be useful for developing small molecules by selectively targeting the Aurora family kinases.

Spatial Compartmentalization Specializes the Function of Aurora A and Aurora B.

Aurora kinase A and B share great similarity in sequences, structures, and phosphorylation motif, yet they show different localizations and play distinct crucial roles. The factors that determine such differences are largely unknown. Here we targeted Aurora A to the localization of Aurora B and found that Aurora A phosphorylates the substrate of Aurora B and substitutes its function in spindle checkpoint. In return, the centrosome targeting of Aurora B substitutes the function of Aurora A in the mitotic entry. Expressing the chimera proteins of the Auroras with exchanged N termini in cells indicates that the divergent N termini are also important for their spatiotemporal localizations and functions. Collectively, we demonstrate that functional divergence of Aurora kinases is determined by spatial compartmentalization, and their divergent N termini also contribute to their spatial and functional differentiation.

An update on the pharmacokinetics and pharmacodynamics of alisertib, a selective Aurora kinase A inhibitor.

Human Aurora kinases, including Aurora kinase A (AURKA), B (AURKB), and C (AURKC), play an essential role in mitotic events such as monitoring of the mitotic checkpoint, creation of bipolar mitotic spindle and alignment of centrosomes on it, also regulating centrosome separation, bio-orientation of chromosomes and cytokinesis. AURKA and AURKB are key regulators of mitosis and centrosome via polymerizing microfilaments and controlling chromatid segregation. In particular, AURKA plays critical roles in the regulation of mitotic entry, centrosome function, bipolar spindle assembly, and chromosome segregation. AURKA has been found to be overexpressed in various solid and haematological cancers and has been linked with poor prognosis. Its important role in cancer initiation, growth, and metastasis has brought the focus to search for potent and selective AURKA inhibitors for cancer treatment. MLN8237, also known as alisertib, is one selective AURKA inhibitor that has shown remarkable anticancer effects in preclinical studies. Alisertib exhibits favourable pharmacokinetic properties. Alisertib has generally showed good partial response rates of 4-52% and good safety profiles in Phase I and II trials when it is solely administered as well as combined with cytotoxic chemotherapeutic drugs. Recently, the multicentre, randomized Phase III study of alisertib in patients with relapsed or refractory peripheral T-cell lymphoma has been discontinued due to unsatisfactory efficacy. The low risk of side effects, accessibility, and effectiveness of alisertib makes it a new promising anticancer therapy and further mechanistic and clinical studies are warranted.

Benzimidazole-Based Quinazolines: In Vitro Evaluation, Quantitative Structure-Activity Relationship, and Molecular Modeling as Aurora Kinase Inhibitors.

A series of benzimidazole-based quinazoline derivatives with different substitutions of primary and secondary amines at the C2 position (1-12) were evaluated for their Aurora kinase inhibitory activities. All compounds except for 3 and 6 showed good activity against Aurora kinase inhibitors, with IC50 values in the range of 0.035-0.532 µM. The ligand efficiency (LE) of the compounds with Aurora A kinase was also determined. The structure-activity relationship and the quantitative structure-activity relationship revealed that the Aurora inhibitory activities of these derivatives primarily depend on the different substitutions of the amine present at the C2 position of the quinazoline core. Molecular docking studies in the active binding site also provided theoretical support for the experimental biological data acquired. The current study identifies a novel class of Aurora kinase inhibitors, which can further be used for the treatment of cancer.

Abrogation of AuroraA-TPX2 by novel natural inhibitors: molecular dynamics-based mechanistic analysis.

INTRODUCTION: Cancer is characterized by uncontrolled cell growth and genetic instabilities. The human Aurora-A kinase protein plays a crucial role in spindle assembly during mitosis and is activated by another candidate oncogene, targeting protein for Xklp2 (TPX2). It has been proposed that dissociation of Aurora A-TPX2 complex leads to disruption of mitotic spindle apparatus, thereby preventing cell division and further tumor growth. MATERIALS AND METHODS: A large natural compound library was docked against the active site of Aurora A-TPX2 complex. The protein-ligand complexes were subjected to molecular dynamics simulation to ascertain their binding stability. The drug properties of the compounds were analyzed to observe their drug-like properties. RESULTS: The virtual screening of natural compound library yielded two high scoring compounds, the first compound CTOM [ZINC ID: 38143674] (Glide score: -9.49) was stable for 17 ns while the second TTOM (Glide score: -9.07) was stable for 15 ns. While CTOM interacted with His280, Thr288 of Aurora A and Tyr34, Lys38 of TPX2, TTOM interacted with Arg285 and Arg286 in addition to the residues involved with CTOM. CONCLUSIONS: We report two natural compounds as potential drugs leads for the disruption of this complex. These ligands show a preferable docking score and have many drugs like properties within in the range of 95% of known drugs. The study provides evidence that CTOM and TTOM can efficiently inhibit the TPX2-mediated activation of Aurora A. Thus, it paves way for an elaborate investigation and establishes the importance of computational approaches as time- and cost-effective techniques.

Identification of ligand efficient, fragment-like hits from an HTS library: structure-based virtual screening and docking investigations of 2H- and 3H-pyrazolo tautomers for Aurora kinase A selectivity.

Furanopyrimidine 1 (IC50 = 273 nM, LE = 0.36, LELP = 10.28) was recently identified by high-throughput screening (HTS) of an in-house library (125,000 compounds) as an Aurora kinase inhibitor. Structure-based hit optimization resulted in lead molecules with in vivo efficacy in a mouse tumour xenograft model, but no oral bioavailability. This is attributed to "molecular obesity", a common problem during hit to lead evolution during which degradation of important molecular properties such as molecular weight (MW) and lipophilicity occurs. This could be effectively tackled by the right choice of hit compounds for optimization. In this regard, ligand efficiency (LE) and ligand efficiency dependent lipophilicity (LELP) indices are more often used to choose fragment-like hits for optimization. To identify hits with appropriate LE, we used a MW cut-off <250, and pyrazole structure to filter HTS library. Next, structure-based virtual screening using software (Libdock and Glide) in the Aurora A crystal structure (PDB ID: 3E5A) was carried out, and the top scoring 18 compounds tested for Aurora A enzyme inhibition. This resulted in the identification of a novel tetrahydro-pyrazolo-isoquinoline hit 7 (IC50 = 852 nM, LE = 0.44, LELP = 8.36) with fragment-like properties suitable for further hit optimization. Moreover, hit 7 was found to be selective for Aurora A (Aurora B IC50 = 35,150 nM) and the possible reasons for selectivity investigated by docking two tautomeric forms (2H- and 3H-pyrazole) of 7 in Auroras A and B (PDB ID: 4AF3) crystal structures. This docking study shows that the major 3H-pyrazole tautomer of 7 binds in Aurora A stronger than in Aurora B.

Benzo[e]pyrimido[5,4-b][1,4]diazepin-6(11H)-one derivatives as Aurora A kinase inhibitors: LQTA-QSAR analysis and detailed systematic validation of the developed model.

Aurora kinases are sub-divided into Aurora A, Aurora B, and Aurora C kinases that are considered as prospective targets for a new class of anticancer drugs. In this work, a 4-D-QSAR model using an LQTA-QSAR approach with previously reported 31 derivatives of benzo[e]pyrimido[5,4 -b][1,4]diazepin -6(11H)-one as potent Aurora kinase A inhibitors has been created. Instead of single conformation, the conformational ensemble profile generated for each ligand by using trajectories and topology information retrieved from molecular dynamics simulations from GROMACS package were aligned and used for the calculation of intermolecular interaction energies at each grid point. The descriptors generated on the basis of these Coulomb and Lennard-Jones potentials as independent variables were used to perform a PLS analysis using biological activity as dependent variable. A good predictive model was generated with nine field descriptors and five latent variables. The model showed [Formula: see text]; [Formula: see text] and [Formula: see text]. This model was further validated systematically by using different validation parameters. This 4D-QSAR model gave valuable information to recognize features essential to adapt and develop novel potential Aurora kinase inhibitors.