Hierarchical Modeling of Activation Mechanisms in the ABL and EGFR Kinase Domains: Thermodynamic and Mechanistic Catalysts of Kinase Activation by Cancer Mutations

Structural and functional studies of the ABL and EGFR kinase domains have recently suggested a common mechanism of activation by cancer-causing mutations. However, dynamics and mechanistic aspects of kinase activation by cancer mutations that stimulate conformational transitions and thermodynamic stabilization of the constitutively active kinase form remain elusive. We present a large-scale computational investigation of activation mechanisms in the ABL and EGFR kinase domains by a panel of clinically important cancer mutants ABL-T315I, ABL-L387M, EGFR-T790M, and EGFR-L858R. We have also simulated the activating effect of the gatekeeper mutation on conformational dynamics and allosteric interactions in functional states of the ABL-SH2-SH3 regulatory complexes. A comprehensive analysis was conducted using a hierarchy of computational approaches that included homology modeling, molecular dynamics simulations, protein stability analysis, targeted molecular dynamics, and molecular docking. Collectively, the results of this study have revealed thermodynamic and mechanistic catalysts of kinase activation by major cancer-causing mutations in the ABL and EGFR kinase domains. By using multiple crystallographic states of ABL and EGFR, computer simulations have allowed one to map dynamics of conformational fluctuations and transitions in the normal (wild-type) and oncogenic kinase forms. A proposed multi-stage mechanistic model of activation involves a series of cooperative transitions between different conformational states, including assembly of the hydrophobic spine, the formation of the Src-like intermediate structure, and a cooperative breakage and formation of characteristic salt bridges, which signify transition to the active kinase form. We suggest that molecular mechanisms of activation by cancer mutations could mimic the activation process of the normal kinase, yet exploiting conserved structural catalysts to accelerate a conformational transition and the enhanced stabilization of the active kinase form. The results of this study reconcile current experimental data with insights from theoretical approaches, pointing to general mechanistic aspects of activating transitions in protein kinases

Conformational Determinants of Phosphotyrosine Peptides Complexed with the Src SH2 Domain

The inhibition of specific SH2 domain mediated protein-protein interactions as an effective chemotherapeutic approach in the treatment of diseases remains a challenge. That different conformations of peptide-ligands are preferred by different SH2 domains is an underappreciated observation from the structural analysis of phosphotyrosine peptide binding to SH2 domains that may aid in future drug design. To explore the nature of ligand binding, we use simulated annealing (SA) to sample the conformational space of phosphotyrosine-containing peptides complexed with the Src SH2 domain. While in good agreement with the crystallographic and NMR studies of high-affinity phosphopeptide-SH2 domain complexes, the results suggest that the structural basis for phopsphopeptide- Src SH2 interactions is more complex than the “two-pronged plug two-hole socket” model. A systematic study of peptides of type pYEEX, where pY is phosphotyrosine and X is a hydrophobic residue, indicates that these peptides can assume two conformations, one extended and one helical, representing the balance between the interaction of residue X with the hydrophobic hole on the surface of the Src SH2 domain, and its contribution to the inherent tendency of the two glutamic acids to form an α-helix. In contrast, a β-turn conformation, almost identical to that observed in the crystal structure of pYVNV bound to the Grb2 SH2 domain, predominates for pYXNX peptides, even in the presence of isoleucine at the third position. While peptide binding affinities, as measured by fluorescence polarization, correlate with the relative proportion of extended peptide conformation, these results suggest a model where all three residues C-terminal to the phosphotyrosine determine the conformation of the bound phosphopeptide. The information obtained in this work can be used in the design of specific SH2 domain inhibitors

Structural and Functional Characterization of Mature Forms of Metalloprotease E495 from Arctic Sea-Ice Bacterium Pseudoalteromonas sp. SM495

