Electrostatic interactions of peptides flanking the tyrosine kinase domain in the epidermal growth factor receptor provides a model for intracellular dimerization and autophosphorylation

Proteins with calmodulin-like domains: structures and functional roles

The appearance of modular proteins is a widespread phenomenon during the evolution of proteins. The combinatorial arrangement of different functional and/or structural domains within a single polypeptide chain yields a wide variety of activities and regulatory properties to the modular proteins. In this review, we will discuss proteins, that in addition to their catalytic, transport, structure, localization or adaptor functions, also have segments resembling the helix-loop-helix EF-hand motifs found in Ca2+-binding proteins, such as calmodulin (CaM). These segments are denoted CaM-like domains (CaM-LDs) and play a regulatory role, making these CaM-like proteins sensitive to Ca2+ transients within the cell, and hence are able to transduce the Ca2+ signal leading to specific cellular responses. Importantly, this arrangement allows to this group of proteins direct regulation independent of other Ca2+-sensitive sensor/transducer proteins, such as CaM. In addition, this review also covers CaM-binding proteins, in which their CaM-binding site (CBS), in the absence of CaM, is proposed to interact with other segments of the same protein denoted CaM-like binding site (CLBS). CLBS are important regulatory motifs, acting either by keeping these CaM-binding proteins inactive in the absence of CaM, enhancing the stability of protein complexes and/or facilitating their dimerization via CBS/CLBS interaction. The existence of proteins containing CaM-LDs or CLBSs substantially adds to the enormous versatility and complexity of Ca2+/CaM signaling.

Screening of nucleotide variations in genomic sequences encoding charged protein regions in the human genome

Background

Studying genetic variation distribution in proteins containing charged regions, called charge clusters (CCs), is of great interest to unravel their functional role. Charge clusters are 20 to 75 residue segments with high net positive charge, high net negative charge, or high total charge relative to the overall charge composition of the protein. We previously developed a bioinformatics tool (FCCP) to detect charge clusters in proteomes and scanned the human proteome for the occurrence of CCs. In this paper we investigate the genetic variations in the human proteins harbouring CCs.

Results

We studied the coding regions of 317 positively charged clusters and 1020 negatively charged ones previously detected in human proteins. Results revealed that coding parts of CCs are richer in sequence variants than their corresponding genes, full mRNAs, and exonic + intronic sequences and that these variants are predominately rare (Minor allele frequency < 0.005). Furthermore, variants occurring in the coding parts of positively charged regions of proteins are more often pathogenic than those occurring in negatively charged ones. Classification of variants according to their types showed that substitution is the major type followed by Indels (Insertions-deletions). Concerning substitutions, it was found that within clusters of both charges, the charged amino acids were the greatest loser groups whereas polar residues were the greatest gainers.

Conclusions

Our findings highlight the prominent features of the human charged regions from the DNA up to the protein sequence which might provide potential clues to improve the current understanding of those charged regions and their implication in the emergence of diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4000-3) contains supplementary material, which is available to authorized users.

Screening for clusters of charge in human virus proteomes

Background

The identification of charge clusters (runs of charged residues) in proteins and their mapping within the protein structure sequence is an important step toward a comprehensive analysis of how these particular motifs mediate, via electrostatic interactions, various molecular processes such as protein sorting, translocation, docking, orientation and binding to DNA and to other proteins. Few algorithms that specifically identify these charge clusters have been designed and described in the literature. In this study, 197 distinctive human viral proteomes were screened for the occurrence of charge clusters (CC) using a new computational approach.

Results

