EGFR: tale of the C-terminal tail

The combined action of the intracellular regions regulates FGFR2 kinase activity

Receptor tyrosine kinases (RTKs) are typically activated through a precise sequence of intracellular phosphorylation events starting with a tyrosine residue on the activation loop (A-loop) of the kinase domain (KD). From this point the mono-phosphorylated enzyme is active, but subject to stringent regulatory mechanisms which can vary dramatically across the different RTKs. In the absence of extracellular stimulation, fibroblast growth factor receptor 2 (FGFR2) exists in the mono-phosphorylated state in which catalytic activity is regulated to allow rapid response upon ligand binding, whilst restricting ligand-independent activation. Failure of this regulation is responsible for pathologic outcomes including cancer. Here we reveal the molecular mechanistic detail of KD control based on combinatorial interactions of the juxtamembrane (JM) and the C-terminal tail (CT) regions of the receptor. JM stabilizes the asymmetric dimeric KD required for substrate phosphorylation, whilst CT binding opposes dimerization, and down-regulates activity. Direct binding between JM and CT delays the recruitment of downstream effector proteins adding a further control step as the receptor proceeds to full activation. Our findings underscore the diversity in mechanisms of RTK oligomerisation and activation.

Structural Insight and Development of EGFR Tyrosine Kinase Inhibitors

Lung cancer has a high prevalence, with a growing number of new cases and mortality every year. Furthermore, the survival rate of patients with non-small-cell lung carcinoma (NSCLC) is still quite low in the majority of cases. Despite the use of conventional therapy such as tyrosine kinase inhibitor for Epidermal Growth Factor Receptor (EGFR), which is highly expressed in most NSCLC cases, there was still no substantial improvement in patient survival. This is due to the drug’s ineffectiveness and high rate of resistance among individuals with mutant EGFR. Therefore, the development of new inhibitors is urgently needed. Understanding the EGFR structure, including its kinase domain and other parts of the protein, and its activation mechanism can accelerate the discovery of novel compounds targeting this protein. This study described the structure of the extracellular, transmembrane, and intracellular domains of EGFR. This was carried out along with identifying the binding pose of commercially available inhibitors in the ATP-binding and allosteric sites, thereby clarifying the research gaps that can be filled. The binding mechanism of inhibitors that have been used clinically was also explained, thereby aiding the structure-based development of new drugs.

Structure and activity of an anti-epidermal growth factor receptor antibody without galactose-α-1,3-galactose residues

Epidermal growth factor receptor (EGFR) is a transmembrane glycoprotein, which has been considered as one of the key targets for cancer therapy. However, currently approved therapeutic anti-EGFR antibody may cause the hypersensitivity reaction induced by galactose-α-1,3-galactose (α-Gal) structure, which is inevitable in insect cell expression system. In this study, the Chinese hamster ovary cell line was used to produce a monoclonal antibody containing simplified glycosylation patterns (code: AB01). And cetuximab was used as a control. The two antibodies were highly similar in molecular weight, secondary structure, binding affinity and endocytosis behavior, whereas the glycotypes are extremely distinct. The flow cytometry assay suggested that AB01 induced cell cycle arrest in G1, thus inhibit cell proliferation. Moreover, both cetuximab and AB01 showed similar sensitivity for all tested cell lines in this research. In conclusion, AB01 could be a potential anti-EGFR drug candidate for cancer therapy.

Phase 2 study of afatinib among patients with recurrent and/or metastatic esophageal squamous cell carcinoma

BACKGROUND

The objective of the current study was to investigate the clinical activity of, safety of, and predictive biomarkers for afatinib, an irreversible pan-ErbB kinase inhibitor, in patients with recurrent and/or metastatic esophageal squamous cell carcinoma (R/M-ESCC).

METHODS

Patients with R/M-ESCC that was refractory to platinum-based chemotherapy were enrolled in the current multicenter, single-arm, phase 2 study and received afatinib at a dose of 40 mg/day. The primary endpoint was the objective response rate. Secondary endpoints included progression-free survival, overall survival, the disease control rate, and the safety profile. To identify predictive biomarkers, single-nucleotide variations, short insertions/deletions, and somatic copy number alterations were assessed using whole-exome sequencing and their associations with clinical outcomes were analyzed.

RESULTS

