Molecular dynamics simulations reveal that Tyr-317 phosphorylation reduces Shc binding affinity for phosphotyrosyl residues of epidermal growth factor receptor

Evolutionary Couplings and Molecular Dynamic Simulations Highlight Details of GPCRs Heterodimers' Interfaces

A growing body of evidence suggests that only a few amino acids ("hot-spots") at the interface contribute most of the binding energy in transient protein-protein interactions. However, experimental protocols to identify these hot-spots are highly labor-intensive and expensive. Computational methods, including evolutionary couplings, have been proposed to predict the hot-spots, but they generally fail to provide details of the interacting amino acids. Here we showed that unbiased evolutionary methods followed by biased molecular dynamic simulations could achieve this goal and reveal critical elements of protein complexes. We applied the methodology to selected G-protein coupled receptors (GPCRs), known for their therapeutic properties. We used the structure-prior-assisted direct coupling analysis (SP-DCA) to predict the binding interfaces of A2aR/D2R, CB1R/D2R, A2aR/CB1R, 5HT2AR/D2R, and 5-HT2AR/mGluR2 receptor heterodimers, which all agreed with published data. In order to highlight details of the interactions, we performed molecular dynamic (MD) simulations using the newly developed AWSEM energy model. We found that these receptors interact primarily through critical residues at the C and N terminal domains and the third intracellular loop (ICL3). The MD simulations showed that these residues are energetically necessary for dimerization and revealed their native conformational state. We subsequently applied the methodology to the 5-HT2AR/5-HTR4R heterodimer, given its implication in drug addiction and neurodegenerative pathologies such as Alzheimer's disease (AD). Further, the SP-DCA analysis showed that 5-HT2AR and 5-HTR4R heterodimerize through the C-terminal domain of 5-HT2AR and ICL3 of 5-HT4R. However, elucidating the details of GPCR interactions would accelerate the discovery of druggable sites and improve our knowledge of the etiology of common diseases, including AD.

A survey of multiscale modeling: Foundations, historical milestones, current status, and future prospects

Research problems in the domains of physical, engineering, biological sciences often span multiple time and length scales, owing to the complexity of information transfer underlying mechanisms. Multiscale modeling (MSM) and high-performance computing (HPC) have emerged as indispensable tools for tackling such complex problems. We review the foundations, historical developments, and current paradigms in MSM. A paradigm shift in MSM implementations is being fueled by the rapid advances and emerging paradigms in HPC at the dawn of exascale computing. Moreover, amidst the explosion of data science, engineering, and medicine, machine learning (ML) integrated with MSM is poised to enhance the capabilities of standard MSM approaches significantly, particularly in the face of increasing problem complexity. The potential to blend MSM, HPC, and ML presents opportunities for unbound innovation and promises to represent the future of MSM and explainable ML that will likely define the fields in the 21st century.

Phosphorylation of threonine residues on Shc promotes ligand binding and mediates crosstalk between MAPK and Akt pathways in breast cancer cells

Scaffold proteins play important roles in regulating signalling network fidelity, the absence of which is often the basis for diseases such as cancer. In the present work, we show that the prototypical scaffold protein Shc is phosphorylated by the extracellular signal-regulated kinase, Erk. In addition, Shc threonine phosphorylation is specifically up-regulated in two selected triple-negative breast cancer (TNBC) cell lines. To explore how Erk-mediated threonine phosphorylation on Shc might play a role in the dysregulation of signalling events, we investigated how Shc affects pathways downstream of EGF receptor. Using an in vitro model and biophysical analysis, we show that Shc threonine phosphorylation is responsible for elevated Akt and Erk signalling, potentially through the recruitment of the 14-3-3 ζ and Pin-1 proteins.

Fortilin binds IRE1α and prevents ER stress from signaling apoptotic cell death

