Flexibility and charge asymmetry in the activation loop of Src tyrosine kinases

Decoupling the dynamic mechanism revealed by FGFR2 mutation-induced population shift

The fibroblast growth factor receptor 2 (FGFR2) is a key component in cellular signaling networks, and its dysfunctional activation has been implicated in various diseases including cancer and developmental disorders. Mutations at the activation loop (A-loop) have been suggested to trigger an increased basal kinase activity. However, the molecular mechanism underlying this highly dynamic process has not been fully understood due to the limitation of static structural information. Here, we conducted multiple, large-scale Gaussian accelerated molecular dynamics simulations of five (K659E, K659N, K659M, K659Q, and K659T) FGFR2 mutants at the A-loop, and comprehensively analyzed the dynamic molecular basis of FGFR2 activation. The results quantified the population shift of each system, revealing that all mutants had a higher proportion of active-like states. Using Markov state models, we extracted the representative structure of different conformational states and identified key residues related to the increased kinase activity. Furthermore, community network analysis showed enhanced information connections in the mutants, highlighting the long-range allosteric communication between the A-loop and the hinge region. Our findings may provide insights into the dynamic mechanism for FGFR2 dysfunctional activation and allosteric drug discovery.Communicated by Ramaswamy H. Sarma.

Oxidative stress promotes fibrosis in systemic sclerosis through stabilization of a kinase-phosphatase complex

Systemic sclerosis (SSc) is a fibrotic autoimmune disease characterized by pathogenic activation of fibroblasts enhanced by local oxidative stress. The tyrosine phosphatase PTP4A1 was identified as a critical promoter of TGF-β signaling in SSc. Oxidative stress is known to functionally inactivate tyrosine phosphatases. Here, we assessed whether oxidation of PTP4A1 modulates its profibrotic action and found that PTP4A1 forms a complex with the kinase SRC in scleroderma fibroblasts, but surprisingly, oxidative stress enhanced rather than reduced PTP4A1’s association with SRC and its profibrotic action. Through structural assessment of the oxo-PTP4A1-SRC complex, we unraveled an unexpected mechanism whereby oxidation of a tyrosine phosphatase promotes its function through modification of its protein complex. Considering the importance of oxidative stress in the pathogenesis of SSc and fibrosis, our findings suggest routes for leveraging PTP4A1 oxidation as a potential strategy for developing antifibrotic agents.

TYROSINE KINASES: COMPLEX MOLECULAR SYSTEMS CHALLENGING COMPUTATIONAL METHODOLOGIES

Classical molecular dynamics (MD) simulations based on atomic models play an increasingly important role in a wide range of applications in physics, biology, and chemistry. Nonetheless, generating genuine knowledge about biological systems using MD simulations remains challenging. Protein tyrosine kinases are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Due to the large conformational changes and long timescales involved in their function, these kinases present particularly challenging problems to modern computational and theoretical frameworks aimed at elucidating the dynamics of complex biomolecular systems. Markov state models have achieved limited success in tackling the broader conformational ensemble and biased methods are often employed to examine specific long timescale events. Recent advances in machine learning continue to push the limitations of current methodologies and provide notable improvements when integrated with the existing frameworks. A broad perspective is drawn from a critical review of recent studies.

Computational studies of anaplastic lymphoma kinase mutations reveal common mechanisms of oncogenic activation

High-risk tumors are genomically heterogeneous, harboring gene amplifications and mutations. The activation status of mutated proteins in cancer can profoundly impact disease progression, patient response, and drug sensitivity. Yet, outside of a few hotspot mutations, functional studies of clinically observed mutations are not commonly pursued. We report a combined experimental profiling and computational analysis of the effects of clinically observed and “test” mutations in the kinase domain of anaplastic lymphoma kinase (ALK), a known oncogenic driver in pediatric neuroblastoma. We find that the activation status of the mutated protein is a good indicator of the transforming ability in NIH 3T3 cells. We also report biophysical as well as data-driven models with predictive power to profile these mutant kinases in silico.

Kinases play important roles in diverse cellular processes, including signaling, differentiation, proliferation, and metabolism. They are frequently mutated in cancer and are the targets of a large number of specific inhibitors. Surveys of cancer genome atlases reveal that kinase domains, which consist of 300 amino acids, can harbor numerous (150 to 200) single-point mutations across different patients in the same disease. This preponderance of mutations—some activating, some silent—in a known target protein make clinical decisions for enrolling patients in drug trials challenging since the relevance of the target and its drug sensitivity often depend on the mutational status in a given patient. We show through computational studies using molecular dynamics (MD) as well as enhanced sampling simulations that the experimentally determined activation status of a mutated kinase can be predicted effectively by identifying a hydrogen bonding fingerprint in the activation loop and the αC-helix regions, despite the fact that mutations in cancer patients occur throughout the kinase domain. In our study, we find that the predictive power of MD is superior to a purely data-driven machine learning model involving biochemical features that we implemented, even though MD utilized far fewer features (in fact, just one) in an unsupervised setting. Moreover, the MD results provide key insights into convergent mechanisms of activation, primarily involving differential stabilization of a hydrogen bond network that engages residues of the activation loop and αC-helix in the active-like conformation (in >70% of the mutations studied, regardless of the location of the mutation).

