Impact of protein and small molecule interactions on kinase conformations.

Protein kinases act as central molecular switches in the control of cellular functions. Alterations in the regulation and function of protein kinases may provoke diseases including cancer. In this study we investigate the conformational states of such disease-associated kinases using the high sensitivity of the kinase conformation (KinCon) reporter system. We first track BRAF kinase activity conformational changes upon melanoma drug binding. Second, we also use the KinCon reporter technology to examine the impact of regulatory protein interactions on LKB1 kinase tumor suppressor functions. Third, we explore the conformational dynamics of RIP kinases in response to TNF pathway activation and small molecule interactions. Finally, we show that CDK4/6 interactions with regulatory proteins alter conformations which remain unaffected in the presence of clinically applied inhibitors. Apart from its predictive value, the KinCon technology helps to identify cellular factors that impact drug efficacies. The understanding of the structural dynamics of full-length protein kinases when interacting with small molecule inhibitors or regulatory proteins is crucial for designing more effective therapeutic strategies.

Use of Structural Alerts for Reactive Metabolites in the Application SpotRM.

Reactive metabolite (RM) formation is widely accepted as playing a crucial role in causing idiosyncratic adverse drug reactions (IADRs), where the liver is most affected. An important goal of drug design is to avoid selection of drug candidates giving rise to RMs and therefore risk causing problems later on involving IADRs. The simplest, initial approach is to avoid test structures that have substructures known or strongly suspected to be associated with IADRs. However, as is evident from the many case reports of IADRs, in most cases a clear association with any (bio)chemical mechanism is lacking, which makes it hard to establish any structure-toxicity relationship. Separate studies of RM formation, in vitro and in vivo, have led to likely evidence and to establishing many structural alerts (SAs) that can be used for fast selection/deselection of planned test compounds. As a background to a discussion of the concept, 25 kinase inhibitor drugs with known problems of hepatotoxicity were probed against a set of SAs contained in the application SpotRM. A clear majority of the probed drugs show liabilities as evident by being flagged by more than one of the fairly established types of SAs. At the same time, no clear SAs were found in three drugs, which is discussed in the broader context of usefulness and selection tactics of SAs in drug design.

Elucidating the Selective Mechanism of Drugs Targeting Cyclin-Dependent Kinases with Integrated MetaD-US Simulation.

Cyclin-dependent kinases (CDKs), including CDK12 and CDK13, play crucial roles in regulating the cell cycle and RNA polymerase II activity, making them vital targets for cancer therapies. SR4835 is a selective inhibitor of CDK12/13, showing significant potential for treating triple-negative breast cancer. To elucidate the selective mechanism of SR4835 among three CDKs (CDK13/12/9), we developed an innovative enhanced sampling method, integrated well-tempered metadynamics-umbrella sampling (IMUS). IMUS synergistically combines the comprehensive pathway exploration capability of well-tempered metadynamics (WT-MetaD) with the precise free energy calculation capability of umbrella sampling, enabling the efficient and accurate characterization of drug-target interactions. The accurate calculation of binding free energy and the detailed analysis of the kinetic mechanism of the drug-target interaction using IMUS successfully elucidate the drug selectivity mechanism targeting the three CDKs, showing that the selectivity is primarily arising from differences in the stability of H-bonds within the Hinge region of the kinases and the interaction patterns during the protein-ligand recognition process. These findings also underscore the utility of IMUS in efficiently and accurately capturing drug-target interaction processes with clear mechanisms.

Identification of Potent ADCK3 Inhibitors through Structure-Based Virtual Screening.

ADCK3 is a member of the UbiB family of atypical protein kinases in humans, with homologues in archaea, bacteria, and eukaryotes. In lieu of protein kinase activity, ADCK3 plays a role in the biosynthesis of coenzyme Q10 (CoQ10), and inactivating mutations can cause a CoQ10 deficiency and ataxia. However, the exact functions of ADCK3 are still unclear, and small-molecule inhibitors could be useful as chemical probes to elucidate its molecular mechanisms. In this study, we applied structure-based virtual screening (VS) to discover a novel chemical series of ADCK3 inhibitors. Through extensive structural analysis of the active-site residues, we developed a pharmacophore model and applied it to a large-scale VS. Out of ∼170,000 compounds virtually screened, 800 top-ranking candidate compounds were selected and tested in both ADCK3 and p38 biochemical assays for hit validation. In total, 129 compounds were confirmed as ADCK3 inhibitors, and among them, 114 compounds are selective against p38, which was used as a counter-target. Molecular dynamics (MD) simulations were then conducted to predict the binding modes of the most potent compounds within the ADCK3 active site. Through metadynamics analysis, we successfully detected the key amino acid residues that govern intermolecular interactions. The findings provided in this study can serve as a promising starting point for drug development.

