Conformational changes in protein loops and helices induced by post-translational phosphorylation

A Tutorial Review of Labeling Methods in Mass Spectrometry-Based Quantitative Proteomics

Recent advancements in mass spectrometry (MS) have revolutionized quantitative proteomics, with multiplex isotope labeling emerging as a key strategy for enhancing accuracy, precision, and throughput. This tutorial review offers a comprehensive overview of multiplex isotope labeling techniques, including precursor-based, mass defect-based, reporter ion-based, and hybrid labeling methods. It details their fundamental principles, advantages, and inherent limitations along with strategies to mitigate the limitation of ratio-distortion. This review will also cover the applications and latest progress in these labeling techniques across various domains, including cancer biomarker discovery, neuroproteomics, post-translational modification analysis, cross-linking MS, and single-cell proteomics. This Review aims to provide guidance for researchers on selecting appropriate methods for their specific goals while also highlighting the potential future directions in this rapidly evolving field.

phosaa14SB and phosaa19SB: Updated Amber Force Field Parameters for Phosphorylated Amino Acids

Phosphorylated amino acids are involved in many cell regulatory networks; proteins containing these post-translational modifications are widely studied both experimentally and computationally. Simulations are used to investigate a wide range of structural and dynamic properties of biomolecules, such as ligand binding, enzyme-reaction mechanisms, and protein folding. However, the development of force field parameters for the simulation of proteins containing phosphorylated amino acids using the Amber program has not kept pace with the development of parameters for standard amino acids, and it is challenging to model these modified amino acids with accuracy comparable to proteins containing only standard amino acids. In particular, the popular ff14SB and ff19SB models do not contain parameters for phosphorylated amino acids. Here, the dihedral parameters for the side chains of the most common phosphorylated amino acids are trained against reference data from QM calculations adopting the ff14SB approach, followed by validation against experimental data. Library files and corresponding parameter files are provided, with versions that are compatible with both ff14SB and ff19SB.

Interaction of Soybean (Glycine max (L.) Merr.) Class II ACBPs with MPK2 and SAPK2 Kinases: New Insights into the Regulatory Mechanisms of Plant ACBPs

Plant acyl-CoA-binding proteins (ACBPs) function in plant development and stress responses, with some ACBPs interacting with protein partners. This study tested the interaction between two Class II GmACBPs (Glycine max ACBPs) and seven kinases, using yeast two-hybrid (Y2H) assays and bimolecular fluorescence complementation (BiFC). The results revealed that both GmACBP3.1 and GmACBP4.1 interact with two soybean kinases, a mitogen-activated protein kinase MPK2, and a serine/threonine-protein kinase SAPK2, highlighting the significance of the ankyrin-repeat (ANK) domain in facilitating protein-protein interactions. Moreover, an in vitro kinase assay and subsequent Phos-tag SDS-PAGE determined that GmMPK2 and GmSAPK2 possess the ability to phosphorylate Class II GmACBPs. Additionally, the kinase-specific phosphosites for Class II GmACBPs were predicted using databases. The HDOCK server was also utilized to predict the binding models of Class II GmACBPs with these two kinases, and the results indicated that the affected residues were located in the ANK region of Class II GmACBPs in both docking models, aligning with the findings of the Y2H and BiFC experiments. This is the first report describing the interaction between Class II GmACBPs and kinases, suggesting that Class II GmACBPs have potential as phospho-proteins that impact signaling pathways.

Design patterns of biological cells

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

Attenphos: General Phosphorylation Site Prediction Model Based on Attention Mechanism

Phosphorylation site prediction has important application value in the field of bioinformatics. It can act as an important reference and help with protein function research, protein structure research, and drug discovery. So, it is of great significance to propose scientific and effective calculation methods to accurately predict phosphorylation sites. In this study, we propose a new method, Attenphos, based on the self-attention mechanism for predicting general phosphorylation sites in proteins. The method not only captures the long-range dependence information of proteins but also better represents the correlation between amino acids through feature vector encoding transformation. Attenphos takes advantage of the one-dimensional convolutional layer to reduce the number of model parameters, improve model efficiency and prediction accuracy, and enhance model generalization. Comparisons between our method and existing state-of-the-art prediction tools were made using balanced datasets from human proteins and unbalanced datasets from mouse proteins. We performed prediction comparisons using independent test sets. The results showed that Attenphos demonstrated the best overall performance in the prediction of Serine (S), Threonine (T), and Tyrosine (Y) sites on both balanced and unbalanced datasets. Compared to current state-of-the-art methods, Attenphos has significantly higher prediction accuracy. This proves the potential of Attenphos in accelerating the identification and functional analysis of protein phosphorylation sites and provides new tools and ideas for biological research and drug discovery.

Computational approaches to investigate fluoride binding, selectivity and transport across the membrane

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

Analysis of Structural Changes in the Protein near the Phosphorylation Site

Modification of the protein after synthesis (PTM) often affects protein function as supported by numerous studies. However, there is no consensus about the degree of structural protein changes after modification. For phosphorylation of serine, threonine, and tyrosine, which is a common PTM in the biology of living organisms, we consider topical issues related to changes in the geometric parameters of a protein (Rg, RMSD, Cα displacement, SASA). The effect of phosphorylation on protein geometry was studied both for the whole protein and at the local level (i.e., in different neighborhoods of the modification site). Heterogeneity in the degree of protein structural changes after phosphorylation was revealed, which allowed for us to isolate a group of proteins having pronounced local structural changes in the neighborhoods of up to 15 amino acid residues from the modification site. This is a comparative study of protein structural changes in neighborhoods of 3-15 amino acid residues from the modified site. Amino acid phosphorylation in proteins with pronounced local changes caused switching from the inactive functional state to the active one.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

The Conformation of the Intrinsically Disordered N-Terminal Region of Barrier-to-Autointegration Factor (BAF) is Regulated by pH and Phosphorylation

