Structure and ion-dependent folding of k-junctions.

k-Junctions are elaborated forms of kink turns with an additional helix on the nonbulged strand, thus forming a three-way helical junction. Two were originally identified in the structures of Arabidopsis and Escherichia coli thiamine pyrophosphate (TPP) riboswitches, and another called DUF-3268 was tentatively identified from sequence information. In this work we show that the Arabidopsis and E. coli riboswitch k-junctions fold in response to the addition of magnesium or sodium ions, and that atomic mutations that should disrupt key hydrogen bonding interactions greatly impair folding. Using X-ray crystallography, we have determined the structure of the DUF-3268 RNA and thus confirmed that it is a k-junction. It also folds upon the addition of metal ions, though requiring a 40-fold lower concentration of either divalent or monovalent ions. The key difference between the DUF-3268 and riboswitch k-junctions is the lack of nucleotides inserted between G1b and A2b in the former. We show that this insertion is primarily responsible for the difference in folding properties. Finally, we show that the DUF-3268 can functionally substitute for the k-junction in the E. coli TPP riboswitch such that the chimera can bind the TPP ligand, although less avidly.

Identification of KIAA1018/FAN1, a DNA repair nuclease recruited to DNA damage by monoubiquitinated FANCD2.

DNA interstrand crosslinks (ICLs) are highly toxic because they block the progression of replisomes. The Fanconi Anemia (FA) proteins, encoded by genes that are mutated in FA, are important for repair of ICLs. The FA core complex catalyzes the monoubiquitination of FANCD2, and this event is essential for several steps of ICL repair. However, how monoubiquitination of FANCD2 promotes ICL repair at the molecular level is unknown. Here, we describe a highly conserved protein, KIAA1018/MTMR15/FAN1, that interacts with, and is recruited to sites of DNA damage by, the monoubiquitinated form of FANCD2. FAN1 exhibits endonuclease activity toward 5' flaps and has 5' exonuclease activity, and these activities are mediated by an ancient VRR_nuc domain. Depletion of FAN1 from human cells causes hypersensitivity to ICLs, defects in ICL repair, and genome instability. These data at least partly explain how ubiquitination of FANCD2 promotes DNA repair.

Crystal structures of the NAD+-II riboswitch reveal two distinct ligand-binding pockets.

We present crystal structures of a new NAD+-binding riboswitch termed NAD+-II, bound to nicotinamide mononucleotide (NMN), nicotinamide adenine dinucleotide (NAD+) and nicotinamide riboside (NR). The RNA structure comprises a number of structural features including three helices, one of which forms a triple helix by interacting with an A5 strand in its minor-groove, and another formed from a long-range pseudoknot. The core of the structure (centrally located and coaxial with the triplex and the pseudoknot) includes two consecutive quadruple base interactions. Unusually the riboswitch binds two molecules of ligand, bound at distinct, non-overlapping sites in the RNA. Binding occurs primarily through the nicotinamide moiety of each ligand, held by specific hydrogen bonding and stacking interactions with the pyridyl ring. The mode of binding is the same for NMN, NR and the nicotinamide moiety of NAD+. In addition, when NAD+ is bound into one site it adopts an elongated conformation such that its diphosphate linker occupies a groove on the surface of the RNA, following which the adenine portion inserts into a pocket and makes specific hydrogen bonding interactions. Thus the NAD+-II riboswitch is distinct from the NAD+-I riboswitch in that it binds two molecules of ligand at separate sites, and that binding occurs principally through the nicotinamide moiety.

Biochemical and mechanistic analysis of the cleavage of branched DNA by human ANKLE1.

ANKLE1 is a nuclease that provides a final opportunity to process unresolved junctions in DNA that would otherwise create chromosomal linkages blocking cell division. It is a GIY-YIG nuclease. We have expressed an active domain of human ANKLE1 containing the GIY-YIG nuclease domain in bacteria, that is monomeric in solution and when bound to a DNA Y-junction, and unilaterally cleaves a cruciform junction. Using an AlphaFold model of the enzyme we identify the key active residues, and show that mutation of each leads to impairment of activity. There are two components in the catalytic mechanism. Cleavage rate is pH dependent, corresponding to a pKa of 6.9, suggesting an involvement of the conserved histidine in proton transfer. The reaction rate depends on the nature of the divalent cation, likely bound by glutamate and asparagine side chains, and is log-linear with the metal ion pKa. We propose that the reaction is subject to general acid-base catalysis, using a combination of tyrosine and histidine acting as general base and water directly coordinated to the metal ion as general acid. The reaction is temperature dependent; activation energy Ea = 37 kcal mol-1, suggesting that cleavage is coupled to opening of DNA in the transition state.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level.

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibility and the role of co-transcriptional folding in these processes.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology.

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4 kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.