Identification of the Capsid Binding Site in the Herpes Simplex Virus 1 Nuclear Egress Complex and Its Role in Viral Primary Envelopment and Replication

During nuclear egress of nascent progeny herpesvirus nucleocapsids, the nucleocapsids acquire a primary envelope by budding through the inner nuclear membrane of infected cells into the perinuclear space between the inner and outer nuclear membranes. Herpes simplex virus 1 (HSV-1) U_L34 and U_L31 proteins form a nuclear egress complex (NEC) and play critical roles in this budding process, designated primary envelopment. To clarify the role of NEC binding to progeny nucleocapsids in HSV-1 primary envelopment, we established an assay system for HSV-1 NEC binding to nucleocapsids and capsid proteins in vitro Using this assay system, we showed that HSV-1 NEC bound to nucleocapsids and to capsid protein U_L25 but not to the other capsid proteins tested (i.e., VP5, VP23, and U_L17) and that HSV-1 NEC binding of nucleocapsids was mediated by the interaction of NEC with U_L25. U_L31 residues arginine-281 (R281) and aspartic acid-282 (D282) were required for efficient NEC binding to nucleocapsids and U_L25. We also showed that alanine substitution of U_L31 R281 and D282 reduced HSV-1 replication, caused aberrant accumulation of capsids in the nucleus, and induced an accumulation of empty vesicles that were similar in size and morphology to primary envelopes in the perinuclear space. These results suggested that NEC binding via U_L31 R281 and D282 to nucleocapsids, and probably to U_L25 in the nucleocapsids, has an important role in HSV-1 replication by promoting the incorporation of nucleocapsids into vesicles during primary envelopment.IMPORTANCE Binding of HSV-1 NEC to nucleocapsids has been thought to promote nucleocapsid budding at the inner nuclear membrane and subsequent incorporation of nucleocapsids into vesicles during nuclear egress of nucleocapsids. However, data to directly support this hypothesis have not been reported thus far. In this study, we have present data showing that two amino acids in the membrane-distal face of the HSV-1 NEC, which contains the putative capsid binding site based on the solved NEC structure, were in fact required for efficient NEC binding to nucleocapsids and for efficient incorporation of nucleocapsids into vesicles during primary envelopment. This is the first report showing direct linkage between NEC binding to nucleocapsids and an increase in nucleocapsid incorporation into vesicles during herpesvirus primary envelopment.

Computational algorithms for in silico profiling of activating mutations in cancer

Methods to catalog and computationally assess the mutational landscape of proteins in human cancers are desirable. One approach is to adapt evolutionary or data-driven methods developed for predicting whether a single-nucleotide polymorphism (SNP) is deleterious to protein structure and function. In cases where understanding the mechanism of protein activation and regulation is desired, an alternative approach is to employ structure-based computational approaches to predict the effects of point mutations. Through a case study of mutations in kinase domains of three proteins, namely, the anaplastic lymphoma kinase (ALK) in pediatric neuroblastoma patients, serine/threonine-protein kinase B-Raf (BRAF) in melanoma patients, and erythroblastic oncogene B 2 (ErbB2 or HER2) in breast cancer patients, we compare the two approaches above. We find that the structure-based method is most appropriate for developing a binary classification of several different mutations, especially infrequently occurring ones, concerning the activation status of the given target protein. This approach is especially useful if the effects of mutations on the interactions of inhibitors with the target proteins are being sought. However, many patients will present with mutations spread across different target proteins, making structure-based models computationally demanding to implement and execute. In this situation, data-driven methods-including those based on machine learning techniques and evolutionary methods-are most appropriate for recognizing and illuminate mutational patterns. We show, however, that, in the present status of the field, the two methods have very different accuracies and confidence values, and hence, the optimal choice of their deployment is context-dependent.

Structural versatility that serves the function of the HRD motif in the catalytic loop of protein tyrosine kinase, Src

Site-directed mutagenesis is a traditional approach for structure-function analysis of protein tyrosine kinases, and it requires the generation, expression, purification, and analysis of each mutant enzyme. In this study, we report a versatile high throughput bacterial screening system that can identify functional kinase mutants by immunological detection of tyrosine phosphorylation. Two key features of this screening system are noteworthy. First, instead of blotting bacterial colonies directly from Agar plates to nitrocellulose membrane, the colonies were cultured in 96-well plates, and then spotted in duplicate onto the membrane with appropriate controls. This made the screening much more reliable compared with direct colony blotting transfer. A second feature is the parallel use of a protein tyrosine phosphatase (PTP)-expressing host and a non-PTP-expressing host. Because high activity Src mutants are toxic to the host, the PTP system allowed the identification of Src mutants with high activity, while the non-PTP system identified Src mutants with low activity. This approach was applied to Src mutant libraries randomized in the highly conserved HRD motif in the catalytic loop, and revealed that structurally diverse residues can replace the His and Arg residues, while the Asp residue is irreplaceable for catalytic activity.