The endoplasmic reticulum, the cytoplasmic organelle that matures a massive amount of nascent secretory polypeptides, is particularly sensitive to stress. Endoplasmic reticulum stress causes unfolded proteins to populate the organelle, eliciting the unfolded protein response. During the unfolded protein response, GRP78—an endoplasmic reticulum master stress regulator—detaches from three endoplasmic reticulum stress sensors (IRE1α, PERK, and ATF6) and allows them to activate the apoptotic signaling pathway. Fortilin, a pro-survival molecule, is known to inhibit apoptosis by binding and inhibiting p53, but its role in endoplasmic reticulum stress-induced apoptosis remains unknown. Here, we report that fortilin directly interacts with the cytoplasmic domain of IRE1α, inhibits both kinase and endoribonuclease (RNase) activities of the stress sensor, and protects cells against apoptotic cell death at both cellular and whole animal levels. Our data support a role of fortilin in the unfolded protein response and its potential participation in human diseases caused by unfolded protein response.

IRE1α is an ER stress sensor, whose activity induces apoptosis. Here, the authors report that fortilin, a pro-survival factor, with yet unknown roles in ER stress, interacts with active IRE1α, inhibits both its kinase end RNase activities, and protects cells from apoptosis both in vitro and in vivo.

Understanding the molecular basis of substrate binding specificity of PTB domains

Protein-protein interactions mediated by phosphotyrosine binding (PTB) domains play a crucial role in various cellular processes. In order to understand the structural basis of substrate recognition by PTB domains, multiple explicit solvent atomistic simulations of 100ns duration have been carried out on 6 PTB-peptide complexes with known binding affinities. MM/PBSA binding energy values calculated from these MD trajectories and residue based statistical pair potential score show good correlation with the experimental dissociation constants. Our analysis also shows that the modeled structures of PTB domains can be used to develop less compute intensive residue level statistical pair potential based approaches for predicting interaction partners of PTB domains.

All-atom molecular dynamics analysis of multi-peptide systems reproduces peptide solubility in line with experimental observations

In order to investigate the contribution of individual amino acids to protein and peptide solubility, we carried out 100 ns molecular dynamics (MD) simulations of 10⁶ ų cubic boxes containing ~3 × 10⁴ water molecules and 27 tetra-peptides regularly positioned at 23 Å from each other and composed of a single amino acid type for all natural amino acids but cysteine and glycine. The calculations were performed using Amber with a standard force field on a special purpose MDGRAPE-3 computer, without introducing any “artificial” hydrophobic interactions. Tetra-peptides composed of I, V, L, M, N, Q, F, W, Y, and H formed large amorphous clusters, and those containing A, P, S, and T formed smaller ones. Tetra-peptides made of D, E, K, and R did not cluster at all. These observations correlated well with experimental solubility tendencies as well as hydrophobicity scales with correlation coefficients of 0.5 to > 0.9. Repulsive Coulomb interactions were dominant in ensuring high solubility, whereas both Coulomb and van der Waals (vdW) energies contributed to the aggregations of low solubility amino acids. Overall, this very first all-atom molecular dynamics simulation of a multi-peptide system appears to reproduce the basic properties of peptide solubility, essentially in line with experimental observations.

In silico prediction and characterization of protein post-translational modifications

This review outlines the computational approaches and procedures for predicting post translational modification (PTM)-induced changes in protein conformation and their influence on protein function(s), the latter being assessed as differential affinity in interaction with either low (ligands for receptors or transporters, substrates for enzymes) or high molecular mass molecules (proteins or nucleic acids in supramolecular assemblies). The scope for an in silico approach is discussed against a summary of the in vitro evidence on the structural and functional outcome of protein PTM.

Physicochemical mechanisms of protein regulation by phosphorylation

Phosphorylation offers a dynamic way to regulate protein activity and subcellular localization, which is achieved through its reversibility and fast kinetics. Adding or removing a dianionic phosphate group somewhere on a protein often changes the protein’s structural properties, its stability and dynamics. Moreover, the majority of signaling pathways involve an extensive set of protein–protein interactions, and phosphorylation can be used to regulate and modulate protein–protein binding. Losses of phosphorylation sites, as a result of disease mutations, might disrupt protein binding and deregulate signal transduction. In this paper we focus on the effects of phosphorylation on protein stability, dynamics, and binding. We describe several physico-chemical mechanisms of protein regulation through phosphorylation and pay particular attention to phosphorylation in protein complexes and phosphorylation in the context of disorder–order and order–disorder transitions. Finally we assess the role of multiple phosphorylation sites in a protein molecule, their possible cooperativity and function.