Molecular Modeling of Single- and Double-Hydrocarbon-Stapled Coiled-Coil Inhibitors against Bcr-Abl: Toward a Treatment Strategy for CML.

The chimeric oncoprotein Bcr-Abl is the causative agent of virtually all chronic myeloid leukemias and a subset of acute lymphoblastic leukemias. As a result of the so-called Philadelphia chromosome translocation t(9;22), Bcr-Abl manifests as a constitutively active tyrosine kinase, which promotes leukemogenesis by activation of cell cycle signaling pathways. Constitutive and oncogenic activation is mediated by an N-terminal coiled-coil oligomerization domain in Bcr (Bcr-CC), presenting a therapeutic target for inhibition of Bcr-Abl activity toward the treatment of Bcr-Abl+ leukemias. Previously, we demonstrated that a rationally designed Bcr-CC mutant, CCmut3, exerts a dominant negative effect upon Bcr-Abl activity by preferential oligomerization with Bcr-CC. Moreover, we have shown that conjugation to a leukemia-specific cell-penetrating peptide (CPP-CCmut3) improves intracellular delivery and activity. However, our full-length CPP-CCmut3 construct (81 aa) is encumbered by an intrinsically high degree of conformational variability and susceptibility to proteolytic degradation relative to traditional small-molecule therapeutics. Here, we iterate a new generation of CCmut3 inhibitors against Bcr-CC-mediated Bcr-Abl assembly designed to address these constraints through incorporation of all-hydrocarbon staples spanning i and i + 7 positions in α-helix 2 (CPP-CCmut3-st). We utilize computational modeling and biomolecular simulation to evaluate single- and double-stapled CCmut3 candidates in silico for dynamics and binding energetics. We further model a truncated system characterized by the deletion of α-helix 1 and the flexible loop linker, which are known to impart high conformational variability. To study the impact of the N-terminal cyclic CPP toward model stability and inhibitor activity, we also model the full-length and truncated systems devoid of the CPP, with a cyclized CPP, and with an open-configuration CPP, for a total of six systems that comprise our library. From this library, we present lead-stapled peptide candidates to be synthesized and evaluated experimentally as our next iteration of inhibitors against Bcr-Abl.

Exploring the binding characteristics of bovine serum albumin with CDK4/6 inhibitors Ribociclib: Multi-spectral analysis and molecular simulation studies.

Ribociclib (RIB), a tyrosine kinase inhibitor, exhibits promising antitumor efficacy and controlled toxicity in HR+/HER2- breast cancer patients, which is closely related to the binding with plasma proteins. This study utilized a combination of spectroscopic techniques including UV spectroscopy, fluorescence spectroscopy, and circular dichroism (CD) as well as molecular docking and molecular dynamic simulation to clarify the binding mechanism between bovine serum albumin (BSA) and RIB. The findings demonstrated that RIB produced a 1:1 stoichiometric complex with BSA, which quenched BSA's fluorescence in the manner of the static quenching mechanism. Site labelling experiments pinpointed Site III on BSA as the primary binding site for RIB, a finding validated by molecular docking. Van der Waals forces and hydrogen bonding interactions as key drivers in the formation of RIB-BSA complexes, a conclusion supported by molecular docking. Molecular simulation studies suggested that the insertion of RIB into the hydrophobic cavity (Site III) of BSA induced subtle conformational changes in the BSA protein, and CD measurements confirmed alterations in BSA secondary structure content. Synchronous and three-dimensional fluorescence spectroscopy further demonstrated that RIB decreased the hydrophobicity of the microenvironment surrounding tyrosine and tryptophan residues. These findings offer valuable insights into the pharmacokinetics and structural modifications of RIB.

Dynamics Playing a Key Role in the Covalent Binding of Inhibitors to Focal Adhesion Kinase.

