Double Hairpin Research Articles

Structure of One-Phase Free Boundaries in the Plane

We study classical solutions to the one-phase free boundary problem in which the free boundary consists of smooth curves and the components of the positive phase are simply-connected. We show that if two components of the free boundary are close, then the solution locally resembles an entire solution discovered by Hauswirth, Helein and Pacard, whose free boundary has the shape of a double hairpin. Our results are analogous to theorems of Colding and Minicozzi characterizing embedded minimal annuli, and a direct connection between our theorems and theirs can be made using a correspondence due to Traizet.

Folding a stable RNA pseudoknot through rearrangement of two hairpin structures

Folding messenger RNA into specific structures is a common regulatory mechanism involved in translation. In Escherichia coli, the operator of the rpsO gene transcript folds into a pseudoknot or double-hairpin conformation. S15, the gene product, binds only to the pseudoknot, thereby repressing its own synthesis when it is present in excess in the cell. The two RNA conformations have been proposed to exist in equilibrium. However, it remained unclear how structural changes can be achieved between these two topologically distinct conformations. We used optical tweezers to study the structural dynamics and rearrangements of the rpsO operator RNA at the single-molecule level. We discovered that the two RNA structures can be interchanged spontaneously and the pseudoknot can exist in conformations that exhibit various levels of stability. Conversion from the double hairpin to a pseudoknot through potential hairpin–hairpin interactions favoured the high-stability conformation. By contrast, mutations that blocked the formation of a hairpin typically resulted in alternative low-stability pseudoknots. These results demonstrate that specific tertiary interactions of RNA can be established and modulated based on the interactions and rearrangements between secondary structural components. Our findings provide new insight into the RNA folding pathway that leads to a regulatory conformation for target protein binding.

Conservation and Functional Element Discovery in 20 Angiosperm Plant Genomes

Here, we describe the construction of a phylogenetically deep, whole-genome alignment of 20 flowering plants, along with an analysis of plant genome conservation. Each included angiosperm genome was aligned to a reference genome, Arabidopsis thaliana, using the LASTZ/MULTIZ paradigm and tools from the University of California-Santa Cruz Genome Browser source code. In addition to the multiple alignment, we created a local genome browser displaying multiple tracks of newly generated genome annotation, as well as annotation sourced from published data of other research groups. An investigation into A. thaliana gene features present in the aligned A. lyrata genome revealed better conservation of start codons, stop codons, and splice sites within our alignments (51% of features from A. thaliana conserved without interruption in A. lyrata) when compared with previous publicly available plant pairwise alignments (34% of features conserved). The detailed view of conservation across angiosperms revealed not only high coding-sequence conservation but also a large set of previously uncharacterized intergenic conservation. From this, we annotated the collection of conserved features, revealing dozens of putative noncoding RNAs, including some with recorded small RNA expression. Comparing conservation between kingdoms revealed a faster decay of vertebrate genome features when compared with angiosperm genomes. Finally, conserved sequences were searched for folding RNA features, including but not limited to noncoding RNA (ncRNA) genes. Among these, we highlight a double hairpin in the 5'-untranslated region (5'-UTR) of the PRIN2 gene and a putative ncRNA with homology targeting the LAF3 protein.

Crystal structure of Δ(185-243)ApoA-I suggests a mechanistic framework for the protein adaptation to the changing lipid load in good cholesterol: from flatland to sphereland via double belt, belt buckle, double hairpin and trefoil/tetrafoil.

Apolipoprotein A-I (apoA-I) is the major protein of plasma high-density lipoproteins (HDLs), macromolecular assemblies of proteins and lipids that remove cell cholesterol and protect against atherosclerosis. HDL heterogeneity, large size (7.7-12 nm), and ability to exchange proteins have prevented high-resolution structural analysis. Low-resolution studies showed that two apoA-I molecules form an antiparallel α-helical "double belt" around an HDL particle. The atomic-resolution structure of the C-terminal truncated lipid-free Δ(185-243)apoA-I, determined recently by Mei and Atkinson, provides unprecedented new insights into HDL structure-function. It allows us to propose a molecular mechanism for the adaptation of the full-length protein to increasing lipid load during cholesterol transport. ApoA-I conformations on small, midsize, and large HDLs are proposed based on the tandem α-helical repeats and the crystal structure of Δ(185-243)apoA-I and are validated by comparison with extensive biophysical data reported by many groups. In our models, the central half of the double belt ("constant" segment 66-184) is structurally conserved while the N- and C-terminal half ("variable" segments 1-65 and 185-243) rearranges upon HDL growth. This includes incremental unhinging of the N-terminal bundle around two flexible regions containing G39 and G65 to elongate the belt, along with concerted swing motion of the double belt around G65-P66 and G185-G186 hinges that are aligned on various-size particles, to confer two-dimensional surface curvature to spherical HDLs. The proposed conformational ensemble integrates and improves several existing HDL models. It helps provide a structural framework necessary to understand functional interactions with over 60 other HDL-associated proteins and, ultimately, improve the cardioprotective function of HDL.

