A biophysical framework for double-drugging kinases

Advances in uncovering the mechanisms of macromolecular conformational entropy

During protein folding, proteins transition from a disordered polymer into a globular structure, markedly decreasing their conformational degrees of freedom, leading to a substantial reduction in entropy. Nonetheless, folded proteins retain substantial entropy as they fluctuate between the conformations that make up their native state. This residual entropy contributes to crucial functions like binding and catalysis, supported by growing evidence primarily from NMR and simulation studies. Here, we propose three major ways that macromolecules use conformational entropy to perform their functions; first, prepaying entropic cost through ordering of the ground state; second, spatially redistributing entropy, in which a decrease in entropy in one area is reciprocated by an increase in entropy elsewhere; third, populating catalytically competent ensembles, in which conformational entropy within the enzymatic scaffold aids in lowering transition state barriers. We also provide our perspective on how solving the current challenge of structurally defining the ensembles encoding conformational entropy will lead to new possibilities for controlling binding, catalysis and allostery.

The Evolving Landscape of Protein Allostery: From Computational and Experimental Perspectives

Protein allostery is a fundamental biological regulatory mechanism that allows communication between distant locations within a protein, modifying its function in response to signals. Experimental techniques, such as NMR spectroscopy and cryo-electron microscopy (cryo-EM), are critical validation tools for computational predictions and provide valuable insights into dynamic conformational changes. Combining these approaches has greatly improved our understanding of classical conformational allostery and complex dynamic coupling mechanisms. Recent advances in machine learning and enhanced sampling methods have broadened the scope of allostery research, identifying cryptic allosteric sites and directing new drug discovery approaches. Despite progress, bridging static structural data with dynamic functional states remains challenging. This review underscores the importance of combining experimental and computational approaches to comprehensively understand protein allostery and its diverse applications in biology and medicine.

Nanohydrogel of Curcumin/Berberine Co-Crystals Induces Apoptosis via Dual Covalent/Noncovalent Inhibition of Caspases in Endometrial Cancer Cell Lines: The Synergy Between Pharmacokinetics and Pharmacodynamics

Endometrial cancer remains a significant therapeutic challenge due to drug resistance and heterogeneity. This study leverages the synergistic potential of curcumin (CUR) and berberine (BBR) co-crystals encapsulated in a nanohydrogel to address these challenges through a pharmacokinetically and pharmacodynamically targeted therapeutic strategy. The nanohydrogel formulation significantly improves the solubility, stability, and bioavailability of CUR/BBR co-crystals, optimizing their therapeutic delivery and sustained release under physiological and tumor microenvironment conditions. On the other hand, the dual inhibitory mechanism of CUR and BBR, with CUR covalently binding to the active site of caspase-3 and BBR non-covalently targeting the allosteric site, achieves enhanced apoptotic activity by disrupting both the catalytic and conformational functions of caspase-3. In vitro cytotoxicity assays demonstrate remarkable efficacy of the CUR/BBR nanohydrogel, achieving an IC50 of 12.36 μg/mL against HEC-59 endometrial cancer cells, significantly outperforming the individual components and the standard drug Camptothecin (IC50: 17.27 μg/mL). Caspase-3/7 assays confirm enhanced apoptosis induction for the nanohydrogel formulation compared to co-crystals alone and Camptothecin. Molecular dynamics simulations and binding free energy analyses further validate the synergistic interaction of CUR and BBR in their dual binding mode. This study introduces a novel therapeutic approach by enhancing drug delivery and dual targeting mechanisms, demonstrating the potential of CUR-BBR nanohydrogel as a robust therapy for EC. This strategy offers a promising platform for addressing drug resistance and improving outcomes in endometrial cancer therapy.

Novel treatment strategies for chronic myeloid leukemia

ABSTRACT

