BCR::ABL1 Kinase N-lobe Mutants Confer Moderate to High Degrees of Resistance to Asciminib.

Secondary kinase domain mutations in BCR::ABL1 represent the most common cause of resistance to tyrosine kinase inhibitor (TKI) therapy in chronic myeloid leukemia patients. The first five approved BCR::ABL1 TKIs target the 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 dynamics 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.

Structural mechanism of a drug-binding process involving a large conformational change of the protein target.

Proteins often undergo large conformational changes when binding small molecules, but atomic-level descriptions of such events have been elusive. Here, we report unguided molecular dynamics simulations of Abl kinase binding to the cancer drug imatinib. In the simulations, imatinib first selectively engages Abl kinase in its autoinhibitory conformation. Consistent with inferences drawn from previous experimental studies, imatinib then induces a large conformational change of the protein to reach a bound complex that closely resembles published crystal structures. Moreover, the simulations reveal a surprising local structural instability in the C-terminal lobe of Abl kinase during binding. The unstable region includes a number of residues that, when mutated, confer imatinib resistance by an unknown mechanism. Based on the simulations, NMR spectra, hydrogen-deuterium exchange measurements, and thermostability measurements and estimates, we suggest that these mutations confer imatinib resistance by exacerbating structural instability in the C-terminal lobe, rendering the imatinib-bound state energetically unfavorable.

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations.

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often "cryptic" without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein-ligand systems.

Targetable HER3 functions driving tumorigenic signaling in HER2-amplified cancers.

Effective inactivation of the HER2-HER3 tumor driver has remained elusive because of the challenging attributes of the pseudokinase HER3. We report a structure-function study of constitutive HER2-HER3 signaling to identify opportunities for targeting. The allosteric activation of the HER2 kinase domain (KD) by the HER3 KD is required for tumorigenic signaling and can potentially be targeted by allosteric inhibitors. ATP binding within the catalytically inactive HER3 KD provides structural rigidity that is important for signaling, but this is mimicked, not opposed, by small molecule ATP analogs, reported here in a bosutinib-bound crystal structure. Mutational disruption of ATP binding and molecular dynamics simulation of the apo KD of HER3 identify a conformational coupling of the ATP pocket with a hydrophobic AP-2 pocket, analogous to EGFR, that is critical for tumorigenic signaling and feasible for targeting. The value of these potential target sites is confirmed in tumor growth assays using gene replacement techniques.

A structural model of a Ras-Raf signalosome.

The protein K-Ras functions as a molecular switch in signaling pathways regulating cell growth. In the human mitogen-activated protein kinase (MAPK) pathway, which is implicated in many cancers, multiple K-Ras proteins are thought to assemble at the cell membrane with Ras effector proteins from the Raf family. Here we propose an atomistic structural model for such an assembly. Our starting point was an asymmetric guanosine triphosphate-mediated K-Ras dimer model, which we generated using unbiased molecular dynamics simulations and verified with mutagenesis experiments. Adding further K-Ras monomers in a head-to-tail fashion led to a compact helical assembly, a model we validated using electron microscopy and cell-based experiments. This assembly stabilizes K-Ras in its active state and presents composite interfaces to facilitate Raf binding. Guided by existing experimental data, we then positioned C-Raf, the downstream kinase MEK1 and accessory proteins (Galectin-3 and 14-3-3σ) on and around the helical assembly. The resulting Ras-Raf signalosome model offers an explanation for a large body of data on MAPK signaling.

Structural Basis of AZD9291 Selectivity for EGFR T790M.

AZD9291 (Osimertinib) is highly effective in treating EGFR-mutated non-small-cell lung cancers (NSCLCs) with T790M-mediated drug resistance. Despite the remarkable success of AZD9291, its binding pose with EGFR T790M remains unclear. Here, we report unbiased, atomic-level molecular dynamics (MD) simulations in which spontaneous association of AZD9291 with EGFR kinases having WT and T790M mutant gatekeepers was observed. Simulation-generated structural models suggest that the binding pose of AZD9291 with T790M differs from its binding pose with the WT, and that AZD9291 interacts extensively with the gatekeeper residue (Met 790) in T790M but not with Thr 790 in the WT, which explains why AZD9291 binds T790M with higher affinity. The MD simulation-generated models were confirmed by experimentally determined EGFR/T790M complex crystal structures. This work may facilitate the rational design of drugs that can overcome resistance mutations to AZD9291, and more generally it suggests the potential of using unbiased MD simulation to elucidate small-molecule binding poses.

The architecture of EGFR's basal complexes reveals autoinhibition mechanisms in dimers and oligomers.