Barrier-to-Autointegration Factor (BAF) is a highly conserved DNA binding protein important for genome integrity. Its localization and function are regulated through phosphorylation. Previously reported structures of BAF suggested that it is fully ordered, but our recent NMR analysis revealed that its N-terminal region is flexible in solution and that S4/T3 di-phosphorylation by VRK1 reduces this flexibility. Here, molecular dynamics (MD) simulation was used to unveil the conformational ensembles accessible to the N-terminal region of BAF either unphosphorylated, mono-phosphorylated on S4 or di-phosphorylated on S4/T3 (pBAF) and to reveal the interactions that contribute to define these ensembles. We show that the intrinsic flexibility observed in the N-terminal region of BAF is reduced by S4 phosphorylation and to a larger extent by S4/T3 di-phosphorylation. Thanks to the atomic description offered by MD supported by the NMR study of several BAF mutants, we identified the dynamic network of salt bridge interactions responsible for the conformational restriction involving pS4 and pT3 with residues located in helix α1 and α6. Using MD, we showed that the flexibility in the N-terminal region of BAF depends on the ionic strength and on the pH. We show that the presence of two negative charges of the phosphoryl groups is required for a substantial decrease in flexibility in pBAF. Using MD supported by NMR, we also showed that H7 deprotonation reduces the flexibility in the N-terminal region of BAF. Thus, the conformation of the intrinsically disordered N-terminal region of BAF is highly tunable, likely related to its diverse functions.

Remodeling of algal photosystem I through phosphorylation

Photosystem I (PSI) with its associated light-harvesting system is the most important generator of reducing power in photosynthesis. The PSI core complex is highly conserved, whereas peripheral subunits as well as light-harvesting proteins (LHCI) reveal a dynamic plasticity. Moreover, in green alga, PSI-LHCI complexes are found as monomers, dimers, and state transition complexes, where two LHCII trimers are associated. Herein, we show light-dependent phosphorylation of PSI subunits PsaG and PsaH as well as Lhca6. Potential consequences of the dynamic phosphorylation of PsaG and PsaH are structurally analyzed and discussed in regard to the formation of the monomeric, dimeric, and LHCII-associated PSI-LHCI complexes.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Cysteine Oxidation in Proteins: Structure, Biophysics, and Simulation

Cysteine side chains can exist in distinct oxidation states depending on the pH and redox potential of the environment, and cysteine oxidation plays important yet complex regulatory roles. Compared with the effects of post-translational modifications such as phosphorylation, the effects of oxidation of cysteine to sulfenic, sulfinic, and sulfonic acid on protein structure and function remain relatively poorly characterized. We present an analysis of the role of cysteine reactivity as a regulatory factor in proteins, emphasizing the interplay between electrostatics and redox potential as key determinants of the resulting oxidation state. A review of current computational approaches suggests underdeveloped areas of research for studying cysteine reactivity through molecular simulations.

Structural analysis of mammalian protein phosphorylation at a proteome level

Phosphorylation is an essential post-translational modification for almost all cellular processes. Several global phosphoproteomics analyses have revealed phosphorylation profiles under different conditions. Beyond identification of phospho-sites, protein structures add another layer of information about their functionality. In this study, we systematically characterize phospho-sites based on their 3D locations in the protein and establish a location map for phospho-sites. More than 250,000 phospho-sites have been analyzed, of which 8,686 sites match at least one structure and are stratified based on their respective 3D positions. Core phospho-sites possess two distinct groups based on their dynamicity. Dynamic core phosphorylations are significantly more functional compared with static ones. The dynamic core and the interface phospho-sites are the most functional among all 3D phosphorylation groups. Our analysis provides global characterization and stratification of phospho-sites from a structural perspective that can be utilized for predicting functional relevance and filtering out false positives in phosphoproteomic studies.

Optimized structure of monoubiquitinated FANCD2 (human) at Lys 561: a theoretical approach

Fanconi anaemia pathway repairs inter-strand cross linking damage (ICL) of the DNA. Monoubiquitination of FANCD2 and FANCI is very crucial for ICL repairing. In this work we have tried to understand the monoubiquitinated FANCD2 structure, which facilitates the FANCD2 for binding the damage part of the chromatin. Crystal structure of the monoubiquitinated FANCD2 alone is not available, therefore we have modelled the optimized structure of the human monoubiquitinated (Lys 561) FANCD2. As there is no suitable software or web server we have developed a method for building up monoubiquitinated product and validated on simplest monoubiquitinated protein, diubiquitin. We have predicted the structure of human monoubiquitinated FANCD2 by using our method and studied the interaction with DNA by docking studies. Molecular Dynamics (MD) simulation has been used to understand the stability of the structure. Large structural differences have been observed between FANCD2 and monoubiquitinated FANCD2. Docking studies with DNA suggest that the binding site varies for the FANCD2 and monoubiquitinated FANCD2.Communicated by Ramaswamy H. Sarma.

A novel cell-penetrating peptide targeting calpain-cleavage of PSD-95 induced by excitotoxicity improves neurological outcome after stroke

Postsynaptic density protein-95 (PSD-95) is a multidomain protein critical to the assembly of signaling complexes at excitatory synapses, required for neuronal survival and function. However, calpain-processing challenges PSD-95 function after overactivation of excitatory glutamate receptors (excitotoxicity) in stroke, a leading cause of death, disability and dementia in need of efficient pharmacological treatments. A promising strategy is neuroprotection of the infarct penumbra, a potentially recoverable area, by promotion of survival signaling. Interference of PSD-95 processing induced by excitotoxicity might thus be a therapeutic target for stroke and other excitotoxicity-associated pathologies.

Methods: The nature and stability of PSD-95 calpain-fragments was analyzed using in vitro assays or excitotoxic conditions induced in rat primary neuronal cultures or a mouse model of stroke. We then sequenced PSD-95 cleavage-sites and rationally designed three cell-penetrating peptides (CPPs) containing these sequences. The peptides effects on PSD-95 stability and neuronal viability were investigated in the cultured neurons, subjected to acute or chronic excitotoxicity. We also analyzed the effect of one of these peptides in the mouse model of stroke by measuring infarct size and evaluating motor coordination and balance.