The functional Hfq-binding module of bacterial sRNAs consists of a double or single hairpin preceded by a U-rich sequence and followed by a 3′ poly(U) tail

Hfq-dependent sRNAs contain, at least, an mRNA base-pairing region, an Hfq-binding site, and a Rho-independent terminator. Recently, we found that the terminator poly(U) of Escherichia coli sRNAs is essential for Hfq binding and therefore for riboregulation. In this study, we tried to identify additional components within Hfq-binding sRNAs required for efficient Hfq binding by using SgrS as a model. We demonstrate by mutational and biochemical studies that an internal hairpin and an immediately upstream U-rich sequence also are required for efficient Hfq binding. We propose that the functional Hfq-binding module of SgrS consists of an internal hairpin preceded by a U-rich sequence and a Rho-independent terminator with a long poly(U) tail. We also show that the Rho-independent terminator alone can act as a functional Hfq-binding module when it is preceded by an internal U-rich sequence. The 3' region of most known sRNAs share the features corresponding to either a double- or single-hairpin-type Hfq-binding module. We also demonstrate that increasing the spacing between the base-pairing region and the Hfq-binding module reduces or impairs the silencing ability. These findings allowed us to design synthetic Hfq-binding sRNAs to target desired mRNAs.

The effect of context on the folding of β-hairpins

Small β-hairpin peptides have been widely used as models for the folding of β-sheets. But how applicable is the folding of such models to β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding simulations of three such proteins that contain the WW domain double hairpin β-sheet motif: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain. We compare the behavior of the free motif in solution and in the context of proteins of different size and architecture. Both Csp proteins lost contacts between the double-hairpin motif and the protein core as the first step of unfolding and proceeded to unfold with loss of the third β-strand, similar to the isolated WW domain. The glucose permease IIA domain is a larger protein and the contacts between the motif and the core were not lost as quickly. Instead the unfolding pathway of glucose permease IIA followed a different pathway with β1 pulling away from the sheet first. Interestingly, when the double hairpin motif was excised from the glucose permease IIA domain and simulated in isolation in water it unfolded by the same pathway as the WW domain, indicating that it is tertiary interactions with the protein that alter the motif's unfolding not a sequence dependent effect on its intrinsic unfolding behavior. With respect to the unfolding of the hairpins, there was no consistent order to the loss of hydrogen bonds between the β-strands in the hairpins in any of the systems. Our results show that while the folding behavior of the isolated WW domain is generally consistent with the double hairpin motif's behavior in the cold shock proteins, it is not the case for the glucose permease IIA domain. So, one must be cautious in extrapolating findings from model systems to larger more complicated proteins where tertiary interactions can overwhelm intrinsic behavior.

A conserved RNA pseudoknot in a putative molecular switch domain of the 3′-untranslated region of coronaviruses is only marginally stable