The Src module: an ancient scaffold in the evolution of cytoplasmic tyrosine kinases

Tyrosine kinases were first discovered as the protein products of viral oncogenes. We now know that this large family of metazoan enzymes includes nearly one hundred structurally diverse members. Tyrosine kinases are broadly classified into two groups: the transmembrane receptor tyrosine kinases, which sense extracellular stimuli, and the cytoplasmic tyrosine kinases, which contain modular ligand-binding domains and propagate intracellular signals. Several families of cytoplasmic tyrosine kinases have in common a core architecture, the "Src module," composed of a Src-homology 3 (SH3) domain, a Src-homology 2 (SH2) domain, and a kinase domain. Each of these families is defined by additional elaborations on this core architecture. Structural, functional, and evolutionary studies have revealed a unifying set of principles underlying the activity and regulation of tyrosine kinases built on the Src module. The discovery of these conserved properties has shaped our knowledge of the workings of protein kinases in general, and it has had important implications for our understanding of kinase dysregulation in disease and the development of effective kinase-targeted therapies.

miR-26a desensitizes non-small cell lung cancer cells to tyrosine kinase inhibitors by targeting PTPN13

Epidermal growth factor receptor (EGFR)-targeted tyrosine kinase inhibitors (TKIs) have emerged as first-line drugs for non-small cell lung cancers (NSCLCs). However, the resistance to TKIs represents the key limitation for their therapeutic efficacy. We found that miR-26a was upregulated in gefitinib-refractory NSCLCs; miR-26a is downstream of EGFR signaling and directly targets and silences protein tyrosine phosphatase non-receptor type 13 (PTPN13) to maintain the activation of Src, a dephosphorylation substrate of PTPN13, thus reinforcing EGFR pathway in a regulatory circuit. miR-26a inhibition significantly improved NSCLC responses to gefitinib. These data revealed a novel mechanism of NSCLC resistance to TKI treatment.

Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases

Protein kinases are enzymes that catalyze the covalent transfer of the γ-phosphate of an adenosine triphosphate (ATP) molecule onto a tyrosine, serine, threonine, or histidine residue in the substrate and thus send a chemical signal to networks of downstream proteins. They are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Unregulated protein kinase activity is often associated with a wide range of diseases, therefore making protein kinases major therapeutic targets. A prototypical system of central interest to understand the regulation of kinase activity is provided by tyrosine kinase c-Src, which belongs to the family of Src-related non-receptor tyrosine kinases (SFKs). Although the broad picture of autoinhibition via the regulatory domains and via the phosphorylation of the C-terminal tail is well characterized from a structural point of view, a detailed mechanistic understanding at the atomic-level is lacking. Advanced computational methods based on all-atom molecular dynamics (MD) simulations are employed to advance our understanding of tyrosine kinase activation. The computational studies suggest that the isolated kinase domain (KD) is energetically most favorable in the inactive conformation when the activation loop (A-loop) of the KD is not phosphorylated. The KD makes transient visits to a catalytically competent active-like conformation. The process of bimolecular trans-autophosphorylation of the A-loop eventually locks the KD in the active state. Activating point mutations may act by slightly increasing the population of the active-like conformation, enhancing the availability of the A-loop to be phosphorylated. The Src-homology 2 (SH2) and Src-homology 3 (SH3) regulatory domains, depending upon their configuration, either promote the inactive or the active state of the kinase domain. In addition to the roles played by the SH3, SH2, and KD, the Src-homology 4-Unique domain (SH4-U) region also serves as a key moderator of substrate specificity and kinase function. Thus, a fundamental understanding of the conformational propensity of the SH4-U region and how this affects the association to the membrane surface are likely to lead to the discovery of new intermediate states and alternate strategies for inhibition of kinase activity for drug discovery. The existence of a multitude of KD conformations poses a great challenge aimed at the design of specific inhibitors. One promising computational strategy to explore the conformational flexibility of the KD is to construct Markov state models from aggregated MD data.

The Activation of c-Src Tyrosine Kinase: Conformational Transition Pathway and Free Energy Landscape

Tyrosine kinases are important cellular signaling allosteric enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Their activity must be tightly controlled, and malfunction can lead to a variety of diseases, particularly cancer. The nonreceptor tyrosine kinase c-Src, a prototypical model system and a representative member of the Src-family, functions as complex multidomain allosteric molecular switches comprising SH2 and SH3 domains modulating the activity of the catalytic domain. The broad picture of self-inhibition of c-Src via the SH2 and SH3 regulatory domains is well characterized from a structural point of view, but a detailed molecular mechanism understanding is nonetheless still lacking. Here, we use advanced computational methods based on all-atom molecular dynamics simulations with explicit solvent to advance our understanding of kinase activation. To elucidate the mechanism of regulation and self-inhibition, we have computed the pathway and the free energy landscapes for the "inactive-to-active" conformational transition of c-Src for different configurations of the SH2 and SH3 domains. Using the isolated c-Src catalytic domain as a baseline for comparison, it is observed that the SH2 and SH3 domains, depending upon their bound orientation, promote either the inactive or active state of the catalytic domain. The regulatory structural information from the SH2-SH3 tandem is allosterically transmitted via the N-terminal linker of the catalytic domain. Analysis of the conformational transition pathways also illustrates the importance of the conserved tryptophan 260 in activating c-Src, and reveals a series of concerted events during the activation process.

