Dynamic Asymmetry and the Role of the Conserved Active-Site Thiol in Rabbit Muscle Creatine Kinase

Symmetric and asymmetric crystal structures of the apo and transition state analogue forms, respectively, of the dimeric rabbit muscle creatine kinase have invoked an induced fit explanation for asymmetry between the two subunits and their active sites. However, previously reported thiol reactivity studies at the dual active-site cysteine 283 residues suggest a more latent asymmetry between the two subunits. The role of that highly conserved active-site cysteine has also not been clearly determined. In this work, the SH vibrations of Cys283 were observed in the unmodified MM isoform enzyme via Raman scattering, and then one and both Cys283 residues in the same dimeric enzyme were modified to covalently attach a cyano group that reports on the active-site environment via its infrared CN stretching absorption band while maintaining the catalytic activity of the enzyme. Unmodified and Cys283-modified enzymes were investigated in the apo and transition state analogue forms of the enzyme. The narrow and invariant SH vibrational bands report a homogeneous environment for the unmodified active-site cysteines, indicating that their thiols are hydrogen bonded to the same H-bond acceptor in the presence and absence of the substrate. The SH peak persists at all physiologically relevant pHs, indicating that Cys283 is protonated at all pHs relevant to enzymatic activity. Molecular dynamics simulations identify the S-H hydrogen bond acceptor as a single, long-resident water molecule and suggest that the role of the conserved yet catalytically unnecessary thiol may be to dynamically rigidify that part of the active site through specific H-bonding to water. The asymmetric and broad CN stretching bands from the CN-modified Cys283 suggest an asymmetric structure in the apo form of the enzyme in which there is a dynamic exchange between spectral subpopulations associated with water-exposed and water-excluded probe environments. Molecular dynamics simulations indicate a homogeneous orientation of the SCN probe group in the active site and thus rule out a local conformational explanation at the residue level for the multipopulation CN stretching bands. The homogeneous simulated SCN orientation suggests strongly that a more global asymmetry between the two subunits is the cause of the CN probes broad and asymmetric infrared line shape. Together, these spectral observations localized at the active-site cysteines indicate an intrinsic, dynamic asymmetry between the two subunits that exists already in the apo form of the dimeric creatine kinase enzyme, rather than being induced by the substrate. Biochemical and methodological consequences of these conclusions are considered.