The ShcD signaling adaptor facilitates ligand-independent phosphorylation of the EGF receptor

Proto-oncogenic Src homology and collagen (Shc) proteins have been considered archetypal adaptors of epidermal growth factor receptor (EGFR)-mediated signaling. We report that in addition to its role as an EGFR-binding partner and Grb2 platform, ShcD acts noncanonically to promote phosphorylation of select EGFR residues. Unexpectedly, Y1068, Y1148, and Y1173 are subject to ShcD-induced, cell-autonomous hyperphosphorylation in the absence of external stimuli. This response is not elicited by other Shc proteins and requires the intrinsic EGFR kinase, as well as the ShcD phosphotyrosine-binding (PTB) domain. Assessments of Erk, Akt, phospholipase C 1γ, and FAK pathways reveal no apparent distal signaling targets of ShcD. Nevertheless, the capacity of cultured cells to repopulate a wounded monolayer is markedly accelerated by ShcD in an EGFR kinase-dependent manner. Furthermore, detection of overexpressed ShcD coincident with EGFR phosphorylation in human gliomas suggests a clinical application for these findings. We thus demonstrate unique and relevant synergy between ShcD and EGFR that is unprecedented among signaling adaptors.

Engineering a single ubiquitin ligase for the selective degradation of all activated ErbB receptor tyrosine kinases

Interrogating specific cellular activities often entails the dissection of posttranslational modifications or functional redundancy conferred by protein families, which demands more sophisticated research tools than simply eliminating a specific gene product by gene targeting or RNA interference. We have developed a novel methodology that involves engineering a single SCF(βTrCP)-based ubiquitin ligase that is capable of not only simultaneously targeting the entire family of ErbB receptor tyrosine kinases for ubiquitination and degradation, but also selectively recruiting only activated ErbBs. The engineered SCF(βTrCP) ubiquitin ligase effectively blocked ErbB signaling and attenuated oncogenicity in breast cancer cells, yet had little effect on the survival and growth of non-cancerous breast epithelial cells. Therefore, engineering ubiquitin ligases offers a simple research tool to dissect the specific traits of tumorigenic protein families, and provides a rapid and feasible means to expand the dimensionality of drug discovery by assessing protein families or posttranslational modifications as potential drug targets.

Teaching an old dogma new tricks: twenty years of Shc adaptor signalling

Shc (Src homology and collagen homology) proteins are considered prototypical signalling adaptors in mammalian cells. Consisting of four unique members, ShcA, B, C and D, and multiple splice isoforms, the family is represented in nearly every cell type in the body, where it engages in an array of fundamental processes to transduce environmental stimuli. Two decades of investigation have begun to illuminate the mechanisms of the flagship ShcA protein, whereas much remains to be learned about the newest discovery, ShcD. It is clear, however, that the distinctive modular architecture of Shc proteins, their promiscuous phosphotyrosine-based interactions with a multitude of membrane receptors, involvement in central cascades including MAPK (mitogen-activated protein kinase) and Akt, and unconventional contributions to oxidative stress and apoptosis all require intricate regulation, and underlie diverse physiological function. From early cardiovascular development and neuronal differentiation to lifespan determination and tumorigenesis, Shc adaptors have proven to be more ubiquitous, versatile and dynamic than their structures alone suggest.

Influence of 63Ser phosphorylation and dephosphorylation on the structure of the stathmin helical nucleation sequence: a molecular dynamics study