Covalent kinase inhibitors (CKIs) have recently garnered considerable attention, yet the rational design of CKIs continues to pose a great challenge. In the discovery of CKIs targeting focal adhesion kinase (FAK), it has been observed that the chemical structure of the linkers plays a key role in achieving covalent targeting of FAK. However, the mechanism behind the observation remains elusive. In this work, we employ a comprehensive suite of advanced computational methods to investigate the mechanism of CKIs covalently targeting FAK. We reveal that the linker of an inhibitor influences the contacts between the warhead and residue(s) and the residence time in active conformation, thereby dictating the inhibitor's capability to bind covalently to FAK. This study reflects the complexity of CKI design and underscores the importance of considering the dynamic interactions and residence times for the successful development of covalent drugs.

Dual-Probe Activity-Based Protein Profiling Reveals Site-Specific Differences in Protein Binding of EGFR-Directed Drugs.

Comparative, dose-dependent analysis of interactions between small molecule drugs and their targets, as well as off-target interactions, in complex proteomes is crucial for selecting optimal drug candidates. The affinity of small molecules for targeted proteins is largely dictated by interactions between amino acid side chains and these drugs. Thus, studying drug-protein interactions at an amino acid resolution provides a comprehensive understanding of the drug selectivity and efficacy. In this study, we further refined the site-specific activity-based protein profiling strategy (ABPP), PhosID-ABPP, on a timsTOF HT mass spectrometer. This refinement enables dual dose-dependent competition of inhibitors within a single cellular proteome. Here, a comparative analysis of two activity-based probes (ABPs), developed to selectively target the epidermal growth factor receptor (EGFR), namely, PF-06672131 (PF131) and PF-6422899 (PF899), facilitated the simultaneous identification of ABP-specific binding sites at a proteome-wide scale within a cellular proteome. Dose-dependent probe-binding preferences for proteinaceous cysteines, even at low nanomolar ABP concentrations, could be revealed. Notably, in addition to the intrinsic affinity of the electrophilic probes for specific sites in targeted proteins, the observed labeling intensity is influenced by several other factors. These include the efficiency of cellular uptake, the stability of the probes, and their intracellular distribution. While both ABPs showed comparable labeling efficiency for EGFR, PF131 had a broader off-target reactivity profile. In contrast, PF899 exhibited a higher labeling efficiency for the ERBB2 receptor and bound to catalytic cysteines in several other enzymes, which is likely to disrupt their catalytic activity. Notably, PF131 effectively labeled ADP/ATP translocase proteins at a concentration of just 1 nm, and we found this affected ATP transport. Analysis of the effect of PF131 and its parent inhibitor Afatinib on murine translocase SLC25A4 (ANT1)-mediated ATP transport strongly indicated that PF131 (10 µM) partially blocked ATP transport. Afatinib was less efficient at inhibiting ATP transport by SLC25A4 than PF131, and the reduction of ATP transport by Afatinib was not significant. Follow-up analysis is required to evaluate the affinity of these inhibitors for ADP/ATP translocase SLC25A4 in more detail. Additionally, the analysis of different binding sites within the EGF receptor and the voltage-dependent anion channel 2 revealed secondary binding sites of both probes and provided insights into the binding poses of inhibitors on these proteins. Insights from the PhosID-ABPP analysis of these two ABPs serve as a valuable resource for understanding drug on- and off-target engagement in a dose- and site-specific manner.

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.

Decoding Inhibitor Egression from Wild-Type and G2019S Mutant LRRK2 Kinase: Insights into Unbinding Mechanisms for Precision Drug Design in Parkinson's Disease.

Leucine-rich repeat kinase 2 (LRRK2) remains a viable target for drug development since the discovery of the association of its mutations with Parkinson's disease (PD). G2019S (in the kinase domain) is the most common mutation for LRRK2-based PD. Though various types of inhibitors have been developed for the kinase domain to reduce the effect of the mutation, understanding the working of these inhibitors at the molecular level is still ongoing. This study focused on the exploration of the dissociation mechanism (pathways) of inhibitors from (WT and G2019S) LRRK2 kinase (using homology model CHK1 kinase), which is one of the crucial aspects in drug discovery. Here, two ATP-competitive type I inhibitors, PF-06447475 and MLi-2 (Comp1 and Comp2 ), and one non-ATP-competitive type II inhibitor, rebastinib (Comp3), were considered for this investigation. To study the unbinding process, random accelerated molecular dynamics simulations were performed. The binding free energies of the three inhibitors for different egression paths were determined using umbrella sampling. This work found four major egression pathways that were adopted by the inhibitors Comp1 (path1, path2, and path3), Comp2 (path1, path2 and path3), and Comp3 (path3 and path4). Also, the mechanism of unbinding for each path and key residues involved in unbinding were explored. Mutation was not observed to impact the preference of the particular egression pathways for both LRRK2-Comp1 and -Comp2 systems. However, the findings suggested that the size of the inhibitor molecules might have an effect on the preference of the egression pathways. The binding energy and residence time of the inhibitors followed a similar trend to experimental observations. The findings of this work might provide insight into designing more potent inhibitors for the G2019S LRRK2 kinase.