Our current understanding of epidermal growth factor receptor (EGFR) autoinhibition is based on X-ray structural data of monomer and dimer receptor fragments and does not explain how mutations achieve ligand-independent phosphorylation. Using a repertoire of imaging technologies and simulations we reveal an extracellular head-to-head interaction through which ligand-free receptor polymer chains of various lengths assemble. The architecture of the head-to-head interaction prevents kinase-mediated dimerisation. The latter, afforded by mutation or intracellular treatments, splits the autoinhibited head-to-head polymers to form stalk-to-stalk flexible non-extended dimers structurally coupled across the plasma membrane to active asymmetric tyrosine kinase dimers, and extended dimers coupled to inactive symmetric kinase dimers. Contrary to the previously proposed main autoinhibitory function of the inactive symmetric kinase dimer, our data suggest that only dysregulated species bear populations of symmetric and asymmetric kinase dimers that coexist in equilibrium at the plasma membrane under the modulation of the C-terminal domain.

EGFR oligomerization organizes kinase-active dimers into competent signalling platforms.

Epidermal growth factor receptor (EGFR) signalling is activated by ligand-induced receptor dimerization. Notably, ligand binding also induces EGFR oligomerization, but the structures and functions of the oligomers are poorly understood. Here, we use fluorophore localization imaging with photobleaching to probe the structure of EGFR oligomers. We find that at physiological epidermal growth factor (EGF) concentrations, EGFR assembles into oligomers, as indicated by pairwise distances of receptor-bound fluorophore-conjugated EGF ligands. The pairwise ligand distances correspond well with the predictions of our structural model of the oligomers constructed from molecular dynamics simulations. The model suggests that oligomerization is mediated extracellularly by unoccupied ligand-binding sites and that oligomerization organizes kinase-active dimers in ways optimal for auto-phosphorylation in trans between neighbouring dimers. We argue that ligand-induced oligomerization is essential to the regulation of EGFR signalling.

Structural analysis of the EGFR/HER3 heterodimer reveals the molecular basis for activating HER3 mutations.

The human epidermal growth factor receptor (HER) tyrosine kinases homo- and heterodimerize to activate downstream signaling pathways. HER3 is a catalytically impaired member of the HER family that contributes to the development of several human malignancies and is mutated in a subset of cancers. HER3 signaling depends on heterodimerization with a catalytically active partner, in particular epidermal growth factor receptor (EGFR) (the founding family member, also known as HER1) or HER2. The activity of homodimeric complexes of catalytically active HER family members depends on allosteric activation between the two kinase domains. To determine the structural basis for HER3 signaling through heterodimerization with a catalytically active HER family member, we solved the crystal structure of the heterodimeric complex formed by the isolated kinase domains of EGFR and HER3. The structure showed HER3 as an allosteric activator of EGFR and revealed a conserved role of the allosteric mechanism in activation of HER family members through heterodimerization. To understand the effects of cancer-associated HER3 mutations at the molecular level, we solved the structures of two kinase domains of HER3 mutants, each in a heterodimeric complex with the kinase domain of EGFR. These structures, combined with biochemical analysis and molecular dynamics simulations, indicated that the cancer-associated HER3 mutations enhanced the allosteric activator function of HER3 by redesigning local interactions at the dimerization interface.

How IGF-1 activates its receptor.

The type I insulin-like growth factor receptor (IGF1R) is involved in growth and survival of normal and neoplastic cells. A ligand-dependent conformational change is thought to regulate IGF1R activity, but the nature of this change is unclear. We point out an underappreciated dimer in the crystal structure of the related Insulin Receptor (IR) with Insulin bound that allows direct comparison with unliganded IR and suggests a mechanism by which ligand regulates IR/IGF1R activity. We test this mechanism in a series of biochemical and biophysical assays and find the IGF1R ectodomain maintains an autoinhibited state in which the TMs are held apart. Ligand binding releases this constraint, allowing TM association and unleashing an intrinsic propensity of the intracellular regions to autophosphorylate. Enzymatic studies of full-length and kinase-containing fragments show phosphorylated IGF1R is fully active independent of ligand and the extracellular-TM regions. The key step triggered by ligand binding is thus autophosphorylation.

Membrane interaction of bound ligands contributes to the negative binding cooperativity of the EGF receptor.

The epidermal growth factor receptor (EGFR) plays a key role in regulating cell proliferation, migration, and differentiation, and aberrant EGFR signaling is implicated in a variety of cancers. EGFR signaling is triggered by extracellular ligand binding, which promotes EGFR dimerization and activation. Ligand-binding measurements are consistent with a negatively cooperative model in which the ligand-binding affinity at either binding site in an EGFR dimer is weaker when the other site is occupied by a ligand. This cooperativity is widely believed to be central to the effects of ligand concentration on EGFR-mediated intracellular signaling. Although the extracellular portion of the human EGFR dimer has been resolved crystallographically, the crystal structures do not reveal the structural origin of this negative cooperativity, which has remained unclear. Here we report the results of molecular dynamics simulations suggesting that asymmetrical interactions of the two binding sites with the membrane may be responsible (perhaps along with other factors) for this negative cooperativity. In particular, in our simulations the extracellular domains of an EGFR dimer spontaneously lay down on the membrane in an orientation in which favorable membrane contacts were made with one of the bound ligands, but could not be made with the other. Similar interactions were observed when EGFR was glycosylated, as it is in vivo.

