The hyperthermophilic Ssh10b from Sulfolobus shibatae is a member of the Sac10b family, which binds RNA in vivo as a physiological substrate, and it has been postulated to play a key role in chromosomal organization in Archaea. Even though the crystal structure of Ssh10b-RNA was resolved successively by X-ray diffraction (Protein Data Bank [PDB] code: 3WBM), the detailed dynamic characteristics of Ssh10b-RNA are still unclear. In this study, molecular dynamics (MDs) simulations at 6 temperatures (300, 350, 375, 400, 450, and 500 K) and molecular mechanics Generalized-Born surface area (MM-GB/SA) free energy calculations were performed to investigate the mechanism of how Ssh10b protects and stabilizes RNA. The simulation results indicate that RNA is stabilized by Ssh10b when the temperature rises up to 375 K. RNA is found to undergo conformational transition between A-RNA and A 0 -RNA when Ssh10b binds to RNA at 3 different temperatures (300, 350, and 375 K). Salt bridges, hydrogen bonds and hydrophobic interactions are observed, and some residues have significant impact on the structural stability of the complex. This study increases our understanding of the dynamics and interaction mechanism of hyperthermophilic proteins and RNA at the atomic level, and offers a model for studying the structural biology of hyperthermophilic proteins and RNA.
K E Y W O R D Shyperthermophilic proteins, molecular dynamics, MM-GB/SA, RNA, Ssh10b