Among 49 enrolled patients, the objective response rate and disease control rate were 14.3% and 73.3%, respectively. With a median follow-up of 6.6 months, the median progression-free survival and overall survival were 3.4 months and 6.3 months, respectively. Treatment-related adverse events were noted to have occurred in 33 patients (67.3%), with the majority being of grade 1 to 2 (adverse events were graded and recorded based on the National Cancer Institute Common Terminology Criteria for Adverse Events [version 4.03]). Whole-exome sequencing demonstrated that the ESCC genomes of patients who demonstrated a response to afatinib were enriched with genomic alterations of TP53 and epidermal growth factor receptor (EGFR). As a predictive marker, a score derived from TP53 disruptive mutations and EGFR amplifications and/or missense mutations demonstrated a significant association with the response to afatinib. The score based on the mutational status of EGFR and TP53 achieved a performance of an area under the curve of 0.86 in predicting the sensitivity of afatinib.

CONCLUSIONS

The results of the current study demonstrated that afatinib can confer modest clinical benefits with manageable toxicity in patients with platinum-resistant R/M-ESCC. Identification of TP53 alterations and EGFR amplifications may serve as predictive markers with which to identify patients with R/M-ESCC who may benefit from afatinib.

LAY SUMMARY

Esophageal squamous cell carcinoma (ESCC) is a type of cancer with a dismal prognosis and very limited treatment options. The clinical efficacy of afatinib was evaluated in patients with recurrent and/or metastatic ESCC, with adverse events demonstrating the modest efficacy with manageable toxicity of this irreversible, pan-ErbB kinase inhibitor. Whole-exome sequencing analysis of 41 cases of ESCC further revealed that the patients harboring epidermal growth factor receptor (EGFR) amplifications and disruptive TP53 mutations are more likely to benefit from treatment with afatinib. The results of the current study have highlighted the clinical value of EGFR and TP53 as predictive biomarkers of platinum-resistant recurrent and/or metastatic ESCC for afatinib sensitivity.

Preparation of a phosphotyrosine-protein standard for use in semiquantitative western blotting with enhanced chemiluminescence

Protein tyrosine phosphorylation is key to activation of receptor tyrosine kinases (RTK) that drive development of some cancers. One challenge of RTK-targeted therapy is identification of those tumors that express non-mutated but activated RTKs. Phosphotyrosine (pTyr) RTK levels should be more predictive of the latter than expressed total protein. Western blotting (WB) with a pTyr antibody and enhanced chemiluminescence (ECL) detection is sufficiently sensitive to detect pTyr-RTKs in human tumor homogenates. Presentation of results by comparing WB images, however, is wanting. Here we describe the preparation of a new pTyr-protein standard, pTyr-ALK48-SB (pA), derived from a commercial anaplastic lymphoma kinase (ALK) recombinant fragment, and its use to quantify pTyr-epidermal growth factor receptor (pTyr-EGFR) in commercial A431 cell lysates. Linearity of one-dimensional (1D) WB plots of pA band density versus load as well as its lower level of detection (0.1 ng, 2 fmole) were determined for standardized conditions. Adding pA to two lots of A431 cell lysates with high and low pTyr-EGFR allowed normalization and quantification of the latter by expressing results as density ratios for both 1D and 2D WB. This approach is semi-quantitative because unknown RTKs may be outside the linear range of detection. Semiquantitative ratios are an improvement over comparisons of images without a reference standard and facilitate comparisons between samples.

Molecular dynamics, MM/PBSA and in vitro validation of a novel quinazoline-based EGFR tyrosine kinase inhibitor identified using structure-based in silico screening

EGFR-TK has been a target strongly associated with the development of NSCLCs. A structure-based virtual screening campaign was launched against EGFR-TK by virtual screening a 3D library of 167 commercially available small molecules downloaded from ChemBridge Corporation. The virtual screen identified 12 virtual hit molecules, which were biologically evaluated against an EGFR-TK inhibitor-sensitive NSCLC cell line, A549. A quinazoline-based molecule 1, was most active and displayed ∼58% cytotoxicity at 20 μM single dose. The mode of cell death suggests molecule 1 induced apoptosis, which is characteristic of EGFR-TK pathway inhibition. A 50 ns MD simulation was conducted on three different systems: free EGFR-TK, molecule 1 complexed to EGFR-TK, and the positive control, lapatinib, complexed to EGFR-TK. The MD simulations showed increase in stabilisation of the EGFR-TK structure for the complexed systems, i.e., lower RMSDs and RMSFs for complexed EGFR-TK structures compared to the free EGFR-TK system. The binding affinities were estimated using MM/PBSA in the last 10 ns of the MD simulation that revealed comparable binding free energies between molecule 1 and lapatinib, ΔG_bind = -25.0 and -23.9 kcal/mol, respectively. Per residue binding free energy decomposition studies revealed non-polar interactions contributed mostly to the binding free energies. Residues Leu718, Arg841 and Phe856 were predicted to contribute most to the binding free energies for molecule 1.