Results: Calpain cleaves three interdomain linker regions in PSD-95 and produces stable fragments corresponding to previously described PSD-95 supramodules (PDZ1-2 and P-S-G) as well as a truncated form SH3-GK. Peptide TP95₄₁₄, containing the cleavage site in the PDZ3-SH3 linker, is able to interfere PSD-95 downregulation and reduces neuronal death by excitotoxicity. Additionally, TP95₄₁₄ is delivered to mice cortex and, in a severe model of permanent ischemia, significantly improves the neurological outcome after brain damage.

Conclusions: Interference of excitotoxicity-induced PSD-95-processing with specific CPPs constitutes a novel and promising therapeutic approach for stroke treatment.

RssB-mediated σS Activation is Regulated by a Two-Tier Mechanism via Phosphorylation and Adaptor Protein - IraD

Regulation of bacterial stress responding σS is a sophisticated process and mediated by multiple interacting partners. Controlled proteolysis of σS is regulated by RssB which maintains minimal level of σS during exponential growth but then elevates σS level while facing stresses. Bacteria developed different strategies to regulate activity of RssB, including phosphorylation of itself and production of anti-adaptors. However, the function of phosphorylation is controversial and the mechanism of anti-adaptors preventing RssB-σS interaction remains elusive. Here, we demonstrated the impact of phosphorylation on the activity of RssB and built the RssB-σS complex model. Importantly, we showed that the phosphorylation site - D58 is at the interface of RssB-σS complex. Hence, mutation or phosphorylation of D58 would weaken the interaction of RssB with σS. We found that the anti-adaptor protein IraD has higher affinity than σS to RssB and its binding interface on RssB overlaps with that for σS. And IraD-RssB complex is preferred over RssB-σS in solution, regardless of the phosphorylation state of RssB. Our study suggests that RssB possesses a two-tier mechanism for regulating σS. First, phosphorylation of RssB provides a moderate and reversible tempering of its activity, followed by a specific and robust inhibition via the anti-adaptor interaction.

Profiling of human lymphocytes reveals a specific network of protein kinases modulated by endurance training status

To date, the effects of endurance exercise training on lymphocyte physiology at the kinome level are largely unknown. Therefore, the present study used a highly sensitive peptide-based kinase activity profiling approach to investigate if the basal activity of tyrosine (Tyr) and serine/threonine (Ser/Thr) kinases of human lymphocytes is affected by the aerobic endurance training status. Results revealed that the activity of various tyrosine kinases of the FGFR family and ZAP70 was increased, whereas the activity of multiple Ser/Thr kinases such as IKKα, CaMK4, PKAα, PKCα+δ (among others) was decreased in lymphocytes of endurance trained athletes (ET). Moreover, functional associations between several differentially regulated kinases in ET-derived lymphocytes were demonstrated by phylogenetic mapping and network analysis. Especially, Ser/Thr kinases of the AGC-kinase (protein kinase A, G, and C) family represent exercise-sensitive key components within the lymphocytes kinase network that may mediate the long-term effects of endurance training. Furthermore, KEGG (Kyoto Encyclopedia of Genes and Genomes) and Reactome pathway analysis indicate that Ras as well as intracellular signaling by second messengers were found to be enriched in the ET individuals. Overall, our data suggest that endurance exercise training improves the adaptive immune competence by modulating the activity of multiple protein kinases in human lymphocytes.

Serotonin and serotonin reuptake inhibitors alter placental aromatase

Serotonin reuptake inhibitors (SRIs) are currently the main molecules prescribed to pregnant women that suffer from depression. Placental cells are exposed to SRIs via maternal blood, and we have previously shown that SRIs alter feto-placental steroidogenesis in an in vitro co-culture model. More specifically, serotonin (5-HT) regulates the estrogen biosynthetic enzyme aromatase (cytochrome P450 19; CYP19), which is disrupted by fluoxetine and its active metabolite norfluoxetine in BeWo choriocarcinoma cells. Based on molecular simulations, the present study illustrates that the SRIs fluoxetine, norfluoxetine, paroxetine, sertraline, citalopram and venlafaxine exhibit binding affinity for the active-site pocket of CYP19, suggesting potential competitive inhibition. Using BeWo cells and primary villous trophoblast cells isolated from normal term placentas, we compared the effects of the SRIs on CYP19 activity. We observed that paroxetine and sertraline induce aromatase activity in BeWo cells, while venlafaxine, fluoxetine, paroxetine and sertraline decrease aromatase activity in primary villous trophoblast. The effects of paroxetine and sertraline in primary villous trophoblasts were observed at the lower doses tested. We also showed that 5-HT and the 5-HT_2A receptor agonist 2,5-dimethoxy-4-iodoamphetamine (DOI) induced CYP19 activity. An increase in phosphorylation of serine and tyrosine and a decrease in threonine phosphorylation of CYP19 was also associated with DOI treatment. Our results contribute to better understanding how 5-HT and SRIs interact with CYP19 and may affect estrogen production. Moreover, this study suggests that alteration of placental 5-HT levels due to depression and/or SRI treatment during pregnancy may be associated with disruption of placental estrogen production.

Distinct Primary Nucleation of Polymorphic Aβ Dimers Yields to Distinguished Fibrillation Pathways

Polymorphic Aβ dimers are the smallest toxic species that play a role in the pathology of Alzheimer's disease. There is great interest in understanding the malfunctions that yield to these toxic species and in providing insights into the molecular mechanisms of the primary nucleation. Herein, we present a first work that demonstrates two distant edges states of Aβ dimers. The first is the so-called "random coil" state dimer that mimics the primary seeding/nucleation that is far from a fibrillation state. The second is the "fibril-like" state dimer that is structurally in close proximity to the fibril, a well-organized state into a fibril-like structure. We show for the first time that a conformational change of one monomer within the dimer impedes primary nucleation, while less fluctuations and relatively large number of interactions in nucleation domains induce the primary nucleation to produce toxic stable species. Overall, the current study exhibits a diversity of primary nucleation in each dimer state, suggesting distinct molecular mechanisms of fibril formation. The conformations of the early stage Aβ dimers that were achieved may provide crucial data for designing inhibitors to impede the primary nucleation.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Structural basis of βTrCP1-associated GLI3 processing