The 3'-untranslated region (UTR) of the group 2 coronavirus mouse hepatitis virus (MHV) genome contains a predicted bulged stem-loop (designated P0ab), a conserved cis-acting pseudoknot (PK), and a more distal stem-loop (designated P2). Base-pairing to create the pseudoknot-forming stem (P1(pk)) is mutually exclusive with formation of stem P0a at the base of the bulged stem-loop; as a result, the two structures cannot be present simultaneously. Herein, we use thermodynamic methods to evaluate the ability of individual subdomains of the 3' UTR to adopt a pseudoknotted conformation. We find that an RNA capable of forming only the predicted PK (58 nt; 3' nucleotides 241-185) adopts the P2 stem-loop with little evidence for P1(pk) pairing in 0.1 M KCl and the absence of Mg(2+); as Mg(2+) or 1 M KCl is added, a new thermal unfolding transition is induced and assignable to P1(pk) pairing. The P1(pk) helix is only marginally stable, ΔG(25) ≈ 1.2 ± 0.3 kcal/mol (5.0 mM Mg(2+), 100 mM K(+)), and unfolded at 37°C. Similar findings characterize an RNA 5' extended through the P0b helix only (89 nt; 294-185). In contrast, an RNA capable of forming either the P0a helix or the pseudoknot (97 nt; 301-185) forms no P1(pk) helix. Thermal unfolding simulations are fully consistent with these experimental findings. These data reveal that the PK forms weakly and only when the competing double-hairpin structure cannot form; in the UTR RNA, the double hairpin is the predominant conformer under all solution conditions.

Presence of Divalent Cation Is Not Mandatory for the Formation of Intramolecular Purine-Motif Triplex Containing Human c-jun Protooncogene Target

Modulation of endogenous gene function, through sequence-specific recognition of double helical DNA via oligonucleotide-directed triplex formation, is a promising approach. Compared to the formation of pyrimidine motif triplexes, which require relatively low pH, purine motif appears to be the most gifted for their stability under physiological conditions. Our previous work has demonstrated formation of magnesium-ion dependent highly stable intermolecular triplexes using a purine third strand of varied lengths, at the purine•pyrimidine (Pu•Py) targets of SIV/HIV-2 (vpx) genes (Svinarchuk, F., Monnot, M., Merle, A., Malvy, C., and Fermandjian, S. (1995) Nucleic Acids Res. 23, 3831-3836). Herein, we show that a designed intramolecular version of the 11-bp core sequence of the said targets, which also constitutes an integral, short, and symmetrical segment (G(2)AG(5)AG(2))•(C(2)TC(5)TC(2)) of human c-jun protooncogene forms a stable triplex, even in the absence of magnesium. The sequence d-C(2)TC(5)TC(2)T(5)G(2)AG(5)AG(2)T(5)G(2)AG(5)AG(2) (I-Pu) folds back twice onto itself to form an intramolecular triple helix via a double hairpin formation. The design ensures that the orientation of the intact third strand is antiparallel with respect to the oligopurine strand of the duplex. The triple helix formation has been revealed by non-denaturating gel assays, UV-thermal denaturation, and circular dichroism (CD) spectroscopy. The monophasic melting curve, recorded in the presence of sodium, represented the dissociation of intramolecular triplex to single strand in one step; however, the addition of magnesium bestowed thermal stability to the triplex. Formation of intramolecular triple helix at neutral pH in sodium, with or without magnesium cations, was also confirmed by gel electrophoresis. The triplex, mediated by sodium alone, destabilizes in the presence of 5'-C(2)TC(5)TC(2)-3', an oligonucleotide complementary to the 3'-oligopurine segments of I-Pu, whereas in the presence of magnesium the triplex remained impervious. CD spectra showed the signatures of triplex structure with A-like DNA conformation. We suggest that the possible formation of pH and magnesium-independent purine-motif triplexes at genomic Pu•Py sequences may be pertinent to gene regulation.

Erratum to: The Presence in tRNA Molecule Sequences of the Double Hairpin, an Evolutionary Stage Through Which the Origin of this Molecule is Thought to have Passed

Erratum to: The Presence in tRNA Molecule Sequences of the Double Hairpin, an Evolutionary Stage Through Which the Origin of this Molecule is Thought to have Passed

The Presence in tRNA Molecule Sequences of the Double Hairpin, an Evolutionary Stage Through Which the Origin of this Molecule is Thought to have Passed

We analyse 6,810 tRNAs, calculating the free energy of the corresponding double hairpin and 'cigar' secondary structures, for which we find a high thermodynamic and statistical significance. We also analyse these tRNAs for similarity and complementarity of their 5' and 3' halves or segments of them in intra-and inter-molecular comparisons. We find very clear signs that the two halves of tRNAs had an evident evolutionary relationship, although it is not totally clear whether this was a relationship of homology or complementarity between the 5' and 3' halves of tRNAs, even if there is strong evidence in favour of the homology hypothesis. Overall, these data favour models for the origin of the tRNA molecule postulating that a duplication event involving a hairpin structure as a precursor was involved in the origin of this molecule. Moreover, we interpret these results and favour the hypothesis that sees the assembly of two hairpin structures sharing a homology relationship as the intermediate evolutionary stage preceding the appearance of the cloverleaf structure of tRNA.

Synthesis, gene-silencing activity and nuclease resistance of 3′–3′-linked double short hairpin RNA

To improve the nuclease resistance of siRNA while reducing its induction of an innate immune response and maintaining its biological activity for possible therapeutic application, we designed and synthesized a series of double short hairpin RNAs (dshRNAs). Each dshRNA consisted of two identical short hairpin RNAs (shRNAs) linked at their 3′ ends by glycerol. The dshRNAs were synthesized on a glycerol-derivatized solid support from amidites with 2-cyanoethoxymethyl (CEM) as the 2′-hydroxyl protecting group. Synthesis was carried out in a single run on a DNA/RNA synthesizer, without the need for enzymatic ligation. The dshRNAs showed structure-dependent gene-silencing activity at the protein level, and dshRNAs in which the 3′ end of the two sense regions were linked showed especially high activity. Inclusion of 2′-O-methyluridine residues in the loop region was associated with 1.6- to 2.4-fold lower induction of interferon-α than was siRNA, without loss of gene-silencing activity. dshRNA also showed higher exonuclease resistance than siRNA or canonical shRNA. Our studies provide a new approach to gene silencing based on the concept of linking the 3′ end of the sense regions of two shRNA molecules to form a double shRNA.

Structure of a Conserved Retroviral RNA Packaging Element by NMR Spectroscopy and Cryo-Electron Tomography

The 5′-untranslated regions of all gammaretroviruses contain a conserved “double-hairpin motif” (Ψ CD) that is required for genome packaging. Both hairpins (SL-C and SL-D) contain GACG tetraloops that, in isolated RNAs, are capable of forming “kissing” interactions stabilized by two intermolecular G–C base pairs. We have determined the three-dimensional structure of the double hairpin from the Moloney murine leukemia virus ([Ψ CD] 2, 132 nt, 42.8 kDa) using a 2H-edited NMR-spectroscopy-based approach. This approach enabled the detection of 1H– 1H dipolar interactions that were not observed in previous studies of isolated SL-C and SL-D hairpin RNAs using traditional 1H– 1H correlated and 1H– 13C-edited NMR methods. The hairpins participate in intermolecular cross-kissing interactions (SL-C to SL-D′ and SLC′ to SL-D) and stack in an end-to-end manner (SL-C to SL-D and SL-C′ to SL-D′) that gives rise to an elongated overall shape (ca 95 Å × 45 Å × 25 Å). The global structure was confirmed by cryo-electron tomography (cryo-ET), making [Ψ CD] 2 simultaneously the smallest RNA to be structurally characterized to date by cryo-ET and among the largest to be determined by NMR. Our findings suggest that, in addition to promoting dimerization, [Ψ CD] 2 functions as a scaffold that helps initiate virus assembly by exposing a cluster of conserved UCUG elements for binding to the cognate nucleocapsid domains of assembling viral Gag proteins.

Deriving four functional anti-HIV siRNAs from a single Pol III-generated transcript comprising two adjacent long hairpin RNA precursors.

Several different approaches exist to generate expressed RNA interference (RNAi) precursors for multiple target inhibition, a strategy referred to as combinatorial (co)RNAi. One such approach makes use of RNA Pol III-expressed long hairpin RNAs (lhRNAs), which are processed by Dicer to generate multiple unique short interfering siRNA effectors. However, because of inefficient intracellular Dicer processing, lhRNA duplexes have been limited to generating two independent effective siRNA species. In this study, we describe a novel strategy whereby four separate anti-HIV siRNAs were generated from a single RNA Pol III-expressed transcript. Two optimized lhRNAs, each comprising two active anti-HIV siRNAs, were placed in tandem to form a double long hairpin (dlhRNA) expression cassette, which encodes four unique and effective siRNA sequences. Processing of the 3′ position lhRNA was more variable but effective multiple processing was possible by manipulating the order of the siRNA-encoding sequences. Importantly, unlike shRNAs, Pol III-expressed dlhRNAs did not compete with endogenous and exogenous microRNAs to disrupt the RNAi pathway. The versatility of expressed lhRNAs is greatly expanded and we provide a mechanism for generating transcripts with modular lhRNAs motifs that contribute to improved coRNAi.

Targeting DNA Hairpin Loops with their Partially Complementary Strands

Antisense, antigene and siRNA strategies are currently used to control the expression of genes. To this end, our laboratory is currently mimicking the targeting of mRNA by targeting DNA hairpin loops with their partially complementary strands. We use a combination of differential scanning calorimetry (DSC) and temperature-dependent UV spectroscopy to investigate the reaction of a variety of DNA hairpin loops (single end loops, dumbbell, three-way junction with two loops and a hairpin with a bulge of 5 nucleotides) with single strands that are complementary to the bases in the loop and to one strand of their stem. The resulting reaction products form duplexes with dangling ends, nicks or with a displaced strand. We determine standard thermodynamic profiles for the unfolding of the reactants (hairpin loop) and products (duplex) of each reaction. The DSC and UV melting curves show monophasic transitions for the unfolding of all DNA single hairpin loops (reactants) and biphasic transitions for the unfolding of the double hairpin loops (reactants) and duplex products. The resulting unfolding data is then used to create thermodynamic (Hess) cycles that correspond to each targeting reaction. All eight targeting reactions investigated yielded favorable free energy contributions that were enthalpy driven. These favorable heat contributions result from the formation of base-pair stacks involving the unpaired bases of the loops, indicating that each single strand was able to disrupt the hairpin loop structure. Supported by Grant MCB-0616005 from NSF.

Computational analysis and determination of a highly conserved surface exposed segment in H5N1 avian flu and H1N1 swine flu neuraminidase

BackgroundCatalytic activity of influenza neuraminidase (NA) facilitates elution of progeny virions from infected cells and prevents their self-aggregation mediated by the catalytic site located in the body region. Research on the active site of the molecule has led to development of effective inhibitors like oseltamivir, zanamivir etc, but the high rate of mutation and interspecies reassortment in viral sequences and the recent reports of oseltamivir resistant strains underlines the importance of determining additional target sites for developing future antiviral compounds. In a recent computational study of 173 H5N1 NA gene sequences we had identified a 50-base highly conserved region in 3'-terminal end of the NA gene.ResultsWe extend the graphical and numerical analyses to a larger number of H5N1 NA sequences (514) and H1N1 swine flu sequences (425) accessed from GenBank. We use a 2D graphical representation model for the gene sequences and a Graphical Sliding Window Method (GSWM) for protein sequences scanning the sequences as a block of 16 amino acids at a time. Using a protein sequence descriptor defined in our model, the protein sliding scan method allowed us to compare the different strains for block level variability, which showed significant statistical correlation to average solvent accessibility of the residue blocks; single amino acid position variability results in no correlation, indicating the impact of stretch variability in chemical environment. Close to the C-terminal end the GSWM showed less descriptor-variability with increased average solvent accessibility (ASA) that is also supported by conserved predicted secondary structure of 3' terminal RNA and visual evidence from 3D crystallographic structure.ConclusionThe identified terminal segment, strongly conserved in both RNA and protein sequences, is especially significant as it is surface exposed and structural chemistry reveals the probable role of this stretch in tetrameric stabilization. It could also participate in other biological processes associated with conserved surface residues. A RNA double hairpin secondary structure found in this segment in a majority of the H5N1 strains also supports this observation. In this paper we propose this conserved region as a probable site for designing inhibitors for broad-spectrum pandemic control of flu viruses with similar NA structure.

Importance of Protein Context on the Unfolding Pathways of β-hairpins

Small β-hairpin peptides as well as three-stranded WW domains have been used as models for the folding of β-hairpins in larger proteins. Previous studies of the FBP28 WW domain proposed that side chain contacts between residues in the strands and not the precise order of backbone hydrogen bond formation guide β-hairpin folding. But how applicable is the folding of model systems, such as FBP28 WW domain, to the folding of β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding molecular dynamics simulations of three proteins that share a double hairpin motif structurally similar to WW domains: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain.

Families of H/ACA ncRNA molecules in Trypanosomatids

H/ACA molecules are small nucleolar RNAs (snoRNAs) that guide the conversion of a uridine into a pseudouridine. This modification is crucial for ribosomal RNA (rRNA) function. In addition, H/ACA RNA also function in ribosomal RNA processing. Unlike the double hairpin structure of H/ACA molecules in other organisms, H/ACA molecules in Trypanosomes have a single hairpin structure. Here we describe the repertoire of H/ACA molecules of the 46 published of Trypanosoma brucei, as well as identify the orthologous counterparts in other Trypanosomatids enlarging the H/ACA collection in this family to over 300.

Diverse modes of 5′‐[4‐(aminoiminomethyl)phenyl]‐[2,2′‐bifuran]‐5‐carboximidamide (DB832) interaction with multi‐stranded DNA structures

The modes of binding of 5'-[4-(aminoiminomethyl)phenyl]-[2,2'-Bifuran]-5-carboximidamide (DB832) to multi-stranded DNAs: human telomere quadruplex, monomolecular R-triplex, pyr/pur/pyr triplex consisting of 12 T*(T x A) triplets, and DNA double helical hairpin were studied. The optical adsorption of the ligand was used for monitoring the binding and for determination of the association constants and the numbers of binding sites. CD spectra of DB832 complexes with the oligonucleotides and the data on the energy transfer from DNA bases to the bound DB832 assisted in elucidating the binding modes. The affinity of DB832 to the studied multi-stranded DNAs was found to be greater (K(ass) approximately 10(7)M(-1)) than to the duplex DNA (K(ass) approximately 2 x 10(5)M(-1)). A considerable stabilizing effect of DB832 binding on R-triplex conformation was detected. The nature of the ligand tight binding differed for the studied multi-stranded DNA depending on their specific conformational features: recombination-type R-triplex demonstrated the highest affinity for DB832 groove binding, while pyr/pur/pyr TTA triplex favored DB832 intercalation at the end stacking contacts and the human telomere quadruplex d[AG(3)(T(2)AG(3))(3)] accommodated the ligand in a capping mode. Additionally, the pyr/pur/pyr TTA triplex and d[AG(3)(T(2)AG(3))(3)] quadruplex bound DB832 into their grooves, though with a markedly lesser affinity. DB832 may be useful for discrimination of the multi-sranded DNA conformations and for R-triplex stabilization.

