Abstract
Ryanodine receptor type 1 (RyR1) is a Ca(2+)-release channel central to skeletal muscle excitation-contraction (EC) coupling. RyR1's cryo-EM structures reveal a zinc-finger motif positioned within the cytoplasmic C-terminal domain (CTD). Yet, owing to limitations in cryo-EM resolution, RyR1 structures lack precision in detailing the metal coordination structure, prompting the need for an accurate model. In this study, we employed molecular dynamics (MD) simulations and the density functional theory (DFT) method to refine the binding characteristics of Zn(2+) in the zinc-finger site of the RyR1 channel. Our findings also highlight substantial conformational changes in simulations conducted in the absence of Zn(2+). Notably, we observed a loss of contact at the interface between protein domains proximal to the zinc-finger site, indicating a crucial role of Zn(2+) in maintaining structural integrity and interdomain interactions within RyR1. Furthermore, this study provides valuable insights into the modulation of ATP, Ca(2+), and caffeine binding, shedding light on the intricate relationship between Zn(2+) coordination and the dynamic behavior of RyR1. Our integrative approach combining MD simulations and DFT calculations enhances our understanding of the molecular mechanisms governing ligand binding in RyR1.