Molecular basis for pseudokinase-dependent autoinhibition of JAK2 tyrosine kinase.

Janus kinase-2 (JAK2) mediates signaling by various cytokines, including erythropoietin and growth hormone. JAK2 possesses tandem pseudokinase and tyrosine-kinase domains. Mutations in the pseudokinase domain are causally linked to myeloproliferative neoplasms (MPNs) in humans. The structure of the JAK2 tandem kinase domains is unknown, and therefore the molecular bases for pseudokinase-mediated autoinhibition and pathogenic activation remain obscure. Using molecular dynamics simulations of protein-protein docking, we produced a structural model for the autoinhibitory interaction between the JAK2 pseudokinase and kinase domains. A striking feature of our model, which is supported by mutagenesis experiments, is that nearly all of the disease mutations map to the domain interface. The simulations indicate that the kinase domain is stabilized in an inactive state by the pseudokinase domain, and they offer a molecular rationale for the hyperactivity of V617F, the predominant JAK2 MPN mutation.

Her2 activation mechanism reflects evolutionary preservation of asymmetric ectodomain dimers in the human EGFR family.

The receptor tyrosine kinase Her2, an intensely pursued drug target, differs from other members of the EGFR family in that it does not bind EGF-like ligands, relying instead on heterodimerization with other (ligand-bound) EGFR-family receptors for activation. The structural basis for Her2 heterodimerization, however, remains poorly understood. The unexpected recent finding of asymmetric ectodomain dimer structures of Drosophila EGFR (dEGFR) suggests a possible structural basis for Her2 heterodimerization, but all available structures for dimers of human EGFR family ectodomains are symmetric. Here, we report results from long-timescale molecular dynamics simulations indicating that a single ligand is necessary and sufficient to stabilize the ectodomain interface of Her2 heterodimers, which assume an asymmetric conformation similar to that of dEGFR dimers. This structural parallelism suggests a dimerization mechanism that has been conserved in the evolution of the EGFR family from Drosophila to human. DOI:http://dx.doi.org/10.7554/eLife.00708.001.

Transitions to catalytically inactive conformations in EGFR kinase.

The epidermal growth factor receptor (EGFR) is a key protein in cellular signaling, and its kinase domain (EGFR kinase) is an intensely pursued target of small-molecule drugs. Although both catalytically active and inactive conformations of EGFR kinase have been resolved crystallographically, experimental characterization of the transitions between these conformations remains difficult. Using unbiased, all-atom molecular dynamics simulations, we observed EGFR kinase spontaneously transition from the active to the so-called "Src-like inactive" conformation by way of two sets of intermediate conformations: One corresponds to a previously identified locally disordered state and the other to previously undescribed "extended" conformations, marked by the opening of the ATP-binding site between the two lobes of the kinase domain. We also simulated the protonation-dependent transition of EGFR kinase to another ["Asp-Phe-Gly-out" ("DFG-out")] inactive conformation and observed similar intermediate conformations. A key element observed in the simulated transitions is local unfolding, or "cracking," which supports a prediction of energy landscape theory. We used hydrogen-deuterium (H/D) exchange measurements to corroborate our simulations and found that the simulated intermediate conformations correlate better with the H/D exchange data than existing active or inactive EGFR kinase crystal structures. The intermediate conformations revealed by our simulations of the transition process differ significantly from the existing crystal structures and may provide unique possibilities for structure-based drug discovery.

Conformational coupling across the plasma membrane in activation of the EGF receptor.

How the epidermal growth factor receptor (EGFR) activates is incompletely understood. The intracellular portion of the receptor is intrinsically active in solution, and to study its regulation, we measured autophosphorylation as a function of EGFR surface density in cells. Without EGF, intact EGFR escapes inhibition only at high surface densities. Although the transmembrane helix and the intracellular module together suffice for constitutive activity even at low densities, the intracellular module is inactivated when tethered on its own to the plasma membrane, and fluorescence cross-correlation shows that it fails to dimerize. NMR and functional data indicate that activation requires an N-terminal interaction between the transmembrane helices, which promotes an antiparallel interaction between juxtamembrane segments and release of inhibition by the membrane. We conclude that EGF binding removes steric constraints in the extracellular module, promoting activation through N-terminal association of the transmembrane helices.

