Related Experiment Video
Updated: Jan 9, 2026

Identification of Cyclin-dependent Kinase 1 Specific Phosphorylation Sites by an In Vitro Kinase Assay
Published on: May 3, 2018
Sticky Salts: Overbinding of Monovalent Cations to Phosphorylations in All-Atom Force Fields
Jules Marien1, Julie Puyo-Fourtine1, Chantal Prévost1
1Université Paris Cité, CNRS, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005, Paris, France.
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 n-phosphate collaborations (nP-collabs), where bulk cations form stable bridges between several phosphoresidues in all-atom molecular dynamic simulations. nP-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 nP-collabs by evaluating the strength of the cation/phosphate interaction through osmotic coefficient (ϕ) calculations on the model 2Na+HPO42- and 2K+HPO42- 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 nP-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 nP-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.
More Related Videos
05:37Rapid in-silico Battery Electrolyte Electrochemical Reaction Generation using 3T-VASP Multi-Scale Energy Minimization
Published on: August 22, 2025
08:54Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid
Published on: January 25, 2020
Related Concept Videos
Phosphorylation
During phosphorylation, protein kinases transfer the terminal phosphate group of ATP to specific amino acid side chains of substrate proteins. Serine, threonine, and tyrosine are the most commonly...
Crystal Field Theory - Octahedral Complexes
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
Protein Kinases and Phosphatases
Protein kinases
Many proteins in the cell are regulated by phosphorylation, the addition of a phosphate group. A family of enzymes called kinases...
Protein Kinases and Phosphatases
Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than the dxy,...
Complexation Equilibria: Overview
The equilibrium constant of the complexation reaction is represented as the formation constant...