🤖 AI Summary
Existing machine learning coarse-grained models struggle to accurately predict temperature-dependent thermodynamic properties of proteins, such as heat capacity, across different temperatures. This work proposes a thermodynamics-informed, temperature-transferable modeling framework that, for the first time, explicitly decomposes energy and entropy contributions within the coarse-grained potential and embeds thermodynamic constraints to ensure physical consistency. The approach enables post-hoc temperature correction without retraining, facilitating both interpolation and extrapolation across temperatures. Validated on the Chignolin protein using 250 microseconds of molecular dynamics data spanning 300–400 K, the method successfully reproduces the temperature dependence of atomistic free energy landscapes and accurately predicts heat capacity.
📝 Abstract
Coarse-grained (CG) molecular simulations offer an efficient alternative to atomistic molecular dynamics to study large and complex biological systems. The accuracy of CG simulations has been increased dramatically by the introduction of machine-learned coarse-grained (MLCG) models. However, these models are typically designed to be used at a single thermodynamic point, lack temperature transferability, and can not be used to predict temperature dependent quantities like the heat capacity. Here we introduce a thermodynamically informed, temperature-transferable MLCG framework for proteins that explicitly decomposes the CG potential of mean force (PMF) into its energetic and entropic components. The model architecture enforces an exact thermodynamic relation between the energetic and entropic components of the PMF and guarantees physically consistent extrapolation and interpolation across temperature regimes. We validate this framework on an extensive dataset spanning a total of 250 $μ$s of molecular dynamics simulations across five temperatures between 300 K and 400 K for the Chignolin protein, and demonstrate that it reproduces the temperature dependency of the reference atomistic free energy surfaces, correcting the temperature-unaware baselines. Furthermore, we show that it is possible to apply an inexpensive, post-hoc temperature-dependent correction that does not require retraining the MLCG potential, accurately recovering the atomistic heat capacity at different temperatures. Overall, this work provides a physically grounded pathway toward thermodynamically transferable MLCG simulations of complex biomolecular systems.