Three hundred and seventy three CC have been identified within the 2549 viral protein sequences screened. The number of protein sequences that are CC-free is 2176 (85.3 %) while 150 and 180 proteins contained positive charge (PCC) and negative charge clusters (NCC), respectively. The NCCs (211 detected) were more prevalent than PCC (162). PCC-containing proteins are significantly longer than those having NCCs (p = 2.10^-16). The most prevalent virus families having PCC and NCC were Herpesviridae followed by Papillomaviridae. However, the single-strand RNA group has in average three times more NCC than PCC. According to the functional domain classification, a significant difference in distribution was observed between PCC and NCC (p = 2. 10⁻⁸) with the occurrence of NCCs being more frequent in C-terminal region while PCC more often fall within functional domains. Only 29 proteins sequences contained both NCC and PCC. Moreover, 101 NCC were conserved in 84 proteins while only 62 PCC were conserved in 60 protein sequences. To understand the mechanism by which the membrane translocation functionalities are embedded in viral proteins, we screened our PCC for sequences corresponding to cell-penetrating peptides (CPPs) using two online databases: CellPPd and CPPpred. We found that all our PCCs, having length varying from 7 to 30 amino-acids were predicted as CPPs. Experimental validation is required to improve our understanding of the role of these PCCs in viral infection process.

Conclusions

Screening distinctive cluster charges in viral proteomes suggested a functional role of these protein regions and might provide potential clues to improve the current understanding of viral diseases in order to tailor better preventive and therapeutic approaches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3086-3) contains supplementary material, which is available to authorized users.

Exploring proteome-wide occurrence of clusters of charged residues in eukaryotes

Clusters of charged residues are one of the key features of protein primary structure since they have been associated to important functions of proteins. Here, we present a proteome wide scan for the occurrence of Charge Clusters in Protein sequences using a new search tool (FCCP) based on a score-based methodology. The FCCP was run to search charge clusters in seven eukaryotic proteomes: Arabidopsis thaliana, Caenorhabditis elegans, Danio rerio, Drosophila melanogaster, Homo sapiens, Mus musculus, and Saccharomyces cerevisiae. We found that negative charge clusters (NCCs) are three to four times more frequent than positive charge clusters (PCCs). The Drosophila proteome is on average the most charged, whereas the human proteome is the least charged. Only 3 to 8% of the studied protein sequences have negative charge clusters, while 1.6 to 3% having PCCs and only 0.07 to 0.6% have both types of clusters. NCCs are localized predominantly in the N-terminal and C-terminal domains, while PCCs tend to be localized within the functional domains of the protein sequences. Furthermore, the gene ontology classification revealed that the protein sequences with negative and PCCs are mainly binding proteins.

A premature termination of human epidermal growth factor receptor transcription in Escherichia coli

Our success in producing an active epidermal growth factor receptor (EGFR) tyrosine kinase in Escherichia coli encouraged us to express the full-length receptor in the same host. Despite its large size, we were successful at producing the full-length EGFR protein fused to glutathione S-transferase (GST) that was detected by Western blot analysis. Moreover, we obtained a majoritarian truncated GST-EGFR form detectable by gel electrophoresis and Western blot. This truncated protein was purified and confirmed by MALDI-TOF/TOF analysis to belong to the N-terminal extracellular region of the EGFR fused to GST. Northern blot analysis showed two transcripts suggesting the occurrence of a transcriptional arrest.

Mutations in the polybasic juxtamembrane sequence of both plasma membrane- and endoplasmic reticulum-localized epidermal growth factor receptors confer ligand-independent cell transformation

Deregulation of ErbB receptor-tyrosine kinases is a hallmark of many human cancers. Conserved in the ErbB family is a cluster of basic amino acid residues in the cytoplasmic juxtamembrane region. We found that charge-silencing mutagenesis within this juxtamembrane region of the epidermal growth factor receptor (EGFR) results in the generation of a mutant receptor (EGFR Mut R1-6) that spontaneously transforms NIH 3T3 cells in a ligand-independent manner. A similar mutant with one additional basic residue, EGFR Mut R1-5, fails to exhibit ligand-independent transformation. The capacity of EGFR Mut R1-6 to mediate this transformation is maintained when this mutant is retained in the endoplasmic reticulum via a single point mutation, L393H, which we describe. We show that EGFR Mut R1-6 with or without L393H exhibits enhanced basal tyrosine phosphorylation when ectopically expressed, and the ligand-independent transforming activity of EGFR Mut R1-6 is sensitive to inhibition of EGFR kinase activity and is particularly dependent on PI3K and mTOR activity. Similar to EGFR Mut R1-6/L393H in NIH 3T3 cells, EGFR variant type III, a highly oncogenic mutant form of EGFR linked to human brain cancers, confers transforming activity while it is wholly endoplasmic reticulum-retained in U87 cells. Our findings highlight the importance of the polybasic juxtamembrane sequence in regulating the oncogenic potential of EGFR signaling.