Computational biophysical characterization of the effect of gatekeeper mutations on the binding of ponatinib to the FGFR kinase.

Fibroblast Growth Factor Receptor (FGFR) is connected to numerous downstream signalling cascades regulating cellular behavior. Any dysregulation leads to a plethora of illnesses, including cancer. Therapeutics are available, but drug resistance driven by gatekeeper mutation impedes the treatment. Ponatinib is an FDA-approved drug against BCR-ABL kinase and has shown effective results against FGFR-mediated carcinogenesis. Herein, we undertake molecular dynamics simulation-based analysis on ponatinib against all the FGFR isoforms having Val to Met gatekeeper mutations. The results suggest that ponatinib is a potent and selective inhibitor for FGFR1, FGFR2, and FGFR4 gatekeeper mutations. The extensive electrostatic and van der Waals interaction network accounts for its high potency. The FGFR3_VM mutation has shown resistance towards ponatinib, which is supported by their lesser binding affinity than wild-type complexes. The disengaged molecular brake and engaged hydrophobic spine were believed to be the driving factors for weak protein-ligand interaction. Taken together, the inhibitory and structural characteristics exhibited by ponatinib may aid in thwarting resistance based on Val-to-Met gatekeeper mutations at an earlier stage of treatment and advance the design and development of other inhibitors targeted at FGFRs harboring gatekeeper mutations.

The Metabolism and Disposition of Brepocitinib in Humans and Characterization of the Formation Mechanism of an Aminopyridine Metabolite.

Brepocitinib is an oral once-daily Janus kinase 1 and Tyrosine kinase 2 selective inhibitor currently in development for the treatment of several autoimmune disorders. Mass balance and metabolic profiles were determined using accelerator mass spectrometry in six healthy male participants following a single oral 60 mg dose of 14C-brepocitinib (∼300 nCi). The average mass balance recovery was 96.7% ± 6.3%, with the majority of dose (88.0% ± 8.0%) recovered in urine and 8.7% ± 2.1% of the dose recovered in feces. Absorption of brepocitinib was rapid, with maximal plasma concentrations of total radioactivity and brepocitinib achieved within 0.5 hours after dosing. Circulating radioactivity consisted primarily of brepocitinib (47.8%) and metabolite M1 (37.1%) derived from hydroxylation at the C5' position of the pyrazole ring. Fractional contributions to metabolism via cytochrome P450 enzymes were determined to be 0.77 for CYP3A4/5 and 0.14 for CYP1A2 based on phenotyping studies in human liver microsomes. However, additional clinical studies are required to understand the potential contribution of CYP1A1. Approximately 83% of the dose was eliminated as N-methylpyrazolyl oxidative metabolites, with 52.1% of the dose excreted as M1 alone. Notably, M1 was not observed as a circulating metabolite in earlier metabolic profiling of human plasma from a multiple ascending dose study with unlabeled brepocitinib. Mechanistic studies revealed that M1 was highly unstable in human plasma and phosphate buffer, undergoing chemical oxidation leading to loss of the 5-hydroxy-1-methylpyrazole moiety and formation of aminopyrimidine cleavage product M2. Time-dependent inhibition and trapping studies with M1 yielded insights into the mechanism of this unusual and unexpected instability. SIGNIFICANCE STATEMENT: This study provides a detailed understanding of the disposition and metabolism of brepocitinib, a JAK1/TYK2 inhibitor for atopic dermatitis, in humans as well as characterization of clearance pathways and pharmacokinetics of brepocitinib and its metabolites.

Allosteric coupling asymmetry mediates paradoxical activation of BRAF by type II inhibitors.