Controlled ubiquitin-mediated protein degradation is essential for various cellular processes. GLI family regulates the transcriptional events of the sonic hedgehog pathway genes that are implicated in almost one fourth of human tumors. GLI3 phosphorylation by Ser/Thr kinases is a primary factor for their transcriptional activity that incurs the formation of both GLI3 repressor and activator forms. GLI3 processing is triggered in an ubiquitin-dependent manner via SCF^βTrCP1 complex; however, structural characterization, mode of action based on sequence of phosphorylation signatures and induced conformational readjustments remain elusive. Here, through structural analysis and molecular dynamics simulation assays, we explored comparative binding pattern of GLI3 phosphopeptides against βTrCP1. A comprehensive and thorough analysis demarcated GLI3 presence in the binding cleft shared by inter-bladed binding grooves of β-propeller. Our results revealed the involvement of all seven WD40 repeats of βTrCP1 in GLI3 interaction. Conversely, GLI3 phosphorylation pattern at primary protein kinase A (PKA) sites and secondary casein kinase 1 (CK1) or glycogen synthase kinase 3 (GSK3) sites was carefully evaluated. Our results indicated that GLI3 processing depends on the 19 phosphorylation sites (849, 852, 855, 856, 860, 861, 864, 865, 868, 872, 873, 876, 877, 880, 899, 903, 906, 907 and 910 positions) by a cascade of PKA, GSK3β and CSKI kinases. The presence of a sequential phosphorylation in the binding induction of GLI3 and βTrCP1 may be a hallmark to authenticate GLI3 processing. We speculate that mechanistic information of the individual residual contributions through structure-guided approaches may be pivotal for the rational design of specific and more potent inhibitors against activated GLI3 with a special emphasis on the anticancer activity.

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.

KCC2 membrane diffusion tunes neuronal chloride homeostasis

Neuronal Cl^- homeostasis is regulated by the activity of two cation chloride co-transporters (CCCs), the K⁺-Cl^- cotransporter KCC2 and the Na⁺-K⁺-Cl^- cotransporter NKCC1, which are primarily extruding and importing chloride in neurons, respectively. Several neurological and psychiatric disorders including epilepsy, neuropathic pain, schizophrenia and autism are associated with altered neuronal chloride (Cl^-) homeostasis. A current view is that the accumulation of intracellular Cl^- in neurons as a result of KCC2 down-regulation and/or NKCC1 up-regulation may weaken inhibitory GABA signaling and thereby promote the development of pathological activities. CCC activity is determined mainly by their level of expression in the plasma membrane. Furthermore, CCCs undergo "diffusion-trapping" in the membrane, a mechanism that is rapidly adjusted by activity-dependent post-translational modifications i.e. phosphorylation/dephosphorylation of key serine and threonine residues. This represents probably the most rapid cellular mechanism for adapting CCC function to changes in neuronal activity. Therefore, interfering with these mechanisms may help restoring Cl^- homeostasis and inhibition under pathological conditions. This article is part of the special issue entitled 'Mobility and trafficking of neuronal membrane proteins'.

How phosphorylation influences E1 subunit pyruvate dehydrogenase: A computational study

Pyruvate (PYR) dehydrogenase complex (PDC) is an enzymatic system that plays a crucial role in cellular metabolism as it controls the entry of carbon into the Krebs cycle. From a structural point of view, PDC is formed by three different subunits (E1, E2 and E3) capable of catalyzing the three reaction steps necessary for the full conversion of pyruvate to acetyl-CoA. Recent investigations pointed out the crucial role of this enzyme in the replication and survival of specific cancer cell lines, renewing the interest of the scientific community. Here, we report the results of our molecular dynamics studies on the mechanism by which posttranslational modifications, in particular the phosphorylation of three serine residues (Ser-264-α, Ser-271-α, and Ser-203-α), influence the enzymatic function of the protein. Our results support the hypothesis that the phosphorylation of Ser-264-α and Ser-271-α leads to (1) a perturbation of the catalytic site structure and dynamics and, especially in the case of Ser-264-α, to (2) a reduction in the affinity of E1 for the substrate. Additionally, an analysis of the channels connecting the external environment with the catalytic site indicates that the inhibitory effect should not be due to the occlusion of the access/egress pathways to/from the active site.

Molecular Dynamics Simulation for All

The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.

Phospho-regulation of intrinsically disordered proteins for actin assembly and endocytosis

Actin filament assembly contributes to the endocytic pathway pleiotropically, with active roles in clathrin-dependent and clathrin-independent endocytosis as well as subsequent endosomal trafficking. Endocytosis comprises a series of dynamic events, including the initiation of membrane curvature, bud invagination, vesicle abscission and subsequent vesicular transport. The ultimate success of endocytosis requires the coordinated activities of proteins that trigger actin polymerization, recruit actin-binding proteins (ABPs) and organize endocytic proteins (EPs) that promote membrane curvature through molecular crowding or scaffolding mechanisms. A particularly interesting phenomenon is that multiple EPs and ABPs contain a substantial percentage of intrinsically disordered regions (IDRs), which can contribute to protein coacervation and phase separation. In addition, intrinsically disordered proteins (IDPs) frequently contain sites for post-translational modifications (PTMs) such as phosphorylation, and these modifications exhibit a high preference for IDR residues [Groban ES et al. (2006) PLoS Comput Biol 2, e32]. PTMs are implicated in regulating protein function by modulating the protein conformation, protein-protein interactions and the transition between order and disorder states of IDPs. The molecular mechanisms by which IDRs of ABPs and EPs fine-tune actin assembly and endocytosis remain mostly unexplored and elusive. In this review, we analyze protein sequences of budding yeast EPs and ABPs, and discuss the potential underlying mechanisms for regulating endocytosis and actin assembly through the emerging concept of IDR-mediated protein multivalency, coacervation, and phase transition, with an emphasis on the phospho-regulation of IDRs. Finally, we summarize the current understanding of how these mechanisms coordinate actin cytoskeleton assembly and membrane curvature formation during endocytosis in budding yeast.