Starting with imatinib, tyrosine kinase inhibitors (TKIs) have turned chronic myeloid leukemia (CML) from a lethal blood cancer into a chronic condition. As patients with access to advanced CML care have an almost normal life expectancy, there is a perception that CML is a problem of the past, and one should direct research resources elsewhere. However, a closer look at the current CML landscape reveals a more nuanced picture. Most patients still require life-long TKI therapy to avoid recurrence of active CML. Chronic TKI toxicity and the high costs of the well-tolerated agents remain challenging. Progression to blast phase still occurs, particularly in socioeconomically disadvantaged parts of the world, where high-risk CML at diagnosis is common. Here, we review the prospects of further improving TKIs to achieve optimal suppression of BCR::ABL1 kinase activity, the potential of combining different classes of TKIs, and the current state of BCR::ABL1 degraders. We cover combination therapy approaches to address TKI resistance in the setting of residual leukemia and in advanced CML. Despite the unprecedented success of TKIs in CML, more work is needed to truly finish the job, and we hope to stimulate innovative research aiming to achieve this goal.

Discovery of N-(2-Acetamidobenzo[d]thiazol-6-yl)-2-phenoxyacetamide Derivatives as Novel Potential BCR-ABL1 Inhibitors Through Structure-Based Virtual Screening

BCR-ABL1 kinase is a critical driver of chronic myeloid leukemia (CML) pathophysiology. The approval of allosteric inhibitor asciminib brings new hope for overcoming drug resistance caused by mutations in the ATP-binding site. To expand the chemical diversity of BCR-ABL1 kinase inhibitors with positive anti-tumor effect with asciminib, structure-based virtual screening and molecular dynamics simulations were employed to discover novel scaffolds. This approach led to the identification of a series of N-(2-acetamidobenzo[d]thiazol-6-yl)-2-phenoxyacetamide derivatives as new BCR-ABL1 inhibitors. The most potent compound, 10m, demonstrated inhibition of BCR-ABL-dependent signaling and showed an anti-tumor effect against K562 cells, with an IC50 value of 0.98 μM. Compound 10m displayed powerful synergistic anti-proliferation and pro-apoptotic effects when combined with asciminib, highlighting its potential as a promising lead for the development of potential BCR-ABL inhibitors.

An SH3-binding allosteric modulator stabilizes the global conformation of the AML-associated Src-family kinase, Hck

While ATP-site inhibitors for protein-tyrosine kinases are often effective drugs, their clinical utility can be limited by off-target activity and acquired resistance mutations due to the conserved nature of the ATP-binding site. However, combining ATP-site and allosteric kinase inhibitors can overcome these shortcomings in a double-drugging framework. Here we explored the allosteric effects of two pyrimidine diamines, PDA1 and PDA2, on the conformational dynamics and activity of the Src-family tyrosine kinase Hck, a promising drug target for acute myeloid leukemia. Using 1H-15N HSQC NMR, we mapped the binding site for both analogs to the SH3 domain. Despite the shared binding site, PDA1 and PDA2 had opposing effects on near-full-length Hck dynamics by hydrogen-deuterium exchange mass spectrometry, with PDA1 stabilizing and PDA2 disrupting the overall kinase conformation. Kinase activity assays were consistent with these observations, with PDA2 enhancing kinase activity while PDA1 was without effect. Molecular dynamics simulations predicted selective bridging of the kinase domain N-lobe and SH3 domain by PDA1, a mechanism of allosteric stabilization supported by site-directed mutagenesis of N-lobe contact sites. Cellular thermal shift assays confirmed SH3 domain-dependent interaction of PDA1 with WT Hck in myeloid leukemia cells and with a kinase domain gatekeeper mutant (T338M). These results identify PDA1 as a starting point for Src-family kinase allosteric inhibitor development that may work in concert with ATP-site inhibitors to suppress the evolution of resistance.

Latent allosteric control of protein interactions by ATP-competitive kinase inhibitors

Protein kinase inhibitors designed to compete with ATP as a primary mode of action turn out to have considerable effects that go beyond their interference of nucleotide binding. New research shows how kinase activation and sometimes noncatalytic functions of protein kinases can be controlled by allosteric properties of kinase inhibitors, communicating perturbations from the active site to distal regulatory regions.