The type II class of RAF inhibitors currently in clinical trials paradoxically activate BRAF at subsaturating concentrations. Activation is mediated by induction of BRAF dimers, but why activation rather than inhibition occurs remains unclear. Using biophysical methods tracking BRAF dimerization and conformation, we built an allosteric model of inhibitor-induced dimerization that resolves the allosteric contributions of inhibitor binding to the two active sites of the dimer, revealing key differences between type I and type II RAF inhibitors. For type II inhibitors the allosteric coupling between inhibitor binding and BRAF dimerization is distributed asymmetrically across the two dimer binding sites, with binding to the first site dominating the allostery. This asymmetry results in efficient and selective induction of dimers with one inhibited and one catalytically active subunit. Our allosteric models quantitatively account for paradoxical activation data measured for 11 RAF inhibitors. Unlike type II inhibitors, type I inhibitors lack allosteric asymmetry and do not activate BRAF homodimers. Finally, NMR data reveal that BRAF homodimers are dynamically asymmetric with only one of the subunits locked in the active αC-in state. This provides a structural mechanism for how binding of only a single αC-in inhibitor molecule can induce potent BRAF dimerization and activation.

Structural and molecular insights from dual inhibitors of EGFR and VEGFR2 as a strategy to improve the efficacy of cancer therapy.

Epidermal growth factor receptor (EGFR) and vascular endothelial growth factor 2 (VEGFR2) are known as valid targets for cancer therapy. Overexpression of EGFR induces uncontrolled cell proliferation and VEGF expression triggering angiogenesis via VEGFR2 signaling. On the other hand, VEGF expression independent of EGFR signaling is already known as one of the mechanisms of resistance to anti-EGFR therapy. Therefore, drugs that act as dual inhibitors of EGFR and VEGFR2 can be a solution to the problem of drug resistance and increase the effectiveness of therapy. In this review, we summarize the relationship between EGFR and VEGFR2 signal transduction in promoting cancer growth and how their kinase domain structures can affect the selectivity of an inhibitor as the basis for designing dual inhibitors. In addition, several recent studies on the development of dual EGFR and VEGFR2 inhibitors involving docking simulations were highlighted in this paper to provide some references such as pharmacophore features of inhibitors and key residues for further research, especially in computer-aided drug design.

PKIB, a Novel Target for Cancer Therapy.

The serine-threonine kinase protein kinase A (PKA) is a cyclic AMP (cAMP)-dependent intracellular protein with multiple roles in cellular biology including metabolic and transcription regulation functions. The cAMP-dependent protein kinase inhibitor ß (PKIB) is one of three known endogenous protein kinase inhibitors of PKA. The role of PKIB is not yet fully understood. Hormonal signaling is correlated with increased PKIB expression through genetic regulation, and increasing PKIB expression is associated with decreased cancer patient prognosis. Additionally, PKIB impacts cancer cell behavior through two mechanisms; the first is the nuclear modulation of transcriptional activation and the second is the regulation of oncogenic AKT signaling. The limited research into PKIB indicates the oncogenic potential of PKIB in various cancers. However, some studies suggest a role of PKIB in non-cancerous disease states. This review aims to summarize the current literature and background of PKIB regarding cancer and related issues. In particular, we will focus on cancer development and therapeutic possibilities, which are of paramount interest in PKIB oncology research.

Guided Docking as a Data Generation Approach Facilitates Structure-Based Machine Learning on Kinases.

Drug discovery pipelines nowadays rely on machine learning models to explore and evaluate large chemical spaces. While including 3D structural information is considered beneficial, structural models are hindered by the availability of protein-ligand complex structures. Exemplified for kinase drug discovery, we address this issue by generating kinase-ligand complex data using template docking for the kinase compound subset of available ChEMBL assay data. To evaluate the benefit of the created complex data, we use it to train a structure-based E(3)-invariant graph neural network. Our evaluation shows that binding affinities can be predicted with significantly higher precision by models that take synthetic binding poses into account compared to ligand- or drug-target interaction models alone.

The Role of Intramolecular Reactions and Chemical Degradation in the Apparent Biotransformation Pathways of a Series of SYK Inhibitors.