Phosphorylation is an important mechanism regulating protein-protein interactions involving intrinsically disordered protein regions. Stathmin, an archetypical example of an intrinsically disordered protein, is a key regulator of microtubule dynamics in which phosphorylation of 63Ser within the helical nucleation sequence strongly down-regulates the tubulin binding and microtubule destabilizing activities of the protein. Experimental studies on a peptide encompassing the 19-residue helical nucleation sequence of stathmin (residues 55-73) indicate that phosphorylation of 63Ser destabilizes the peptide's secondary structure by disrupting the salt bridges supporting its helical conformation. In order to investigate this hypothesis at atomic resolution, we performed molecular dynamics simulations of nonphosphorylated and phosphorylated stathmin-[55-73] at room temperature and pressure, neutral pH, and explicit solvation using the recently released GROMOS force field 54A7. In the simulations of nonphosphorylated stathmin-[55-73] emerged salt bridges associated with helical configurations. In the simulations of 63Ser phosphorylated stathmin-[55-73] these configurations dispersed and were replaced by a proliferation of salt bridges yielding disordered configurations. The transformation of the salt bridges was accompanied by emergence of numerous interactions between main and side chains, involving notably the oxygen atoms of the phosphorylated 63Ser. The loss of helical structure induced by phosphorylation is reversible, however, as a final simulation showed. The results extend the hypothesis of salt bridge derangement suggested by experimental observations of the stathmin nucleation sequence, providing new insights into regulation of intrinsically disordered protein systems mediated by phosphorylation.

Experimental identification and theoretical analysis of a thermally stabilized green fluorescent protein variant

In this study, we aim to relate experimentally measured macroscopic properties to dynamic and structural changes as calculated by molecular dynamics (MD) simulations. We performed the analysis on four GFP (green fluorescent protein) variants, which have amino acid replacements or insertion in a flexible region on the protein surface and which resulted from a previous protein splicing reaction optimization experiment. The variants are a reference GFP (CEGFP), GFP-N144C, GFP-N144C/Y145F, and a GFP with five residues inserted between Y145 and N146 (GFP-5ins). As a result, we identified a single Y145F mutation that increased the thermal stability of GFP-N144C/Y145F by 3-4 °C. Because circular dichroism measurements indicated that the overall GFP β-barrel fold was maintained in all variants, we presumed that the fluorescence activity and thermal stability related to local changes that could be detected by standard MD simulations. The 60 ns MD simulations indicated that the Y145's hydroxyl group, which is straight and buried in the crystal structure, was bent avoiding the hydrophobic core during the simulation in both CEGFP and GFP-N144C. This local strain was relieved in GFP-N144C/Y145F, where the tyrosine's hydroxyl group was replaced with the F145 hydrophobic aliphatic carbon. F145 remained indeed buried during the simulation maintaining local compactness, which presumably reflected the improved thermal stability of GFP-N144C/Y145F. Furthermore, the analysis of internal water molecules localized within the GFP's β-barrel suggested that a change in the local hydrogen bonding pattern around the chromophore correlated with a strong fluorescence activity decrease in GFP-5ins. Although relating experimental observation with calculated molecular features proved to be delicate, this study suggested that some microscopic features could be useful reporters for redesigning GFPs and other proteins. The newly identified GFP-N144C/Y145F was among the most stable GFP variant and demonstrates the potential of such computer-aided design.

Structural systems biology and multiscale signaling models

We review current advances in experimental as well as computational modeling and simulation approaches to structural systems biology, whose overall aim is to build quantitative models of signaling networks while retaining the crucial elements of molecular specificity. We briefly discuss the current and emerging experimental and computational methods, particularly focusing on hybrid and multiscale methods, and highlight several applications in cell signaling with quantitative and predictive capabilities. The scope of such models range from delineating protein-protein interactions to describing clinical implications.

An adaptive coarse graining method for signal transduction in three dimensions