Cocrystallization of the Src-Family Kinase Hck with the ATP-Site Inhibitor A-419259 Stabilizes an Extended Activation Loop Conformation

Hematopoietic cell kinase (Hck) is a member of the Src kinase family and is a promising drug target in myeloid leukemias. Here, we report the crystal structure of human Hck in complex with the pyrrolopyrimidine inhibitor A-419259, determined at a resolution of 1.8 Å. This structure reveals the complete Hck active site in the presence of A-419259, including the αC-helix, the DFG motif, and the activation loop. A-419259 binds at the ATP-site of Hck and induces an overall closed conformation of the kinase with the regulatory SH3 and SH2 domains bound intramolecularly to their respective internal ligands. A-419259 stabilizes the DFG-in/αC-helix-out conformation observed previously with Hck and the pyrazolopyrimidine inhibitor PP1 (PDB: 1QCF). However, the activation loop conformations are distinct, with PP1 inducing a folded loop structure with the tyrosine autophosphorylation site (Tyr416) pointing into the ATP binding site, while A-419259 stabilizes an extended loop conformation with Tyr416 facing out into the solvent. Autophosphorylation also induces activation loop extension and significantly reduces the Hck sensitivity to PP1 but not A-419259. In cancer cells where Hck is constitutively active, the extended autophosphorylation loop may render Hck more sensitive to inhibitors like A-419259 which prefer this kinase conformation. More generally, these results provide additional insight into targeted kinase inhibitor design and how conformational preferences of inhibitors may impact selectivity and potency.

Predicting Mutation-Induced Allosteric Changes in Structures and Conformational Ensembles of the ABL Kinase Using AlphaFold2 Adaptations with Alanine Sequence Scanning

Despite the success of AlphaFold2 approaches in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and have been challenged to accurately capture the effects of single point mutations that induced significant structural changes. We examined several implementations of AlphaFold2 methods to predict conformational ensembles for state-switching mutants of the ABL kinase. The results revealed that a combination of randomized alanine sequence masking with shallow multiple sequence alignment subsampling can significantly expand the conformational diversity of the predicted structural ensembles and capture shifts in populations of the active and inactive ABL states. Consistent with the NMR experiments, the predicted conformational ensembles for M309L/L320I and M309L/H415P ABL mutants that perturb the regulatory spine networks featured the increased population of the fully closed inactive state. The proposed adaptation of AlphaFold can reproduce the experimentally observed mutation-induced redistributions in the relative populations of the active and inactive ABL states and capture the effects of regulatory mutations on allosteric structural rearrangements of the kinase domain. The ensemble-based network analysis complemented AlphaFold predictions by revealing allosteric hotspots that correspond to state-switching mutational sites which may explain the global effect of regulatory mutations on structural changes between the ABL states. This study suggested that attention-based learning of long-range dependencies between sequence positions in homologous folds and deciphering patterns of allosteric interactions may further augment the predictive abilities of AlphaFold methods for modeling of alternative protein sates, conformational ensembles and mutation-induced structural transformations.

BCR::ABL1 kinase N-lobe mutants confer moderate to high degrees of resistance to asciminib

ABSTRACT

Secondary kinase domain mutations in BCR::ABL1 represent the most common cause of resistance to tyrosine kinase inhibitor (TKI) therapy in patients with chronic myeloid leukemia. The first 5 approved BCR::ABL1 TKIs target the adenosine triphosphate (ATP)-binding pocket. Mutations confer resistance to these ATP-competitive TKIs and those approved for other malignancies by decreasing TKI affinity and/or increasing ATP affinity. Asciminib, the first highly active allosteric TKI approved for any malignancy, targets an allosteric regulatory pocket in the BCR::ABL1 kinase C-lobe. As a non-ATP-competitive inhibitor, the activity of asciminib is predicted to be impervious to increases in ATP affinity. Here, we report several known mutations that confer resistance to ATP-competitive TKIs in the BCR::ABL1 kinase N-lobe that are distant from the asciminib binding pocket yet unexpectedly confer in vitro resistance to asciminib. Among these is BCR::ABL1 M244V, which confers clinical resistance even to escalated asciminib doses. We demonstrate that BCR::ABL1 M244V does not impair asciminib binding, thereby invoking a novel mechanism of resistance. Molecular dynamic simulations of the M244V substitution implicate stabilization of an active kinase conformation through impact on the α-C helix as a mechanism of resistance. These N-lobe mutations may compromise the clinical activity of ongoing combination studies of asciminib with ATP-competitive TKIs.