Architecture and membrane interactions of the EGF receptor.

Dimerization-driven activation of the intracellular kinase domains of the epidermal growth factor receptor (EGFR) upon extracellular ligand binding is crucial to cellular pathways regulating proliferation, migration, and differentiation. Inactive EGFR can exist as both monomers and dimers, suggesting that the mechanism regulating EGFR activity may be subtle. The membrane itself may play a role but creates substantial difficulties for structural studies. Our molecular dynamics simulations of membrane-embedded EGFR suggest that, in ligand-bound dimers, the extracellular domains assume conformations favoring dimerization of the transmembrane helices near their N termini, dimerization of the juxtamembrane segments, and formation of asymmetric (active) kinase dimers. In ligand-free dimers, by holding apart the N termini of the transmembrane helices, the extracellular domains instead favor C-terminal dimerization of the transmembrane helices, juxtamembrane segment dissociation and membrane burial, and formation of symmetric (inactive) kinase dimers. Electrostatic interactions of EGFR's intracellular module with the membrane are critical in maintaining this coupling.

Crystal structures of the JAK2 pseudokinase domain and the pathogenic mutant V617F.

The protein tyrosine kinase JAK2 mediates signaling through numerous cytokine receptors. JAK2 possesses a pseudokinase domain (JH2) and a tyrosine kinase domain (JH1). Through unknown mechanisms, JH2 regulates the catalytic activity of JH1, and hyperactivating mutations in the JH2 region of human JAK2 cause myeloproliferative neoplasms (MPNs). We showed previously that JAK2 JH2 is, in fact, catalytically active. Here we present crystal structures of human JAK2 JH2, including both wild type and the most prevalent MPN mutant, V617F. The structures reveal that JH2 adopts the fold of a prototypical protein kinase but binds Mg-ATP noncanonically. The structural and biochemical data indicate that the V617F mutation rigidifies α-helix C in the N lobe of JH2, facilitating trans-phosphorylation of JH1. The crystal structures of JH2 afford new opportunities for the design of novel JAK2 therapeutics targeting MPNs.

Regulation of RAS oncogenicity by acetylation.

Members of the RAS small GTPase family regulate cellular responses to extracellular stimuli by mediating the flux through downstream signal transduction cascades. RAS activity is strongly dependent on its subcellular localization and its nucleotide-binding status, both of which are modulated by posttranslational modification. We have determined that RAS is posttranslationally acetylated on lysine 104. Molecular dynamics simulations suggested that this modification affects the conformational stability of the Switch II domain, which is critical for the ability of RAS to interact with guanine nucleotide exchange factors. Consistent with this model, an acetylation-mimetic mutation in K-RAS4B suppressed guanine nucleotide exchange factor-induced nucleotide exchange and inhibited in vitro transforming activity. These data suggest that lysine acetylation is a negative regulatory modification on RAS. Because mutations in RAS family members are extremely common in cancer, modulation of RAS acetylation may constitute a therapeutic approach.

Oncogenic mutations counteract intrinsic disorder in the EGFR kinase and promote receptor dimerization.

The mutation and overexpression of the epidermal growth factor receptor (EGFR) are associated with the development of a variety of cancers, making this prototypical dimerization-activated receptor tyrosine kinase a prominent target of cancer drugs. Using long-timescale molecular dynamics simulations, we find that the N lobe dimerization interface of the wild-type EGFR kinase domain is intrinsically disordered and that it becomes ordered only upon dimerization. Our simulations suggest, moreover, that some cancer-linked mutations distal to the dimerization interface, particularly the widespread L834R mutation (also referred to as L858R), facilitate EGFR dimerization by suppressing this local disorder. Corroborating these findings, our biophysical experiments and kinase enzymatic assays indicate that the L834R mutation causes abnormally high activity primarily by promoting EGFR dimerization rather than by allowing activation without dimerization. We also find that phosphorylation of EGFR kinase domain at Tyr845 may suppress the intrinsic disorder, suggesting a molecular mechanism for autonomous EGFR signaling.

How does a drug molecule find its target binding site?

Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, detailed experimental characterization of the process by which such binding occurs has proven challenging. We conducted relatively long, unguided molecular dynamics simulations in which a ligand (the cancer drug dasatinib or the kinase inhibitor PP1) was initially placed at a random location within a box that also contained a protein (Src kinase) to which that ligand was known to bind. In several of these simulations, the ligand correctly identified its target binding site, forming a complex virtually identical to the crystallographically determined bound structure. The simulated trajectories provide a continuous, atomic-level view of the entire binding process, revealing persistent and noteworthy intermediate conformations and shedding light on the role of water molecules. The technique we employed, which does not assume any prior knowledge of the binding site's location, may prove particularly useful in the development of allosteric inhibitors that target previously undiscovered binding sites.