Detection of rare and novel EGFR mutations in NSCLC patients: Implications for treatment-decision

OBJECTIVES

Mutations in the gene that encodes epidermal growth factor receptor (EGFR) are biomarkers that predict how non-small cell lung cancer (NSCLC) patients respond to EGFR-targeted therapies collectively known as tyrosine kinase inhibitors (TKIs). Thus, EGFR genotyping provides crucial information for treatment decision. Both Sanger sequencing and real-time PCR methodologies are used for EGFR genotyping. However, methods based on real-time PCR have limitations, as they may not detect rare or novel mutations. The aim of this study was to determine the prevalence of rare mutations in the tyrosine kinase domain (exons 18-21) of the EGFR gene not targeted by the most frequently used real-time PCR approaches, i.e., the cobas® EGFR Mutation Test, and the Idylla™ EGFR Mutation Assay.

METHODS

A total of 1228 NSCLC patients were screened for mutations in exons 18-21 of the EGFR gene using Sanger sequencing.

RESULTS

We observed that 252 patients (∼20%) had at least one mutation in the EGFR gene, and 38 (∼3%) carried uncommon genetic alterations that would not be identified by the cobas® or the Idylla™ tests. We further found six new single mutations and seven previously unreported compound mutations. Clinical information and patient outcome are presented for these cases.

CONCLUSIONS

This study highlights the value of sequencing-based approaches to identify rare mutations. Our results add to the inventory of known EGFR mutations, thus contributing to improved lung cancer precision treatment.

Type II cGMP-dependent protein kinase phosphorylates EGFR at threonine 669 and thereby inhibits its activation

Our previous study demonstrated that type II cGMP-dependent protein kinase (PKG II) inhibited epidermal growth factor (EGF) induced tyrosine phosphorylation/activation of the EGF receptor (EGFR). This paper was designed to investigate the mechanism of the inhibition of PKG II on EGFR activation. Gastric cancer cells HGC-27 and AGS were infected with an adenoviral vector encoding the cDNA of PKG II (Ad-PKG II) to overexpress PKG II and treated with 8-(4-chlorophenylthio) guanosine-3',5'-cyclic monophosphate (8-pCPT-cGMP) to activate the kinase. Co-immunoprecipitation (Co-IP) and bimolecular fluorescence complementation (BiFC) assay were performed to detect the interaction between PKG II and EGFR. Western blotting, mass spectrometry (MS) and site mutagenesis were performed to detect the PKG II-specific phosphorylation site on EGFR. The results showed that in living COS-7 cells, which were infected with Ad-PKG II and treated with 8-pCPT-cGMP, there was an interaction between PKG II and EGFR. The results also showed that PKG II caused threonine 669 (T669) phosphorylation of EGFR in HGC-27 and AGS cells infected with Ad-PKG II and treated with 8-pCPT-cGMP, and then inhibited the activation of EGFR. When T669 of EGFR was mutated to alanine, the inhibitory effect of PKG II on the activation of EGFR was eradicated. These findings suggested a PKG II-specific phosphorylation site on EGFR, and might be beneficial to illuminate the anti-tumor role of PKG II.

Tracing the origin and evolution of pseudokinases across the tree of life

