Conformational regulation of the EGFR kinase core by the juxtamembrane and C-terminal tail: A molecular dynamics study

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.

Expanding the Disorder-Function Paradigm in the C-Terminal Tails of Erbbs

ErbBs are receptor tyrosine kinases involved not only in development, but also in a wide variety of diseases, particularly cancer. Their extracellular, transmembrane, juxtamembrane, and kinase folded domains were described extensively over the past 20 years, structurally and functionally. However, their whole C-terminal tails (CTs) following the kinase domain were only described at atomic resolution in the last 4 years. They were shown to be intrinsically disordered. The CTs are known to be tyrosine-phosphorylated when the activated homo- or hetero-dimers of ErbBs are formed. Their phosphorylation triggers interaction with phosphotyrosine binding (PTB) or Src Homology 2 (SH2) domains and activates several signaling pathways controling cellular motility, proliferation, adhesion, and apoptosis. Beyond this passive role of phosphorylated domain and site display for partners, recent structural and function studies unveiled active roles in regulation of phosphorylation and interaction: the CT regulates activity of the kinase domain; different phosphorylation states have different compaction levels, potentially modulating the succession of phosphorylation events; and prolines have an important role in structure, dynamics, and possibly regulatory interactions. Here, we review both the canonical role of the disordered CT domains of ErbBs as phosphotyrosine display domains and the recent findings that expand the known range of their regulation functions linked to specific structural and dynamic features.

Anti-EGFR Binding Nanobody Delivery System to Improve the Diagnosis and Treatment of Solid Tumours

BACKGROUND

Epidermal Growth Factor Receptor (EGFR) and members of its homologous protein family mediate transmembrane signal transduction by binding to a specific ligand, which leads to regulated cell growth, differentiation, proliferation and metastasis. With the development and application of Genetically Engineered Antibodies (GEAs), Nanobodies (Nbs) constitute a new research hot spot in many diseases. A Nb is characterized by its low molecular weight, deep tissue penetration, good solubility and high antigen-binding affinity, the anti-EGFR Nbs are of significance for the diagnosis and treatment of EGFR-positive tumours.

OBJECTIVE

This review aims to provide a comprehensive overview of the information about the molecular structure of EGFR and its transmembrane signal transduction mechanism, and discuss the anti-EGFR-Nbs influence on the diagnosis and treatment of solid tumours.

METHODS

Data were obtained from PubMed, Embase and Web of Science. All patents are searched from the following websites: the World Intellectual Property Organization (WIPO®), the United States Patent Trademark Office (USPTO®) and Google Patents.

RESULTS

EGFR is a key target for regulating transmembrane signaling. The anti-EGFR-Nbs for targeted drugs could effectively improve the diagnosis and treatment of solid tumours.

CONCLUSION

EGFR plays a role in transmembrane signal transduction. The Nbs, especially anti- EGFR-Nbs, have shown effectiveness in the diagnosis and treatment of solid tumours. How to increase the affinity of Nb and reduce its immunogenicity remain a great challenge.

Unveiling the Structural Insights into the Selective Inhibition of Protein Kinase D1

BACKGROUND

Although protein kinase D1 (PKD1) has been proved to be an efficient target for anticancer drug development, lack of structural details and substrate binding mechanisms are the main obstacles for the development of selective inhibitors with therapeutic benefits.

OBJECTIVE

The present study described the in silico dynamics behaviors of PKD1 in binding with selective and non-selective inhibitors and revealed the critical binding site residues for the selective kinase inhibition.

METHODS

Here, the three dimensional model of PKD1 was initially constructed by homology modeling along with binding site characterization to explore the non-conserved residues. Subsequently, two known inhibitors were docked to the catalytic site and the detailed ligand binding mechanisms and post binding dyanmics were investigated by molecular dynamics simulation and binding free energy calculations.

RESULTS

According to the binding site analysis, PKD1 serves several non-conserved residues in the G-loop, hinge and catalytic subunits. Among them, the residues including Leu662, His663, and Asp665 from hinge region made polar interactions with selective PKD1 inhibitor in docking simulation, which were further validated by the molecular dynamics simulation. Both inhibitors strongly influenced the structural dynamics of PKD1 and their computed binding free energies were in accordance with experimental bioactivity data.

CONCLUSION

The identified non-conserved residues likely to play critical role on molecular reorganization and inhibitor selectivity. Taken together, this study explained the molecular basis of PKD1 specific inhibition, which may help to design new selective inhibitors for better therapies to overcome cancer and PKD1 dysregulated disorders.

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.

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.

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

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

Cancer-associated arginine-to-histidine mutations confer a gain in pH sensing to mutant proteins