Ribocentre: a database of ribozymes.

Ribozymes are excellent systems in which to study 'sequence - structure - function' relationships in RNA molecules. Understanding these relationships may greatly help structural modeling and design of functional RNA structures and some functional structural modules could be repurposed in molecular design. At present, there is no comprehensive database summarising all the natural ribozyme families. We have therefore created Ribocentre, a database that collects together sequence, structure and mechanistic data on 21 ribozyme families. This includes available information on timelines, sequence families, secondary and tertiary structures, catalytic mechanisms, applications of the ribozymes together with key publications. The database is publicly available at https://www.ribocentre.org.

Structure and mechanism of a methyltransferase ribozyme.

Known ribozymes in contemporary biology perform a limited range of chemical catalysis, but in vitro selection has generated species that catalyze a broader range of chemistry; yet, there have been few structural and mechanistic studies of selected ribozymes. A ribozyme has recently been selected that can catalyze a site-specific methyl transfer reaction. We have solved the crystal structure of this ribozyme at a resolution of 2.3 Å, showing how the RNA folds to generate a very specific binding site for the methyl donor substrate. The structure immediately suggests a catalytic mechanism involving a combination of proximity and orientation and nucleobase-mediated general acid catalysis. The mechanism is supported by the pH dependence of the rate of catalysis. A selected methyltransferase ribozyme can thus use a relatively sophisticated catalytic mechanism, broadening the range of known RNA-catalyzed chemistry.

Structure and folding of four putative kink turns identified in structured RNA species in a test of structural prediction rules.

k-Turns are widespread key architectural elements that occur in many classes of RNA molecules. We have shown previously that their folding properties (whether or not they fold into their tightly kinked structure on addition of metal ions) and conformation depend on their local sequence, and we have elucidated a series of rules for prediction of these properties from sequence. In this work, we have expanded the rules for prediction of folding properties, and then applied the full set to predict the folding and conformation of four probable k-turns we have identified amongst 224 structured RNA species found in bacterial intergenenic regions by the Breaker lab (1). We have analyzed the ion-dependence of folding of the four k-turns using fluorescence resonance energy transfer, and determined the conformation of two of them using X-ray crystallography. We find that the experimental data fully conform to both the predicted folding and conformational properties. We conclude that our folding rules are robust, and can be applied to new k-turns of unknown characteristics with confidence.

Structure and ligand binding of the ADP-binding domain of the NAD+ riboswitch.

The nadA motif is the first known NAD+-dependent riboswitch, comprising two similar tandem bulged stem-loop structures. We have determined the structure of the 5' domain 1 of the riboswitch. It has three coaxial helical segments, separated by an ACANCCCC bulge and by an internal loop, with a tertiary contact between them that includes two C:G base pairs. We have determined the structure with a number of ligands related to NADH, but in each case only the ADP moiety is observed. The adenosine adopts an anti conformation, forms multiple hydrogen bonds across the width of the sugar edge of the penultimate C:G base pair of the helix preceding the bulge, and the observed contacts have been confirmed by mutagenesis and calorimetry. Two divalent metal ions play a key structural role at the narrow neck of the bulge. One makes direct bonding contacts to the diphosphate moiety, locking it into position. Thus the nucleobase, ribose, and phosphate groups of the ADP moiety are all specifically recognized by the RNA. The NAD+ riboswitch is modular. Domain 1 is an ADP binding domain that may be ancient and could potentially be used in combination with other ligand binding motifs such as CoA.

Functional organization of box C/D RNA-guided RNA methyltransferase.

Box C/D RNA protein complexes (RNPs) catalyze site-specific 2'-O-methylation of RNA with specificity determined by guide RNAs. In eukaryotic C/D RNP, the paralogous Nop58 and Nop56 proteins specifically associate with terminal C/D and internal C'/D' motifs of guide RNAs, respectively. We have reconstituted active C/D RNPs with recombinant proteins of the thermophilic yeast Chaetomium thermophilum. Nop58 and Nop56 could not distinguish between the two C/D motifs in the reconstituted enzyme, suggesting that the assembly specificity is imposed by trans-acting factors in vivo. The two C/D motifs are functionally independent and halfmer C/D RNAs can also guide site-specific methylation. Extensive pairing between C/D RNA and substrate is inhibitory to modification for both yeast and archaeal C/D RNPs. N6-methylated adenine at box D/D' interferes with the function of the coupled guide. Our data show that all C/D RNPs share the same functional organization and mechanism of action and provide insight into the assembly specificity of eukaryotic C/D RNPs.

Crystal structure and ligand-induced folding of the SAM/SAH riboswitch.