In vitro metabolism studies of the spleen tyrosine kinase inhibitors AZ-A and AZ-B identified four unusual metabolites. M1 (mass-to-charge ratio 411) was formed by both molecules and was common to several analogs (AZ-C to AZ-H) sharing the same core structure, appearing to derive from the complete loss of a pendent 3,4-diaminotetrahydropyran ring and pyrazole ring cleavage resulting in a nonobvious metabolite. M2-M4 were formed by AZ-A and a subset of the other compounds only and apparently resulted from a sequential loss of H2 from parent. Initial attempts to isolate M3 for identification were unsuccessful due to sample degradation, and it was subsequently found that M2 and M3 underwent sequential chemical degradation steps to M4. M4 was successfully isolated and shown by mass spectrometry and NMR spectroscopy to be a tricyclic species incorporating the pyrazole and the 3,4-diaminotetrahydropyran groups. We propose that this arises from an intramolecular reaction between the primary amine on the tetrahydropyran and a putative epoxide intermediate on the adjacent pyrazole ring, evidence for which was generated in a ß-mercaptoethanol-trapping experiment. The loss of the tetrahydropyran moiety observed in M1 was found to be enhanced in an analog that was unable to undergo the intramolecular reaction step, leading us to propose two possible reaction pathways originating from the reactive intermediate. Ultimately, we conclude that the apparently complex and unusual metabolism of this series of compounds likely resulted from a single metabolic activation step forming an epoxide intermediate, which subsequently underwent intramolecular rearrangement and/or chemical degradation to form the final observed products. SIGNIFICANCE STATEMENT: The current work provides an unusual biotransformation example showing the potential for intramolecular reactions and chemical degradation to give the appearance of complex metabolism arising from a single primary route of metabolism.

Unexpected Noncovalent Off-Target Activity of Clinical BTK Inhibitors Leads to Discovery of a Dual NUDT5/14 Antagonist.

Cofactor mimicry represents an attractive strategy for the development of enzyme inhibitors but can lead to off-target effects due to the evolutionary conservation of binding sites across the proteome. Here, we uncover the ADP-ribose (ADPr) hydrolase NUDT5 as an unexpected, noncovalent, off-target of clinical BTK inhibitors. Using a combination of biochemical, biophysical, and intact cell NanoBRET assays as well as X-ray crystallography, we confirm catalytic inhibition and cellular target engagement of NUDT5 and reveal an unusual binding mode that is independent of the reactive acrylamide warhead. Further investigation of the prototypical BTK inhibitor ibrutinib also revealed potent inhibition of the largely unstudied NUDIX hydrolase family member NUDT14. By exploring structure-activity relationships (SARs) around the core scaffold, we identify a potent, noncovalent, and cell-active dual NUDT5/14 inhibitor. Cocrystallization experiments yielded new insights into the NUDT14 hydrolase active site architecture and inhibitor binding, thus providing a basis for future chemical probe design.

Low-Dose Sorafenib Promotes Cancer Stem Cell Expansion and Accelerated Tumor Progression in Soft Tissue Sarcomas.

The cancer stem cell (CSC) hypothesis postulates that heterogeneous human cancers harbor a population of stem-like cells which are resistant to cytotoxic therapies, thus providing a reservoir of relapse following conventional therapies like chemotherapy and radiation (RT). CSCs have been observed in multiple human cancers, and their presence has been correlated with worse clinical outcomes. Here, we sought to evaluate the impact of drug dosing of the multi-tyrosine kinase inhibitor, sorafenib, on CSC and non-CSCs in soft tissue sarcoma (STS) models, hypothesizing differential effects of sorafenib based on dose and target cell population. In vitro, human cancer cell lines and primary STS from surgical specimens were exposed to escalating doses of sorafenib to determine cell viability and expression of CSC marker aldehyde dehydrogenase (ALDH). In vivo, ALDHbright CSCs were isolated, exposed to sorafenib, and xenograft growth and survival analyses were performed. We observed that sarcoma CSCs appear to paradoxically respond to the tyrosine kinase inhibitor sorafenib at low doses with increased proliferation and stem-like function of CSCs, whereas anti-viability effects dominated at higher doses. Importantly, STS patients receiving neoadjuvant sorafenib and RT on a clinical trial (NCT00864032) showed increased CSCs post therapy, and higher ALDH scores post therapy were associated with worse metastasis-free survival. These data suggest that low-dose sorafenib may promote the CSC phenotype in STS with clinically significant effects, including increased tumor growth and higher rates of metastasis formation in sarcoma patients.

Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases.

Understanding kinase-inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A-1) of the insulin-like growth factor-I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine-arginine-aspartic acid (HRD) and aspartic acid-phenylalanine-glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C-helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.