The intracellular pH (pHi) of most cancers is constitutively higher than that of normal cells and enhances proliferation and cell survival. We found that increased pHi enabled the tumorigenic behaviors caused by somatic arginine-to-histidine mutations, which are frequent in cancer and confer pH sensing not seen with wild-type proteins. Experimentally raising the pHi increased the activity of R776H mutant epidermal growth factor receptor (EGFR-R776H), thereby increasing proliferation and causing transformation in fibroblasts. An Arg-to-Gly mutation did not confer these effects. Molecular dynamics simulations of EGFR suggested that decreased protonation of His⁷⁷⁶ at high pH causes conformational changes in the αC helix that may stabilize the active form of the kinase. An Arg-to-His, but not Arg-to-Lys, mutation in the transcription factor p53 (p53-R273H) decreased its transcriptional activity and attenuated the DNA damage response in fibroblasts and breast cancer cells with high pHi. Lowering pHi attenuated the tumorigenic effects of both EGFR-R776H and p53-R273H. Our data suggest that some somatic mutations may confer a fitness advantage to the higher pHi of cancer cells.

Molecular Mechanism of Nucleotide-Dependent Allosteric Regulation in AMP-Activated Protein Kinase

The AMP-activated protein kinase (AMPK), a central enzyme in the regulation of energy homeostasis, is an important drug target for type 2 diabetes, obesity, and cancer. Binding of adenosine nucleotides to the regulatory γ-subunit tightly regulates the activity of this enzyme. Though recent crystal structures of AMPK have provided important insights into the allosteric activation of AMPK, molecular details of the regulatory mechanism of AMPK activation is still elusive. Here, we have performed extensive all-atom molecular dynamics (MD) simulations and shown that the kinase domain (KD) and γ-subunit come closer resulting in a more compact heterotrimeric AMPK complex in AMP-bound state compared to the ATP-bound state. The binding of ATP at site 3 of regulatory γ-subunit allosterically inhibits AMPK by destabilizing different regulatory regions of α-subunit: the autoinhibitory domain, the linker region, and the activation loop of the kinase core. The catalytically important residues experience a change in mechanical stress, and major rearrangements in community structure derived from residue-residue interaction energy-based network are observed in KD and α-linker region upon binding of different nucleotides. Our results also highlight the role of conserved charged residues forming an ionic network near the site 3 of γ-subunit in allosteric communications.

Dynamic Allostery Mediated by a Conserved Tryptophan in the Tec Family Kinases

Bruton’s tyrosine kinase (Btk) is a Tec family non-receptor tyrosine kinase that plays a critical role in immune signaling and is associated with the immunological disorder X-linked agammaglobulinemia (XLA). Our previous findings showed that the Tec kinases are allosterically activated by the adjacent N-terminal linker. A single tryptophan residue in the N-terminal 17-residue linker mediates allosteric activation, and its mutation to alanine leads to the complete loss of activity. Guided by hydrogen/deuterium exchange mass spectrometry results, we have employed Molecular Dynamics simulations, Principal Component Analysis, Community Analysis and measures of node centrality to understand the details of how a single tryptophan mediates allostery in Btk. A specific tryptophan side chain rotamer promotes the functional dynamic allostery by inducing coordinated motions that spread across the kinase domain. Either a shift in the rotamer population, or a loss of the tryptophan side chain by mutation, drastically changes the coordinated motions and dynamically isolates catalytically important regions of the kinase domain. This work also identifies a new set of residues in the Btk kinase domain with high node centrality values indicating their importance in transmission of dynamics essential for kinase activation. Structurally, these node residues appear in both lobes of the kinase domain. In the N-lobe, high centrality residues wrap around the ATP binding pocket connecting previously described Catalytic-spine residues. In the C-lobe, two high centrality node residues connect the base of the R- and C-spines on the αF-helix. We suggest that the bridging residues that connect the catalytic and regulatory architecture within the kinase domain may be a crucial element in transmitting information about regulatory spine assembly to the catalytic machinery of the catalytic spine and active site.

Bruton’s tyrosine kinase (Btk) belongs to the Tec family of protein tyrosine kinases, and plays a crucial role in the signaling pathway in B-cells. Alteration of Btk activity results in the serious immunological disorder, X-linked agammaglobulinemia. Btk is a multi-domain protein and the activity of the kinase domain is regulated by the adjacent non-catalytic domains, which mediate their effect by means of a conserved tryptophan residue. In this work, we have investigated the mechanism of regulation by this tryptophan residue, W395, in the linker preceding the Btk kinase domain. Using hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations we identify structural elements within the kinase domain that are required for function by transmitting the allosteric effects of W395. Molecular Dynamics simulations further guided us to delineate the kinase domain into dynamically correlated sets of residues using community analysis, thereby identifying the important communication nodes that connect the various elements of the kinase domain required for function. The analyses performed indicate clearly how the W395A mutant changes the communication pathway required for function.