Conformational Selection and Induced Fit Mechanisms in the Binding of an Anticancer Drug to the c-Src Kinase

Understanding the conformational changes associated with the binding of small ligands to their biological targets is a fascinating and meaningful question in chemistry, biology and drug discovery. One of the most studied and important is the so-called "DFG-flip" of tyrosine kinases. The conserved three amino-acid DFG motif undergoes an "in to out" movement resulting in a particular inactive conformation to which "type II" kinase inhibitors, such as the anti-cancer drug Imatinib, bind. Despite many studies, the details of this prototypical conformational change are still debated. Here we combine various NMR experiments and surface plasmon resonance with enhanced sampling molecular dynamics simulations to shed light into the conformational dynamics associated with the binding of Imatinib to the proto-oncogene c-Src. We find that both conformational selection and induced fit play a role in the binding mechanism, reconciling opposing views held in the literature. Moreover, an external binding pose and local unfolding (cracking) of the aG helix are observed.

Computational study of the W260A activating mutant of Src tyrosine kinase

Tyrosine kinases are enzymes playing a critical role in cellular signaling. Mutations causing increased in kinase activity are often associated with cancer and various pathologies. One example in Src tyrosine kinases is offered by the substitution of the highly conserved tryptophan 260 by an alanine (W260A), which has been shown to cause an increase in activity. Here, molecular dynamics simulations based on atomic models are carried out to characterize the conformational changes in the linker region and the catalytic (kinase) domain of Src kinase to elucidate the impact of the W260A mutation. Umbrella sampling calculations show that the conformation of the linker observed in the assembled down-regulated state of the kinase is most favored when the kinase domain is in the inactive state, whereas the conformation of the linker observed in the re-assembled up-regulated state of the kinase is favored when the kinase domain is in the unphosphorylated active-like state. The calculations further indicate that there are only small differences between the WT and W260A mutant. In both cases, the intermediates states are very similar and the down-regulated inactive conformation is the most stable state. However, the calculations also show that the free energy cost to reach the unphosphorylated active-like conformation is slightly smaller for the W260A mutant compared with WT. A simple kinetic model is developed and submitted to a Bayesian Monte Carlo analysis to illustrate how such small differences can contribute to accelerate the trans-autophosphorylation reaction and yield a large increase in the activity of the mutant as observed experimentally.

Advanced ensemble modelling of flexible macromolecules using X-ray solution scattering

Dynamic ensembles of macromolecules mediate essential processes in biology. Understanding the mechanisms driving the function and molecular interactions of 'unstructured' and flexible molecules requires alternative approaches to those traditionally employed in structural biology. Small-angle X-ray scattering (SAXS) is an established method for structural characterization of biological macromolecules in solution, and is directly applicable to the study of flexible systems such as intrinsically disordered proteins and multi-domain proteins with unstructured regions. The Ensemble Optimization Method (EOM) [Bernadó et al. (2007 ▶). J. Am. Chem. Soc. 129, 5656-5664] was the first approach introducing the concept of ensemble fitting of the SAXS data from flexible systems. In this approach, a large pool of macromolecules covering the available conformational space is generated and a sub-ensemble of conformers coexisting in solution is selected guided by the fit to the experimental SAXS data. This paper presents a series of new developments and advancements to the method, including significantly enhanced functionality and also quantitative metrics for the characterization of the results. Building on the original concept of ensemble optimization, the algorithms for pool generation have been redesigned to allow for the construction of partially or completely symmetric oligomeric models, and the selection procedure was improved to refine the size of the ensemble. Quantitative measures of the flexibility of the system studied, based on the characteristic integral parameters of the selected ensemble, are introduced. These improvements are implemented in the new EOM version 2.0, and the capabilities as well as inherent limitations of the ensemble approach in SAXS, and of EOM 2.0 in particular, are discussed.

Autoinhibition of Bruton's tyrosine kinase (Btk) and activation by soluble inositol hexakisphosphate

Bruton's tyrosine kinase (Btk), a Tec-family tyrosine kinase, is essential for B-cell function. We present crystallographic and biochemical analyses of Btk, which together reveal molecular details of its autoinhibition and activation. Autoinhibited Btk adopts a compact conformation like that of inactive c-Src and c-Abl. A lipid-binding PH-TH module, unique to Tec kinases, acts in conjunction with the SH2 and SH3 domains to stabilize the inactive conformation. In addition to the expected activation of Btk by membranes containing phosphatidylinositol triphosphate (PIP3), we found that inositol hexakisphosphate (IP6), a soluble signaling molecule found in both animal and plant cells, also activates Btk. This activation is a consequence of a transient PH-TH dimerization induced by IP6, which promotes transphosphorylation of the kinase domains. Sequence comparisons with other Tec-family kinases suggest that activation by IP6 is unique to Btk.

