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; 64:6570-6582. [PMID: 39092904 DOI: 10.1021/acs.jcim.4c00742] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [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
|
Huang C, Pettitt BM. Parameter Dependence of the Solubility Limit for Disodium Phosphate. J Phys Chem B 2023; 127:8690-8696. [PMID: 37788371 PMCID: PMC10913426 DOI: 10.1021/acs.jpcb.3c05343] [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] [Indexed: 10/05/2023]
Abstract
The solubility limit was calculated for supersaturated solutions of disodium phosphate in water as a function of the sodium-oxygen Lennard-Jones radius parameter Rmin. We found that changes in the sodium-oxygen Rmin were clearly exponentially related to the concentration of the solubility limit. Starting from standard force fields more suited to nucleic acids and phospholipids, only relatively small changes were required to achieve the experimentally known solubility limit. Simultaneously, we found that it was possible to achieve the solubility limit and the osmotic pressure with the same model parameters. Based on transferability, the adjusted Rmin parameter can be used to more accurately model phosphorylated proteins.
Collapse
Affiliation(s)
- Caleb Huang
- Biochemistry and Molecular Biology, University of Texas Medical Branch, Galveston, Texas 77555-0304, United States
| | - B Montgomery Pettitt
- Biochemistry and Molecular Biology, University of Texas Medical Branch, Galveston, Texas 77555-0304, United States
| |
Collapse
|
5
|
Šinko G. Modeling of a near-attack conformation of oxime in phosphorylated acetylcholinesterase via a reactivation product, a phosphorylated oxime. Chem Biol Interact 2023; 383:110656. [PMID: 37579936 DOI: 10.1016/j.cbi.2023.110656] [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] [Received: 02/09/2023] [Revised: 07/26/2023] [Accepted: 08/09/2023] [Indexed: 08/16/2023]
Abstract
At the present, only four antidotes are in use in therapy for poisoning by organophosphorus compounds: 2-PAM, HI-6, obidoxime and trimedoxime. Numerous compounds have been designed and synthetized to be more effective reactivators than those currently in use. Many of those new compounds fail at the enzyme level because interactions formed within the AChE active site are not favourable ones that lead to a successful reactivation. The approach in which the modeling of a phosphorylated oxime (POX), a product of successful reactivation in the AChE active site, may be a way to better understand the role of active site residues during the process of formation of the Michaelis type of complex between an enzyme and oxime. After reactivation, a change in phosphorus stereochemistry occurs leading to a different spatial arrangement of attached substituents, now including an oxime. To study interactions between the AChE oxyanion hole and a phosphorylated oxime, an S203G mutant was used to avoid the steric hindrance caused by the catalytic serine. In this way, the POX could be positioned close to the oxyanion hole. In the final step, the oxime without a phosphoester moiety was transferred into the phosphorylated AChE and molecular dynamics was used to test the stability of the near-attack conformation of the oxime near the phosphorylated serine.
Collapse
Affiliation(s)
- Goran Šinko
- Institute for Medical Research and Occupational Health, Ksaverska cesta 2, HR-10 000, Zagreb, Croatia.
| |
Collapse
|
6
|
Pandey AK, Ganguly HK, Sinha SK, Daniels KE, Yap GPA, Patel S, Zondlo NJ. An Inherent Difference between Serine and Threonine Phosphorylation: Phosphothreonine Strongly Prefers a Highly Ordered, Compact, Cyclic Conformation. ACS Chem Biol 2023; 18:1938-1958. [PMID: 37595155 DOI: 10.1021/acschembio.3c00068] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/20/2023]
Abstract
Phosphorylation and dephosphorylation of proteins by kinases and phosphatases are central to cellular responses and function. The structural effects of serine and threonine phosphorylation were examined in peptides and in proteins, by circular dichroism, NMR spectroscopy, bioinformatics analysis of the PDB, small-molecule X-ray crystallography, and computational investigations. Phosphorylation of both serine and threonine residues induces substantial conformational restriction in their physiologically more important dianionic forms. Threonine exhibits a particularly strong disorder-to-order transition upon phosphorylation, with dianionic phosphothreonine preferentially adopting a cyclic conformation with restricted ϕ (ϕ ∼ -60°) stabilized by three noncovalent interactions: a strong intraresidue phosphate-amide hydrogen bond, an n → π* interaction between consecutive carbonyls, and an n → σ* interaction between the phosphate Oγ lone pair and the antibonding orbital of C-Hβ that restricts the χ2 side-chain conformation. Proline is unique among the canonical amino acids for its covalent cyclization on the backbone. Phosphothreonine can mimic proline's backbone cyclization via noncovalent interactions. The preferred torsions of dianionic phosphothreonine are ϕ,ψ = polyproline II helix > α-helix (ϕ ∼ -60°); χ1 = g-; χ2 ∼ +115° (eclipsed C-H/O-P bonds). This structural signature is observed in diverse proteins, including in the activation loops of protein kinases and in protein-protein interactions. In total, these results suggest a structural basis for the differential use and evolution of threonine versus serine phosphorylation sites in proteins, with serine phosphorylation typically inducing smaller, rheostat-like changes, versus threonine phosphorylation promoting larger, step function-like switches, in proteins.
Collapse
Affiliation(s)
- Anil K Pandey
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Himal K Ganguly
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Sudipta Kumar Sinha
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
- Department of Chemistry, Indian Institute of Technology Ropar, Nangal Road, Rupnagar 140001, India
| | - Kelly E Daniels
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Glenn P A Yap
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Sandeep Patel
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Neal J Zondlo
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| |
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
|
Cool AM, Lindert S. Umbrella Sampling Simulations Measure Switch Peptide Binding and Hydrophobic Patch Opening Free Energies in Cardiac Troponin. J Chem Inf Model 2022; 62:5666-5674. [PMID: 36283742 PMCID: PMC9712266 DOI: 10.1021/acs.jcim.2c00508] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The cardiac troponin (cTn) complex is an important regulatory protein in heart contraction. Upon binding of Ca2+, cTn undergoes a conformational shift that allows the troponin I switch peptide (cTnISP) to be released from the actin filament and bind to the troponin C hydrophobic patch (cTnCHP). Mutations and modifications to this complex can change its sensitivity to Ca2+ and alter the energetics of the transition from the Ca2+-unbound, cTnISP-unbound form to the Ca2+-bound, cTnISP-bound form. We utilized targeted molecular dynamics (TMD) to obtain a trajectory of this transition pathway, followed by umbrella sampling to estimate the free energy associated with the cTnISP-cTnCHP binding and the cTnCHP opening events for wild-type (WT) cTn. We were able to reproduce experimental values for the cTnISP-cTnCHP binding event and obtain cTnCHP opening free energies in agreement with previous computational measurements of smaller cTnC systems. This excellent agreement for WT cTn demonstrated the strength of computational methods in studying the dynamics and energetics of the cTn complex. We then introduced mutations to the cTn complex that cause cardiomyopathy or alter its Ca2+ sensitivity and observed a general decrease in the free energy of opening the cTnCHP. For these same mutations, we observed no general trend in the effect on the cTnISP-cTnCHP binding event. Our method sets the stage for future computational studies on this system that predict the consequences of yet uncharacterized mutations on cTn dynamics and energetics.
Collapse
Affiliation(s)
- Austin M. Cool
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, 43210
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, 43210
| |
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
|
Künzel N, Helms V. How Peptides Bind to PSD-95/Discs-Large/ZO-1 Domains. J Chem Theory Comput 2022; 18:3845-3859. [PMID: 35608157 DOI: 10.1021/acs.jctc.1c01140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
PSD-95/discs-large/ZO-1 (PDZ) domains form a large family of adaptor proteins that bind to the C-terminal tails of their binding partner proteins. Via extensive molecular dynamics simulations and alchemical free energy calculations, we characterized the binding modi of phosphorylated and unphosphorylated EQVSAV peptides and of a EQVEAV phosphate mimic to the hPTP1E PDZ2 and MAGI1 PDZ1 domains. The simulations reproduced the well-known binding characteristics such as tight coordination of the peptidic carboxyl tail and pronounced hydrogen bonding between the peptide backbone and the backbone atoms of a β-sheet in PDZ. Overall, coordination by hPTP1E PDZ2 appeared tighter than by MAGI1 PDZ1. Simulations of wild-type PDZ and arginine mutants suggest that contacts with Arg79/85 in hPTP1E/MAGI1 are more important for the EQVEAV peptide than for EQVSAV. Alchemical free energy calculations and PaCS-MD simulations could well reproduce the difference in binding free energy between unphosphorylated EQVSAV and EQVEAV peptides and the absolute binding free energy of EQVSAV. However, likely due to small force field inaccuracies, the simulations erroneously favored binding of the phosphorylated peptide instead of its unphosphorylated counterpart, which is in contrast to the experiment.
Collapse
Affiliation(s)
- Nicolas Künzel
- Center for Bioinformatics, Saarland University, P.O. Box 15 11 50, D-66041 Saarbrücken, Germany
| | - Volkhard Helms
- Center for Bioinformatics, Saarland University, P.O. Box 15 11 50, D-66041 Saarbrücken, Germany
| |
Collapse
|
11
|
Rani L, Mallajosyula SS. Site-Specific Stabilization and Destabilization of α Helical Peptides upon Phosphorylation and O-GlcNAcylation. J Phys Chem B 2021; 125:13444-13459. [PMID: 34870441 DOI: 10.1021/acs.jpcb.1c09419] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Helices (α-helix) are the most common type of secondary structure motif present in proteins. In this study, we have investigated the structural influence of phosphorylation and O-GlcNAcylation, common intracellular post-translational modifications (PTMs), on the α-helical conformation. The simulation studies were performed on the Baldwin model α-helical peptide sequence (Ac-AKAAAAKAAAAKAA-NH2). The Baldwin sequences were chosen due to the availability of site-specific experimental post-translational data for cross-validation with the simulations. The influence of PTMs was examined across the span of the α-helix, namely, at the N-terminus, position 10 (interior region), and the C-terminus for both serine and threonine residues placed at these positions. Molecular dynamics (MD) simulations revealed that phosphorylation and O-GlcNAcylation at the N-terminus lead to the stabilization of the helical conformation. PTMs in the interior or the C-terminus were found to disrupt helicity, with the disruption being more pronounced for PTMs in the interior region, in accordance with experimental studies. It was found that phosphorylation-derived destabilization was mainly due to the formation of an intraresidue HN-PO32- electrostatic interaction and interactions between the phosphate group and the side chain of adjacent lysine residues (NH3···PO32-). Hydrophobic and steric clashes were the main causes of destabilization in the case of O-GlcNAcylation. The structural disruptions were found to be more pronounced for PTM at the threonine site when compared to the serine site. The salt-bridge-dependent stability of the α-helix was found to be highly position specific, an i → i + 4 interaction stabilizing the helix, with other placements leading to the destabilization of the helix.
Collapse
Affiliation(s)
- Lata Rani
- Department of Chemistry, Indian Institute of Technology Gandhinagar, Palaj, Gandhinagar 382355, Gujarat, India
| | - Sairam S Mallajosyula
- Department of Chemistry, Indian Institute of Technology Gandhinagar, Palaj, Gandhinagar 382355, Gujarat, India
| |
Collapse
|
12
|
Peng W, Wang T, Liang XR, Yang YS, Wang QZ, Cheng HF, Peng YK, Ding F. Characterizing the potentially neuronal acetylcholinesterase reactivity toward chiral pyraclofos: Enantioselective insights from spectroscopy, in silico docking, molecular dynamics simulation and per-residue energy decomposition studies. J Mol Graph Model 2021; 110:108069. [PMID: 34773872 DOI: 10.1016/j.jmgm.2021.108069] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2021] [Revised: 10/27/2021] [Accepted: 10/28/2021] [Indexed: 10/20/2022]
Abstract
Chiral organophosphorus agents are distributed ubiquitously in the environment, but the neuroactivity of these asymmetric chemicals to humans remains uncertain. This scenario was to explore the stereoselective neurobiological response of human acetylcholinesterase (AChE) to chiral pyraclofos at the enantiomeric scale, and then decipher the microscopic basis of enantioselective neurotoxicity of pyraclofos enantiomers. The results indicated that (R)-/(S)-pyraclofos can form the bioconjugates with AChE with a stoichiometric ratio of 1:1, but the neuronal affinity of (R)-pyraclofos (K = 6.31 × 104 M-1) with AChE was larger than that of (S)-pyraclofos (K = 1.86 × 104 M-1), and significant enantioselectivity was existed in the biochemical reaction. The modes of neurobiological action revealed that pyraclofos enantiomers were situated at the substrate binding domain, and the strength of the overall noncovalent bonds between (S)-pyraclofos and the residues was weaker than that of (R)-pyraclofos, resulting in the high inhibitory effect of (R)-pyraclofos toward the activity of AChE. Dynamic enantioselective biointeractions illustrated that the intervention of inherent conformational flexibility in the AChE-(R)-pyraclofos was greater than that of the AChE-(S)-pyraclofos, which arises from the big spatial displacement and the conformational flip of the binding domain composed of the residues Thr-64~Asn-89, Gly-122~Asp-134, and Thr-436~Tyr-449. Energy decomposition exhibited that the Gibbs free energies of the AChE-(R)-/(S)-pyraclofos were ΔG° = -37.4/-30.2 kJ mol-1, respectively, and the disparity comes from the electrostatic energy during the stereoselective neurochemical reactions. Quantitative conformational analysis further confirmed the atomic-scale computational chemistry conclusions, and the perturbation of (S)-pyraclofos on the AChE's ordered conformation was lower than that of (R)-pyraclofos, which is germane to the interaction energies of the crucial residues, e.g. Tyr-124, Tyr-337, Asp-74, Trp-86, and Tyr-119. Evidently, this attempt will contribute mechanistic information to uncovering the neurobiological effects of chiral organophosphates on the body.
Collapse
Affiliation(s)
- Wei Peng
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China; State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China
| | - Tao Wang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China
| | - Xiang-Rong Liang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China
| | - Yu-Sen Yang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China
| | - Qi-Zhao Wang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China
| | - Hong-Fei Cheng
- School of Earth Science and Resources, Chang'an University, Xi'an, 710054, China
| | - Yu-Kui Peng
- Xining Center for Agricultural Product Quality and Safety Testing, Xining, 810016, China
| | - Fei Ding
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China; Department of Agricultural Chemistry, Qingdao Agricultural University, Qingdao, 266109, China.
| |
Collapse
|
13
|
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
|
14
|
Perdikari TM, Jovic N, Dignon GL, Kim YC, Fawzi NL, Mittal J. A predictive coarse-grained model for position-specific effects of post-translational modifications. Biophys J 2021; 120:1187-1197. [PMID: 33582133 DOI: 10.1016/j.bpj.2021.01.034] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Revised: 12/29/2020] [Accepted: 01/19/2021] [Indexed: 02/06/2023] Open
Abstract
Biomolecules undergo liquid-liquid phase separation (LLPS), resulting in the formation of multicomponent protein-RNA membraneless organelles in cells. However, the physiological and pathological role of post-translational modifications (PTMs) on the biophysics of phase behavior is only beginning to be probed. To study the effect of PTMs on LLPS in silico, we extend our transferable coarse-grained model of intrinsically disordered proteins to include phosphorylated and acetylated amino acids. Using the parameters for modified amino acids available for fixed-charge atomistic force fields, we parameterize the size and atomistic hydropathy of the coarse-grained-modified amino acid beads and, hence, the interactions between the modified and natural amino acids. We then elucidate how the number and position of phosphorylated and acetylated residues alter the protein's single-chain compactness and its propensity to phase separate. We show that both the number and the position of phosphorylated threonines/serines or acetylated lysines can serve as a molecular on/off switch for phase separation in the well-studied disordered regions of Fused in Sarcoma (FUS) and DDX3X, respectively. We also compare modified residues to their commonly used PTM mimics for their impact on chain properties. Importantly, we show that the model can predict and capture experimentally measured differences in the phase behavior for position-specific modifications, showing that the position of modifications can dictate phase separation. In sum, this model will be useful for studying LLPS of post-translationally modified intrinsically disordered proteins and predicting how modifications control phase behavior with position-specific resolution.
Collapse
Affiliation(s)
| | - Nina Jovic
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania
| | - Gregory L Dignon
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania; Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York
| | - Young C Kim
- Center for Materials Physics and Technology, U.S. Naval Research Laboratory, Washington D.C
| | - Nicolas L Fawzi
- Department of Molecular Pharmacology, Physiology, and Biotechnology; Robert J. and Nancy D. Carney Institute for Brain Science, Brown University, Providence, Rhode Island.
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania.
| |
Collapse
|
15
|
Oliveira NFB, Pires IDS, Machuqueiro M. Improved GROMOS 54A7 Charge Sets for Phosphorylated Tyr, Ser, and Thr to Deal with pH-Dependent Binding Phenomena. J Chem Theory Comput 2020; 16:6368-6376. [PMID: 32809819 DOI: 10.1021/acs.jctc.0c00529] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Phosphorylation is a ubiquitous post-translational modification in proteins, and the phosphate group is present constitutively or transiently in most biological building blocks. These phosphorylated biomolecules are involved in many high-affinity binding/unbinding events that rely predominantly on electrostatic interactions. To build accurate models of these molecules, we need an improved description of the atomic partial charges for all relevant protonation states. In this work, we showed that the commonly used protocols to derive atomic partial charges using well-solvated molecules are inadequate to model the protonation equilibria in binding events. We introduced a protocol based on PB/MC calculations with a single representative conformation (of both protonation states) and used the resulting pKa estimations to help manually curate the atomic partial charges. The final charge set, which is fully compatible with the GROMOS 54A7 force field, proved to be very effective in modeling the protonation equilibrium in different phosphorylated peptides in the free (tetrapeptides, pentapeptides, and pY1021) and protein-complexed forms (pY1021/PLC-γ1 complex). This was particularly important in the case of the pY1021 bound to the SH2 domain of PLC-γ1, where only our curated charge set captured the correct protonation equilibrium at the neutral to slightly acidic pH range. The binding/unbinding phenomena in that pH range are biologically relevant, and to improve our models, we need to go beyond the commonly used protocols and obtain revised force field parameters for these molecules.
Collapse
Affiliation(s)
- Nuno F B Oliveira
- BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| | - Inês D S Pires
- BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| | - Miguel Machuqueiro
- BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| |
Collapse
|
16
|
Rieloff E, Skepö M. Phosphorylation of a Disordered Peptide-Structural Effects and Force Field Inconsistencies. J Chem Theory Comput 2020; 16:1924-1935. [PMID: 32050065 DOI: 10.1021/acs.jctc.9b01190] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Phosphorylation is one of the most abundant types of post-translational modifications of intrinsically disordered proteins (IDPs). This study examines the conformational changes in the 15-residue-long N-terminal fragment of the IDP statherin upon phosphorylation, using computer simulations with two different force fields: AMBER ff99SB-ILDN and CHARMM36m. The results from the simulations are compared with experimental small-angle X-ray scattering (SAXS) and circular dichroism data. In the unphosphorylated state, the two force fields are in excellent agreement regarding global structural properties such as size and shape. However, they exhibit some differences in the extent and type of the secondary structure. In the phosphorylated state, neither of the force fields performs well compared to the experimental data. Both force fields show a compaction of the peptide upon phosphorylation, greater than what is seen in SAXS experiments, although they differ in the local structure. While the CHARMM force field increases the fraction of bends in the peptide as a response to strong interactions between the phosphorylated residues and arginines, the AMBER force field shows an increase of the helical content in the N-terminal part of the peptide, where the phosphorylated residues reside, in better agreement with circular dichroism results.
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
|
17
|
Sieradzan AK, Bogunia M, Mech P, Ganzynkowicz R, Giełdoń A, Liwo A, Makowski M. Introduction of Phosphorylated Residues into the UNRES Coarse-Grained Model: Toward Modeling of Signaling Processes. J Phys Chem B 2019; 123:5721-5729. [PMID: 31194908 DOI: 10.1021/acs.jpcb.9b03799] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Phosphorylated proteins take part in many signaling pathways and play a key role in homeostasis regulation. The all-atom force fields enable us to study the systems containing phosphorylated proteins, but they are limited to short time scales. In this paper, we report the extension of the physics-based coarse-grained UNRES force field to treat systems with phosphorylated amino-acid residues. To derive the respective potentials, appropriate physics-based analytical expressions were fitted to the potentials of mean force of systems modeling phosphorylated amino-acid residues computed in our previous work and implemented in UNRES. The extended UNRES performed well in ab initio simulations of two miniproteins containing phosphorylated residues, strongly suggesting that realistic large-scale simulations of processes involving phosphorylated proteins, especially signaling processes, are now possible.
Collapse
Affiliation(s)
- Adam K Sieradzan
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Małgorzata Bogunia
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Paulina Mech
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Robert Ganzynkowicz
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Artur Giełdoń
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Adam Liwo
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Mariusz Makowski
- Faculty of Chemistry , University of Gdańsk , ul. Wita Stwosza 63 , 80-308 Gdańsk , Poland
| |
Collapse
|