Mechanistic Insights into R776H Mediated Activation of Epidermal Growth Factor Receptor Kinase

The epidermal growth factor receptor (EGFR) kinase is activated by a variety of mutations in human cancers. R776H is one such recurrent mutation (R752H in another numbering system) in the αC-β4 loop of the tyrosine kinase domain that activates EGFR in the absence of the activating EGF ligand. However, the mechanistic details of how R776H contributes to kinase activation are not well understood. Here using cell-based cotransfection assays, we show that the R776H mutation activates EGFR in a dimerization-dependent manner by preferentially adopting the acceptor position in the asymmetric dimer. The acceptor function, but not the donor function, is enhanced for the R776H mutant, supporting the "superacceptor" hypothesis proposed for oncogenic mutations in EGFR. We also find that phosphorylation of monomeric EGFR is increased by R776H mutation, providing insights into EGFR lateral phosphorylation and oligomerization. On the basis of molecular modeling and molecular dynamics simulation, we propose a model in which loss of key autoinhibitory αC-helix capping interaction and alteration of coconserved cis regulatory interactions between the kinase domain and the flanking regulatory segments contribute to mutational activation. Since the R776 equivalent position is mutated in ErbB2 and ErbB4, our studies have implications for understanding kinase mutational activation in other ErbB family members as well.

Effect of phosphorylation on EGFR dimer stability probed by single-molecule dynamics and FRET/FLIM

Deregulation of epidermal growth factor receptor (EGFR) signaling has been correlated with the development of a variety of human carcinomas. EGF-induced receptor dimerization and consequent trans- auto-phosphorylation are among the earliest events in signal transduction. Binding of EGF is thought to induce a conformational change that consequently unfolds an ectodomain loop required for dimerization indirectly. It may also induce important allosteric changes in the cytoplasmic domain. Despite extensive knowledge on the physiological activation of EGFR, the effect of targeted therapies on receptor conformation is not known and this particular aspect of receptor function, which can potentially be influenced by drug treatment, may in part explain the heterogeneous clinical response among cancer patients. Here, we used Förster resonance energy transfer/fluorescence lifetime imaging microscopy (FRET/FLIM) combined with two-color single-molecule tracking to study the effect of ATP-competitive small molecule tyrosine kinase inhibitors (TKIs) and phosphatase-based manipulation of EGFR phosphorylation on live cells. The distribution of dimer on-times was fitted to a monoexponential to extract dimer off-rates (koff). Our data show that pretreatment with gefitinib (active conformation binder) stabilizes the EGFR ligand-bound homodimer. Overexpression of EGFR-specific DEP-1 phosphatase was also found to have a stabilizing effect on the homodimer. No significant difference in the koff of the dimer could be detected when an anti-EGFR antibody (425 Snap single-chain variable fragment) that allows for dimerization of ligand-bound receptors, but not phosphorylation, was used. These results suggest that both the conformation of the extracellular domain and phosphorylation status of the receptor are involved in modulating the stability of the dimer. The relative fractions of these two EGFR subpopulations (interacting versus free) were obtained by a fractional-intensity analysis of ensemble FRET/FLIM images. Our combined imaging approach showed that both the fraction and affinity (surrogate of conformation at a single-molecule level) increased after gefitinib pretreatment or DEP-1 phosphatase overexpression. Using an EGFR mutation (I706Q, V948R) that perturbs the ability of EGFR to dimerize intracellularly, we showed that a modest drug-induced increase in the fraction/stability of the EGFR homodimer may have a significant biological impact on the tumor cell's proliferation potential.

Structure-based network analysis of activation mechanisms in the ErbB family of receptor tyrosine kinases: the regulatory spine residues are global mediators of structural stability and allosteric interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.

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.

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.

Ligand binding and dynamics of the monomeric epidermal growth factor receptor ectodomain

The ectodomain of the human epidermal growth factor receptor (hEGFR) controls input to several cell signalling networks via binding with extracellular growth factors. To gain insight into the dynamics and ligand binding of the ectodomain, the hEGFR monomer was subjected to molecular dynamics simulation. The monomer was found to be substantially more flexible than the ectodomain dimer studied previously. Simulations where the endogeneous ligand EGF binds to either Subdomain I or Subdomain III, or where hEGFR is unbound, show significant differences in dynamics. The molecular mechanics Poisson-Boltzmann surface area method has been used to derive relative free energies of ligand binding, and we find that the ligand is capable of binding either subdomain with a slight preference for III. Alanine-scanning calculations for the effect of selected ligand mutants on binding reproduce the trends of affinity measurements. Taken together, these results emphasize the possible role of the ectodomain monomer in the initial step of ligand binding, and add details to the static picture obtained from crystal structures.