Integration of a Randomized Sequence Scanning Approach in AlphaFold2 and Local Frustration Profiling of Conformational States Enable Interpretable Atomistic Characterization of Conformational Ensembles and Detection of Hidden Allosteric States in the ABL1 Protein Kinase

Despite the success of AlphaFold methods in predicting single protein structures, these methods showed intrinsic limitations in the characterization of multiple functional conformations of allosteric proteins. The recent NMR-based structural determination of the unbound ABL kinase in the active state and discovery of the inactive low-populated functional conformations that are unique for ABL kinase present an ideal challenge for the AlphaFold2 approaches. In the current study, we employ several adaptations of the AlphaFold2 methodology to predict protein conformational ensembles and allosteric states of the ABL kinase including randomized alanine sequence scanning combined with the multiple sequence alignment subsampling proposed in this study. We show that the proposed new AlphaFold2 adaptation combined with local frustration profiling of conformational states enables accurate prediction of the protein kinase structures and conformational ensembles, also offering a robust approach for interpretable characterization of the AlphaFold2 predictions and detection of hidden allosteric states. We found that the large high frustration residue clusters are uniquely characteristic of the low-populated, fully inactive ABL form and can define energetically frustrated cracking sites of conformational transitions, presenting difficult targets for AlphaFold2. The results of this study uncovered previously unappreciated fundamental connections between local frustration profiles of the functional allosteric states and the ability of AlphaFold2 methods to predict protein structural ensembles of the active and inactive states. This study showed that integration of the randomized sequence scanning adaptation of AlphaFold2 with a robust landscape-based analysis allows for interpretable atomistic predictions and characterization of protein conformational ensembles, providing a physical basis for the successes and limitations of current AlphaFold2 methods in detecting functional allosteric states that play a significant role in protein kinase regulation.

Prediction of Conformational Ensembles and Structural Effects of State-Switching Allosteric Mutants in the Protein Kinases Using Comparative Analysis of AlphaFold2 Adaptations with Sequence Masking and Shallow Subsampling

Despite the success of AlphaFold2 approaches in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and have been challenged to accurately capture of the effects of single point mutations that induced significant structural changes. We systematically examined several implementations of AlphaFold2 methods to predict conformational ensembles for state-switching mutants of the ABL kinase. The results revealed that a combination of randomized alanine sequence masking with shallow multiple sequence alignment subsampling can significantly expand the conformational diversity of the predicted structural ensembles and capture shifts in populations of the active and inactive ABL states. Consistent with the NMR experiments, the predicted conformational ensembles for M309L/L320I and M309L/H415P ABL mutants that perturb the regulatory spine networks featured the increased population of the fully closed inactive state. On the other hand, the predicted conformational ensembles for the G269E/M309L/T334I and M309L/L320I/T334I triple ABL mutants that share activating T334I gate-keeper substitution are dominated by the active ABL form. The proposed adaptation of AlphaFold can reproduce the experimentally observed mutation-induced redistributions in the relative populations of the active and inactive ABL states and capture the effects of regulatory mutations on allosteric structural rearrangements of the kinase domain. The ensemble-based network analysis complemented AlphaFold predictions by revealing allosteric mediating centers that often directly correspond to state-switching mutational sites or reside in their immediate local structural proximity, which may explain the global effect of regulatory mutations on structural changes between the ABL states. This study suggested that attention-based learning of long-range dependencies between sequence positions in homologous folds and deciphering patterns of allosteric interactions may further augment the predictive abilities of AlphaFold methods for modeling of alternative protein sates, conformational ensembles and mutation-induced structural transformations.