Dissecting Structure-Encoded Determinants of Allosteric Cross-Talk between Post-Translational Modification Sites in the Hsp90 Chaperones

Post-translational modifications (PTMs) represent an important regulatory instrument that modulates structure, dynamics and function of proteins. The large number of PTM sites in the Hsp90 proteins that are scattered throughout different domains indicated that synchronization of multiple PTMs through a combinatorial code can be invoked as an important mechanism to orchestrate diverse chaperone functions and recognize multiple client proteins. In this study, we have combined structural and coevolutionary analysis with molecular simulations and perturbation response scanning analysis of the Hsp90 structures to characterize functional role of PTM sites in allosteric regulation. The results reveal a small group of conserved PTMs that act as global mediators of collective dynamics and allosteric communications in the Hsp90 structures, while the majority of flexible PTM sites serve as sensors and carriers of the allosteric structural changes. This study provides a comprehensive structural, dynamic and network analysis of PTM sites across Hsp90 proteins, identifying specific role of regulatory PTM hotspots in the allosteric mechanism of the Hsp90 cycle. We argue that plasticity of a combinatorial PTM code in the Hsp90 may be enacted through allosteric coupling between effector and sensor PTM residues, which would allow for timely response to structural requirements of multiple modified enzymes.

GABAA receptor dependent synaptic inhibition rapidly tunes KCC2 activity via the Cl--sensitive WNK1 kinase

The K+-Cl- co-transporter KCC2 (SLC12A5) tunes the efficacy of GABAA receptor-mediated transmission by regulating the intraneuronal chloride concentration [Cl-]i. KCC2 undergoes activity-dependent regulation in both physiological and pathological conditions. The regulation of KCC2 by synaptic excitation is well documented; however, whether the transporter is regulated by synaptic inhibition is unknown. Here we report a mechanism of KCC2 regulation by GABAA receptor (GABAAR)-mediated transmission in mature hippocampal neurons. Enhancing GABAAR-mediated inhibition confines KCC2 to the plasma membrane, while antagonizing inhibition reduces KCC2 surface expression by increasing the lateral diffusion and endocytosis of the transporter. This mechanism utilizes Cl- as an intracellular secondary messenger and is dependent on phosphorylation of KCC2 at threonines 906 and 1007 by the Cl--sensing kinase WNK1. We propose this mechanism contributes to the homeostasis of synaptic inhibition by rapidly adjusting neuronal [Cl-]i to GABAAR activity.

Real-Time Multiplex Kinase Phosphorylation Sensors in Living Cells

Phosphorylation is an important post-translational modification implicated in cellular signaling and regulation. However, current methods to study protein phosphorylation by various kinases lack spatiotemporal resolution or the ability to simultaneously observe in real time the activity of multiple kinases in live cells. We present a peptide biosensor strategy with time correlated single photon counting-fluorescence lifetime imaging (TCSPC-FLIM) to interrogate the spatial and temporal dynamics of VEGFR-2 and AKT phosphorylation activity in real time in live 2D and 3D cell culture models at single cell resolution. By recording the increase in fluorescence lifetime due to a change in the solvatochromic environment of the sensor upon phosphorylation, we demonstrate that spatiotemporal maps of protein kinase activity can be obtained. Our results suggest that fluorescence lifetime imaging of peptide biosensors can be effectively and specifically used to monitor and quantify phosphorylation of multiple kinases in live cells.

Regulation of Cx37 channel and growth-suppressive properties by phosphorylation

Growth suppression mediated by connexin 37 (Cx37; also known as GJA4) requires interaction between its C-terminus and functional pore-forming domain. Using rat insulinoma cells, we show that Cx37 induces cell death and cell cycle arrest, and slowed cell cycling. Whether differential phosphorylation might regulate intramolecular interactions, and consequently the growth-suppressive phenotype, is unknown. Protein kinase C inhibition increased the open state probability of low-conductance gap junction channels (GJChs) and reduced GJCh closed state probability. Substituting alanine at serine residues 275, 302 and 328 eliminated Cx37-induced cell death, supported proliferation and reduced the GJCh closed state probability. With additional alanine for serine substitutions at residues 285, 319, 321 and 325, Cx37-induced cell death was eliminated and the growth arrest period prolonged, and GJCh closed state probability was restored. With aspartate substitution at these seven sites, apoptosis was induced and the open state probability of large conductance GJChs (and hemichannels) was increased. These data suggest that differential phosphorylation of the C-terminus regulates channel conformation and, thereby, cell cycle progression and cell survival.

Regulation of the Equilibrium between Closed and Open Conformations of Annexin A2 by N-Terminal Phosphorylation and S100A4-Binding

Annexin A2 (ANXA2) has a versatile role in membrane-associated functions including membrane aggregation, endo- and exocytosis, and it is regulated by post-translational modifications and protein-protein interactions through the unstructured N-terminal domain (NTD). Our sequence analysis revealed a short motif responsible for clamping the NTD to the C-terminal core domain (CTD). Structural studies indicated that the flexibility of the NTD and CTD are interrelated and oppositely regulated by Tyr24 phosphorylation and Ser26Glu phosphomimicking mutation. The crystal structure of the ANXA2-S100A4 complex showed that asymmetric binding of S100A4 induces dislocation of the NTD from the CTD and, similar to the Ser26Glu mutation, unmasks the concave side of ANXA2. In contrast, pTyr24 anchors the NTD to the CTD and hampers the membrane-bridging function. This inhibition can be restored by S100A4 and S100A10 binding. Based on our results we provide a structural model for regulation of ANXA2-mediated membrane aggregation by NTD phosphorylation and S100 binding.