Increased serum level of epidermal growth factor receptor in liver cancer patients and its association with exposure to arsenic

Arsenic is a human carcinogen and can activate epidermal growth factor receptor (EGFR) in human cell lines. As EGFR is associated with the occurrence of cancers, we conducted a study to evaluate whether serum EGFR may increase in liver cancer patients, particularly in those with exposure to arsenic. We recruited 100 patients of liver cancer and 100 age- and sex-matched controls in Taiwan and determined EGFR levels in sera by enzyme-linked immunosorbent assay. The patients had higher EGFR levels (668.1 vs. 243.1 fmol/mL, p<0.01), and after adjusting for hepatitis B and C, they still had an average EGFR level 406.1 fmol/mL higher than that of the controls (p<0.01). When we compared 22 patients residing in an endemic area of arsenic intoxication to 22 age- and sex-matched patients residing outside the area, we found that patients from the endemic area had higher EGFR levels (882.8 vs. 511.6 fmol/mL, p = 0.04). We concluded that EGFR is over-expressed in patients of liver cancer, particularly in those with exposure to arsenic, and therefore, serum EGFR level is not only a potential biomarker of liver cancer, but also a potential biomarker of cancers associated with arsenic exposure.

Regulation of the ligand-dependent activation of the epidermal growth factor receptor by calmodulin

Calmodulin (CaM) is the major component of calcium signaling pathways mediating the action of various effectors. Transient increases in the intracellular calcium level triggered by a variety of stimuli lead to the formation of Ca(2+)/CaM complexes, which interact with and activate target proteins. In the present study the role of Ca(2+)/CaM in the regulation of the ligand-dependent activation of the epidermal growth factor receptor (EGFR) has been examined in living cells. We show that addition of different cell permeable CaM antagonists to cultured cells or loading cells with a Ca(2+) chelator inhibited ligand-dependent EGFR auto(trans)phosphorylation. This occurred also in the presence of inhibitors of protein kinase C, CaM-dependent protein kinase II and calcineurin, which are known Ca(2+)- and/or Ca(2+)/CaM-dependent EGFR regulators, pointing to a direct effect of Ca(2+)/CaM on the receptor. Furthermore, we demonstrate that down-regulation of CaM in conditional CaM knock out cells stably transfected with the human EGFR decreased its ligand-dependent phosphorylation. Substitution of six basic amino acid residues within the CaM-binding domain (CaM-BD) of the EGFR by alanine resulted in a decreased phosphorylation of the receptor and of its downstream substrate phospholipase Cγ1. These results support the hypothesis that Ca(2+)/CaM regulates the EGFR activity by directly interacting with the CaM-BD of the receptor located at its cytosolic juxtamembrane region.

Ligand-based receptor tyrosine kinase partial agonists: New paradigm for cancer drug discovery?

INTRODUCTION: Receptor tyrosine kinases (RTKs) are validated targets for oncology drug discovery and several RTK antagonists have been approved for the treatment of human malignancies. Nonetheless, the discovery and development of RTK antagonists has lagged behind the discovery and development of agents that target G-protein coupled receptors. In part, this is because it has been difficult to discover analogs of naturally-occurring RTK agonists that function as antagonists. AREAS COVERED: Here we describe ligands of ErbB receptors that function as partial agonists for these receptors, thereby enabling these ligands to antagonize the activity of full agonists for these receptors. We provide insights into the mechanisms by which these ligands function as antagonists. We discuss how information concerning these mechanisms can be translated into screens for novel small molecule- and antibody-based antagonists of ErbB receptors and how such antagonists hold great potential as targeted cancer chemotherapeutics. EXPERT OPINION: While there have been a number of important key findings into this field, the identification of the structural basis of ligand functional specificity is still of the greatest importance. While it is true that, with some notable exceptions, peptide hormones and growth factors have not proven to be good platforms for oncology drug discovery; addressing the fundamental issues of antagonistic partial agonists for receptor tyrosine kinases has the potential to steer oncology drug discovery in new directions. Mechanism based approaches are now emerging to enable the discovery of RTK partial agonists that may antagonize both agonist-dependent and -independent RTK signaling and may hold tremendous promise as targeted cancer chemotherapeutics.

