Sticky Salts: Overbinding of Monovalent Cations to Phosphorylations in All-Atom Force Fields

Phosphorylation is a major post-translational modification that is involved in the regulation of the dynamics and function of intrinsically disordered proteins (IDPs). We recently characterized a phenomenon, which we termed <i>n</i>-phosphate collaborations (<i>n</i>P-collabs), where bulk cations form stable bridges between several phosphoresidues in all-atom molecular dynamic simulations. <i>n</i>P-collabs were found to be sensitive to the combination of force fields and cation types. Here, we attempt to assess the physical relevance of these <i>n</i>P-collabs by evaluating the strength of the cation/phosphate interaction through osmotic coefficient (ϕ) calculations on the model 2Na<sup>+</sup>HPO<sub>4</sub><sup>2-</sup> and 2K<sup>+</sup>HPO<sub>4</sub><sup>2-</sup> salts, using different classical force fields for phosphorylations. All force fields were found to overestimate the strength of the interaction to various degrees. We thus designed new parameters for CHARMM36m and AmberFF99SB-ILDN using the Electronic Continuum Correction (ECC) approach, which provide remarkable agreement for ϕ values for both cation types and over a range of concentrations. We provide a preliminary test of these ECC parameters for phosphorylations by simulating the 7-fold-phosphorylated rhodopsin peptide 7PP and comparing secondary chemical shifts to experimental data. Conformational ensembles resulting from the ECC-derived phosphorylated force fields display both qualitative and quantitative improvements with regard to full-charge force fields. We thus conclude that long-lasting <i>n</i>P-collabs are artifacts for classical force fields born from the lack of explicit polarization, and propose a possible computational strategy for the extensive parametrization of phosphorylations. The presence of long-lived <i>n</i>P-collabs in simulations produced using classical force fields is therefore a serious concern for the accurate modeling of multiphosphorylated peptides and IDPs, which are at the center of research questions regarding neurodegenerative diseases such as Alzheimer's or Parkinson's.

Authors

Institutions

Publication Details

Journal
Journal of Chemical Information and Modeling
Published
2025-12-05
DOI
https://doi.org/10.1021/acs.jcim.5c02154
Primary Topic
Photoreceptor and optogenetics research
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Sticky Salts: Overbinding of Monovalent Cations to Phosphorylations in All-Atom Force Fields

Chantal Prévost, Sophie Sacquin‐Mora, Jules Marien, Julie Puyo-Fourtine et al.
Journal of Chemical Information and Modeling
Photoreceptor and optogenetics research
article

Sticky Salts: Overbinding of Monovalent Cations to Phosphorylations in All-Atom Force Fields

Chantal Prévost, Sophie Sacquin‐Mora, Jules Marien, Julie Puyo-Fourtine, Élise Duboué-Dijon
article en

Abstract

Phosphorylation is a major post-translational modification that is involved in the regulation of the dynamics and function of intrinsically disordered proteins (IDPs). We recently characterized a phenomenon, which we termed <i>n</i>-phosphate collaborations (<i>n</i>P-collabs), where bulk cations form stable bridges between several phosphoresidues in all-atom molecular dynamic simulations. <i>n</i>P-collabs were found to be sensitive to the combination of force fields and cation types. Here, we attempt to assess the physical relevance of these <i>n</i>P-collabs by evaluating the strength of the cation/phosphate interaction through osmotic coefficient (ϕ) calculations on the model 2Na<sup>+</sup>HPO<sub>4</sub><sup>2-</sup> and 2K<sup>+</sup>HPO<sub>4</sub><sup>2-</sup> salts, using different classical force fields for phosphorylations. All force fields were found to overestimate the strength of the interaction to various degrees. We thus designed new parameters for CHARMM36m and AmberFF99SB-ILDN using the Electronic Continuum Correction (ECC) approach, which provide remarkable agreement for ϕ values for both cation types and over a range of concentrations. We provide a preliminary test of these ECC parameters for phosphorylations by simulating the 7-fold-phosphorylated rhodopsin peptide 7PP and comparing secondary chemical shifts to experimental data. Conformational ensembles resulting from the ECC-derived phosphorylated force fields display both qualitative and quantitative improvements with regard to full-charge force fields. We thus conclude that long-lasting <i>n</i>P-collabs are artifacts for classical force fields born from the lack of explicit polarization, and propose a possible computational strategy for the extensive parametrization of phosphorylations. The presence of long-lived <i>n</i>P-collabs in simulations produced using classical force fields is therefore a serious concern for the accurate modeling of multiphosphorylated peptides and IDPs, which are at the center of research questions regarding neurodegenerative diseases such as Alzheimer's or Parkinson's.

Journal of Chemical Information and ModelingVol. 65(24)
Centre National de la Recherche Scientifique (FR), Université Paris Cité (FR), Laboratoire de Biochimie Théorique (FR)
Openalex Percentile: Top 58%
Photoreceptor and optogenetics research
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.