1
|
Marien J, Prévost C, Sacquin-Mora S. nP-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat. J Chem Inf Model 2024. [PMID: 39092904 DOI: 10.1021/acs.jcim.4c00742] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Tau is an intrinsically disordered (IDP) microtubule-associated protein (MAP) that plays a key part in microtubule assembly and organization. The function of tau can be regulated by multiple phosphorylation sites. These post-translational modifications are known to decrease the binding affinity of tau for microtubules, and abnormal tau phosphorylation patterns are involved in Alzheimer's disease. Using all-atom molecular dynamics simulations, we compared the conformational landscapes explored by the tau R2 repeat domain (which comprises a strong tubulin binding site) in its native state and with multiple phosphorylations on the S285, S289, and S293 residues, with four different standard force field (FF)/water model combinations. We find that the different parameters used for the phosphate groups (which can be more or less flexible) in these FFs and the specific interactions between bulk cations and water lead to the formation of a specific type of counterion bridge, termed nP-collab (for nphosphate collaboration, with n being an integer), where counterions form stable structures binding with two or three phosphate groups simultaneously. The resulting effect of nP-collabs on the tau-R2 conformational space differs when using sodium or potassium cations and is likely to impact the peptide overall dynamics and how this MAP interacts with tubulins. We also investigated the effect of phosphoresidue spacing and ionic concentration by modeling polyalanine peptides containing two phosphoserines located one-six residues apart. Three new metrics specifically tailored for IDPs (proteic Menger curvature, local curvature, and local flexibility) were introduced, which allow us to fully characterize the impact of nP-collabs on the dynamics of disordered peptides at the residue level.
Collapse
Affiliation(s)
- Jules Marien
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Chantal Prévost
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Sophie Sacquin-Mora
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
2
|
Bickel D, Vranken W. Effects of Phosphorylation on Protein Backbone Dynamics and Conformational Preferences. J Chem Theory Comput 2024; 20:4998-5011. [PMID: 38830621 PMCID: PMC11210476 DOI: 10.1021/acs.jctc.4c00206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Revised: 05/22/2024] [Accepted: 05/22/2024] [Indexed: 06/05/2024]
Abstract
Phosphorylations are the most common and extensively studied post-translational modification (PTM) of proteins in eukaryotes. They constitute a major regulatory mechanism, modulating protein function, protein-protein interactions, as well as subcellular localization. Phosphorylation sites are preferably located in intrinsically disordered regions and have been shown to trigger structural rearrangements and order-to-disorder transitions. They can therefore have a significant effect on protein backbone dynamics or conformation, but only sparse experimental data are available. To obtain a more general description of how and when phosphorylations have a significant effect on protein behavior, molecular dynamics (MD) currently provides the only suitable framework to study these effects at a large scale in atomistic detail. This study develops a systematic MD simulation framework to explore the influence of phosphorylations on the local backbone dynamics and conformational propensities of proteins. Through a series of glycine-backbone peptides, we studied the effects of amino acid residues including the three most common phosphorylations (Ser, Thr, and Tyr), on local backbone dynamics and conformational propensities. We further extended our study to investigate the interactions of all such residues between position i to positions i + 1, i + 2, i + 3, and i + 4 in such peptides. The final data set comprises structural ensembles for 3393 sequences with more than 1 μs of sampling for each ensemble. To validate the relevance of the results, the structural and conformational properties extracted from the MD simulations are compared to NMR data from the Biological Magnetic Resonance Data Bank. The systematic nature of this study enables the projection of the gained knowledge onto any phosphorylation site in the proteome and provides a general framework for the study of further PTMs. The full data set is publicly available, as a training and reference set.
Collapse
Affiliation(s)
- David Bickel
- Interuniversity
Institute of Bioinformatics in Brussels, 1050 Brussels, Belgium
- Structural
Biology Brussels, Vrije Universiteit Brussels, 1050 Brussels, Belgium
| | - Wim Vranken
- Interuniversity
Institute of Bioinformatics in Brussels, 1050 Brussels, Belgium
- Structural
Biology Brussels, Vrije Universiteit Brussels, 1050 Brussels, Belgium
| |
Collapse
|
3
|
Usher ET, Fossat MJ, Holehouse AS. Phosphorylation of disordered proteins tunes local and global intramolecular interactions. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.10.598315. [PMID: 38915510 PMCID: PMC11195077 DOI: 10.1101/2024.06.10.598315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared to published findings. We simulated short (< 20 residues) and long (> 80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multi-phosphorylated IDRs. Our results support and expand on prior observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.
Collapse
Affiliation(s)
- Emery T. Usher
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, USA
- Center for Biomolecular Condensates (CBC), Washington University in St. Louis, St. Louis, MO, USA
| | - Martin J. Fossat
- Department of Biological Physics, Max Planck Institute of Immunobiology and Epigenetics, Freiburg im Breisgau, Germany
| | - Alex S. Holehouse
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, USA
- Center for Biomolecular Condensates (CBC), Washington University in St. Louis, St. Louis, MO, USA
| |
Collapse
|
4
|
Barbhuiya TK, Jayarathna DK, Gilmour R, Smet-Nocca C, Gandhi NS. Characterization of Posttranslationally Modified PHF-1 Tau Peptides Using Gaussian Accelerated Molecular Dynamics Simulation. Methods Mol Biol 2024; 2754:3-31. [PMID: 38512658 DOI: 10.1007/978-1-0716-3629-9_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
The microtubule-associated protein, Tau, is an intrinsically disordered protein that plays a crucial role in neurodegenerative diseases like Alzheimer's disease. The posttranslational modifications across the Tau protein domains are involved in regulating Tau protein's function and disease onset. Of the various posttranslational modifications at Ser, Thr, and Tyr sites, O-GlcNAcylation and phosphorylation are the most critical ones, playing a vital role in Tau aggregation and tauopathies. To understand the function, it is essential to characterize the structural changes associated with Tau modification. Previous experimental studies have focused on high-resolution nuclear magnetic resonance techniques to structurally characterize the effect of phosphorylation, O-GlcNAcylation, and combination of both PTMs on Tau conformation in small peptides centered on the PHF-1 epitope from amino acid 392 to 411. The structural characterization using atomistic molecular dynamics simulation of such disordered peptides requires long simulation time, proper sampling method, and utilization of appropriate force fields for accurate determination of conformational ensembles, resembling the experimental data. This chapter details the protocol for the structural characterization of modified Tau peptides using the CHARMM36m force field and enhanced sampling methods like Gaussian accelerated molecular dynamics (GaMD) simulation. We have focused on a detailed explanation of the GaMD method and analyses of molecular dynamics trajectories to explain the relationship between two modifications, phospho- and glyco-, at C-terminus of Tau protein and its stable conformation over the longer simulation timeframes. The analyses involve energetics reweighting, clustering of simulation trajectories, and characterization of secondary structure using circular dichroism data from the simulation. The reader can utilize this protocol to investigate the structures of complex proteins, especially the disordered ones.
Collapse
Affiliation(s)
- Tabassum Khair Barbhuiya
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia
- Cancer and Ageing Research Program, Translational Research Institute, Woolloongabba, QLD, Australia
| | - Dulari K Jayarathna
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia
| | - Raechelle Gilmour
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia
| | - Caroline Smet-Nocca
- University of Lille, Inserm, CHU Lille, Institut Pasteur de Lille, U1167 - RID-AGE - Risk Factors and Molecular Determinants of Aging-Related Diseases, Lille, France
- CNRS EMR9002 Integrative Structural Biology, Lille, France
| | - Neha S Gandhi
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia.
- Cancer and Ageing Research Program, Translational Research Institute, Woolloongabba, QLD, Australia.
| |
Collapse
|
5
|
Amith W, Dutagaci B. Complex Conformational Space of the RNA Polymerase II C-Terminal Domain upon Phosphorylation. J Phys Chem B 2023; 127:9223-9235. [PMID: 37870995 PMCID: PMC10626582 DOI: 10.1021/acs.jpcb.3c02655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Revised: 10/03/2023] [Indexed: 10/25/2023]
Abstract
Intrinsically disordered proteins (IDPs) have been closely studied during the past decade due to their importance in many biological processes. The disordered nature of this group of proteins makes it difficult to observe its full span of the conformational space using either experimental or computational studies. In this article, we explored the conformational space of the C-terminal domain (CTD) of RNA polymerase II (Pol II), which is also an intrinsically disordered low complexity domain, using enhanced sampling methods. We provided a detailed conformational analysis of model systems of CTD with different lengths; first with the last 44 residues of the human CTD sequence and finally the CTD model with 2-heptapeptide repeating units. We then investigated the effects of phosphorylation on CTD conformations by performing simulations at different phosphorylated states. We obtained broad conformational spaces in nonphosphorylated CTD models, and phosphorylation has complex effects on the conformations of the CTD. These complex effects depend on the length of the CTD, spacing between the multiple phosphorylation sites, ion coordination, and interactions with the nearby residues.
Collapse
Affiliation(s)
- Weththasinghage
D. Amith
- Department of Molecular and
Cell Biology, University of California,
Merced, Merced, California 95343, United States
| | - Bercem Dutagaci
- Department of Molecular and
Cell Biology, University of California,
Merced, Merced, California 95343, United States
| |
Collapse
|
6
|
Wang X, Wang Y, Guo M, Wang X, Li Y, Zhang JZH. Assessment of an Electrostatic Energy-Based Charge Model for Modeling the Electrostatic Interactions in Water Solvent. J Chem Theory Comput 2023; 19:6294-6312. [PMID: 37656610 DOI: 10.1021/acs.jctc.3c00467] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/03/2023]
Abstract
The protein force field based on the restrained electrostatic potential (RESP) charges has limitations in accurately describing hydrogen bonding interactions in proteins. To address this issue, we propose an alternative approach called the electrostatic energy-based charges (EEC) model, which shows improved performance in describing electrostatic interactions (EIs) of hydrogen bonds in proteins. In this study, we further investigate the performance of the EEC model in modeling EIs in water solvent. Our findings demonstrate that the fixed EEC model can effectively reproduce the quantum mechanics/molecular mechanics (QM/MM)-calculated EIs between a water molecule and various water solvent environments. However, to achieve the same level of computational accuracy, the electrostatic potential (ESP) charge model needs to fluctuate according to the electrostatic environment. Our analysis indicates that the requirement for charge adjustments depends on the specific mathematical and physical representation of EIs as a function of the environment for deriving charges. By comparing with widely used empirical water models calibrated to reproduce experimental properties, we confirm that the performance of the EEC model in reproducing QM/MM EIs is similar to that of general purpose TIP4P-like water models such as TIP4P-Ew and TIP4P/2005. When comparing the computed 10,000 distinct EI values within the range of -40 to 0 kcal/mol with the QM/MM results calculated at the MP2/aug-cc-pVQZ/TIP3P level, we noticed that the mean unsigned error (MUE) for the EEC model is merely 0.487 kcal/mol, which is remarkably similar to the MUE values of the TIP4P-Ew (0.63 kcal/mol) and TIP4P/2005 (0.579 kcal/mol) models. However, both the RESP method and the TIP3P model exhibit a tendency to overestimate the EIs, as evidenced by their higher MUE values of 1.761 and 1.293 kcal/mol, respectively. EEC-based molecular dynamics simulations have demonstrated that, when combined with appropriate van der Waals parameters, the EEC model can closely reproduce oxygen-oxygen radial distribution function and density of water, showing a remarkable similarity to the well-established TIP4P-like empirical water models. Our results demonstrate that the EEC model has the potential to build force fields with comparable accuracy to more sophisticated empirical TIP4P-like water models.
Collapse
Affiliation(s)
- Xianwei Wang
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Yiying Wang
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Man Guo
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Xuechao Wang
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Yang Li
- College of Information Science and Engineering, Shandong Agricultural University, Tai'an, Shandong 271018, China
| | - John Z H Zhang
- Shenzhen Institute of Synthetic Biology, Faculty of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong 518055, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
7
|
Song G, Zhong B, Zhang B, Rehman AU, Chen HF. Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins. J Chem Inf Model 2023; 63:1602-1614. [PMID: 36800279 DOI: 10.1021/acs.jcim.3c00112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Phosphorylation of proteins plays an important regulatory role at almost all levels of cellular organization. Molecular dynamics (MD) simulation is a promising tool to reveal the mechanism of how phosphorylation regulates many key biological processes at the atomistic level. MD simulation accuracy depends on force field precision, while the current force fields for phospho-amino acids have resulted in notable inconsistency with experimental data. Here, a new force field parameter (named FB18CMAP) is generated by fitting against quantum mechanics (QM) energy in aqueous solution with φ/ψ dihedral potential-energy surfaces optimized using CMAP parameters. MD simulations of phosphorylated dipeptides, intrinsically disordered proteins (IDPs), and ordered (folded) proteins show that FB18CMAP can mimic NMR observables and structural characteristics of phosphorylated dipeptides and proteins more accurately than the FB18 force field. These findings suggest that FB18CMAP performs well in both the simulation of ordered and disordered states of phosphorylated proteins.
Collapse
Affiliation(s)
- Ge Song
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Bozitao Zhong
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Bo Zhang
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ashfaq Ur Rehman
- Departments of Molecular Biology and Biochemistry, University of California, Irvine, California 92697, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China.,Shanghai Center for Bioinformation Technology, Shanghai 200240, China
| |
Collapse
|
8
|
How phosphorylation impacts intrinsically disordered proteins and their function. Essays Biochem 2022; 66:901-913. [PMID: 36350035 PMCID: PMC9760426 DOI: 10.1042/ebc20220060] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 10/14/2022] [Accepted: 10/17/2022] [Indexed: 11/10/2022]
Abstract
Phosphorylation is the most common post-translational modification (PTM) in eukaryotes, occurring particularly frequently in intrinsically disordered proteins (IDPs). These proteins are highly flexible and dynamic by nature. Thus, it is intriguing that the addition of a single phosphoryl group to a disordered chain can impact its function so dramatically. Furthermore, as many IDPs carry multiple phosphorylation sites, the number of possible states increases, enabling larger complexities and novel mechanisms. Although a chemically simple and well-understood process, the impact of phosphorylation on the conformational ensemble and molecular function of IDPs, not to mention biological output, is highly complex and diverse. Since the discovery of the first phosphorylation site in proteins 75 years ago, we have come to a much better understanding of how this PTM works, but with the diversity of IDPs and their capacity for carrying multiple phosphoryl groups, the complexity grows. In this Essay, we highlight some of the basic effects of IDP phosphorylation, allowing it to serve as starting point when embarking on studies into this topic. We further describe how recent complex cases of multisite phosphorylation of IDPs have been instrumental in widening our view on the effect of protein phosphorylation. Finally, we put forward perspectives on the phosphorylation of IDPs, both in relation to disease and in context of other PTMs; areas where deep insight remains to be uncovered.
Collapse
|
9
|
Balanced Force Field ff03CMAP Improving the Dynamics Conformation Sampling of Phosphorylation Site. Int J Mol Sci 2022; 23:ijms231911285. [PMID: 36232586 PMCID: PMC9569523 DOI: 10.3390/ijms231911285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 09/08/2022] [Accepted: 09/11/2022] [Indexed: 11/30/2022] Open
Abstract
Phosphorylation plays a key role in plant biology, such as the accumulation of plant cells to form the observed proteome. Statistical analysis found that many phosphorylation sites are located in disordered regions. However, current force fields are mainly trained for structural proteins, which might not have the capacity to perfectly capture the dynamic conformation of the phosphorylated proteins. Therefore, we evaluated the performance of ff03CMAP, a balanced force field between structural and disordered proteins, for the sampling of the phosphorylated proteins. The test results of 11 different phosphorylated systems, including dipeptides, disordered proteins, folded proteins, and their complex, indicate that the ff03CMAP force field can better sample the conformations of phosphorylation sites for disordered proteins and disordered regions than ff03. For the solvent model, the results strongly suggest that the ff03CMAP force field with the TIP4PD water model is the best combination for the conformer sampling. Additional tests of CHARMM36m and FB18 force fields on two phosphorylated systems suggest that the overall performance of ff03CMAP is similar to that of FB18 and better than that of CHARMM36m. These results can help other researchers to choose suitable force field and solvent models to investigate the dynamic properties of phosphorylation proteins.
Collapse
|
10
|
Sluchanko NN. Recent advances in structural studies of 14-3-3 protein complexes. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2022; 130:289-324. [PMID: 35534110 DOI: 10.1016/bs.apcsb.2021.12.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
Being phosphopeptide-binding hubs, 14-3-3 proteins coordinate multiple cellular processes in eukaryotes, including the regulation of apoptosis, cell cycle, ion channels trafficking, transcription, signal transduction, and hormone biosynthesis. Forming constitutive α-helical dimers, 14-3-3 proteins predominantly recognize specifically phosphorylated Ser/Thr sites within their partners; this generally stabilizes phosphotarget conformation and affects its activity, intracellular distribution, dephosphorylation, degradation and interactions with other proteins. Not surprisingly, 14-3-3 complexes are involved in the development of a range of diseases and are considered promising drug targets. The wide interactome of 14-3-3 proteins encompasses hundreds of different phosphoproteins, for many of which the interaction is well-documented in vitro and in vivo but lack the structural data that would help better understand underlying regulatory mechanisms and develop new drugs. Despite obtaining structural information on 14-3-3 complexes is still lagging behind the research of 14-3-3 interactions on a proteome-wide scale, recent works provided some advances, including methodological improvements and accumulation of new interesting structural data, that are discussed in this review.
Collapse
Affiliation(s)
- Nikolai N Sluchanko
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology" of the Russian Academy of Sciences, Moscow, Russian Federation.
| |
Collapse
|
11
|
Li L, Wang N, Xiong Y, Guo G, Zhu M, Gu Y. Transcription Factor FOSL1 Enhances Drug Resistance of Breast Cancer Through DUSP7-Mediated Dephosphorylation of PEA15. Mol Cancer Res 2021; 20:515-526. [PMID: 34907034 DOI: 10.1158/1541-7786.mcr-21-0658] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Revised: 11/02/2021] [Accepted: 12/07/2021] [Indexed: 11/16/2022]
Abstract
Breast cancer (BC) represents one of the commonest and deadliest malignancies in women. However, drug resistance has always been a major obstacle to cancer treatment. Transcription factors have been reported to have close association with drug resistance of tumors. Recently, by analyzing the data from Gene Expression Omnibus (GEO) database (id: GSE76540), we found that transcription factor FOSL1 was significantly up-regulated in the transcriptome of doxorubicin-resistant BC cells compared with that in sensitive parental cells. Therefore, we aim to explore the regulatory mechanism of FOSL1 in affecting the drug resistance of BC cells. FOSL1 expression in doxorubicin-resistant BC cells was firstly examined through RT-qPCR, and then its influence on the drug resistance of BC cells was explored through a series of in vitro and in vivo mechanism assays. Results showed that FOSL1 promoted the drug resistance of BC cells to doxorubicin both in intro and in vivo. It positively regulated the transcription of DUSP7 in BC doxorubicin-resistant cells and DUSP7 also enhanced the drug resistance of BC cells. Furthermore, FOSL1 promoted the dephosphorylation of PEA15 through DUSP7. In conclusion, it was verified that FOSL1 promoted the drug resistance in breast cancer through DUSP7-mediated dephosphorylation of PEA15. Implications: These initial findings suggest that the FOSL1/DUSP7/PEA15 pathway may provide a theoretical guidance for BC treatment.
Collapse
Affiliation(s)
- Lin Li
- First Affiliated Hospital of Zhengzhou University
| | - Nan Wang
- First Affiliated Hospital of Zhengzhou University
| | - Youyi Xiong
- First Affiliated Hospital of Zhengzhou University
| | | | - Mingzhi Zhu
- First Affiliated Hospital of Zhengzhou University
| | - Yuanting Gu
- Department of Breast Surgery, First Affiliated Hospital of Zhengzhou University
| |
Collapse
|
12
|
Stoppelman JP, Ng TT, Nerenberg PS, Wang LP. Development and Validation of AMBER-FB15-Compatible Force Field Parameters for Phosphorylated Amino Acids. J Phys Chem B 2021; 125:11927-11942. [PMID: 34668708 DOI: 10.1021/acs.jpcb.1c07547] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Phosphorylation of select amino acid residues is one of the most common biological mechanisms for regulating protein structures and functions. While computational modeling can be used to explore the detailed structural changes associated with phosphorylation, most molecular mechanics force fields developed for the simulation of phosphoproteins have been noted to be inconsistent with experimental data. In this work, we parameterize force fields for the phosphorylated forms of the amino acids serine, threonine, and tyrosine using the ForceBalance software package with the goal of improving agreement with experiments for these residues. Our optimized force field, denoted as FB18, is parameterized using high-quality ab initio potential energy scans and is designed to be fully compatible with the AMBER-FB15 protein force field. When utilized in MD simulations together with the TIP3P-FB water model, we find that FB18 consistently enhances the prediction of experimental quantities such as 3J NMR couplings and intramolecular hydrogen-bonding propensities in comparison to previously published models. As was reported with AMBER-FB15, we also see improved agreement with the reference QM calculations in regions at and away from local minima. We thus believe that the FB18 parameter set provides a promising route for the further investigation of the varied effects of protein phosphorylation.
Collapse
Affiliation(s)
- John P Stoppelman
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Tracey T Ng
- Department of Physics & Astronomy, California State University, Los Angeles, California 90032, United States
| | - Paul S Nerenberg
- Department of Physics & Astronomy, California State University, Los Angeles, California 90032, United States.,Department of Biological Sciences, California State University, Los Angeles, California 90032, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, California 95616, United States
| |
Collapse
|
13
|
Rieloff E, Skepö M. The Effect of Multisite Phosphorylation on the Conformational Properties of Intrinsically Disordered Proteins. Int J Mol Sci 2021; 22:11058. [PMID: 34681718 PMCID: PMC8541499 DOI: 10.3390/ijms222011058] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 10/07/2021] [Accepted: 10/08/2021] [Indexed: 02/06/2023] Open
Abstract
Intrinsically disordered proteins are involved in many biological processes such as signaling, regulation, and recognition. A common strategy to regulate their function is through phosphorylation, as it can induce changes in conformation, dynamics, and interactions with binding partners. Although phosphorylated intrinsically disordered proteins have received increased attention in recent years, a full understanding of the conformational and structural implications of phosphorylation has not yet been achieved. Here, we present all-atom molecular dynamics simulations of five disordered peptides originated from tau, statherin, and β-casein, in both phosphorylated and non-phosphorylated state, to compare changes in global dimensions and structural elements, in an attempt to gain more insight into the controlling factors. The changes are in qualitative agreement with experimental data, and we observe that the net charge is not enough to predict the impact of phosphorylation on the global dimensions. Instead, the distribution of phosphorylated and positively charged residues throughout the sequence has great impact due to the formation of salt bridges. In statherin, a preference for arginine-phosphoserine interaction over arginine-tyrosine accounts for a global expansion, despite a local contraction of the phosphorylated region, which implies that also non-charged residues can influence the effect of phosphorylation.
Collapse
Affiliation(s)
- Ellen Rieloff
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden;
| | - Marie Skepö
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden;
- LINXS—Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
14
|
Nassar R, Dignon GL, Razban RM, Dill KA. The Protein Folding Problem: The Role of Theory. J Mol Biol 2021; 433:167126. [PMID: 34224747 PMCID: PMC8547331 DOI: 10.1016/j.jmb.2021.167126] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Revised: 06/21/2021] [Accepted: 06/26/2021] [Indexed: 10/20/2022]
Abstract
The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.
Collapse
Affiliation(s)
- Roy Nassar
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA; Department of Chemistry, Stony Brook University, Stony Brook, NY, USA
| | - Gregory L Dignon
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Rostam M Razban
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA; Department of Chemistry, Stony Brook University, Stony Brook, NY, USA; Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY, USA.
| |
Collapse
|
15
|
Rieloff E, Skepö M. Molecular Dynamics Simulations of Phosphorylated Intrinsically Disordered Proteins: A Force Field Comparison. Int J Mol Sci 2021; 22:10174. [PMID: 34576338 PMCID: PMC8470740 DOI: 10.3390/ijms221810174] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2021] [Revised: 09/15/2021] [Accepted: 09/16/2021] [Indexed: 12/21/2022] Open
Abstract
Phosphorylation is a common post-translational modification among intrinsically disordered proteins and regions, which helps regulate function by changing the protein conformations, dynamics, and interactions with binding partners. To fully comprehend the effects of phosphorylation, computer simulations are a helpful tool, although they are dependent on the accuracy of the force field used. Here, we compared the conformational ensembles produced by Amber ff99SB-ILDN+TIP4P-D and CHARMM36m, for four phosphorylated disordered peptides ranging in length from 14-43 residues. CHARMM36m consistently produced more compact conformations with a higher content of bends, mainly due to more stable salt bridges. Based on comparisons with experimental size estimates for the shortest and longest peptide, CHARMM36m appeared to overestimate the compactness. The difference between the force fields was largest for the peptide showing the greatest separation between positively charged and phosphorylated residues, in line with the importance of charge distribution. For this peptide, the conformational ensemble did not change significantly upon increasing the ionic strength from 0 mM to 150 mM, despite a reduction of the salt-bridging probability in the CHARMM36m simulations, implying that salt concentration has negligible effects in this study.
Collapse
Affiliation(s)
- Ellen Rieloff
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden;
| | - Marie Skepö
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden;
- LINXS—Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
16
|
How multisite phosphorylation impacts the conformations of intrinsically disordered proteins. PLoS Comput Biol 2021; 17:e1008939. [PMID: 33945530 PMCID: PMC8148376 DOI: 10.1371/journal.pcbi.1008939] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2020] [Revised: 05/25/2021] [Accepted: 04/06/2021] [Indexed: 12/15/2022] Open
Abstract
Phosphorylation of intrinsically disordered proteins (IDPs) can produce changes in structural and dynamical properties and thereby mediate critical biological functions. How phosphorylation effects intrinsically disordered proteins has been studied for an increasing number of IDPs, but a systematic understanding is still lacking. Here, we compare the collapse propensity of four disordered proteins, Ash1, the C-terminal domain of RNA polymerase (CTD2’), the cytosolic domain of E-Cadherin, and a fragment of the p130Cas, in unphosphorylated and phosphorylated forms using extensive all-atom molecular dynamics (MD) simulations. We find all proteins to show V-shape changes in their collapse propensity upon multi-site phosphorylation according to their initial net charge: phosphorylation expands neutral or overall negatively charged IDPs and shrinks positively charged IDPs. However, force fields including those tailored towards and commonly used for IDPs overestimate these changes. We find quantitative agreement of MD results with SAXS and NMR data for Ash1 and CTD2’ only when attenuating protein electrostatic interactions by using a higher salt concentration (e.g. 350 mM), highlighting the overstabilization of salt bridges in current force fields. We show that phosphorylation of IDPs also has a strong impact on the solvation of the protein, a factor that in addition to the actual collapse or expansion of the IDP should be considered when analyzing SAXS data. Compared to the overall mild change in global IDP dimension, the exposure of active sites can change significantly upon phosphorylation, underlining the large susceptibility of IDP ensembles to regulation through post-translational modifications. Intrinsically disordered proteins (IDPs) are a class of proteins that lack secondary and tertiary structures and instead explore a broad conformational ensemble. Their functions, from transcriptional regulation to signal transmission, are tightly regulated. IDPs are subject of extensive reversible post-translational modifications (PTMs), such as phosphorylation, methylation and glycosylation. Among these PTMs, phosphorylation is one of the most common and important PTMs. However, the mechanism of how phosphorylation affects the conformations and functions of IDPs remains unclear. To answer this question, we have performed extensive all-atom molecular dynamics simulations for four representative IDPs: Ash1, E-Cadherin, CTD2’ and p130Cas in their unphosphorylated and phosphorylated forms. Our results showed that all IDPs undergo a mild change upon multi-site phosphorylation, which is V-shaped: phosphorylation moderately expands neutral or overall negatively charged IDPs and shrinks positively charged IDPs. More importantly, in two of these IDPs, only two biologically relevant phosphorylation sites suffice to render the adjacent negatively charged active site significantly more exposed to the environment, which implies a higher probability to interact with other binding partners.
Collapse
|
17
|
Argudo PG, Giner-Casares JJ. Folding and self-assembly of short intrinsically disordered peptides and protein regions. NANOSCALE ADVANCES 2021; 3:1789-1812. [PMID: 36133101 PMCID: PMC9417027 DOI: 10.1039/d0na00941e] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Accepted: 01/17/2021] [Indexed: 05/15/2023]
Abstract
Proteins and peptide fragments are highly relevant building blocks in self-assembly for nanostructures with plenty of applications. Intrinsically disordered proteins (IDPs) and protein regions (IDRs) are defined by the absence of a well-defined secondary structure, yet IDPs/IDRs show a significant biological activity. Experimental techniques and computational modelling procedures for the characterization of IDPs/IDRs are discussed. Directed self-assembly of IDPs/IDRs allows reaching a large variety of nanostructures. Hybrid materials based on the derivatives of IDPs/IDRs show a promising performance as alternative biocides and nanodrugs. Cell mimicking, in vivo compartmentalization, and bone regeneration are demonstrated for IDPs/IDRs in biotechnological applications. The exciting possibilities of IDPs/IDRs in nanotechnology with relevant biological applications are shown.
Collapse
Affiliation(s)
- Pablo G Argudo
- Université de Bordeaux, CNRS, Bordeaux INP, LCPO 16 Avenue Pey-Berland 33600 Pessac France
| | - Juan J Giner-Casares
- Departamento de Química Física y T. Aplicada, Instituto Universitario de Nanoquímica IUNAN, Facultad de Ciencias, Universidad de Córdoba (UCO) Campus de Rabanales, Ed. Marie Curie E-14071 Córdoba Spain
| |
Collapse
|
18
|
Wang W. Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins. Phys Chem Chem Phys 2021; 23:777-784. [PMID: 33355572 DOI: 10.1039/d0cp05818a] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.
Collapse
Affiliation(s)
- Wenning Wang
- Department of Chemistry, Multiscale Research Institute of Complex Systems and Institute of Biomedical Sciences, Fudan University, Shanghai 200438, China.
| |
Collapse
|
19
|
Fagerberg E, Månsson LK, Lenton S, Skepö M. The Effects of Chain Length on the Structural Properties of Intrinsically Disordered Proteins in Concentrated Solutions. J Phys Chem B 2020; 124:11843-11853. [PMID: 33337879 PMCID: PMC7872433 DOI: 10.1021/acs.jpcb.0c09635] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Intrinsically disordered proteins (IDP) are proteins that sample
a heterogeneous ensemble of conformers in solution. An estimated 25–30%
of all eukaryotic proteins belong to this class. In vivo, IDPs function under conditions that are highly crowded by other
biological macromolecules. Previous research has highlighted that
the presence of crowding agents can influence the conformational ensemble
sampled by IDPs, resulting in either compaction or expansion. The
effects of self-crowding of the disordered protein Histatin 5 has,
in an earlier study, been found to have limited influence on the conformational
ensemble. In this study, it is examined whether the short chain length
of Histatin 5 can explain the limited effects of crowding observed,
by introducing (Histatin 5)2, a tandem repeat of Histatin
5. By utilizing small-angle X-ray scattering, it is shown that the
conformational ensemble is conserved at high protein concentrations,
in resemblance with Histatin 5, although with a lowered protein concentration
at which aggregation arises. Under dilute conditions, atomistic molecular
dynamics and coarse-grained Monte Carlo simulations, as well as an
established scaling law, predicted more extended conformations than
indicated by experimental data, hence implying that (Histatin 5)2 does not behave as a self-avoiding random walk.
Collapse
Affiliation(s)
- Eric Fagerberg
- Theoretical Chemistry, Lund University, P.O. Box 124, Lund SE-221 00, Sweden
| | - Linda K Månsson
- Theoretical Chemistry, Lund University, P.O. Box 124, Lund SE-221 00, Sweden
| | - Samuel Lenton
- Theoretical Chemistry, Lund University, P.O. Box 124, Lund SE-221 00, Sweden.,LINXS - Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, Lund SE-223 70, Sweden
| | - Marie Skepö
- Theoretical Chemistry, Lund University, P.O. Box 124, Lund SE-221 00, Sweden.,LINXS - Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, Lund SE-223 70, Sweden
| |
Collapse
|
20
|
Zheng W, Dignon G, Brown M, Kim YC, Mittal J. Hydropathy Patterning Complements Charge Patterning to Describe Conformational Preferences of Disordered Proteins. J Phys Chem Lett 2020; 11:3408-3415. [PMID: 32227994 PMCID: PMC7450210 DOI: 10.1021/acs.jpclett.0c00288] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Understanding the conformational ensemble of an intrinsically disordered protein (IDP) is of great interest due to its relevance to critical intracellular functions and diseases. It is now well established that the polymer scaling behavior can provide a great deal of information about the conformational properties as well as liquid-liquid phase separation of an IDP. It is, therefore, extremely desirable to be able to predict an IDP's scaling behavior from the protein sequence itself. The work in this direction so far has focused on highly charged proteins and how charge patterning can perturb their structural properties. As naturally occurring IDPs are composed of a significant fraction of uncharged amino acids, the rules based on charge content and patterning are only partially helpful in solving the problem. Here, we propose a new order parameter, sequence hydropathy decoration, which can provide a near-quantitative understanding of scaling and structural properties of IDPs devoid of charged residues. We combine this with a charge patterning parameter, sequence charge decoration, to obtain a general equation, parametrized from extensive coarse-grained simulation data, for predicting protein dimensions from the sequence. We finally test this equation against available experimental data and find a semiquantitative match in predicting the scaling behavior. We also provide guidance on how to extend this approach to experimental data, which should be feasible in the near future.
Collapse
Affiliation(s)
- Wenwei Zheng
- College of Integrative Sciences and Arts, Arizona State University, Mesa, Arizona 85212, United States
| | - Gregory Dignon
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Matthew Brown
- College of Integrative Sciences and Arts, Arizona State University, Mesa, Arizona 85212, United States
| | - Young C Kim
- Center for Materials Physics and Technology, Naval Research Laboratory, Washington, D.C. 20375, United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|