Palmitoylation of the EGF receptor impairs signal transduction and abolishes high-affinity ligand binding

The intracellular juxtamembrane domain of the EGF receptor has been shown to be involved in the stimulation of the receptor's tyrosine kinase activity. To further explore the function of this portion of the EGF receptor, a consensus site for protein palmitoylation was inserted at the beginning of the juxtamembrane domain of the receptor. The altered EGF receptor incorporated [(3)H]palmitate, demonstrating that it was palmitoylated. Compared to the wild-type EGF receptor, the palmitoylated EGF receptor was significantly impaired in EGF-stimulated receptor autophosphorylation as well as ligand-induced receptor internalization. While both the wild-type and the palmitoylated EGF receptors exhibited a similar propensity to associate with lipid rafts, only the wild-type receptor exited lipid rafts in response to EGF. Binding of [(125)I]EGF to the wild-type EGF receptor showed a curvilinear Scatchard plot with both high- and low-affinity forms of the receptor. By contrast, the palmitoylated receptor exhibited only low-affinity EGF binding. These data suggest that the cytoplasmic juxtamembrane domain is involved not only in the transmission of the proliferative signal generated by ligand binding but also in facilitating the adoption of the high-affinity conformation by the extracellular ligand binding domain.

Investigating the function of three non-synonymous SNPs in EGFR gene: structural modelling and association with breast cancer

Non-synonymous single nucleotide polymorphisms (nsSNPs) represent common genomic variations that alter protein sequence and function. Some nsSNPs affecting conserved amino acids have been reported to be associated with cancer susceptibility. Interestingly, Epidermal Growth Factor Receptor (EGFR) is commonly overexpressed and mutated in many cancers. In this study, we investigated the structural effect of three deleterious nsSNPs: rs17337451 (R962G), rs1140476 (R977C) and rs17290699 (H988P) within EGFR using computational tools. The modelled mutant dimers showed less stability than wild type EGFR dimer. Furthermore, we showed the important role of R962 and H988 residues in the EGFR dimer formation. We also report preliminary experimental data for SNP R977C suggesting that the variant C977 might confer greater risk for breast cancer. These results contribute to an improved understanding of the EGFR dimer stability and provide new elements for understanding the relationship between EGFR and cancer.

A conserved threonine residue in the juxtamembrane domain of the XA21 pattern recognition receptor is critical for kinase autophosphorylation and XA21-mediated immunity

Despite the key role that pattern recognition receptors (PRRs) play in regulating immunity in plants and animals, the mechanism of activation of the associated non-arginine-aspartate (non-RD) kinases is unknown. The rice PRR XA21 recognizes the pathogen-associated molecular pattern, Ax21 (activator of XA21-mediated immunity). Here we show that the XA21 juxtamembrane (JM) domain is required for kinase autophosphorylation. Threonine 705 in the XA21 JM domain is essential for XA21 autophosphorylation in vitro and XA21-mediated innate immunity in vivo. The replacement of Thr(705) by an alanine or glutamic acid abolishes XA21 autophosphorylation and eliminates interactions between XA21 and four XA21-binding proteins in yeast and rice. Although threonine residues analogous to Thr(705) of XA21 are present in the JM domains of most RD and non-RD plant receptor-like kinases, this residue is not required for autophosphorylation of the Arabidopsis RD RLK BRI1 (brassinosteroid insensitive 1). The threonine 705 of XA21 is conserved only in the JM domains of plant RLKs but not in those of fly, human, or mouse suggesting distinct regulatory mechanisms. These results contribute to growing knowledge regarding the mechanism by which non-RD RLKs function in plant.