The spatio-temporal landscape of the plasma membrane regulates activation and signal transduction of membrane bound receptors by restricting their two-dimensional mobility and by inducing receptor clustering. This regulation also extends to complex formation between receptors and adaptor proteins, which are the intermediate signaling molecules involved in cellular signaling that relay the received cues from cell surface to cytoplasm and eventually to the nucleus. Although their investigation poses challenging technical difficulties, there is a crucial need to understand the impact of the receptor diffusivity, clustering, and spatial heterogeneity, and of receptor-adaptor protein complex formation on the cellular signal transduction patterns. Building upon our earlier studies, we have developed an adaptive coarse-grained Monte Carlo method that can be used to investigate the role of diffusion, clustering and membrane corralling on receptor association and receptor-adaptor protein complex formation dynamics in three dimensions. The new Monte Carlo lattice based approach allowed us to introduce spatial resolution on the 2-D plasma membrane and to model the cytoplasm in three-dimensions. Being a multi-resolution approach, our new method makes it possible to represent various parts of the cellular system at different levels of detail and enabled us to utilize the locally homogeneous assumption when justified (e.g., cytoplasmic region away from the cell membrane) and avoid its use when high spatial resolution is needed (e.g., cell membrane and cytoplasmic region near the membrane) while keeping the required computational complexity manageable. Our results have shown that diffusion has a significant impact on receptor-receptor dimerization and receptor-adaptor protein complex formation kinetics. We have observed an "adaptor protein hopping" mechanism where the receptor binding proteins may hop between receptors to form short-lived transient complexes. This increased residence time of the adaptor proteins near cell membrane and their ability to frequently change signaling partners may explain the increase in signaling efficiency when receptors are clustered. We also hypothesize that the adaptor protein hopping mechanism can cause concurrent or sequential activation of multiple signaling pathways, thus leading to crosstalk between diverse biological functions.

Dynamic conformational changes in munc18 prevent syntaxin binding

The Sec1/munc18 protein family is essential for vesicle fusion in eukaryotic cells via binding to SNARE proteins. Protein kinase C modulates these interactions by phosphorylating munc18a thereby reducing its affinity to one of the central SNARE members, syntaxin-1a. The established hypothesis is that the reduced affinity of the phosphorylated munc18a to syntaxin-1a is a result of local electrostatic repulsion between the two proteins, which interferes with their compatibility. The current study challenges this paradigm and offers a novel mechanistic explanation by revealing a syntaxin-non-binding conformation of munc18a that is induced by the phosphomimetic mutations. In the present study, using molecular dynamics simulations, we explored the dynamics of the wild-type munc18a versus phosphomimetic mutant munc18a. We focused on the structural changes that occur in the cavity between domains 3a and 1, which serves as the main syntaxin-binding site. The results of the simulations suggest that the free wild-type munc18a exhibits a dynamic equilibrium between several conformations differing in the size of its cavity (the main syntaxin-binding site). The flexibility of the cavity's size might facilitate the binding or unbinding of syntaxin. In silico insertion of phosphomimetic mutations into the munc18a structure induces the formation of a conformation where the syntaxin-binding area is rigid and blocked as a result of interactions between residues located on both sides of the cavity. Therefore, we suggest that the reduced affinity of the phosphomimetic mutant/phosphorylated munc18a is a result of the closed-cavity conformation, which makes syntaxin binding energetically and sterically unfavorable. The current study demonstrates the potential of phosphoryalation, an essential biological process, to serve as a driving force for dramatic conformational changes of proteins modulating their affinity to target proteins.

Protein phosphorylation plays a significant regulatory role in multi-component systems engaged in signal transduction or coordination of cellular processes, by activating or deactivating proteins. The potential of phosphorylation to induce substantial conformational changes in proteins, thereby changing their affinity to target proteins, has already been shown but the dynamics of the process is not fully elucidated. In the present study, we investigated, by molecular dynamics simulations, the dynamic conformational changes in munc18a, a protein that is crucial for neurotransmitter release and interacts tightly with the SNARE syntaxin-1. We further investigated the conformational changes that occur in munc18a when it is phosphorylated, reducing its affinity to syntaxin-1a. The results of the simulations suggest that there is a conformational flexibility of the syntaxin-unbounded munc18a that allows changes in the shape of the syntaxin-1a binding cavity. In silico insertion of phosphomimetic mutations into munc18a led to a reduction in the flexibility and closure of the syntaxin-binding site. We suggest that the reduced affinity of phosphorylated munc18a to syntaxin-1a stems from the difficulty of syntaxin-1a to bind to the munc18a closed-cavity conformation, induced by the PKC phosphorylation of munc18a.