A phosphorylation-motif for tuneable helix stabilisation in intrinsically disordered proteins - Lessons from the sodium proton exchanger 1 (NHE1)

Intrinsically disordered proteins (IDPs) are involved in many pivotal cellular processes including phosphorylation and signalling. The structural and functional effects of phosphorylation of IDPs remain poorly understood and difficult to predict. Thus, a need exists to identify motifs that confer phosphorylation-dependent perturbation of the local preferences for forming e.g. helical structures as well as motifs that do not. The disordered distal tail of the Na⁺/H⁺ exchanger 1 (NHE1) is six-times phosphorylated (S693, S723, S726, S771, T779, S785) by the mitogen activated protein kinase 2 (MAPK1, ERK2). Using NMR spectroscopy, we found that two out of those six phosphorylation sites had a stabilizing effect on transient helices. One of these was further investigated by circular dichroism and NMR spectroscopy as well as by molecular dynamic simulations, which confirmed the stabilizing effect and resulted in the identification of a short linear motif for helix stabilisation: [S/T]-P-{3}-[R/K] where [S/T] is the phosphorylation-site. By analysing IDP and phosphorylation site databases we found that the motif is significantly enriched around known phosphorylation sites, supporting a potential wider-spread role in phosphorylation-mediated regulation of intrinsically disordered proteins. The identification of such motifs is important for understanding the molecular mechanism of cellular signalling, and is crucial for the development of predictors for the structural effect of phosphorylation; a tool of relevance for understanding disease-promoting mutations that for example interfere with signalling for instance through constitutive active and often cancer-promoting signalling.

Understanding the role of protein interaction motifs in transcriptional regulators: implications for crop improvement

Recently, the conjunction of disciplines such as developmental biology and proteomics enabled the dissection of diverse cellular processes, by analysis of their transcriptional regulatory pathways. In particular, it has been shown that transcription factor interactions play critical roles in the development of many complex traits and control cellular phenotypic plasticity, whereas protein phosphorylation modifications regulate protein activity at the posttranslational level. The present work posits that protein-protein interactions by functional motifs, as well as the phosphorylation state in these sites, are fundamental plant biological phenotype determinants, whose elucidation and understanding will allow manipulation of complex traits, thereby contributing to the design of novel methodologies for molecular breeders and plant physiologists.

Dynamic, structural and thermodynamic basis of insulin-like growth factor 1 kinase allostery mediated by activation loop phosphorylation

Despite the importance of kinases' catalytic activity regulation in cell signaling, detailed mechanisms underlying their activity regulation are poorly understood. Herein, using insulin-like growth factor 1 receptor kinase (IGF-1RK) as a model, the mechanisms of kinase regulation by its activation loop (A-loop) phosphorylation were investigated through molecular dynamics (MD) and alchemical free energy simulations. Analyses of the simulation results and free energy landscapes determined for the entire catalytic cycle of the kinase revealed that A-loop phosphorylation affects each step in the IGF-1RK catalytic cycle, including conformational change, substrate binding/product release and catalytic phosphoryl transfer. Specifically, the conformational equilibrium of the kinase is shifted by 13.2 kcal mol-1 to favor the active conformation after A-loop phosphorylation, which increases substrate binding affinity of the activated kinase. This free energy shift is achieved primarily via destabilization of the inactive conformation. The free energy of the catalytic reaction is also changed by 3.3 kcal mol-1 after the phosphorylation and in the end, facilitates product release. Analyses of MD simulations showed that A-loop phosphorylation produces these energetic effects by perturbing the side chain interactions around each A-loop tyrosine. These interaction changes are propagated to the remainder of the kinase to modify the orientations and dynamics of the αC-helix and A-loop, and together yield the observed free energy changes. Since many protein kinases share similar interactions identified in this work, the mechanisms of kinase allostery and catalysis unraveled here can be applicable to them.

Didehydrophenylalanine, an abundant modification in the beta subunit of plant polygalacturonases

The structure and the activity of proteins are often regulated by transient or stable post- translational modifications (PTM). Different from well-known, abundant modifications such as phosphorylation and glycosylation some modifications are limited to one or a few proteins across a broad range of related species. Although few examples of the latter type are known, the evolutionary conservation of these modifications and the enzymes responsible for their synthesis suggest an important physiological role. Here, the first observation of a new, fold-directing PTM is described. During the analysis of alfalfa cell wall proteins a -2Da mass shift was observed on phenylalanine residues in the repeated tetrapeptide FxxY of the beta-subunit of polygalacturonase. This modular protein is known to be involved in developmental and stress-responsive processes. The presence of this modification was confirmed using in-house and external datasets acquired by different commonly used techniques in proteome studies. Based on these analyses it was found that all identified phenylalanine residues in the sequence FxxY of this protein were modified to α,β-didehydro-Phe (ΔPhe). Besides showing the reproducible identification of ΔPhe in different species arguments that substantiate the fold-determining role of ΔPhe are given.

Use of restrained molecular dynamics to predict the conformations of phosphorylated receiver domains in two-component signaling systems

Two‐component signaling (TCS) is the primary means by which bacteria, as well as certain plants and fungi, respond to external stimuli. Signal transduction involves stimulus‐dependent autophosphorylation of a sensor histidine kinase and phosphoryl transfer to the receiver domain of a downstream response regulator. Phosphorylation acts as an allosteric switch, inducing structural and functional changes in the pathway's components. Due to their transient nature, phosphorylated receiver domains are challenging to characterize structurally. In this work, we provide a methodology for simulating receiver domain phosphorylation to predict conformations that are nearly identical to experimental structures. Using restrained molecular dynamics, phosphorylated conformations of receiver domains can be reliably sampled on nanosecond timescales. These simulations also provide data on conformational dynamics that can be used to identify regions of functional significance related to phosphorylation. We first validated this approach on several well‐characterized receiver domains and then used it to compare the upstream and downstream components of the fungal Sln1 phosphorelay. Our results demonstrate that this technique provides structural insight, obtained in the absence of crystallographic or NMR information, regarding phosphorylation‐induced conformational changes in receiver domains that regulate the output of their associated signaling pathway. To our knowledge, this is the first time such a protocol has been described that can be broadly applied to TCS proteins for predictive purposes. Proteins 2016; 85:155–176. © 2016 Wiley Periodicals, Inc.