Protein phosphorylation by eukaryotic protein kinases (ePKs) is a fundamental mechanism of cell signaling in all organisms. In model vertebrates, ~10% of ePKs are classified as pseudokinases, which have amino acid changes within the catalytic machinery of the kinase domain that distinguish them from their canonical kinase counterparts. However, pseudokinases still regulate various signaling pathways, usually doing so in the absence of their own catalytic output. To investigate the prevalence, evolutionary relationships, and biological diversity of these pseudoenzymes, we performed a comprehensive analysis of putative pseudokinase sequences in available eukaryotic, bacterial, and archaeal proteomes. We found that pseudokinases are present across all domains of life, and we classified nearly 30,000 eukaryotic, 1500 bacterial, and 20 archaeal pseudokinase sequences into 86 pseudokinase families, including ~30 families that were previously unknown. We uncovered a rich variety of pseudokinases with notable expansions not only in animals but also in plants, fungi, and bacteria, where pseudokinases have previously received cursory attention. These expansions are accompanied by domain shuffling, which suggests roles for pseudokinases in plant innate immunity, plant-fungal interactions, and bacterial signaling. Mechanistically, the ancestral kinase fold has diverged in many distinct ways through the enrichment of unique sequence motifs to generate new families of pseudokinases in which the kinase domain is repurposed for noncanonical nucleotide binding or to stabilize unique, inactive kinase conformations. We further provide a collection of annotated pseudokinase sequences in the Protein Kinase Ontology (ProKinO) as a new mineable resource for the signaling community.

Structure and Dynamics of the EGF Receptor as Revealed by Experiments and Simulations and Its Relevance to Non-Small Cell Lung Cancer

The epidermal growth factor receptor (EGFR) is historically the prototypical receptor tyrosine kinase, being the first cloned and the first where the importance of ligand-induced dimer activation was ascertained. However, many years of structure determination has shown that EGFR is not completely understood. One challenge is that the many structure fragments stored at the PDB only provide a partial view because full-length proteins are flexible entities and dynamics play a key role in their functionality. Another challenge is the shortage of high-resolution data on functionally important higher-order complexes. Still, the interest in the structure/function relationships of EGFR remains unabated because of the crucial role played by oncogenic EGFR mutants in driving non-small cell lung cancer (NSCLC). Despite targeted therapies against EGFR setting a milestone in the treatment of this disease, ubiquitous drug resistance inevitably emerges after one year or so of treatment. The magnitude of the challenge has inspired novel strategies. Among these, the combination of multi-disciplinary experiments and molecular dynamic (MD) simulations have been pivotal in revealing the basic nature of EGFR monomers, dimers and multimers, and the structure-function relationships that underpin the mechanisms by which EGFR dysregulation contributes to the onset of NSCLC and resistance to treatment.

Reading Conformational Changes in Proteins with a New Colloidal-Based Interfacial Biosensing System

Many biological events such as mutations or aberrant post-translational modifications can alter the conformation and/or folding stability of proteins and their subsequent biological function, which may trigger the onset of diseases like cancer. Evaluating protein folding is hence crucial for the diagnosis of these diseases. Yet, it is still challenging to detect changes in protein folding, especially if they are subtle, in a simple and highly sensitive manner with the current assays. Herein, we report a new colloidal-based interfacial biosensing approach for qualitative and quantitative profiling of various types of changes in protein folding; from denaturation to variant conformations in native proteins, such as protein activation via mutations or phosphorylation. The approach is based on the direct interfacial interaction of proteins freely available in solution with added tannic-acid-capped gold nanoparticles, to interrogate their folding status in their solubilized form. We found that under the optimized conditions, proteins can modulate colloids solvation according to their folding or conformational status, which can be visualized in a single step, by the naked eye, with minimal protein input requirements (limit of detection of 1 ng/μL). Protein folding detection was achieved regardless of protein topology and size without using conformation-specific antibodies and mutational analysis, which are the most common assays for sensing malfunctioning proteins. The approach showed excellent sensitivity, superior to circular dichroism, for the detection of the very subtle conformational changes induced by activating mutations and phosphorylation in epidermal growth factor receptor (EGFR) and extracellular signal-regulated kinase (ERK) proteins. This enabled their detection even in complex samples derived from lung cancer cells, which contained up to 95% excess of their wild-type forms. A broader clinical translation was shown via monitoring the action of conformation-restoring drugs, such as tyrosine kinase inhibitors, on EGFR conformation and its downstream protein network, using the ERK protein as a surrogate.

Coupled regulation by the juxtamembrane and sterile α motif (SAM) linker is a hallmark of ephrin tyrosine kinase evolution