Mechanism of displacement of a catalytically essential loop from the active site of mammalian fructose-1,6-bisphosphatase

AMP triggers a 15° subunit-pair rotation in fructose-1,6-bisphosphatase (FBPase) from its active R state to its inactive T state. During this transition, a catalytically essential loop (residues 50-72) leaves its active (engaged) conformation. Here, the structures of Ile(10) → Asp FBPase and molecular dynamic simulations reveal factors responsible for loop displacement. The AMP/Mg(2+) and AMP/Zn(2+) complexes of Asp(10) FBPase are in intermediate quaternary conformations (completing 12° of the subunit-pair rotation), but the complex with Zn(2+) provides the first instance of an engaged loop in a near-T quaternary state. The 12° subunit-pair rotation generates close contacts involving the hinges (residues 50-57) and hairpin turns (residues 58-72) of the engaged loops. Additional subunit-pair rotation toward the T state would make such contacts unfavorable, presumably causing displacement of the loop. Targeted molecular dynamics simulations reveal no steric barriers to subunit-pair rotations of up to 14° followed by the displacement of the loop from the active site. Principal component analysis reveals high-amplitude motions that exacerbate steric clashes of engaged loops in the near-T state. The results of the simulations and crystal structures are in agreement: subunit-pair rotations just short of the canonical T state coupled with high-amplitude modes sterically displace the dynamic loop from the active site.

EGFR: tale of the C-terminal tail

The carboxy terminal tail of epidermal growth factor receptor (EGFR) plays a critical role in the regulation of the enzyme activity of the kinase. There is a good structural model for the mechanism by which the C-terminal tail proximal to the kinase domain contributes to the negative regulation of the activity. Its conformation in the active state, conversely, has remained elusive due to its dynamic nature. A recently published structure of EGFR kinase domain shows the conformation of the proximal C-terminal tail in the active kinase. Analysis of this conformational state of the C-terminal tail is presented, and some of the mutagenesis data is revisited.

Structural modelling and dynamics of proteins for insights into drug interactions

Proteins are the workhorses of biomolecules and their function is affected by their structure and their structural rearrangements during ligand entry, ligand binding and protein-protein interactions. Hence, the knowledge of protein structure and, importantly, the dynamic behaviour of the structure are critical for understanding how the protein performs its function. The predictions of the structure and the dynamic behaviour can be performed by combinations of structure modelling and molecular dynamics simulations. The simulations also need to be sensitive to the constraints of the environment in which the protein resides. Standard computational methods now exist in this field to support the experimental effort of solving protein structures. This review presents a comprehensive overview of the basis of the calculations and the well-established computational methods used to generate and understand protein structure and function and the study of their dynamic behaviour with the reference to lung-related targets.

Identification of potent EGFR inhibitors from TCM Database@Taiwan

Overexpression of epidermal growth factor receptor (EGFR) has been associated with cancer. Targeted inhibition of the EGFR pathway has been shown to limit proliferation of cancerous cells. Hence, we employed Traditional Chinese Medicine Database (TCM Database@Taiwan) (http://tcm.cmu.edu.tw) to identify potential EGFR inhibitor. Multiple Linear Regression (MLR), Support Vector Machine (SVM), Comparative Molecular Field Analysis (CoMFA), and Comparative Molecular Similarities Indices Analysis (CoMSIA) models were generated using a training set of EGFR ligands of known inhibitory activities. The top four TCM candidates based on DockScore were 2-O-caffeoyl tartaric acid, Emitine, Rosmaricine, and 2-O-feruloyl tartaric acid, and all had higher binding affinities than the control Iressa®. The TCM candidates had interactions with Asp855, Lys716, and Lys728, all which are residues of the protein kinase binding site. Validated MLR (r² = 0.7858) and SVM (r² = 0.8754) models predicted good bioactivity for the TCM candidates. In addition, the TCM candidates contoured well to the 3D-Quantitative Structure-Activity Relationship (3D-QSAR) map derived from the CoMFA (q² = 0.721, r² = 0.986) and CoMSIA (q² = 0.662, r² = 0.988) models. The steric field, hydrophobic field, and H-bond of the 3D-QSAR map were well matched by each TCM candidate. Molecular docking indicated that all TCM candidates formed H-bonds within the EGFR protein kinase domain. Based on the different structures, H-bonds were formed at either Asp855 or Lys716/Lys728. The compounds remained stable throughout molecular dynamics (MD) simulation. Based on the results of this study, 2-O-caffeoyl tartaric acid, Emitine, Rosmaricine, and 2-O-feruloyl tartaric acid are suggested to be potential EGFR inhibitors.

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.

The energy landscape analysis of cancer mutations in protein kinases

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