E495 is the most abundant protease secreted by the Arctic sea-ice bacterium Pseudoalteromonas sp. SM495. As a thermolysin family metalloprotease, E495 was found to have multiple active forms in the culture of strain SM495. E495-M (containing only the catalytic domain) and E495-M-C1 (containing the catalytic domain and one PPC domain) were two stable mature forms, and E495-M-C1-C2 (containing the catalytic domain and two PPC domains) might be an intermediate. Compared to E495-M, E495-M-C1 had similar affinity and catalytic efficiency to oligopeptides, but higher affinity and catalytic efficiency to proteins. The PPC domains from E495 were expressed as GST-fused proteins. Both of the recombinant PPC domains were shown to have binding ability to proteins C-phycocyanin and casein, and domain PPC1 had higher affinity to C-phycocyanin than domain PPC2. These results indicated that the domain PPC1 in E495-M-C1 could be helpful in binding protein substrate, and therefore, improving the catalytic efficiency. Site-directed mutagenesis on the PPC domains showed that the conserved polar and aromatic residues, D26, D28, Y30, Y/W65, in the PPC domains played key roles in protein binding. Our study may shed light on the mechanism of organic nitrogen degradation in the Arctic sea ice

Allosteric Interactions between the Myristate- and ATP-Site of the Abl Kinase

Abl kinase inhibitors targeting the ATP binding pocket are currently employed as potent anti-leukemogenic agents but drug resistance has become a significant clinical limitation. Recently, a compound that binds to the myristate pocket of Abl (GNF-5) was shown to act cooperatively with nilotinib, an ATP-competitive inhibitor to target the recalcitrant “T315I” gatekeeper mutant of Bcr-Abl. To uncover an explanation for how drug binding at a distance from the kinase active site could lead to inhibition and how inhibitors could combine their effects, hydrogen exchange mass spectrometry (HX MS) was employed to monitor conformational effects in the presence of both dasatinib, a clinically approved ATP-site inhibitor, and GNF-5. While dasatinib binding to wild type Abl clearly influenced Abl conformation, no binding was detected between dasatinib and T315I. GNF-5, however, elicited the same conformational changes in both wild type and T315I, including changes to dynamics within the ATP site located approximately 25 Å from the site of GNF-5 interaction. Simultaneous binding of dasatinib and GNF-5 to T315I caused conformational and/or dynamics changes in Abl such that effects of dasatinib on T315I were the same as when it bound to wild type Abl. These results provide strong biophysical evidence that allosteric interactions play a role in Abl kinase downregulation and that targeting sites outside the ATP binding site can provide an important pharmacological tool to overcome mutations that cause resistance to ATP-competitive inhibitors

The Energy Landscape Analysis of Cancer Mutations in Protein Kinases

The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling may involve a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. The presented study has demonstrated that activation cancer mutations may affect the thermodynamic equilibrium between kinase states by allosterically altering the distribution of locally frustrated sites and increasing the local frustration in the inactive form, while eliminating locally frustrated sites and restoring structural rigidity of the active form. The energy landsape analysis of protein kinases and the proposed role of locally frustrated sites in activation mechanisms may have useful implications for bioinformatics-based screening and detection of functional sites critical for allosteric regulation in complex biomolecular systems

Probing Molecular Mechanisms of the Hsp90 Chaperone: Biophysical Modeling Identifies Key Regulators of Functional Dynamics

Deciphering functional mechanisms of the Hsp90 chaperone machinery is an important objective in cancer biology aiming to facilitate discovery of targeted anti-cancer therapies. Despite significant advances in understanding structure and function of molecular chaperones, organizing molecular principles that control the relationship between conformational diversity and functional mechanisms of the Hsp90 activity lack a sufficient quantitative characterization. We combined molecular dynamics simulations, principal component analysis, the energy landscape model and structure-functional analysis of Hsp90 regulatory interactions to systematically investigate functional dynamics of the molecular chaperone. This approach has identified a network of conserved regions common to the Hsp90 chaperones that could play a universal role in coordinating functional dynamics, principal collective motions and allosteric signaling of Hsp90. We have found that these functional motifs may be utilized by the molecular chaperone machinery to act collectively as central regulators of Hsp90 dynamics and activity, including the inter-domain communications, control of ATP hydrolysis, and protein client binding. These findings have provided support to a long-standing assertion that allosteric regulation and catalysis may have emerged via common evolutionary routes. The interaction networks regulating functional motions of Hsp90 may be determined by the inherent structural architecture of the molecular chaperone. At the same time, the thermodynamics-based “conformational selection” of functional states is likely to be activated based on the nature of the binding partner. This mechanistic model of Hsp90 dynamics and function is consistent with the notion that allosteric networks orchestrating cooperative protein motions can be formed by evolutionary conserved and sparsely connected residue clusters. Hence, allosteric signaling through a small network of distantly connected residue clusters may be a rather general functional requirement encoded across molecular chaperones. The obtained insights may be useful in guiding discovery of allosteric Hsp90 inhibitors targeting protein interfaces with co-chaperones and protein binding clients