Ephrin (Eph) receptor tyrosine kinases have evolutionarily diverged from other tyrosine kinases to respond to specific activation and regulatory signals that require close coupling of kinase catalytic and regulatory functions. However, the evolutionary basis for such functional coupling is not fully understood. We employed an evolutionary systems approach involving statistical mining of large sequence and structural data sets to define the hallmarks of Eph kinase evolution and functional specialization. We found that some of the most distinguishing Eph-specific residues structurally tether the flanking juxtamembrane and sterile α motif (SAM) linker regions to the kinase domain, and substitutions of these residues in EphA3 resulted in faster kinase activation. We report for the first time that the SAM domain linker is functionally coupled to the juxtamembrane through co-conserved residues in the kinase domain and that together these residues provide a structural framework for coupling catalytic and regulatory functions. The unique organization of Eph-specific tethering networks and the identification of other Eph-specific sequence features of unknown functions provide new hypotheses for future functional studies and new clues to disease mutations altering Eph kinase-specific functions.

Biophysical Evidence for Intrinsic Disorder in the C-terminal Tails of the Epidermal Growth Factor Receptor (EGFR) and HER3 Receptor Tyrosine Kinases

The epidermal growth factor receptor (EGFR)/ErbB family of receptor tyrosine kinases includes oncogenes important in the progression of breast and other cancers, and they are targets for many drug development strategies. Each member of the ErbB family possesses a unique, structurally uncharacterized C-terminal tail that plays an important role in autophosphorylation and signal propagation. To determine whether these C-terminal tails are intrinsically disordered regions, we conducted a battery of biophysical experiments on the EGFR and HER3 tails. Using hydrogen/deuterium exchange mass spectrometry, we measured the conformational dynamics of intracellular half constructs and compared the tails with the ordered kinase domains. The C-terminal tails demonstrate more rapid deuterium exchange behavior when compared with the kinase domains. Next, we expressed and purified EGFR and HER3 tail-only constructs. Results from circular dichroism spectroscopy, size exclusion chromatography with multiangle light scattering, dynamic light scattering, analytical ultracentrifugation, and small angle X-ray scattering each provide evidence that the EGFR and HER3 C-terminal tails are intrinsically disordered with extended, non-globular structure in solution. The intrinsic disorder and extended conformation of these tails may be important for their function by increasing the capture radius and reducing the thermodynamic barriers for binding of downstream signaling proteins.

Constitutive asymmetric dimerization drives oncogenic activation of epidermal growth factor receptor carboxyl-terminal deletion mutants

Genomic alterations targeting the Epidermal Growth Factor Receptor (EGFR) gene have been strongly associated with cancer pathogenesis. The clinical effectiveness of EGFR targeted therapies, including small molecules directed against the kinase domain such as gefitinib, erlotinib and afatinib, have been proven successful in treating non-small cell lung cancer patients with tumors harboring EGFR kinase domain mutations. Recent large-scale genomic studies in glioblastoma and lung cancer have identified an additional class of oncogenic mutations caused by the intragenic deletion of carboxy-terminal coding regions. Here, we report that combinations of exonic deletions of exon 25 to 28 lead to the oncogenic activation of EGF receptor in the absence of ligand and consequent cellular transformation, indicating a significant role of C-terminal domain in modulating EGFR activation. Furthermore, we show that the oncogenic activity of the resulting C-terminal deletion mutants are efficiently inhibited by EGFR-targeted drugs including erlotinib, afatinib, dacomitinib as well as cetuximab, expanding the therapeutic rationale of cancer genome-based EGFR targeted approaches. Finally, in vivo and in vitro preclinical studies demonstrate that constitutive asymmetric dimerization in mutant EGFR is a key mechanism for oncogenic activation and tumorigenesis by C-terminal deletion mutants. Therefore, our data provide compelling evidence for oncogenic activation of C-terminal deletion mutants at the molecular level and we propose that C-terminal deletion status of EGFR can be considered as a potential genomic marker for EGFR-targeted therapy.

A tale of the epidermal growth factor receptor: The quest for structural resolution on cells

The challenge of determining the architecture and geometry of oligomers of the epidermal growth factor receptor (EGFR) on the cell surface has been approached using a variety of biochemical and biophysical methods. This review is intended to provide a narrative of how key concepts in the field of EGFR research have evolved over the years, from the origins of the prevalent EGFR signalling dimer hypothesis through to the development and implementation of methods that are now challenging the conventional view. The synergy between X-ray crystallography and cellular fluorescence microscopy has become particularly important, precisely because the results from these two methods diverged and highlighted the complexity of the challenge. We illustrate how developments in super-resolution microscopy are now bridging this gap. Exciting times lie ahead where knowledge of the nature of the complexes can assist with the development of a new generation of anti-cancer drugs.