While most SAM riboswitches strongly discriminate between SAM and SAH, the SAM/SAH riboswitch responds to both ligands with similar apparent affinities. We have determined crystal structures of the SAM/SAH riboswitch bound to SAH, SAM and other variant ligands at high resolution. The riboswitch forms an H-type pseudoknot structure with coaxial alignment of the stem-loop helix (P1) and the pseudoknot helix (PK). An additional three base pairs form at the non-open end of P1, and the ligand is bound at the interface between the P1 extension and the PK helix. The adenine nucleobase is stacked into the helix and forms a trans Hoogsteen-Watson-Crick base pair with a uridine, thus becoming an integral part of the helical structure. The majority of the specific interactions are formed with the adenosine. The methionine or homocysteine chain lies in the groove making a single hydrogen bond, and there is no discrimination between the sulfonium of SAM or the thioether of SAH. Single-molecule FRET analysis reveals that the riboswitch exists in two distinct conformations, and that addition of SAM or SAH shifts the population into a stable state that likely corresponds to the form observed in the crystal. A model for translational regulation is presented whereby in the absence of ligand the riboswitch is largely unfolded, lacking the PK helix so that translation can be initiated at the ribosome binding site. But the presence of ligand stabilizes the folded conformation that includes the PK helix, so occluding the ribosome binding site and thus preventing the initiation of translation.

The role of RNA structure in translational regulation by L7Ae protein in archaea.

A recent study has shown that archaeal L7Ae binds to a putative k-turn structure in the 5'-leader of the mRNA of its structural gene to regulate translation. To function as a regulator, the RNA should be unstructured in the absence of protein, but it should adopt a k-turn-containing stem-loop on binding L7Ae. Sequence analysis of UTR sequences indicates that their k-turn elements will be unable to fold in the absence of L7Ae, and we have demonstrated this experimentally in solution using FRET for the Archaeoglobus fulgidus sequence. We have solved the X-ray crystal structure of the complex of the A. fulgidus RNA bound to its cognate L7Ae protein. The RNA adopts a standard k-turn conformation that is specifically recognized by the L7Ae protein, so stabilizing the stem-loop. In-line probing of the natural-sequence UTR shows that the RNA is unstructured in the absence of L7Ae binding, but folds on binding the protein such that the ribosome binding site is occluded. Thus, L7Ae regulates its own translation by switching the conformation of the RNA to alter accessibility.

Structure-guided design of a high-affinity ligand for a riboswitch.

We have designed structure-based ligands for the guanidine-II riboswitch that bind with enhanced affinity, exploiting the twin binding sites created by loop-loop interaction. We synthesized diguanidine species, comprising two guanidino groups covalently connected by Cn linkers where n = 4 or 5. Calorimetric and fluorescent analysis shows that these ligands bind with a 10-fold higher affinity to the riboswitch compared to guanidine. We determined X-ray crystal structures of the riboswitch bound to the new ligands, showing that the guanidino groups are bound to both nucleobases and backbone within the binding pockets, analogously to guanidine binding. The connecting chain passes through side openings in the binding pocket and traverses the minor groove of the RNA. The combination of the riboswitch loop-loop interaction and our novel ligands has potential applications in chemical biology.

Junction resolving enzymes use multivalency to keep the Holliday junction dynamic.

Holliday junction (HJ) resolution by resolving enzymes is essential for chromosome segregation and recombination-mediated DNA repair. HJs undergo two types of structural dynamics that determine the outcome of recombination: conformer exchange between two isoforms and branch migration. However, it is unknown how the preferred branch point and conformer are achieved between enzyme binding and HJ resolution given the extensive binding interactions seen in static crystal structures. Single-molecule fluorescence resonance energy transfer analysis of resolving enzymes from bacteriophages (T7 endonuclease I), bacteria (RuvC), fungi (GEN1) and humans (hMus81-Eme1) showed that both types of HJ dynamics still occur after enzyme binding. These dimeric enzymes use their multivalent interactions to achieve this, going through a partially dissociated intermediate in which the HJ undergoes nearly unencumbered dynamics. This evolutionarily conserved property of HJ resolving enzymes provides previously unappreciated insight on how junction resolution, conformer exchange and branch migration may be coordinated.

Cooperative control of holliday junction resolution and DNA repair by the SLX1 and MUS81-EME1 nucleases.

Holliday junctions (HJs) are X-shaped DNA structures that arise during homologous recombination, which must be removed to enable chromosome segregation. The SLX1 and MUS81-EME1 nucleases can both process HJs in vitro, and they bind in close proximity on the SLX4 scaffold, hinting at possible cooperation. However, the cellular roles of mammalian SLX1 are not yet known. Here, we use mouse genetics and structure function analysis to investigate SLX1 function. Disrupting the murine Slx1 and Slx4 genes revealed that they are essential for HJ resolution in mitotic cells. Moreover, SLX1 and MUS81-EME1 act together to resolve HJs in a manner that requires tethering to SLX4. We also show that SLX1, like MUS81-EME1, is required for repair of DNA interstrand crosslinks, but this role appears to be independent of HJ cleavage, at least in mouse cells. These findings shed light on HJ resolution in mammals and on maintenance of genome stability.

