Nature CommunicationsBiophysicsPeer-reviewed journal article
Molecular simulations of enzymatic phosphorylation of disordered proteins and their condensates
Emanuele ZippoiD, Dorothee DormanniD, Thomas SpeckiD, Lukas S. StelzliD
Abstract
Understanding the condensation and aggregation of intrinsically disordered proteins in a non-equilibrium environment is crucial for unraveling many biological processes. Active enzymes catalyse many processes by consuming chemical fuels such as ATP. Enzymes called kinases phosphorylate disordered regions of proteins and thus profoundly affect their properties and interactions. Protein phosphorylation is implicated in neurodegenerative diseases and may modulate pathogenesis. However, how protein sequence and molecular recognition of a disordered protein by kinases determine phosphorylation patterns is not understood. In principle, molecular dynamics simulations hold the promise to resolve how phosphorylation affects disordered proteins and their assemblies. In practice, chemically-detailed simulations of enzymatic reactions and the dynamics of enzymes are highly challenging, in particular it is difficult to verify whether implementations of driven simulations are thermodynamically consistent. We can now address this problem with residue-level coarse-grained molecular dynamics simulations, integrating Metropolis Monte Carlo steps to model chemical reactions. Importantly, we show how to verify by Markov-state modeling that the realisation of a non-equilibrium steady state satisfies local-detailed balance. We investigate TDP-43 phosphorylation by the kinase CK1{delta} in simulations, examining patterns of phosphorylation and assessing its preventive role in chain aggregation, which may be a cytoprotective mechanism in neurodegenerative diseases. We find that the degree of residue phosphorylation is determined by sequence preference and charges, rather than by the position in the chain. The phosphorylation frequency is also affected by the phosphorylation patterns, since the interactions between CK1{delta} and TDP-43 actively change after each reaction. For TDP-43, our simulations show condensates dissolution through phosphorylation with kinases binding to the condensates and phosphorylating TDP-43 in the condensates.
- Citations
- 28
- FWCI
- 6.73
- Field percentile
- top 2%
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Kinase-phosphatase cycles drive self-sustained oscillations of protein condensates
Emanuele Zippo, Le Qiao, … Lukas S. Stelzl·bioRxiv·Oct 30, 2026
- A coarse-grained model for simulations of phosphorylated disordered proteins
Arriën Symon Rauh, Gustav Stausbøll Hedemark, … Kresten Lindorff‐Larsen·Biophysical JournalPeer-reviewed·Jul 7, 2025
- Phosphorylation of disordered proteins tunes local and global intramolecular interactions
Emery T. Usher, Martin J. Fossat, Alex S. Holehouse·Biophysical JournalPeer-reviewed·Nov 13, 2024
- Effects of phosphorylation on protein backbone dynamics and conformational preferences
David Bickel, Wim Vranken·Journal of Chemical Theory and ComputationPeer-reviewed·Jun 3, 2024
- Slow diffusion limits phosphorylation in a biomolecular condensate
Nicolás S. González Foutel, Ankush Garg, … Magnus Kjærgaard·bioRxiv·Feb 2, 2026
- Mechanism of MEK1 activation by phosphorylation
Everett L. Jin, Levent Sari, Milo M. Lin·Biophysical JournalPeer-reviewed·Feb 26, 2026