The EGFR family: not so prototypical receptor tyrosine kinases

The epidermal growth factor receptor (EGFR) was among the first receptor tyrosine kinases (RTKs) for which ligand binding was studied and for which the importance of ligand-induced dimerization was established. As a result, EGFR and its relatives have frequently been termed "prototypical" RTKs. Many years of mechanistic studies, however, have revealed that--far from being prototypical--the EGFR family is quite unique. As we discuss in this review, the EGFR family uses a distinctive "receptor-mediated" dimerization mechanism, with ligand binding inducing a dramatic conformational change that exposes a dimerization arm. Intracellular kinase domain regulation in this family is also unique, being driven by allosteric changes induced by asymmetric dimer formation rather than the more typical activation-loop phosphorylation. EGFR family members also distinguish themselves from other RTKs in having an intracellular juxtamembrane (JM) domain that activates (rather than autoinhibits) the receptor and a very large carboxy-terminal tail that contains autophosphorylation sites and serves an autoregulatory function. We discuss recent advances in mechanistic aspects of all of these components of EGFR family members, attempting to integrate them into a view of how RTKs in this important class are regulated at the cell surface.

Coarse-grained molecular simulation of epidermal growth factor receptor protein tyrosine kinase multi-site self-phosphorylation

Upon the ligand-dependent dimerization of the epidermal growth factor receptor (EGFR), the intrinsic protein tyrosine kinase (PTK) activity of one receptor monomer is activated, and the dimeric receptor undergoes self-phosphorylation at any of eight candidate phosphorylation sites (P-sites) in either of the two C-terminal (CT) domains. While the structures of the extracellular ligand binding and intracellular PTK domains are known, that of the ∼225-amino acid CT domain is not, presumably because it is disordered. Receptor phosphorylation on CT domain P-sites is critical in signaling because of the binding of specific signaling effector molecules to individual phosphorylated P-sites. To investigate how the combination of conventional substrate recognition and the unique topological factors involved in the CT domain self-phosphorylation reaction lead to selectivity in P-site phosphorylation, we performed coarse-grained molecular simulations of the P-site/catalytic site binding reactions that precede EGFR self-phosphorylation events. Our results indicate that self-phosphorylation of the dimeric EGFR, although generally believed to occur in trans, may well occur with a similar efficiency in cis, with the P-sites of both receptor monomers being phosphorylated to a similar extent. An exception was the case of the most kinase-proximal P-site-992, the catalytic site binding of which occurred exclusively in cis via an intramolecular reaction. We discovered that the in cis interaction of P-site-992 with the catalytic site was facilitated by a cleft between the N-terminal and C-terminal lobes of the PTK domain that allows the short CT domain sequence tethering P-site-992 to the PTK core to reach the catalytic site. Our work provides several new mechanistic insights into the EGFR self-phosphorylation reaction, and demonstrates the potential of coarse-grained molecular simulation approaches for investigating the complexities of self-phosphorylation in molecules such as EGFR (HER/ErbB) family receptors and growth factor receptor PTKs in general.

The epidermal growth factor receptor (EGFR) is one of a large group of cell surface receptors that allow cells to respond to growth-stimulating signals in their environment. Upon sensing of growth factor, the EGFR is activated, which triggers a signaling cascade leading to the cell nucleus and ultimately initiating cell division. The first event following receptor activation is an intramolecular kinase reaction that results in the introduction of phosphate groups onto several specific amino acids (phosphorylation sites or P-sites) in the tail of the EGFR protein. Thus, the tail of the receptor undergoes self-phosphorylation, which involves conformational motions enabling the various P-sites to access the catalytic site. The structure of the tail of the receptor is unknown, and hence the mechanism of the self-phosphorylation reaction is not well understood. To investigate this mechanism, we generated a structural model of the EGFR protein and performed computer simulations of EGFR P-site/catalytic site binding reactions. These simulations indicated how the distribution of P-sites along the tail of the receptor and restrictions in molecular movements of the tail lead to selectivity in the phosphorylation of the different P-sites. Our simulations yielded unique insights into the mechanism of EGFR self-phosphorylation that have important biological implications.