Hotspot mutations in KIT receptor differentially modulate its allosterically coupled conformational dynamics: impact on activation and drug sensitivity

Receptor tyrosine kinase KIT controls many signal transduction pathways and represents a typical allosterically regulated protein. The mutation-induced deregulation of KIT activity impairs cellular physiological functions and causes serious human diseases. The impact of hotspots mutations (D816H/Y/N/V and V560G/D) localized in crucial regulatory segments, the juxtamembrane region (JMR) and the activation (A-) loop, on KIT internal dynamics was systematically studied by molecular dynamics simulations. The mutational outcomes predicted in silico were correlated with in vitro and in vivo activation rates and drug sensitivities of KIT mutants. The allosteric regulation of KIT in the native and mutated forms is described in terms of communication between the two remote segments, JMR and A-loop. A strong correlation between the communication profile and the structural and dynamical features of KIT in the native and mutated forms was established. Our results provide new insight on the determinants of receptor KIT constitutive activation by mutations and resistance of KIT mutants to inhibitors. Depiction of an intra-molecular component of the communication network constitutes a first step towards an integrated description of vast communication pathways established by KIT in physiopathological contexts.

Receptor tyrosine kinase KIT plays a crucial role in the regulation of cell signaling. This allosterically controlled activity may be affected by gain-of-function mutations that promote the development of several cancers. Identification of the molecular basis of KIT constitutive activation and allosteric regulation has inspired computational study of KIT hotspot mutations. In the present contribution, we investigated the mutation-induced effects on KIT conformational dynamics and intra-protein communication conditionally on the mutation location and the nature of the substituting amino acid. Our data elucidate that all studied mutations stabilize an inactive non-autoinhibited state of KIT over the inactive auto-inhibited state prevalent for the native protein. This shift in the protein conformational landscape promotes KIT constitutive activation. Our in silico analysis established correlations between the structural and dynamical effects induced by oncogenic mutations and the mutants auto-activation rates and drug sensitivities measured in vitro and in vivo. Particularly, the A-loop mutations stabilize the drug-resistant forms, while the JMR mutations may facilitate inhibitors binding to the active site. Cross-correlations established between local and long-range structural and dynamical effects demonstrate the allosteric character of the gain-of-function mutations mode of action.

SRC-family tyrosine kinases in oogenesis, oocyte maturation and fertilization: an evolutionary perspective

The oocyte is a highly specialized cell poised to respond to fertilization with a unique set of actions needed to recognize and incorporate a single sperm, complete meiosis, reprogram maternal and paternal genomes and assemble them into a unique zygotic genome, and finally initiate the mitotic cell cycle. Oocytes accomplish this diverse series of events through an array of signal transduction pathway components that include a characteristic collection of protein tyrosine kinases. The src-family protein kinases (SFKs) figure importantly in this signaling array and oocytes characteristically express certain SFKs at high levels to provide for the unique actions that the oocyte must perform. The SFKs typically exhibit a distinct pattern of subcellular localization in oocytes and perform critical functions in different subcellular compartments at different steps during oocyte maturation and fertilization. While many aspects of SFK signaling are conserved among oocytes from different species, significant differences exist in the extent to which src-family-mediated pathways are used by oocytes from species that fertilize externally vs those which are fertilized internally. The observation that several oocyte functions which require SFK signaling appear to represent common points of failure during assisted reproductive techniques in humans, highlights the importance of these signaling pathways for human reproductive health.

Locking the active conformation of c-Src kinase through the phosphorylation of the activation loop

Molecular dynamics umbrella sampling simulations are used to compare the relative stability of the active conformation of the catalytic domain of c-Src kinase while the tyrosine 416 in the activation loop (A-loop) is either unphosphorylated or phosphorylated. When the A-loop is unphosphorylated, there is considerable flexibility of the kinase. While the active conformation of the kinase is not forbidden and can be visited transiently, it is not the predominant state. This is consistent with the view that c-Src displays some catalytic activity even when the A-loop is unphosphorylated. In contrast, phosphorylation of the A-loop contributes to stabilize several structural features that are critical for catalysis, such as the hydrophobic regulatory spine, the HRD motif, and the electrostatic switch. In summary, the free-energy landscape calculations demonstrate that phosphorylation of tyrosine 416 in the A-loop essentially "locks" the kinase into its catalytically competent conformation.

Cracking the molecular origin of intrinsic tyrosine kinase activity through analysis of pathogenic gain-of-function mutations

