Combining CHARMM36m with OPC water improves accuracy
- 1Max Planck Institute for Multidisciplinary Sciences
Abstract
Atomistic simulations of intrinsically disordered proteins (IDPs) are notoriously sensitive to force field inaccuracies, either regarding the protein or the water model, yielding inaccurate observables such as compactness, secondary structure propensities, or kinetics. The currently most widely used IDP force fields are Amber99sb-disp (A99disp) and Charmm36m (C36m). A99disp includes a new water model and thus optimized both, the protein and the water interactions. In contrast, C36m used the Tip3p water model and optimized only protein interactions. In many cases, C36m+Tip3p underestimates radii of gyration compared to FRET or SAXS experiments. Such overly compact structural ensembles are believed to arise from an imbalance between protein-protein, protein-water, and water-water interactions, which might be due to Tip3p inaccuracies. Here, we aim at re-balancing these interactions by combining C36m with the Optimal Point Charge (OPC) water model. Recently, this C36m+OPC combination showed improved accuracy for the disordered domain of the measles virus nucleoprotein. Here we present a systematic assessment, comparing C36m+OPC, C36m+Tip3p, C36m+Tip4p, C22*, A03ws, A99sb-ws, and A99disp for five IDPs, as well as a subset of those for five globular proteins, a set of disordered AGQ-repeat peptides, and the fast folding miniprotein CLN025. We compared extensive MD simulations (> 8.5 ms) with SAXS, NMR, circular dichroism, photo-induced electron transfer (PET), T-jump infrared spectroscopy, and X-ray crystallography measurements. We found that combining C36m with OPC improved accuracy over C36m+Tip3p for IDP ensembles without compromising its accuracy for globular proteins. While also the kinetics of the AGQ-peptides were more accurate for C36m+OPC, those of CLN025 folding were less accurate. Overall, C36m+OPC showed similar accuracy as A99sb-ws and A99disp, the latter is currently considered among the most accurate protein force fields.
Related
- Large language model consensus substantially improves the cell type annotation accuracy for scRNA-seq data
- Integration of exosomal SNORD46 and AL592064.2 with conventional tumor markers improves diagnostic accuracy for colorectal cancer
- Toward mechanical markers of muscle state exercise: combining ex vivo and in vivo approaches
- Hindlimb Muscle Morphology and Function in a New Atrophy Model Combining Spinal Cord Injury and Cast Immobilization
- Combining AlphaFold 3 with microscale thermophoresis for aptamer optimisation
BibTeX
@misc{kozlowski2026combining,
title = {{Combining CHARMM36m with OPC water improves accuracy}},
author = {N. Kozlowski and G. Nagy and C. F. Ottoway and L. V. Bock and S. M. Vaiana and H. Grubmuller},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.25.754057},
url = {https://www.biorxiv.org/content/10.64898/2026.09.25.754057v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Combining CHARMM36m with OPC water improves accuracy
AU - N. Kozlowski
AU - G. Nagy
AU - C. F. Ottoway
AU - L. V. Bock
AU - S. M. Vaiana
AU - H. Grubmuller
PY - 2026
DA - 2026/09/25
PB - bioRxiv
DO - 10.64898/2026.09.25.754057
UR - https://www.biorxiv.org/content/10.64898/2026.09.25.754057v1
AB - Atomistic simulations of intrinsically disordered proteins (IDPs) are notoriously sensitive to force field inaccuracies, either regarding the protein or the water model, yielding inaccurate observables such as compactness, secondary structure propensities, or kinetics. The currently most widely used IDP force fields are Amber99sb-disp (A99disp) and Charmm36m (C36m). A99disp includes a new water model and thus optimized both, the protein and the water interactions. In contrast, C36m used the Tip3p water model and optimized only protein interactions. In many cases, C36m+Tip3p underestimates radii of gyration compared to FRET or SAXS experiments. Such overly compact structural ensembles are believed to arise from an imbalance between protein-protein, protein-water, and water-water interactions, which might be due to Tip3p inaccuracies. Here, we aim at re-balancing these interactions by combining C36m with the Optimal Point Charge (OPC) water model. Recently, this C36m+OPC combination showed improved accuracy for the disordered domain of the measles virus nucleoprotein. Here we present a systematic assessment, comparing C36m+OPC, C36m+Tip3p, C36m+Tip4p, C22*, A03ws, A99sb-ws, and A99disp for five IDPs, as well as a subset of those for five globular proteins, a set of disordered AGQ-repeat peptides, and the fast folding miniprotein CLN025. We compared extensive MD simulations (> 8.5 ms) with SAXS, NMR, circular dichroism, photo-induced electron transfer (PET), T-jump infrared spectroscopy, and X-ray crystallography measurements. We found that combining C36m with OPC improved accuracy over C36m+Tip3p for IDP ensembles without compromising its accuracy for globular proteins. While also the kinetics of the AGQ-peptides were more accurate for C36m+OPC, those of CLN025 folding were less accurate. Overall, C36m+OPC showed similar accuracy as A99sb-ws and A99disp, the latter is currently considered among the most accurate protein force fields.
ER -