Inference in signal transduction pathways using em algorithm and an implicit algorithm: incomplete data case

We summarize here the Implicit statistical inference approach as an alternative to Bayesian networks and we give an effective iterative algorithm analogous to the Expectation Maximization algorithm to infer signal transduction network when the set of data is incomplete. We proved the convergence of our algorithm that we called Implicit algorithm and we apply it to simulated data for a simplified signal transduction pathway of the EGFR protein.

Applications of computational tools to predict functional SNPs effects in human ErbB genes

Understanding the functions of single nucleotide polymorphisms (SNPs) can greatly help to understand the genetic basis of human complex diseases such as cancer. However, identifying functional SNPs among the huge number of available SNPs is challenging. In this study, we analyzed the genetic variations that can alter the expression and function of ErbB proteins using different computational tools. For noncoding SNP, we found that one SNP located in 59UTR of ErbB1 gene might change protein expression level and two SNPS located in regulatory regions might affect transcriptional regulation of Erbb1 and Erbb4. For coding SNPs we predicted that 25 nonsynonymous SNPs (most of them in ErbB1 gene) might disrupt the protein function among which 22 might alter protein structure. Prediction regarding the potential effect of the SNPs showed that 13 of them located within the tyrosine kinase or the ligand binding domain are likely to be associated with cancer.

Calmodulin-mediated regulation of the epidermal growth factor receptor

In this review, we first describe the mechanisms by which the epidermal growth factor receptor generates a Ca(2+) signal and, subsequently, we compile the available experimental evidence regarding the role that the Ca(2+)/calmodulin complex, formed after the rise in cytosolic free Ca(2+) concentration, exerts on the receptor. We focus not only on the indirect action that Ca(2+)/calmodulin exerts on the epidermal growth factor receptor, as a result of the activation of distinct calmodulin-dependent kinases, but also, and more extensively, on the direct interaction of Ca(2+)/calmodulin with the receptor. We also describe several mechanistic models that could account for the Ca(2+)/calmodulin-mediated regulation of epidermal growth factor receptor activity. The control exerted by calmodulin on distinct epidermal growth factor receptor-mediated cellular functions is also discussed. Finally, the phosphorylation of this Ca(2+) sensor by the epidermal growth factor receptor is highlighted.

The intracellular juxtamembrane domain of the epidermal growth factor (EGF) receptor is responsible for the allosteric regulation of EGF binding

We have previously shown that the binding of epidermal growth factor (EGF) to its receptor can best be described by a model that involves negative cooperativity in an aggregating system (Macdonald, J. L., and Pike, L. J. (2008) Proc. Natl. Acad. Sci. U. S. A. 105, 112-117). However, despite the fact that biochemical analyses indicate that EGF induces dimerization of its receptor, the binding data provided no evidence for positive linkage between EGF binding and dimer assembly. By analyzing the binding of EGF to a number of receptor mutants, we now report that in naive, unphosphorylated EGF receptors, ligand binding is positively linked to receptor dimerization but the linkage is abolished upon autophosphorylation of the receptor. Both phosphorylated and unphosphorylated EGF receptors exhibit negative cooperativity, indicating that mechanistically, cooperativity is distinct from the phenomenon of linkage. Nonetheless, both the positive linkage and the negative cooperativity observed in EGF binding require the presence of the intracellular juxtamembrane domain. This indicates the existence of inside-out signaling in the EGF receptor system. The intracellular juxtamembrane domain has previously been shown to be required for the activation of the EGF receptor tyrosine kinase (Thiel, K. W., and Carpenter, G. (2007) Proc. Natl. Acad. Sci. U. S. A. 104, 19238-19243). Our experiments expand the role of this domain to include the allosteric control of ligand binding by the extracellular domain.

