Blind Prediction of Complex Water and Ion Ensembles Around RNA in CASP16

Proteins: Structure, Function, And Bioinformatics, n/a, 1–22, 2025

Authors: Kretsch Rachael C., Posani Elisa, Baulin Eugene F., Bujnicki Janusz M., Bussi Giovanni, Cheatham Iii Thomas E., Chen Shi-Jie, Elofsson Arne, Farsani Masoud Amiri, Fisher Olivia N., Gromiha M. Michael, Gupta Ayush, Hamada Michiaki, Harini K., Hu Gang, Huang David, Iwakiri Junichi, Jain Anika, Kagaya Yuki, Kihara Daisuke, Kmiecik Sebastian, Krishnan Sowmya Ramaswamy, Kurisaki Ikuo, Languin-Cattoën Olivier, Li Jun, Li Shanshan, Malekzadeh Karim, Nakamura Tsukasa, Ni Wentao, Chandran Nithin, Palo Michael Z., Park Joon Hong, Pilla Smita, Poblete Simón, Pucci Fabrizio, Punuru Pranav, Saha Anouka, Sato Kengo, Srivastava Ambuj, Terashi Genki, Tugolukova Emilia, Verburgt Jacob, Wuyun Qiqige, Zerze Gül H., Zhang Kaiming, Zhang Sicheng, Zheng Wei, Zhou Yuanzhe, Chiu Wah, Case David A., Das Rhiju

Abstract

ABSTRACT Biomolecules rely on water and ions for stable folding, but these interactions are often transient, dynamic, or disordered and thus hidden from experiments and evaluation challenges that represent biomolecules as single, ordered structures. Here, we compare blindly predicted ensembles of water and ion structure to the cryo-EM densities observed around the Tetrahymena ribozyme at 2.2–2.3 Å resolution, collected through target R1260 in the CASP16 competition. Twenty-six groups participated in this solvation “cryo-ensemble” prediction challenge, submitting over 350 million atoms in total, offering the first opportunity to compare blind predictions of dynamic solvent shell ensembles to cryo-EM density. Predicted atomic ensembles were converted to density through local alignment and these densities were compared to the cryo-EM densities using Pearson correlation, Spearman correlation, mutual information, and precision-recall curves. These predictions show that an ensemble representation is able to capture information of transient or dynamic water and ions better than traditional atomic models, but there remains a large accuracy gap to the performance ceiling set by experimental uncertainty. Overall, molecular dynamics approaches best matched the cryo-EM density, with blind predictions from bussilab\_plain\_md, SoutheRNA, bussilab\_replex, coogs2, and coogs3 outperforming the baseline molecular dynamics prediction. This study indicates that simulations of water and ions can be quantitatively evaluated with cryo-EM maps. We propose that further community-wide blind challenges can drive and evaluate progress in modeling water, ions, and other previously hidden components of biomolecular systems.