Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent ({approx}1 µs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex

Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain

Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg2+-ions. Furthermore, a role for Mg2+-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA–adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg2+, but becomes partially pre-organized for ligand binding in the presence of Mg2+. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg2+ is more pronounced

Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain

Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg2+-ions. Furthermore, a role for Mg2+-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA–adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg2+, but becomes partially pre-organized for ligand binding in the presence of Mg2+. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg2+ is more pronounced

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (≈1 μs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex

Recognition of two distinct elements in the RNA substrate by the RNA-binding domain of the T. thermophilus DEAD box helicase Hera

DEAD box helicases catalyze the ATP-dependent destabilization of RNA duplexes. Whereas duplex separation is mediated by the helicase core shared by all members of the family, flanking domains often contribute to binding of the RNA substrate. The Thermus thermophilus DEAD-box helicase Hera (for “heat-resistant RNA-binding ATPase”) contains a C-terminal RNA-binding domain (RBD). We have analyzed RNA binding to the Hera RBD by a combination of mutational analyses, nuclear magnetic resonance and X-ray crystallography, and identify residues on helix α1 and the C-terminus as the main determinants for high-affinity RNA binding. A crystal structure of the RBD in complex with a single-stranded RNA resolves the RNA–protein interactions in the RBD core region around helix α1. Differences in RNA binding to the Hera RBD and to the structurally similar RBD of the Bacillus subtilis DEAD box helicase YxiN illustrate the versatility of RNA recognition motifs as RNA-binding platforms. Comparison of chemical shift perturbation patterns elicited by different RNAs, and the effect of sequence changes in the RNA on binding and unwinding show that the RBD binds a single-stranded RNA region at the core and simultaneously contacts double-stranded RNA through its C-terminal tail. The helicase core then unwinds an adjacent RNA duplex. Overall, the mode of RNA binding by Hera is consistent with a possible function as a general RNA chaperone

Interplay of ‘induced fit’ and preorganization in the ligand induced folding of the aptamer domain of the guanine binding riboswitch

Riboswitches are highly structured elements in the 5′-untranslated regions (5′-UTRs) of messenger RNA that control gene expression by specifically binding to small metabolite molecules. They consist of an aptamer domain responsible for ligand binding and an expression platform. Ligand binding in the aptamer domain leads to conformational changes in the expression platform that result in transcription termination or abolish ribosome binding. The guanine riboswitch binds with high-specificity to guanine and hypoxanthine and is among the smallest riboswitches described so far. The X-ray-structure of its aptamer domain in complex with guanine/hypoxanthine reveals an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base pairing interactions. We analyzed the conformational transitions of the aptamer domain induced by binding of hypoxanthine using high-resolution NMR-spectroscopy in solution. We found that the long-range base pairing interactions are already present in the free RNA and preorganize its global fold. The ligand binding core region is lacking hydrogen bonding interactions and therefore likely to be unstructured in the absence of ligand. Mg(2+)-ions are not essential for ligand binding and do not change the structure of the RNA-ligand complex but stabilize the structure at elevated temperatures. We identified a mutant RNA where the long-range base pairing interactions are disrupted in the free form of the RNA but form upon ligand binding in an Mg(2+)-dependent fashion. The tertiary interaction motif is stable outside the riboswitch context

Recognition of two distinct elements in the RNA substrate by the RNA-binding domain of the T. thermophilus DEAD box helicase Hera

DEAD box helicases catalyze the ATP-dependent destabilization of RNA duplexes. Whereas duplex separation is mediated by the helicase core shared by all members of the family, flanking domains often contribute to binding of the RNA substrate. The Thermus thermophilus DEAD-box helicase Hera (for "heat-resistant RNA-binding ATPase”) contains a C-terminal RNA-binding domain (RBD). We have analyzed RNA binding to the Hera RBD by a combination of mutational analyses, nuclear magnetic resonance and X-ray crystallography, and identify residues on helix α1 and the C-terminus as the main determinants for high-affinity RNA binding. A crystal structure of the RBD in complex with a single-stranded RNA resolves the RNA-protein interactions in the RBD core region around helix α1. Differences in RNA binding to the Hera RBD and to the structurally similar RBD of the Bacillus subtilis DEAD box helicase YxiN illustrate the versatility of RNA recognition motifs as RNA-binding platforms. Comparison of chemical shift perturbation patterns elicited by different RNAs, and the effect of sequence changes in the RNA on binding and unwinding show that the RBD binds a single-stranded RNA region at the core and simultaneously contacts double-stranded RNA through its C-terminal tail. The helicase core then unwinds an adjacent RNA duplex. Overall, the mode of RNA binding by Hera is consistent with a possible function as a general RNA chaperon

The crystal structure of Nep1 reveals an extended SPOUT-class methyltransferase fold and a pre-organized SAM-binding site

Ribosome biogenesis in eukaryotes requires the participation of a large number of ribosome assembly factors. The highly conserved eukaryotic nucleolar protein Nep1 has an essential but unknown function in 18S rRNA processing and ribosome biogenesis. In Saccharomyces cerevisiae the malfunction of a temperature-sensitive Nep1 protein (nep1-1ts) was suppressed by the addition of S-adenosylmethionine (SAM). This suggests the participation of Nep1 in a methyltransferase reaction during ribosome biogenesis. In addition, yeast Nep1 binds to a 6-nt RNA-binding motif also found in 18S rRNA and facilitates the incorporation of ribosomal protein Rps19 during the formation of pre-ribosomes. Here, we present the X-ray structure of the Nep1 homolog from the archaebacterium Methanocaldococcus jannaschii in its free form (2.2 Å resolution) and bound to the S-adenosylmethionine analog S-adenosylhomocysteine (SAH, 2.15 Å resolution) and the antibiotic and general methyltransferase inhibitor sinefungin (2.25 Å resolution). The structure reveals a fold which is very similar to the conserved core fold of the SPOUT-class methyltransferases but contains a novel extension of this common core fold. SAH and sinefungin bind to Nep1 at a preformed binding site that is topologically equivalent to the cofactor-binding site in other SPOUT-class methyltransferases. Therefore, our structures together with previous genetic data suggest that Nep1 is a genuine rRNA methyltransferase

The crystal structure of Nep1 reveals an extended SPOUT-class methyltransferase fold and a pre-organized SAM-binding site

Ribosome biogenesis in eukaryotes requires the participation of a large number of ribosome assembly factors. The highly conserved eukaryotic nucleolar protein Nep1 has an essential but unknown function in 18S rRNA processing and ribosome biogenesis. In Saccharomyces cerevisiae the malfunction of a temperature-sensitive Nep1 protein (nep1-1ts) was suppressed by the addition of S-adenosylmethionine (SAM). This suggests the participation of Nep1 in a methyltransferase reaction during ribosome biogenesis. In addition, yeast Nep1 binds to a 6-nt RNA-binding motif also found in 18S rRNA and facilitates the incorporation of ribosomal protein Rps19 during the formation of pre-ribosomes. Here, we present the X-ray structure of the Nep1 homolog from the archaebacterium Methanocaldococcus jannaschii in its free form (2.2 Å resolution) and bound to the S-adenosylmethionine analog S-adenosylhomocysteine (SAH, 2.15 Å resolution) and the antibiotic and general methyltransferase inhibitor sinefungin (2.25 Å resolution). The structure reveals a fold which is very similar to the conserved core fold of the SPOUT-class methyltransferases but contains a novel extension of this common core fold. SAH and sinefungin bind to Nep1 at a preformed binding site that is topologically equivalent to the cofactor-binding site in other SPOUT-class methyltransferases. Therefore, our structures together with previous genetic data suggest that Nep1 is a genuine rRNA methyltransferase

The ribosome assembly factor Nep1 responsible for Bowen–Conradi syndrome is a pseudouridine-N1-specific methyltransferase

Nep1 (Emg1) is a highly conserved nucleolar protein with an essential function in ribosome biogenesis. A mutation in the human Nep1 homolog causes Bowen–Conradi syndrome—a severe developmental disorder. Structures of Nep1 revealed a dimer with a fold similar to the SPOUT-class of RNA-methyltransferases suggesting that Nep1 acts as a methyltransferase in ribosome biogenesis. The target for this putative methyltransferase activity has not been identified yet. We characterized the RNA-binding specificity of Methanocaldococcus jannaschii Nep1 by fluorescence- and NMR-spectroscopy as well as by yeast three-hybrid screening. Nep1 binds with high affinity to short RNA oligonucleotides corresponding to nt 910–921 of M. jannaschii 16S rRNA through a highly conserved basic surface cleft along the dimer interface. Nep1 only methylates RNAs containing a pseudouridine at a position corresponding to a previously identified hypermodified N1-methyl-N3-(3-amino-3-carboxypropyl) pseudouridine (m1acp3-Ψ) in eukaryotic 18S rRNAs. Analysis of the methylated nucleoside by MALDI-mass spectrometry, HPLC and NMR shows that the methyl group is transferred to the N1 of the pseudouridine. Thus, Nep1 is the first identified example of an N1-specific pseudouridine methyltransferase. This enzymatic activity is also conserved in human Nep1 suggesting that Nep1 is the methyltransferase in the biosynthesis of m1acp3-Ψ in eukaryotic 18S rRNAs