Insight into the recognition patterns of the ErbB receptor family transmembrane domains: heterodimerization models through molecular dynamics search

ErbB receptors undergo a complex interaction network defining hierarchical and competition relationships. Dimerization is driven entirely by receptor-receptor interactions and the transmembrane domains play a role in modulating the specificity and the selection of the partners during signal transduction. To shed light on the role of the GxxxG-like dimerization motifs in the formation of ErbB transmembrane heterodimers, we propose structural models resulting from conformational search method combined with molecular dynamics simulations. Left-handed structures of the transmembrane heterodimers are found preponderant over right-handed structures. All heterotypic heterodimers undergo two modes of association either via the N-terminal motif or the C-terminal motif. The transmembrane domain of ErbB3 impairs this C-terminal motif but also associates with the other partners owing to the presence of Gly residues. The two dimerization modes involve different orientations of the two helices. Thus, a molecular-switch model allowing the transition between the two dimerizing states may apply to the heterodimers and could help interpret receptor competition for the formation of homodimers and heterodimers. The comparison between experimental and theoretical results on the dimerization hierarchy of the transmembrane domains is not straightforward. However, we demonstrate that the intrinsic properties of the transmembrane sequences are an important component in heterodimer formation and that the ErbB2 and ErbB3 transmembrane domains have a strong power for heterodimerization as observed experimentally.

Epidermal growth factor receptor juxtamembrane region regulates allosteric tyrosine kinase activation

Structural studies of the extracellular and tyrosine kinase domains of the epidermal growth factor receptor (ErbB-1) provide considerable insight into facets of the receptor activation mechanism, but the contributions of other regions of ErbB-1 have not been ascertained. This study demonstrates that the intracellular juxtamembrane (JM) region plays a vital role in the kinase activation mechanism. In the experiments described herein, the entire ErbB-1 intracellular domain (ICD) has been expressed in mammalian cells to explore the significance of the JM region in kinase activity. Deletion of the JM region (DeltaJM) results in a severe loss of ICD tyrosine phosphorylation, indicating that this region is required for maximal activity of the tyrosine kinase domain. Coexpression of DeltaJM and dimerization-deficient kinase domain ICD mutants revealed that the JM region is indispensable for allosteric kinase activation and productive monomer interactions within a dimer. Studies with the intact receptor confirmed the role of the JM region in kinase activation. Within the JM region, Thr-654 is a known protein kinase C (PKC) phosphorylation site that modulates kinase activity in the context of the intact ErbB-1 receptor; yet, the mechanism is not known. Whereas a T654A mutation promotes increased ICD tyrosine phosphorylation, the phosphomimetic T654D mutant generates a 50% reduction in ICD tyrosine phosphorylation. Similar to the DeltaJM mutants, the T654D mutant ICD failed to interact with a wild-type monomer. This study reveals an integral role for the intracellular JM region of ErbB-1 in allosteric kinase activation.

Membrane-permeable calmodulin inhibitors (e.g. W-7/W-13) bind to membranes, changing the electrostatic surface potential: dual effect of W-13 on epidermal growth factor receptor activation

Membrane-permeable calmodulin inhibitors, such as the napthalenesulfonamide derivatives W-7/W-13, trifluoperazine, and calmidazolium, are used widely to investigate the role of calcium/calmodulin (Ca2+/CaM) in living cells. If two chemically different inhibitors (e.g. W-7 and trifluoperazine) produce similar effects, investigators often assume the effects are due to CaM inhibition. Zeta potential measurements, however, show that these amphipathic weak bases bind to phospholipid vesicles at the same concentrations as they inhibit Ca2+/CaM; this suggests that they also bind to the inner leaflet of the plasma membrane, reducing its negative electrostatic surface potential. This change will cause electrostatically bound clusters of basic residues on peripheral (e.g. Src and K-Ras4B) and integral (e.g. epidermal growth factor receptor (EGFR)) proteins to translocate from the membrane to the cytoplasm. We measured inhibitor-mediated translocation of a simple basic peptide corresponding to the calmodulin-binding juxtamembrane region of the EGFR on model membranes; W-7/W-13 causes translocation of this peptide from membrane to solution, suggesting that caution must be exercised when interpreting the results obtained with these inhibitors in living cells. We present evidence that they exert dual effects on autophosphorylation of EGFR; W-13 inhibits epidermal growth factor-dependent EGFR autophosphorylation under different experimental conditions, but in the absence of epidermal growth factor, W-13 stimulates autophosphorylation of the receptor in four different cell types. Our interpretation is that the former effect is due to W-13 inhibition of Ca2+/CaM, but the latter results could be due to binding of W-13 to the plasma membrane.