Structure and ligand binding of the glutamine-II riboswitch.

We have determined the structure of the glutamine-II riboswitch ligand binding domain using X-ray crystallography. The structure was solved using a novel combination of homology modeling and molecular replacement. The structure comprises three coaxial helical domains, the central one of which is a pseudoknot with partial triplex character. The major groove of this helix provides the binding site for L-glutamine, which is extensively hydrogen bonded to the RNA. Atomic mutation of the RNA at the ligand binding site leads to loss of binding shown by isothermal titration calorimetry, explaining the specificity of the riboswitch. A metal ion also plays an important role in ligand binding. This is directly bonded to a glutamine carboxylate oxygen atom, and its remaining inner-sphere water molecules make hydrogen bonding interactions with the RNA.

The kink-turn in the structural biology of RNA.

The kink-turn (k-turn) is a widespread structural motif found in functional RNA species. It typically comprises a three-nucleotide bulge followed by tandem trans sugar edge-Hoogsteen G:A base pairs. It introduces a sharp kink into the axis of duplex RNA, juxtaposing the minor grooves. Cross-strand H-bonds form at the interface, accepted by the conserved adenine nucleobases of the G:A basepairs. Alternative acceptors for one of these divides the k-turns into two conformational classes N3 and N1. The base pair that follows the G:A pairs (3b:3n) determines which conformation is adopted by a given k-turn. k-turns often mediate tertiary contacts in folded RNA species and frequently bind proteins. Common k-turn binding proteins include members of the L7Ae family, such as the human 15·5k protein. A recognition helix within these proteins binds in the widened major groove on the outside of the k-turn, that makes specific H-bonds with the conserved guanine nucleobases of the G:A pairs. L7Ae binds with extremely high affinity, and single-molecule data are consistent with folding by conformational selection. The standard, simple k-turn can be elaborated in a variety of ways, that include the complex k-turns and the k-junctions. In free solution in the absence of added metal ions or protein k-turns do not adopt the tightly-kinked conformation. They undergo folding by the binding of proteins, by the formation of tertiary contacts, and some (but not all) will fold on the addition of metal ions. Whether or not folding occurs in the presence of metal ions depends on local sequence, including the 3b:3n position, and the -1b:-1n position (5' to the bulge). In most cases -1b:-1n = C:G, so that the 3b:3n position is critical since it determines both folding properties and conformation. In general, the selection of these sequence matches a given k-turn to its biological requirements. The k-turn structure is now very well understood, to the point at which they can be used as a building block for the formation of RNA nano-objects, including triangles and squares.

Structure and ligand binding of the SAM-V riboswitch.

SAM-V is one of the class of riboswitches that bind S-adenosylmethione, regulating gene expression by controlling translation. We have solved the crystal structure of the metY SAM-V riboswitch bound to its SAM ligand at 2.5 Å resolution. The RNA folds as an H-type pseudoknot, with a major-groove triple helix in which resides the SAM ligand binding site. The bound SAM adopts an elongated conformation aligned with the axis of the triple helix, and is held at either end by hydrogen bonding to the adenine and the amino acid moieties. The central sulfonium cation makes electrostatic interactions with an U:A.U base triple, so conferring specificity. We propose a model in which SAM binding leads to association of the triplex third strand that stabilizes a short helix and occludes the ribosome binding site. Thus the new structure explains both ligand specificity and the mechanism of genetic control.

A monovalent ion in the DNA binding interface of the eukaryotic junction-resolving enzyme GEN1.

GEN1 is a member of the FEN/EXO family of structure-selective nucleases that cleave 1 nt 3' to a variety of branchpoints. For each, the H2TH motif binds a monovalent ion and plays an important role in binding one helical arm of the substrates. We investigate here the importance of this metal ion on substrate specificity and GEN1 structure. In the presence of K+ ions the substrate specificity is wider than in Na+, yet four-way junctions remain the preferred substrate. In a combination of K+ and Mg2+ second strand cleavage is accelerated 17-fold, ensuring bilateral cleavage of the junction. We have solved crystal structures of Chaetomium thermophilum GEN1 with Cs+, K+ and Na+ bound. With bound Cs+ the loop of the H2TH motif extends toward the active site so that D199 coordinates a Mg2+, buttressed by an interaction of the adjacent Y200. With the lighter ions bound the H2TH loop changes conformation and retracts away from the active site. We hypothesize this conformational change might play a role in second strand cleavage acceleration.