The basal (ligand-independent) kinase activity of receptor tyrosine kinases (RTKs) promotes trans-phosphorylation on activation loop tyrosines upon ligand-induced receptor dimerization, thus upregulating intrinsic kinase activity and triggering intracellular signaling. To understand the molecular determinants of intrinsic kinase activity, we used X-ray crystallography and NMR spectroscopy to analyze pathogenic FGF receptor mutants with gradations in gain-of-function activity. These structural analyses revealed a "two-state" dynamic equilibrium model whereby the kinase toggles between an "inhibited," structurally rigid ground state and a more dynamic and heterogeneous active state. The pathogenic mutations have different abilities to shift this equilibrium toward the active state. The increase in the fractional population of FGF receptors in the active state correlates with the degree of gain-of-function activity and clinical severity. Our data demonstrate that the fractional population of RTKs in the active state determines intrinsic kinase activity and underscore how a slight increase in the active population of kinases can have grave consequences for human health.

"DFG-flip" in the insulin receptor kinase is facilitated by a helical intermediate state of the activation loop

We have characterized a large-scale inactive-to-active conformational change in the activation-loop of the insulin receptor kinase domain at the atomistic level via untargeted temperature-accelerated molecular dynamics (TAMD) and free-energy calculations using the string method. TAMD simulations consistently show folding of the A-loop into a helical conformation followed by unfolding to an active conformation, causing the highly conserved DFG-motif (Asp(1150), Phe(1151), and Gly(1152)) to switch from the inactive "D-out/F-in" to the nucleotide-binding-competent "D-in/F-out" conformation. The minimum free-energy path computed from the string method preserves these helical intermediates along the inactive-to-active path, and the thermodynamic free-energy differences are consistent with previous work on various other kinases. The mechanisms revealed by TAMD also suggest that the regulatory spine can be dynamically assembled/disassembled either by DFG-flip or by movement of the αC-helix. Together, these findings both broaden our understanding of kinase activation and point to intermediates as specific therapeutic targets.

Allosteric communication across the native and mutated KIT receptor tyrosine kinase

A fundamental goal in cellular signaling is to understand allosteric communication, the process by which signals originated at one site in a protein propagate dependably to affect remote functional sites. Here, we describe the allosteric regulation of the receptor tyrosine kinase KIT. Our analysis evidenced that communication routes established between the activation loop (A-loop) and the distant juxtamembrane region (JMR) in the native protein were disrupted by the oncogenic mutation D816V positioned in the A-loop. In silico mutagenesis provided a plausible way of restoring the protein communication detected in the native KIT by introducing a counter-balancing second mutation D792E. The communication patterns observed in the native and mutated KIT correlate perfectly with the structural and dynamical features of these proteins. Particularly, a long-distance effect of the D816V mutation manifested as an important structural re-organization of the JMR in the oncogenic mutant was completely vanished in the double mutant D816V/D792E. This detailed characterization of the allosteric communication in the different forms of KIT, native and mutants, was performed by using a modular network representation composed of communication pathways and independent dynamic segments. Such representation permits to enrich a purely mechanistic interaction-based model of protein communication by the introduction of concerted local atomic fluctuations. This method, validated on KIT receptor, may guide a rational modulation of the physiopathological activities of other receptor tyrosine kinases.

The majority of functionally important biological processes are regulated by allosteric communication within individual proteins and across protein complexes. Receptor tyrosine kinases (RTKs) control signal transduction pathways and consequently represent a typical paradigm. The mutation-induced deregulation of RTK activity impairs crucial cellular physiological functions and causes serious human diseases. The present study focuses on the allosteric communication across the three-dimensional structure of the RTK KIT cytoplasmic region. Combining a mechanistic model of information transmission with the analysis of concerted local atomic fluctuations we examined and compared the communication profiles in the native and D816V-mutated proteins. This approach permitted to localize and visualize communication routes in the native KIT and revealed that these routes were disrupted in the mutant D816V. We proposed in silico mutagenesis as a mean to restore the communication detected in the native KIT. Our work sheds light on the allosteric communication in RTKs, a phenomenon playing an essential role in signaling pathways albeit experiments do not provide the atomic details of the path followed in going from one structural element to the other. A rational understanding of the molecular determinants underlying the effects of disease-related kinase mutations may contribute to the improvement of targeted therapies.

αC helix as a switch in the conformational transition of Src/CDK-like kinase domains

One mechanism of regulating the catalytic activity of protein kinases is through conformational transitions. Despite great diversity in the structural changes involved in the transitions, a certain set of changes within the kinase domain (KD) has been observed for many kinases including Src and CDK2. We investigated this conformational transition computationally to identify the topological features that are energetically critical to the transition. Results from both molecular dynamics sampling and transition path optimization highlight the displacement of the αC helix as the major energy barrier, mediating the switch of the KD between the active and down-regulated states. The critical role of the αC helix is noteworthy by providing a rationale for a number of activation and deactivation mechanisms known to occur in cells. We find that kinases with the αC helix displacement exist throughout the kinome, suggesting that this feature may have emerged early in evolution.

Insights into the conformational switching mechanism of the human vascular endothelial growth factor receptor type 2 kinase domain