Loss of Hsp90 association up-regulates Src-dependent ErbB2 activity

The receptor tyrosine kinase ErbB2 plays a crucial role in tumorigenesis. We showed previously that the molecular chaperone Hsp90 protects ErbB2 from proteasome-mediated degradation by binding to a short loop structure in the N-lobe of the kinase domain. Here we show that loss of Hsp90 binding correlates with enhanced ErbB2 kinase activity and its transactivating potential, concomitant with constitutively increased phosphorylation of Tyr877, located in the activation loop of the kinase domain. We show further that Tyr877 phosphorylation is mediated by Src and that it is necessary for the enhanced kinase activity of ErbB2. Finally, computer modeling of the kinase domain suggests a phosphorylation-dependent reorientation of the activation loop, denoting the importance of Tyr877 phosphorylation for ErbB2 activity. These findings suggest that Hsp90 binding to ErbB2 participates in regulation of kinase activity as well as kinase stability.

Phosphorylation of Thr654 but not Thr669 within the juxtamembrane domain of the EGF receptor inhibits calmodulin binding

Calcium-calmodulin (CaM) binding to the epidermal growth factor receptor (EGFR) has been shown to both inhibit and stimulate receptor activity. CaM binds to the intracellular juxtamembrane (JM) domain (Met645-Phe688) of EGFR. Protein kinase C (PKC) mediated phosphorylation of Thr654 occurs within this domain. CaM binding to the JM domain inhibits PKC phosphorylation and conversely PKC mediated phosphorylation of Thr654 or Glu substitution of Thr654 inhibits CaM binding. A second threonine residue (Thr669) within the JM domain is phosphorylated by the mitogen-activated protein kinase (MAPK). Previous results have shown that CaM interferes with EGFR-induced MAPK activation. If and how phosphorylation of Thr669 affects CaM-EGFR interaction is however not known. In the present study we have used surface plasmon resonance (BIAcore) to study the influence of Thr669 phosphorylation on real time interactions between the intracellular juxtamembrane (JM) domain of EGFR and CaM. The EGFR-JM was expressed as GST fusion proteins in Escherichia coli and phosphorylation was mimicked by generating Glu substitutions of either Thr654 or Thr669. Purified proteins were coupled to immobilized anti-GST antibodies at the sensor surface and increasing concentration of CaM was applied. When mutating Thr654 to Glu654 no specific CaM binding could be detected. However, neither single substitutions of Thr669 (Gly669 or Glu669) nor double mutants Gly654/Gly669 or Gly654/Glu669 influenced the binding of CaM to the EGFR-JM. This clearly shows that PKC may regulate EGF-mediated CaM signalling through phosphorylation of Thr654 whereas phosphorylation of Thr669 seems to play a CaM independent regulatory role. The role of both residues in the EGFR-calmodulin interaction was also studied in silico. Our modelling work supports a scenario where Thr654 from the JM domain interacts with Glu120 in the calmodulin molecule. Phosphorylation of Thr654 or Glu654 substitution creates a repulsive electrostatic force that would diminish CaM binding to the JM domain. These results are in line with the Biacore experiments showing a weak binding of the CaM to the JM domain with Thr654 mutated to Glu. Furthermore, these results provide a hypothesis to how CaM binding to EGFR might both positively and negatively interfere with EGFR-activity.