A NMR strategy to unambiguously distinguish nucleic acid hairpin and duplex conformations applied to a Xist RNA A-repeat

All RNA sequences that fold into hairpins possess the intrinsic potential to form intermolecular duplexes because of their high self-complementarity. The thermodynamically more stable duplex conformation is favored under high salt conditions and at high RNA concentrations, posing a challenging problem for structural studies of small RNA hairpin conformations. We developed and applied a novel approach to unambiguously distinguish RNA hairpin and duplex conformations for the structural analysis of a Xist RNA A-repeat. Using a combination of a quantitative HNN-COSY experiment and an optimized double isotope-filtered NOESY experiment we could define the conformation of the 26-mer A-repeat RNA. In contrast to a previous secondary structure prediction of a double hairpin structure, the NMR data show that only the first predicted hairpin is formed, while the second predicted hairpin mediates dimerization of the A-repeat by duplex formation with a second A-repeat. The strategy employed here will be generally applicable to identify and quantify populations of hairpin and duplex conformations and to define RNA folding topology from inter- and intra-molecular base-pairing patterns.

A Biochemically Defined System for Coding Joint Formation in V(D)J Recombination

V(D)J recombination is one of the most complex DNA transactions in biology. The RAG complex makes double-stranded breaks adjacent to signal sequences and creates hairpin coding ends. Here, we find that the kinase activity of the Artemis:DNA-PKcs complex can be activated by hairpin DNA ends in cis, thereby allowing the hairpins to be nicked and then to undergo processing and joining by nonhomologous DNA end joining. Based on these insights, we have reconstituted many aspects of the antigen receptor diversification of V(D)J recombination by using 13 highly purified polypeptides, thereby permitting variable domain exon assembly by using this fully defined system in accord with the 12/23 rule for this process. The features of the recombination sites created by this system include all of the features observed in vivo (nucleolytic resection, P nucleotides, and N nucleotide addition), indicating that most, if not all, of the end modification enzymes have been identified.