Investigating Molecular Mechanisms of Activation and Mutation of the HER2 Receptor Tyrosine Kinase through Computational Modeling and Simulation

Human epidermal growth factor receptor 2 (HER2)/ErbB2 is a receptor tyrosine kinase belonging to the EGFR/ErbB family and is overexpressed in 20-30% of human breast cancers. Since there is a growing effort to develop pharmacological inhibitors of the HER2 kinase for the treatment of breast cancer, it is clinically valuable to rationalize how specific mutations impact the molecular mechanism of receptor activation. Although several crystal structures of the ErbB kinases have been solved, the precise mechanism of HER2 activation remains unknown, and it has been suggested that HER2 is unique in its requirement for phosphorylation of Y877, a key tyrosine residue located in the activation loop (A-loop). In our studies, discussed here, we have investigated the mechanisms that are important in HER2 kinase domain regulation and compared them with the other ErbB family members, namely EGFR and ErbB4, to determine the molecular basis for HER2's unique mode of activation. We apply computational simulation techniques at the atomic level and at the electronic structure (quantum mechanical) level to elucidate details of the mechanisms governing the kinase domains of these ErbB members. Through analysis of our simulation results, we have discovered potential regulatory mechanisms common to EGFR, HER2, and ErbB4, including a tight coupling between the A-loop and catalytic loop that may contribute to alignment of residues required for catalysis in the active kinase. We further postulate an autoinhibitory mechanism whereby the inactive kinase is stabilized through sequestration of catalytic residues. In HER2, we also predict a role for phosphorylated Y877 in bridging a network of hydrogen bonds that fasten the A-loop in its active conformation, suggesting that HER2 may be unique among the ErbB members in requiring A-loop tyrosine phosphorylation for functionality. In EGFR, HER2, and ErbB4, we discuss the possible effects of activating mutations. Delineation of the activation mechanism of HER2 in the context of the other ErbB members is crucial for understanding how the activated kinase might interact with downstream molecules and couple to signaling cascades that promote cancer. Our comparative analysis furthers insight into the mechanics of activation of the HER2 kinase and enables us to predict the effect of an identified insertion mutation on HER2 activation. Further understanding of the mechanism of HER2 kinase activation at the atomic scale and how it couples to downstream signaling at the cellular scale will elucidate predictive molecular phenotypes that may indicate likelihood of response to specific therapies for HER2-mediated cancers.

Molecular systems biology of ErbB1 signaling: bridging the gap through multiscale modeling and high-performance computing

The complexity in intracellular signaling mechanisms relevant for the conquest of many diseases resides at different levels of organization with scales ranging from the subatomic realm relevant to catalytic functions of enzymes to the mesoscopic realm relevant to the cooperative association of molecular assemblies and membrane processes. Consequently, the challenge of representing and quantifying functional or dysfunctional modules within the networks remains due to the current limitations in our understanding of mesoscopic biology, i.e., how the components assemble into functional molecular ensembles. A multiscale approach is necessary to treat a hierarchy of interactions ranging from molecular (nm, ns) to signaling (microm, ms) length and time scales, which necessitates the development and application of specialized modeling tools. Complementary to multiscale experimentation (encompassing structural biology, mechanistic enzymology, cell biology, and single molecule studies) multiscale modeling offers a powerful and quantitative alternative for the study of functional intracellular signaling modules. Here, we describe the application of a multiscale approach to signaling mediated by the ErbB1 receptor which constitutes a network hub for the cell's proliferative, migratory, and survival programs. Through our multiscale model, we mechanistically describe how point-mutations in the ErbB1 receptor can profoundly alter signaling characteristics leading to the onset of oncogenic transformations. Specifically, we describe how the point mutations induce cascading fragility mechanisms at the molecular scale as well as at the scale of the signaling network to preferentially activate the survival factor Akt. We provide a quantitative explanation for how the hallmark of preferential Akt activation in cell-lines harboring the constitutively active mutant ErbB1 receptors causes these cell-lines to be addicted to ErbB1-mediated generation of survival signals. Consequently, inhibition of ErbB1 activity leads to a remarkable therapeutic response in the addicted cell lines.