Molecular Dynamics Analysis of 4E-BP2 Protein Fold Stabilization Induced by Phosphorylation

Protein phosphorylation can affect the interaction with partner proteins but can also induce conformational transitions. In case of the eukaryotic translation initiation factor 4E-binding protein 2 (4E-BP2) threonine (Thr) phosphorylation at two turn motifs results in transition from a disordered to a folded structure. In order to elucidate the stabilizing mechanism we employed comparative molecular dynamics (MD) free energy simulations on the turn motifs indicating that Thr-phosphorylation favors a folded whereas dephosphorylation or substitution by Glu residues destabilizes the turn structure. In multiple unrestrained MD simulations at elevated temperature of the 4E-BP2 domain only the double phosphorylated variant remained close to the folded structure in agreement with experiment. Three surface Arg residues were identified as additional key elements for the tertiary structure stabilization of the whole phosphorylated domain. In addition to the local turn structure double phosphorylation also leads to an overall electrostatic stabilization of the folded form compared to wild type and other investigated variants of 4E-BP2. The principles of phosphorylation mediated fold stabilization identified in the present study may also be helpful for identifying other structural motifs that can be affected by phosphorylation or provide a route to design such motifs.

Phosphorylation of the Brome Mosaic Virus Capsid Regulates the Timing of Viral Infection

The four brome mosaic virus (BMV) RNAs (RNA1 to RNA4) are encapsidated in three distinct virions that have different disassembly rates in infection. The mechanism for the differential release of BMV RNAs from virions is unknown, since 180 copies of the same coat protein (CP) encapsidate each of the BMV genomic RNAs. Using mass spectrometry, we found that the BMV CP contains a complex pattern of posttranslational modifications. Treatment with phosphatase was found to not significantly affect the stability of the virions containing RNA1 but significantly impacted the stability of the virions that encapsidated BMV RNA2 and RNA3/4. Cryo-electron microscopy reconstruction revealed dramatic structural changes in the capsid and the encapsidated RNA. A phosphomimetic mutation in the flexible N-terminal arm of the CP increased BMV RNA replication and virion production. The degree of phosphorylation modulated the interaction of CP with the encapsidated RNA and the release of three of the BMV RNAs. UV cross-linking and immunoprecipitation methods coupled to high-throughput sequencing experiments showed that phosphorylation of the BMV CP can impact binding to RNAs in the virions, including sequences that contain regulatory motifs for BMV RNA gene expression and replication. Phosphatase-treated virions affected the timing of CP expression and viral RNA replication in plants. The degree of phosphorylation decreased when the plant hosts were grown at an elevated temperature. These results show that phosphorylation of the capsid modulates BMV infection.

IMPORTANCE

How icosahedral viruses regulate the release of viral RNA into the host is not well understood. The selective release of viral RNA can regulate the timing of replication and gene expression. Brome mosaic virus (BMV) is an RNA virus, and its three genomic RNAs are encapsidated in separate virions. Through proteomic, structural, and biochemical analyses, this work shows that posttranslational modifications, specifically, phosphorylation, on the capsid protein regulate the capsid-RNA interaction and the stability of the virions and affect viral gene expression. Mutational analysis confirmed that changes in modification affected virion stability and the timing of viral infection. The mechanism for modification of the virion has striking parallels to the mechanism of regulation of chromatin packaging by nucleosomes.

Impact of Ser17 Phosphorylation on the Conformational Dynamics of the Oncoprotein MDM2

MDM2 is an important oncoprotein that downregulates the activity of the tumor suppressor protein p53 via binding of its N-terminal domain to the p53 transactivation domain. The first 24 residues of the MDM2 N-terminal domain form an intrinsically disordered "lid" region that interconverts on a millisecond time scale between "open" and "closed" states in unliganded MDM2. While the former conformational state is expected to facilitate p53 binding, the latter competes in a pseudo-substrate manner with p53 for its binding site. Phosphorylation of serine 17 in the MDM2 lid region is thought to modulate the equilibrium between "open" and "closed" lid states, but contradictory findings on the favored lid conformational state upon phosphorylation have been reported. Here, the nature of the conformational states of MDM2 pSer17 and Ser17Asp variants was addressed by means of enhanced sampling molecular dynamics simulations. Detailed analyses of the computed lid conformational ensembles indicate that both lid variants stabilize a "closed" state, with respect to wild type. Nevertheless, the nature of the closed-state conformational ensembles differs significantly between the pSer17 and Ser17Asp variants. Thus, care should be applied in the interpretation of biochemical experiments that use phosphomimetic variants to model the effects of phosphorylation on the structure and dynamics of this disordered protein region.

Fimbrin phosphorylation by metaphase Cdk1 regulates actin cable dynamics in budding yeast

Actin cables, composed of actin filament bundles nucleated by formins, mediate intracellular transport for cell polarity establishment and maintenance. We previously observed that metaphase cells preferentially promote actin cable assembly through cyclin-dependent kinase 1 (Cdk1) activity. However, the relevant metaphase Cdk1 targets were not known. Here we show that the highly conserved actin filament crosslinking protein fimbrin is a critical Cdk1 target for actin cable assembly regulation in budding yeast. Fimbrin is specifically phosphorylated on threonine 103 by the metaphase cyclin–Cdk1 complex, in vivo and in vitro. On the basis of conformational simulations, we suggest that this phosphorylation stabilizes fimbrin's N-terminal domain, and modulates actin filament binding to regulate actin cable assembly and stability in cells. Overall, this work identifies fimbrin as a key target for cell cycle regulation of actin cable assembly in budding yeast, and suggests an underlying mechanism.

