Deep mutational scanning results of the LINE-1 and OR4K15 ribozymes.

(A) Average relative activity of mutations at each sequence position of the original LINE-1 ribozyme (LINE-1-ori). (B) Average relative activity of mutations at each sequence position of the original OR4K15 ribozyme (OR4K15-ori). (C) (Left) Base pairing maps inferred from LINE-1-ori deep mutation data by CODA (lower triangle) and by CODA in combination with Monte Carlo simulated annealing (CODA+MC) (upper triangle). (Right) Same as (Left) but for the contiguous functional region corresponds to 54-99 of LINE-1-ori (LINE-1-rbz). (D) (Left) Base pairing maps inferred from OR4K15-ori deep mutation data by CODA (lower triangle) and by CODA+MC (upper triangle). (Right) Same as (Left) but for the contiguous functional region which corresponds to 70-116 of OR4K15-ori (OR4K15-rbz).

Secondary structures of the functional regions of the LINE-1 and OR4K15 ribozymes.

(A) Relative activity of all mutations at each sequence position of the functional region of LINE-1 ribozyme (LINE-1-rbz). (B) Comparison between the base pairing maps inferred from CODA analysis of deep mutational scanning of the functional region of LINE-1 ribozyme (LINE-1-rbz) (lower triangle) and that from the corresponding region (54-99 nt) in the deep mutational scanning of the original LINE-1 ribozyme (LINE-1-ori) (upper triangle). (C) Comparison between the base pairing map inferred from LINE-1-rbz deep mutational data by CODA (lower triangle) and that after Monte-Carlo simulated annealing (CODA+MC, upper triangle). (D) Secondary structure model (left: LINE-1-rbz; right: OR4K15-rbz) inferred from the deep mutational scanning result. The red trianglesindicate the cleavage sites of ribozymes. (E) Consensus sequence and secondary structure model (left: LINE-1-rbz; right: OR4K15-rbz) based on the alignment of functional mutants with RA′ 3 0.5. The red triangles indicate the ribozyme cleavage sites. Positions with conservation of 95, 98, and 100% were marked with gray, black and red nucleotides, respectively; positions in which nucleotide identity is less conserved are represented by circles. Green shading denotes predicted base pairs supported by covariation. R and Y denote purine and pyrimidine, respectively.

Homology modelled 3D structure of the two-stranded LINE-1-core and OR4K15-core.

(A) Comparison between the secondary structure of LINE-1-rbz, OR4K15-rbz, four-way (PDB ID: 5Y87) and three-way junctional twister-sister ribozyme (PDB ID: 5T5A) reveals strikingly similar internal loops surrounding the cleavage site (red triangles) with few nucleotide differences at each side of the internal loops L1. (B) A cartoon view of the homology modelled structure of LINE-1-core (left) and OR4K15-core (right). (C) A detailed view of the catalytic core L1 of LINE-1-core in the cartoon representation.

Special features of the homology modelled LINE-1-core structure.

(A) The network of H-bonding interactions involving G7. (B) Intermolecular contacts at the C10-A11 cleavage site in the modelled structure. The scissile phosphate is colored magenta. (C) PAGE analysis of the cleavage activity of wild-type (WT) and mutants C10U, C10G, C10A, and A11U.

(A) (Top) The components of the bimolecular construct of LINE-1-core during cleavage, substrate strand (54-71) and enzyme strand (83-99).

(Bottom) The components of the bimolecular construct of OR4K15-core during cleavage, the enzyme strand (70-84) and substrate strand (101-116). The red arrowhead indicates the cleavage site. (B) PAGE-based cleavage assays of the bimolecular construct with the all-RNA substrate and a substrate analog (dC) wherein a 2′-deoxycytosine is substituted for the cytosine ribonucleotide at position 10 or 40 of the substrate RNA. The substrate RNA was incubated with the enzyme RNA at 37ºC for 1h in the presence (+) or absence (-) of 1mM MgCl2. The single-channel fluorescent images (left) were generated by using UV excitation (302 nm) and 590/110 nm emission on the ChemiDoc MP imaging system. The multi-channel fluorescent images (right) are overlays of two scans. They were generated from the ChemiDoc MP imaging system (Bio-Rad), Fluorescein (excitation: Epi-Blue 460-490 nm, emission: 532/28) for FAM, Cy3 (excitation: Epi-Green 520-545 nm, emission: 602/50) for TAMRA. (C) Cleavage assays in the absence (-) or presence of cobalt hexammine chloride [Co(NH3)6Cl3] or MgCl2 at 5mM for 1h. (D) Cleavage assays at 37ºC for 1h, in the absence (-) or presence (+) of various metal ions. Divalent metal ions (left) were used at a final concentration of 1mM, and monovalent metal ions (right) were used at a final concentration of 1M.

The number of mutants and the coverage of single and double mutations of LINE-1 ribozyme in two deep mutational scanning experiments for the original sequence (LINE-1-ori) and the functional region (LINE-1-rbz), respectively.

Summary of two deep sequencing results for the original sequence (LINE-1-ori) and the functional region (LINE-1-rbz), respectively.

Sequences with similar motifs found in a specialized Rfam database.

Oligonucleotides used in this study.

Mutation rates of two deep mutational scanning results at each nucleotide position of LINE-1 ribozyme.

Mutation rates at each nucleotide position of OR4K15 ribozyme

Distribution of different mutation types in deep sequencing result of LINE-1 ribozyme.

Distribution of different mutation types in the deep sequencing result of OR4K15 ribozyme.

Single mutation-based average RA′/RA distribution of LINE-1-rbz.

(A) Structure descriptor for LINE-1-rbz. (B) Structure descriptor for OR4K15-rbz.

(A) The structure of LINE-1-core in a surface representation, except for the C10-A11 site, which is shown in a stick representation. (B) Continuous stacking of stems P1 and P2. (C) Stacking interactions in loop L1 region.

Base-pairing interactions in loop L1 of LINE-1-core that extend the stem P1 and stem P2.

(A) A detailed view of the interactions in loop L1 of OR4K15-core. (B) Stacking interactions in loop L1 region. (C) The network of H-bonding interactions involving G37. (D) Intermolecular contacts at the C40-A41 cleavage site in the modelled structure. The scissile phosphate is colored magenta.

PAGE-based cleavage assay of LINE-1-core in a 24-hour range.

The dependence of two lantern ribozymes’ rate constants on Mg2+/Mn2+ concentration.

A representative time course for the bimolecular ribozyme construct LINE-1-core under the conditions indicated as used for determining kobs values.

Experimental pipeline of two deep mutational scanning experiments.

(A) The mutation library of the original sequence of LINE-1 ribozyme was generated by error-prone PCR. The RNA-seq library was constructed by using the cDNA synthesized by reverse transcription and template switching. (B) The chemically synthesized doped library was used to construct the mutation library of the contiguous, functional region of LINE-1 ribozyme. RtcB ligase was used to captured 3′ fragment after cleavage, then the ligated product was used to construct the RNA-seq library.