The molecular basis of Abelson kinase regulation by its αI-helix

Abelson tyrosine kinase (Abl) is regulated by the arrangement of its regulatory core, consisting sequentially of the SH3, SH2, and kinase (KD) domains, where an assembled or disassembled core corresponds to low or high kinase activity, respectively. It was recently established that binding of type II ATP site inhibitors, such as imatinib, generates a force from the KD N-lobe onto the SH3 domain and in consequence disassembles the core. Here, we demonstrate that the C-terminal αI-helix exerts an additional force toward the SH2 domain, which correlates both with kinase activity and type II inhibitor-induced disassembly. The αI-helix mutation E528K, which is responsible for the ABL1 malformation syndrome, strongly activates Abl by breaking a salt bridge with the KD C-lobe and thereby increasing the force onto the SH2 domain. In contrast, the allosteric inhibitor asciminib strongly reduces Abl's activity by fixating the αI-helix and reducing the force onto the SH2 domain. These observations are explained by a simple mechanical model of Abl activation involving forces from the KD N-lobe and the αI-helix onto the KD/SH2SH3 interface.

FLEXR GUI: a graphical user interface for multi-conformer modeling of proteins

Proteins are well known 'shapeshifters' which change conformation to function. In crystallography, multiple conformational states are often present within the crystal and the resulting electron-density map. Yet, explicitly incorporating alternative states into models to disentangle multi-conformer ensembles is challenging. We previously reported the tool FLEXR, which, within a few minutes, automatically separates conformational signal from noise and builds the corresponding, often missing, structural features into a multi-conformer model. To make the method widely accessible for routine multi-conformer building as part of the computational toolkit for macromolecular crystallography, we present a graphical user interface (GUI) for FLEXR, designed as a plugin for Coot 1. The GUI implementation seamlessly connects FLEXR models with the existing suite of validation and modeling tools available in Coot. We envision that FLEXR will aid crystallographers by increasing access to a multi-conformer modeling method that will ultimately lead to a better representation of protein conformational heterogeneity in the Protein Data Bank. In turn, deeper insights into the protein conformational landscape may inform biology or provide new opportunities for ligand design. The code is open source and freely available on GitHub at https://github.com/TheFischerLab/FLEXR-GUI.

Interpretable Atomistic Prediction and Functional Analysis of Conformational Ensembles and Allosteric States in Protein Kinases Using AlphaFold2 Adaptation with Randomized Sequence Scanning and Local Frustration Profiling

The groundbreaking achievements of AlphaFold2 (AF2) approaches in protein structure modeling marked a transformative era in structural biology. Despite the success of AF2 tools in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and fold-switching systems. The recent NMR-based structural determination of the unbound ABL kinase in the active state and two inactive low-populated functional conformations that are unique for ABL kinase presents an ideal challenge for AF2 approaches. In the current study we employ several implementations of AF2 methods to predict protein conformational ensembles and allosteric states of the ABL kinase including (a) multiple sequence alignments (MSA) subsampling approach; (b) SPEACH_AF approach in which alanine scanning is performed on generated MSAs; and (c) introduced in this study randomized full sequence mutational scanning for manipulation of sequence variations combined with the MSA subsampling. We show that the proposed AF2 adaptation combined with local frustration mapping of conformational states enable accurate prediction of the ABL active and intermediate structures and conformational ensembles, also offering a robust approach for interpretable characterization of the AF2 predictions and limitations in detecting hidden allosteric states. We found that the large high frustration residue clusters are uniquely characteristic of the low-populated, fully inactive ABL form and can define energetically frustrated cracking sites of conformational transitions, presenting difficult targets for AF2 methods. This study uncovered previously unappreciated, fundamental connections between distinct patterns of local frustration in functional kinase states and AF2 successes/limitations in detecting low-populated frustrated conformations, providing a better understanding of benefits and limitations of current AF2-based adaptations in modeling of conformational ensembles.