Metaphase cells preferentially promote actin cable assembly through cyclin-dependent kinase 1 (Cdk1) activity. Here the authors identify fimbrin as one of the main metaphase Cdk1 targets for cell cycle regulation of actin cable assembly in budding yeast.

Monte Carlo Sampling with Hierarchical Move Sets: POSH Monte Carlo

We present a new Monte Carlo method for sampling rugged energy landscapes that allows for efficient transitions across sparsely distributed local basins. The trial move consists of two steps. The first step is a large initial trial move, and the second step is a Monte Carlo trajectory generated using smaller trial moves. To maintain detailed balance, a reverse transition probability is estimated along a path that differs from the forward path. Since the forward and reverse transitions are different, we name the algorithm POSH (port out, starboard home) Monte Carlo. The process obeys detailed balance to the extent that the transition probabilities are correctly estimated. There is an optimal range of performance for a given energy landscape, which depends on how sparsely the low energy states of the system are distributed. For simple model systems, adequate precision is obtained over a large range of inner steps settings. Side chain sampling of residues in the binding region of progesterone antibody 1dba are studied, and show that significant improvement over a comparable standard protocol can be obtained using POSH sampling. To compare with experimental data, the phosphopeptide Ace-Gly-Ser-pSer-Ser-Nma is also studied, and the resulting NMR observables compare well with experiment. For the biomolecular systems studied, we show that POSH sampling generates precise distributions using the number of inner steps set up to 20.

E2 superfamily of ubiquitin-conjugating enzymes: constitutively active or activated through phosphorylation in the catalytic cleft

Protein phosphorylation is a modification that offers a dynamic and reversible mechanism to regulate the majority of cellular processes. Numerous diseases are associated with aberrant regulation of phosphorylation-induced switches. Phosphorylation is emerging as a mechanism to modulate ubiquitination by regulating key enzymes in this pathway. The molecular mechanisms underpinning how phosphorylation regulates ubiquitinating enzymes, however, are elusive. Here, we show the high conservation of a functional site in E2 ubiquitin-conjugating enzymes. In catalytically active E2s, this site contains aspartate or a phosphorylatable serine and we refer to it as the conserved E2 serine/aspartate (CES/D) site. Molecular simulations of substrate-bound and -unbound forms of wild type, mutant and phosphorylated E2s, provide atomistic insight into the role of the CES/D residue for optimal E2 activity. Both the size and charge of the side group at the site play a central role in aligning the substrate lysine toward E2 catalytic cysteine to control ubiquitination efficiency. The CES/D site contributes to the fingerprint of the E2 superfamily. We propose that E2 enzymes can be divided into constitutively active or regulated families. E2s characterized by an aspartate at the CES/D site signify constitutively active E2s, whereas those containing a serine can be regulated by phosphorylation.

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.

Prediction of Functionally Important Phospho-Regulatory Events in Xenopus laevis Oocytes

The African clawed frog Xenopus laevis is an important model organism for studies in developmental and cell biology, including cell-signaling. However, our knowledge of X. laevis protein post-translational modifications remains scarce. Here, we used a mass spectrometry-based approach to survey the phosphoproteome of this species, compiling a list of 2636 phosphosites. We used structural information and phosphoproteomic data for 13 other species in order to predict functionally important phospho-regulatory events. We found that the degree of conservation of phosphosites across species is predictive of sites with known molecular function. In addition, we predicted kinase-protein interactions for a set of cell-cycle kinases across all species. The degree of conservation of kinase-protein interactions was found to be predictive of functionally relevant regulatory interactions. Finally, using comparative protein structure models, we find that phosphosites within structured domains tend to be located at positions with high conformational flexibility. Our analysis suggests that a small class of phosphosites occurs in positions that have the potential to regulate protein conformation.

Towards the computational design of protein post-translational regulation

Protein post-translational modifications (PTMs) are a fast and versatility mechanism used by the cell to regulate the function of proteins in response to changing conditions. PTMs can alter the activity of proteins by allosteric regulation or by controlling protein interactions, localization and abundance. Recent advances in proteomics have revealed the extent of regulation by PTMs and the different mechanisms used in nature to exert control over protein function via PTMs. These developments can serve as the foundation for the rational design of protein regulation. Here we review the advances in methods to determine the function of PTMs, protein allosteric control and examples of rational design of PTM regulation. These advances create an opportunity to move synthetic biology forward by making use of a level of regulation that is of yet unexplored.

Molecular mechanisms of asymmetric RAF dimer activation

Protein phosphorylation is one of the most common post-translational modifications in cell regulatory mechanisms. Dimerization plays also a crucial role in the kinase activity of many kinases, including RAF, CDK2 (cyclin-dependent kinase 2) and EGFR (epidermal growth factor receptor), with heterodimers often being the most active forms. However, the structural and mechanistic details of how phosphorylation affects the activity of homo- and hetero-dimers are largely unknown. Experimentally, synthesizing protein samples with fully specified and homogeneous phosphorylation states remains a challenge for structural biology and biochemical studies. Typically, multiple changes in phosphorylation lead to activation of the same protein, which makes structural determination methods particularly difficult. It is also not well understood how the occurrence of phosphorylation and dimerization processes synergize to affect kinase activities. In the present article, we review available structural data and discuss how MD simulations can be used to model conformational transitions of RAF kinase dimers, in both their phosphorylated and unphosphorylated forms.

Function of the Herpes Simplex Virus 1 Small Capsid Protein VP26 Is Regulated by Phosphorylation at a Specific Site

Replacement of the herpes simplex virus 1 small capsid protein VP26 phosphorylation site Thr-111 with alanine reduced viral replication and neurovirulence to levels observed with the VP26 null mutation. This mutation reduced VP26 expression and mislocalized VP26 and its binding partner, the major capsid protein VP5, in the nucleus. VP5 mislocalization was also observed with the VP26 null mutation. Thus, we postulate that phosphorylation of VP26 at Thr-111 regulates VP26 function in vitro and in vivo.