Atomic-level insights into metabolite recognition and specificity of the SAM-II riboswitch

Although S-adenosylhomocysteine (SAH), a metabolic by-product of S-adenosylmethionine (SAM), differs from SAM only by a single methyl group and an overall positive charge, SAH binds the SAM-II riboswitch with more than 1000-fold less affinity than SAM. Using atomistic molecular dynamics simulations, we investigated the molecular basis of such high selectivity in ligand recognition by SAM-II riboswitch. The biosynthesis of SAM exclusively generates the (S,S) stereoisomer, and (S,S)-SAM can spontaneously convert to the (R,S) form. We, therefore, also examined the effects of (R,S)-SAM binding to SAM-II and its potential biological function. We find that the unfavorable loss in entropy in SAM-II binding is greater for (S,S)- and (R,S)-SAM than SAH, which is compensated by stabilizing electrostatic interactions with the riboswitch. The positively charged sulfonium moiety on SAM acts as the crucial anchor point responsible for the formation of key ionic interactions as it fits favorably in the negatively charged binding pocket. In contrast, SAH, with its lone pair of electrons on the sulfur, experiences repulsion in the binding pocket of SAM-II and is enthalpically destabilized. In the presence of SAH, similar to the unbound riboswitch, the pseudoknot structure of SAM-II is not completely formed, thus exposing the Shine-Dalgarno sequence. Unlike SAM, this may further facilitate ribosomal assembly and translation initiation. Our analysis of the conformational ensemble sampled by SAM-II in the absence of ligands and when bound to SAM or SAH reveals that ligand binding follows a combination of conformational selection and induced-fit mechanisms.