Human vascular endothelial growth factor receptor type 2 (h-VEFGR2) is a receptor tyrosine kinase involved in the angiogenesis process and regarded as an interesting target for the design of anticancer drugs. Its activation/inactivation mechanism is related to conformational changes in its cytoplasmatic kinase domain, involving first among all the αC-helix in N-lobe and the A-loop in C-lobe. Affinity of inhibitors for the active or inactive kinase form could dictate the open or closed conformation of the A-loop, thus making the different conformations of the kinase domain receptor (KDR) domain different drug targets in drug discovery. In this view, a detailed knowledge of the conformational landscape of KDR domain is of central relevance to rationalize the efficiency and selectivity of kinase inhibitors. Here, molecular dynamics simulations were used to gain insight into the conformational switching activity of the KDR domain and to identify intermediate conformations between the two limiting active and inactive conformations. Specific energy barriers have been selectively removed to induce, and hence highlight at the atomistic level, the regulation mechanism of the A-loop opening. The proposed strategy allowed to repeatedly observe the escape of the KDR domain from the DFG-out free energy basin and to identify rare intermediate conformations between the DFG-out and the DFG-in structures to be employed in a structure-based drug discovery process.

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.

Despite recent progress in computational and experimental studies of dynamic regulation in protein kinases, a mechanistic understanding of long-range communication and mechanisms of mutation-induced signaling controlling kinase activity remains largely qualitative. In this study, we have performed a systematic modeling and analysis of allosteric activation in ABL and EGFR kinases at the increasing level of complexity - from catalytic domain to multi-domain regulatory complexes. The results of this study have revealed organizing structural and mechanistic principles of allosteric signaling in protein kinases. Although activation mechanisms in ABL and EGFR kinases have evolved through acquisition of structurally different regulatory complexes, we have found that long-range interdomain communication between common functional segments (αF-helix and αC-helix) may be important for allosteric activation. The results of study have revealed molecular signatures of activating cancer mutations and have shed the light on general mechanistic aspects of mutation-induced signaling in protein kinases. An advanced understanding and further characterization of molecular signatures of kinase mutations may aid in a better rationalization of mutational effects on clinical outcomes and facilitate molecular-based therapeutic strategies to combat kinase mutation-dependent tumorigenesis.

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.

Protein kinases are involved in a huge amount of cellular processes through phosphorylation, a crucial mechanism in cell signaling, and their misregulation often results in disease. The deactivation of protein tyrosine kinases (PTKs) or their oncogenic activation arises from mutations which affect the protein primary structure and the configuration of the enzymatic site apparently by stabilizing the activation loop (A-loop) extended conformation. Particularly, mutation D816V of receptor tyrosine kinase (RTK) KIT, found in patients with pediatric mastocytosis, acute leukemia or germ cell tumors, can be considered as the archetype of mutation inducing a displacement of the population equilibrium toward the active conformation. We present a comprehensive computational study of the activating mechanism(s) of this mutation. Our multi-approach in silico procedure evidenced a local alteration of the A-loop structure, 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. Our results provided a plausible conception of how the observed departure of JMR from kinase domain in the mutant promotes kinase mutant dimerization without requiring extra-cellular ligand binding. The pocket profiles we obtained suggested putative allosteric binding sites that could be targeted by ligands/modulators that trap the mutated enzyme.

Transglutaminase 2 null macrophages respond to lipopolysaccharide stimulation by elevated proinflammatory cytokine production due to an enhanced αvβ3 integrin-induced Src tyrosine kinase signaling

Transglutaminase 2 (TG2) is a protein crosslinking enzyme with several additional biochemical functions. Loss of TG2 in vivo results in impaired phagocytosis of apoptotic cells and altered proinflammatory cytokine production by macrophages engulfing apoptotic cells leading to autoimmunity. It has been proposed that TG2 acts as an integrin β(3) coreceptor in the engulfment process, while altered proinflammatory cytokine production is related to the lack of latent TGFβ activation by TG2 null macrophages. Here we report that TG2 null macrophages respond to lipopolysaccharide treatment by elevated IL-6 and TNFα production. Though TGFβ has been proposed to act as a feed back regulator of proinflammatory cytokine production in LPS-stimulated macrophages, this phenomenon is not related to the lack of active TGFβ production. Instead, in the absence of TG2 integrin β(3) maintains an elevated basal Src family kinase activity in macrophages, which leads to enhanced phosphorylation and degradation of the IκBα. Low basal levels of IκBα explain the enhanced sensitivity of TG2 null macrophages to signals that regulate NF-κB. Our data suggest that TG2 null macrophages bear a proinflammatory phenotype, which might contribute to the enhanced susceptibility of these mice to develop autoimmunity and atherosclerosis.

Nanofibers Self-assembled from Structural Complementary Borono-decapeptides

