Predicting Saturation Concentrations of Phase-Separating Proteins via Thermodynamic Integration

Phase separation of proteins and nucleic acids into biomolecular condensates contributes to the regulation of cellular compartmentalisation in membrane-less environments. A key parameter controlling the onset of biomolecular condensate formation via liquid—liquid phase separation is the saturation concentration ( C _sat )— the threshold concentration above which condensation takes place. While measuring C _sat for protein solutions in vitro is experimentally accessible, determining this quantity in simulations remains challenging due to the extremely low equilibrium concentrations at which many proteins phase separate. This occurs because the gold standard in simulations consists on combining a residue-resolution coarse-grained model with the Direct Coexistence simulation method, which yields poor estimates of the equilibrium concentrations of the dilute phase due to lack of statistics. In this work, we present two independent thermodynamic integration (TI) schemes which, when combined with Direct Coexistence simulations, enable accurate calculation of saturation concentrations and phase diagrams—facilitating direct comparison with experimental measurements across a wide range of conditions. Our methods, combined with the Mpipi-Recharged residue-resolution coarse-grained model, accurately estimate C _sat for a wide range of intrinsically disordered and multi-domain proteins, including disease-associated RNA- and DNA-binding proteins involved in the formation of stress granules and P granules, as well as engineered mutants of hnRNPA1. Furthermore, we compare our TI methods against a computationally efficient machine-learning predictor trained to estimate saturation concentrations at sub-physiological temperatures. While both approaches yield realistic predictions, explicit molecular dynamics simulations enable the calculation of complete phase diagrams and provide insight into the molecular mechanisms and interactions driving phase-separation. Overall, our approach offers a robust, physically grounded framework for improving and validating coarse-grained models of biomolecular phase behaviour, effectively bridging the gap between simulation and experiment.