Mutation D816V Alters the Internal Structure and Dynamics of c-KIT Receptor Cytoplasmic Region: Implications for Dimerization and Activation Mechanisms

The type III receptor tyrosine kinase (RTK) KIT plays a crucial role in the transmission of cellular signals through phosphorylation events that are associated with a switching of the protein conformation between inactive and active states. D816V KIT mutation is associated with various pathologies including mastocytosis and cancers. D816V-mutated KIT is constitutively active, and resistant to treatment with the anti-cancer drug Imatinib. To elucidate the activating molecular mechanism of this mutation, we applied a multi-approach procedure combining molecular dynamics (MD) simulations, normal modes analysis (NMA) and binding site prediction. Multiple 50-ns MD simulations of wild-type KIT and its mutant D816V were recorded using the inactive auto-inhibited structure of the protein, characteristic of type III RTKs. Computed free energy differences enabled us to quantify the impact of D816V on protein stability in the inactive state. We evidenced a local structural alteration of the activation loop (A-loop) upon mutation, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. A thorough normal mode analysis of several MD conformations led to a plausible molecular rationale to propose that JMR is able to depart its auto-inhibitory position more easily in the mutant than in wild-type KIT and is thus able to promote kinase mutant dimerization without the need for extra-cellular ligand binding. Pocket detection at the surface of NMA-displaced conformations finally revealed that detachment of JMR from the kinase domain in the mutant was sufficient to open an access to the catalytic and substrate binding sites

Sequence and Structure Signatures of Cancer Mutation Hotspots in Protein Kinases

Protein kinases are the most common protein domains implicated in cancer, where somatically acquired mutations are known to be functionally linked to a variety of cancers. Resequencing studies of protein kinase coding regions have emphasized the importance of sequence and structure determinants of cancer-causing kinase mutations in understanding of the mutation-dependent activation process. We have developed an integrated bioinformatics resource, which consolidated and mapped all currently available information on genetic modifications in protein kinase genes with sequence, structure and functional data. The integration of diverse data types provided a convenient framework for kinome-wide study of sequence-based and structure-based signatures of cancer mutations. The database-driven analysis has revealed a differential enrichment of SNPs categories in functional regions of the kinase domain, demonstrating that a significant number of cancer mutations could fall at structurally equivalent positions (mutational hotspots) within the catalytic core. We have also found that structurally conserved mutational hotspots can be shared by multiple kinase genes and are often enriched by cancer driver mutations with high oncogenic activity. Structural modeling and energetic analysis of the mutational hotspots have suggested a common molecular mechanism of kinase activation by cancer mutations, and have allowed to reconcile the experimental data. According to a proposed mechanism, structural effect of kinase mutations with a high oncogenic potential may manifest in a significant destabilization of the autoinhibited kinase form, which is likely to drive tumorigenesis at some level. Structure-based functional annotation and prediction of cancer mutation effects in protein kinases can facilitate an understanding of the mutation-dependent activation process and inform experimental studies exploring molecular pathology of tumorigenesis

3D-QSAR and Molecular Docking Studies of 4-Anilinoquinazoline Derivatives: a Rational Approach to Anticancer Drug Design

The present article is an attempt to formulate the three-dimensional quantitative structure–activity relationship (3D-QSAR) modeling of 4-anilinoquinazoline derivatives having promising anticancer activities inhibiting epidermal growth factor (EGFR) kinase. Molecular field analysis was applied for the generation of steric and electrostatic descriptors based on aligned structures. Partial least-squares (PLS) methodwas applied for QSAR model development considering training and test set approaches. The PLSmodels showed some interesting results in terms of internal and external predictability against EGFR kinase inhibition for such type of anilinoquinazoline derivatives. Steric and electrostatic field effects are discussed in the light of contribution plot generated. Finally, molecular docking analysis was carried out to better understand of the interactions between EGFR target and inhibitors in this series. Hydrophobic and hydrogenbond interactions lead to identification of active binding sites of EGFR protein in the docked complex