A series of structural complementary decapeptides with phenyl boronic acid tails or borono-decapeptides (BPs) were designed and synthesized for supramolecular self-assembly. After dissolving these borono-decapeptides in deionized (DI) water, well-defined nanofibers were formed in BP1 (B(OH)(2) VEKELVKEKL-OH) and BP3 (B(OH)(2) AELELARARL-OH). It was found that the self-assembled borono-decapeptide BP1 and BP3 have a parallel β-sheet conformation in the formed nanofibers. The strategy demonstrated here shows a great prospect in preparation of well-ordered nanofibers via rationally designing the molecular structures of peptides.

New pyrazolo-[3,4-d]-pyrimidine derivative Src kinase inhibitors lead to cell cycle arrest and tumor growth reduction of human medulloblastoma cells

Medulloblastoma is the most common malignant brain tumor in children, and despite improvements in the overall survival rate, it still lacks an effective treatment. Src plays an important role in cancer, and recently high Src activity was documented in medulloblastoma. In this report, we examined the effects of novel pyrazolo-[3,4-d]-pyrimidine derivative Src inhibitors in medulloblastoma. By MTS assay, we showed that the pyrimidine derivatives indicated as S7, S29, and SI163 greatly reduce the growth rate of medulloblastoma cells by inhibiting Src phosphorylation, compared with HT22 non-neoplastic nerve cells. These compounds also halt cells in the G(2)/M phase, and this effect likely occurs through the regulation of cdc2 and CDC25C phosphorylation, as shown by Western blot. Moreover, the exposure to pyrimidine derivatives induces apoptosis, assayed by the supravital propidium iodide assay, through modulation of the apoptotic proteins Bax and Bcl2, and inhibits tumor growth in vivo in a mouse model. Notably, S7, S29, and SI163 show major inhibitory effects on medulloblastoma cell growth compared with the chemotherapeutic agents cisplatin and etoposide. In conclusion, our results suggest that S7, S29, and SI163 could be novel attractive candidates for the treatment of medulloblastoma or tumors characterized by high Src activity.

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.

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.

Mutations in protein kinases are implicated in many cancers, and an important goal of cancer research is to elucidate molecular effects of mutated kinase genes that contribute to tumorigenesis. We present a comprehensive computational study of molecular mechanisms of kinase activation by cancer-causing mutations. Using a battery of computational approaches, we have systematically investigated the effects of clinically important cancer mutants on dynamics of the ABL and EGFR kinase domains and regulatory multi-protein complexes. The results of this study have illuminated common and specific features of the activation mechanism in the normal and oncogenic forms of ABL and EGFR. We have found that mutants with the higher oncogenic activity may cause a partial destabilization of the inactive structure, while simultaneously facilitating activating transitions and the enhanced stabilization of the active conformation. Our results provided useful insights into thermodynamic and mechanistic aspects of the activation mechanism and highlighted the role of structurally distinct conformational states in kinase regulation. Ultimately, molecular signatures of activation mechanisms in the normal and oncogenic states may aid in the correlation of mutational effects with clinical outcomes and facilitate the development of therapeutic strategies to combat kinase mutation-dependent tumorigenesis.

Atomistic view of the conformational activation of Src kinase using the string method with swarms-of-trajectories

The inactive-to-active conformational transition of the catalytic domain of human c-Src tyrosine kinase is characterized using the string method with swarms-of-trajectories with all-atom explicit solvent molecular dynamics simulations. The activation process occurs in two main steps in which the activation loop (A-loop) opens first, followed by the rotation of the alphaC helix. The computed potential of mean force energy along the activation pathway displays a local minimum, which allows the identification of an intermediate state. These results show that the string method with swarms-of-trajectories is an effective technique to characterize complex and slow conformational transitions in large biomolecular systems.

Mapping the conformational transition in Src activation by cumulating the information from multiple molecular dynamics trajectories

The Src-family kinases are allosteric enzymes that play a key role in the regulation of cell growth and proliferation. In response to cellular signals, they undergo large conformational changes to switch between distinct inactive and active states. A computational strategy for characterizing the conformational transition pathway is presented to bridge the inactive and active states of the catalytic domain of Hck. The information from a large number (78) of independent all-atom molecular dynamics trajectories with explicit solvent is combined together to assemble a connectivity map of the conformational transition. Two intermediate states along the activation pathways are identified, and their structural features are characterized. A coarse free-energy landscape is built in terms of the collective motions corresponding to the opening of the activation loop (A-loop) and the rotation of the alphaC helix. This landscape shows that the protein can adopt a multitude of conformations in which the A-loop is partially open, while the alphaC helix remains in the orientation characteristic of the inactive conformation. The complete transition leading to the active conformation requires a concerted movement involving further opening of the A-loop, the relative alignment of N-lobe and C-lobe, and the rotation of the alphaC helix needed to recruit the residues necessary for catalysis in the active site. The analysis leads to a dynamic view of the full-length kinase activation, whereby transitions of the catalytic domain to intermediate configurations with a partially open A-loop are permitted, even while the SH2-SH3 clamp remains fully engaged. These transitions would render Y416 available for the transphosphorylation event that ultimately locks down the active state. The results provide a broad framework for picturing the conformational transitions leading to kinase activation.