Next-generation sequencing has proven an extremely effective technology for molecular counting applications where the number of sequence reads provides a digital readout for RNA-seq, ChIP-seq, Tn-seq and other applications. The extremely large number of sequence reads that can be obtained per run permits the analysis of increasingly complex samples. For lower complexity samples, however, a point of diminishing returns is reached when the number of counts per sequence results in oversampling with no increase in data quality. A solution to making next-generation sequencing as efficient and affordable as possible involves assaying multiple samples in a single run. Here, we report the successful 96-plexing of complex pools of DNA barcoded yeast mutants and show that such ‘Bar-seq’ assessment of these samples is comparable with data provided by barcode microarrays, the current benchmark for this application. The cost reduction and increased throughput permitted by highly multiplexed sequencing will greatly expand the scope of chemogenomics assays and, equally importantly, the approach is suitable for other sequence counting applications that could benefit from massive parallelization.

Site-specific labeling of Escherichia coli ribosomes has allowed application of single-molecule fluorescence spectroscopy and force methods to probe the mechanism of translation. To apply these approaches to eukaryotic translation, eukaryotic ribosomes must be specifically labeled with fluorescent labels and molecular handles. Here, we describe preparation and labeling of the small and large yeast ribosomal subunits. Phylogenetically variable hairpin loops in ribosomal RNA are mutated to allow hybridization of oligonucleotides to mutant ribosomes. We demonstrate specific labeling of the ribosomal subunits, and their use in single-molecule fluorescence and force experiments.

RNA interference (RNAi) has been revolutionary for the specific inhibition of gene expression. However, the application of RNAi has been hampered by the fact that many siRNAs induce dose-dependent unwanted secondary effects. Therefore, new methods to increase inhibition of gene expression with low doses of inhibitors are required. We have tested the combination of RNAi and U1i (U1 small nuclear RNA—snRNA—interference). U1i is based on U1 inhibitors (U1in), U1 snRNA molecules modified to target a pre-mRNA and inhibit its gene expression by blocking nuclear polyadenylation. The combination of RNAi and U1i resulted in stronger inhibition of reporter or endogenous genes than that obtained using either of the techniques alone. The increased inhibition observed is stable over time and allows higher inhibition than the best obtained with either of the inhibitors alone even with decreased doses of the inhibitors. We believe that the combination of RNAi and U1i will be of interest when higher inhibition is required or when potent inhibitors are not available. Also, the combination of these techniques would allow functional inhibition with a decreased dose of inhibitors, avoiding toxicity due to dose-dependent unwanted effects.

Continuous efforts have been made to improve next-generation sequencing methods for increased robustness and for applications on low amounts of starting material. We applied double-stranded library protocols for the Roche 454 platform to avoid the yield-reducing steps associated with single-stranded library preparation, and applied a highly sensitive Taqman MGB-probe-based quantitative polymerase chain reaction (qPCR) method. The MGB-probe qPCR, which can detect as low as 100 copies, was used to quantify the amount of effective library, i.e. molecules that form functional clones in emulsion PCR. We also demonstrate that the distribution of library molecules on capture beads follows a Poisson distribution. Combining the qPCR and Poisson statistics, the labour-intensive and costly titration can be eliminated and trace amounts of starting material such as precious clinical samples, transcriptomes of small tissue samples and metagenomics on low biomass environments is applicable.

RNA editing regulates mitochondrial gene expression in trypanosomatid pathogens by creating functional mRNAs. It is catalyzed by a multi-protein complex (the editosome), and is found to be essential in both insect stage and mammalian blood stream form of Trypanosoma brucei. This particular form of RNA editing is unique to trypanosomatids, and thus provides a suitable drug target in trypanosomatid pathogens. Here, we demonstrate the feasibility of a rapid and sensitive fluorescence-based reporter assay to monitor RNA editing based on ribozyme activity. We could validate our new assay using previously identified inhibitors against the essential RNA editing ligase. The principle advantages of this assay are: (i) the use of non-radioactively labeled materials, (ii) sensitivity afforded by fluorescence instrumentation applicable to high-throughput screening of chemical inhibitors against the essential editosome and (iii) a rapid and convenient ‘mix and measure’ type of assay in low volume with a high signal to noise ratio. This assay should enhance rapid identification and characterization of the editosome inhibitors primarily based on the overall composition of the editosomes from T. brucei. These inhibitors could also be tested against the editosomes from the closely related pathogens including T. cruzi and Leishmania species.

Although identification of active motifs in large random sequence pools is central to RNA in vitro selection, no systematic computational equivalent of this process has yet been developed. We develop a computational approach that combines target pool generation, motif scanning and motif screening using secondary structure analysis for applications to 1012–1014-sequence pools; large pool sizes are made possible using program redesign and supercomputing resources. We use the new protocol to search for aptamer and ribozyme motifs in pools up to experimental pool size (1014 sequences). We show that motif scanning, structure matching and flanking sequence analysis, respectively, reduce the initial sequence pool by 6–8, 1–2 and 1 orders of magnitude, consistent with the rare occurrence of active motifs in random pools. The final yields match the theoretical yields from probability theory for simple motifs and overestimate experimental yields, which constitute lower bounds, for aptamers because screening analyses beyond secondary structure information are not considered systematically. We also show that designed pools using our nucleotide transition probability matrices can produce higher yields for RNA ligase motifs than random pools. Our methods for generating, analyzing and designing large pools can help improve RNA design via simulation of aspects of in vitro selection.

MicroRNAs (miRNAs) are naturally occurring small RNAs that regulate the expression of several genes. MiRNAs’ targeting rules are based on sequence complementarity between their mature products and targeted genes’ mRNAs. Based on our present understanding of those rules, we developed an algorithm to design artificial miRNAs to target simultaneously a set of predetermined genes. To validate in silico our algorithm, we tested different sets of genes known to be targeted by a single miRNA. The algorithm finds the seed of the corresponding miRNA among the solutions, which also include the seeds of new artificial miRNA sequences potentially capable of targeting these genes as well. We also validated the functionality of some artificial miRNAs designed to target simultaneously members of the E2F family. These artificial miRNAs reproduced the effects of E2Fs inhibition in both normal human fibroblasts and prostate cancer cells where they inhibited cell proliferation and induced cellular senescence. We conclude that the current miRNA targeting rules based on the seed sequence work to design multiple-target artificial miRNAs. This approach may find applications in both research and therapeutics.

In vitro reconstitution of a biological complex or process normally involves assembly of multiple individually purified protein components. Here we present a strategy that couples expression and assembly of multiple gene products with functional detection in an in vitro reconstituted protein synthesis system. The strategy potentially allows experimental reconstruction of a multi-component biological complex or process using only DNA templates instead of purified proteins. We applied this strategy to bacterial transcription initiation by co-expressing genes encoding Escherichia coli RNA polymerase subunits and sigma factors in the reconstituted protein synthesis system and by coupling the synthesis and assembly of a functional RNA polymerase holoenzyme with the expression of a reporter gene. Using such a system, we demonstrated sigma-factor-dependent, promoter-specific transcription initiation. Since protein synthesis, complex formation and enzyme catalysis occur in the same in vitro reaction mixture, this reconstruction process resembles natural biosynthetic pathways and avoids time-consuming expression and purification of individual proteins. The strategy can significantly reduce the time normally required by conventional reconstitution methods, allow rapid generation and detection of genetic mutations, and provide an open and designable platform for in vitro study and intervention of complex biological processes.

5-Aza-2'-deoxycytidine (5-aza-dC) is a nucleoside analogue with cytotoxic and DNA demethylating effects. Here we show that 5-aza-dC induces the proteasomal degradation of free (non-chromatin bound) DNMT1 through a mechanism which is dependent on DNA synthesis and the targeting of incorporated 5-aza-dC residues by DNMT1 itself. Thus, 5-aza-dC induces Dnmt1 degradation in wild-type mouse ES cells, but not in Dnmt [3a–/–, 3b–/–] mouse ES cells which express Dnmt1 but lack DNA methylation (<0.7% of CpG methylated) and contain few hemi-methylated CpG sites, these being the preferred substrates for Dnmt1. We suggest that adducts formed between DNMT1 and 5-aza-dC molecules in DNA induce a ubiquitin-E3 ligase activity which preferentially targets free DNMT1 molecules for degradation by the proteasome. The proteasome inhibitor MG132 prevents DNMT1 degradation and reduces hypomethylation induced by 5-aza-dC.

Identification of regulatory elements and their target genes is complicated by the fact that regulatory elements can act over large genomic distances. Identification of long-range acting elements is particularly important in the case of disease genes as mutations in these elements can result in human disease. It is becoming increasingly clear that long-range control of gene expression is facilitated by chromatin looping interactions. These interactions can be detected by chromosome conformation capture (3C). Here, we employed 3C as a discovery tool for identification of long-range regulatory elements that control the cystic fibrosis transmembrane conductance regulator gene, CFTR. We identified four elements in a 460-kb region around the locus that loop specifically to the CFTR promoter exclusively in CFTR expressing cells. The elements are located 20 and 80 kb upstream; and 109 and 203 kb downstream of the CFTR promoter. These elements contain DNase I hypersensitive sites and histone modification patterns characteristic of enhancers. The elements also interact with each other and the latter two activate the CFTR promoter synergistically in reporter assays. Our results reveal novel long-range acting elements that control expression of CFTR and suggest that 3C-based approaches can be used for discovery of novel regulatory elements.

G-quadruplexes are nucleic acid secondary structures for which many biological roles have been proposed but whose existence in vivo has remained elusive. To assess their formation, highly specific G-quadruplex ligands are needed. Here, we tested Phen-DC3 and Phen-DC6, two recently released ligands of the bisquinolinium class. In vitro, both compounds exhibit high affinity for the G4 formed by the human minisatellite CEB1 and inhibit efficiently their unwinding by the yeast Pif1 helicase. In vivo, both compounds rapidly induced recombination-dependent rearrangements of CEB1 inserted in the Saccharomyces cerevisiae genome, but did not affect the stability of other tandem repeats lacking G-quadruplex forming sequences. The rearrangements yielded simple-deletion, double-deletion or complex reshuffling of the polymorphic motif units, mimicking the phenotype of the Pif1 inactivation. Treatment of Pif1-deficient cells with the Phen-DC compounds further increased CEB1 instability, revealing additional G4 formation per cell. In sharp contrast, the commonly used N-methyl-mesoporphyrin IX G-quadruplex ligand did not affect CEB1 stability. Altogether, these results demonstrate that the Phen-DC bisquinolinium compounds are potent molecular tools for probing the formation of G-quadruplexes in vivo, interfere with their processing and elucidate their biological roles.

During meiosis there is an imperative to create sufficient crossovers for homologue segregation. This can be achieved during repair of programmed DNA double-strand breaks (DSBs), which are biased towards using a homologue rather than sister chromatid as a repair template. Various proteins contribute to this bias, one of which is a meiosis specific kinase Mek1. It has been proposed that Mek1 establishes the bias by creating a barrier to sister chromatid repair, as distinct from enforcing strand invasion with the homologue. We looked for evidence that Mek1 positively stimulates strand invasion of the homologue. This was done by analysing repair of DSBs induced by the VMA1-derived endonuclease (VDE) and flanked by directly repeated sequences that can be used for intrachromatid single-strand annealing (SSA). SSA competes with interhomologue strand invasion significantly more successfully when Mek1 function is lost. We suggest the increase in intrachromosomal SSA reflects an opportunistic default repair pathway due to loss of a MEK1 stimulated bias for strand invasion of the homologous chromosome. Making use of an inhibitor sensitive mek1-as1 allele, we found that Mek1 function influences the repair pathway throughout the first4–5 h of meiosis. Perhaps reflecting a particular need to create bias for successful interhomologue events before chromosome pairing is complete.

The RecA family of proteins is essential in homologous recombination, a critical step in DNA repair. Here, we report that a rationally-designed small peptide based on the crystal structure of Escherichia coli RecA–DNA complex can promote homologous recombination through the enhancement of both RecA-mediated strand assimilation and three-strand exchange activity. Among 17 peptides tested, peptide #3 with the amino acid sequence of IRFLTARRR has the most potent activity in promoting the RecA-mediated D-loop formation by ~7.2-fold at 37°C. Other peptides such as IRFLTAKKK and IRLLTARRR also have similar, albeit lower, activities. Therefore, hydrophobicity and poly-positive charges, and the space between them in those small peptides are crucial features for such activities. The enhancement of recombination by these peptides appears to be a general phenomenon as similar results were seen by using different plasmids. Remarkably, peptide #3 alone without RecA can also promote the D-loop formation at elevated temperature. Cell viability assays showed that the peptide elevates mammalian cell resistance to two cytotoxic DNA drugs, cisplatin and doxorubicin. The rescue of viability may result from increased DNA repair efficiency. Such peptides may find future biological applications.

DNA primases catalyze the synthesis of the oligoribonucleotides required for the initiation of lagging strand DNA synthesis. Biochemical studies have elucidated the mechanism for the sequence-specific synthesis of primers. However, the physical interactions of the primase with the DNA template to explain the basis of specificity have not been demonstrated. Using a combination of surface plasmon resonance and biochemical assays, we show that T7 DNA primase has only a slightly higher affinity for DNA containing the primase recognition sequence (5'-TGGTC-3') than for DNA lacking the recognition site. However, this binding is drastically enhanced by the presence of the cognate Nucleoside triphosphates (NTPs), Adenosine triphosphate (ATP) and Cytosine triphosphate (CTP) that are incorporated into the primer, pppACCA. Formation of the dimer, pppAC, the initial step of sequence-specific primer synthesis, is not sufficient for the stable binding. Preformed primers exhibit significantly less selective binding than that observed with ATP and CTP. Alterations in subdomains of the primase result in loss of selective DNA binding. We present a model in which conformational changes induced during primer synthesis facilitate contact between the zinc-binding domain and the polymerase domain.

A strong, negative cis-element located at the first intron +502/+835 (I300) of zebrafish myf5 has been reported. To elucidate the molecular mechanism underlying this repression network, we microinjected zebrafish single-cell embryos with I300 RNA, resulting in the dramatic reduction of luciferase activity driven by the myf5 promoter. Within this I300 segment, we identified an intronic microRNA (miR-In300) located at +609/+632 and found that it was more highly expressed in the older mature somites than those newly formed, which negatively correlated with the distribution of zebrafish myf5 transcripts. We proved that miR-In300 suppressed the transcription of myf5 through abolishing myf5 promoter activity, and we subsequently identified the long isoform of the Dickkopf-3 gene (dkk3) as the target gene of miR-In300. We further found that injection of the dkk3-morpholinos (MOs) resulted in downregulation of myf5 transcripts in somites, whereas co-injection of myf5 mRNA with dkk3-MO1 enabled rescue of the defects induced by dkk3-MO1 alone. Finally, injection of miR-In300-MO enhanced both myf5 transcripts in somites and the level of Dkk3 protein in zebrafish embryos. Based on these findings, we concluded that miR-In300 binds to its target gene dkk3, which inhibits the translation of dkk3 mRNA and, in turn, suppresses zebrafish myf5 promoter activity.

Twisted intercalating nucleic acid (TINA) is a novel intercalator and stabilizer of Hoogsteen type parallel triplex formations (PT). Specific design rules for position of TINA in triplex forming oligonucleotides (TFOs) have not previously been presented. We describe a complete collection of easy and robust design rules based upon more than 2500 melting points (Tm) determined by FRET. To increase the sensitivity of PT, multiple TINAs should be placed with at least 3 nt in-between or preferable one TINA for each half helixturn and/or whole helixturn. We find that Tm of base mismatches on PT is remarkably high (between 7.4 and 15.2°C) compared to antiparallel duplexes (between 3.8 and 9.4°C). The specificity of PT by Tm increases when shorter TFOs and higher pH are chosen. To increase Tms, base mismatches should be placed in the center of the TFO and when feasible, A, C or T to G base mismatches should be avoided. Base mismatches can be neutralized by intercalation of a TINA on each side of the base mismatch and masked by a TINA intercalating direct 3' (preferable) or 5' of it. We predict that TINA stabilized PT will improve the sensitivity and specificity of DNA based clinical diagnostic assays.

BLM, one of the human RecQ helicases, plays a fundamental role in homologous recombination-based error-free DNA repair pathways, which require its translocation and DNA unwinding activities. Although translocation is essential in vivo during DNA repair processes and it provides a framework for more complex activities of helicases, including strand separation and nucleoprotein displacement, its mechanism has not been resolved for any human DNA helicase. Here, we present a quantitative model for the translocation of a monomeric form of BLM along ssDNA. We show that BLM performs translocation at a low adenosine triphosphate (ATP) coupling ratio (1 ATP consumed/1 nucleotide traveled) and moderate processivity (with a mean number of 50 nucleotides traveled in a single run). We also show that the rate-limiting step of the translocation cycle is a transition between two ADP-bound enzyme states. Via opening of the helicase core, this structural change may drive the stepping of BLM along the DNA track by a directed inchworm mechanism. The data also support the conclusion that BLM performs double-stranded DNA unwinding by fully active duplex destabilization.

N1-meA and N3-meC are cytotoxic DNA base methylation lesions that can accumulate in the genomes of various organisms in the presence of SN2 type methylating agents. We report here the structural characterization of these base lesions in duplex DNA using a cross-linked protein–DNA crystallization system. The crystal structure of N1-meA:T pair shows an unambiguous Hoogsteen base pair with a syn conformation adopted by N1-meA, which exhibits significant changes in the opening, roll and twist angles as compared to the normal A:T base pair. Unlike N1-meA, N3-meC does not establish any interaction with the opposite G, but remains partially intrahelical. Also, structurally characterized is the N6-meA base modification that forms a normal base pair with the opposite T in duplex DNA. Structural characterization of these base methylation modifications provides molecular level information on how they affect the overall structure of duplex DNA. In addition, the base pairs containing N1-meA or N3-meC do not share any specific characteristic properties except that both lesions create thermodynamically unstable regions in a duplex DNA, a property that may be explored by the repair proteins to locate these lesions.

Human immunodeficiency virus reverse transcriptase (HIV-RT) binds more stably in binary complexes with RNA–DNA versus DNA–DNA. Current results indicate that only the -2 and -4 RNA nucleotides (-1 hybridized to the 3' recessed DNA base) are required for stable binding to RNA–DNA, and even a single RNA nucleotide conferred significantly greater stability than DNA–DNA. Replacing 2'- hydroxyls on pivotal RNA bases with 2'-O-methyls did not affect stability, indicating that interactions between hydroxyls and RT amino acids do not stabilize binding. RT’s Kd (koff/kon) for DNA–DNA and RNA–DNA were similar, although koff differed almost 40-fold, suggesting a faster kon for DNA–DNA. Avian myeloblastosis and Moloney murine leukemia virus RTs also bound more stably to RNA–DNA, but the difference was less pronounced than with HIV-RT. We propose that the H- versus B-form structures of RNA–DNA and DNA–DNA, respectively, allow the former to conform more easily to HIV-RT’s binding cleft, leading to more stable binding. Biologically, the ability of RT to form a more stable complex on RNA–DNA may aid in degradation of RNA fragments that remain after DNA synthesis.

CCA-adding enzymes are polymerases existing in two distinct enzyme classes that both synthesize the C-C-A triplet at tRNA 3'-ends. Class II enzymes (found in bacteria and eukaryotes) carry a flexible loop in their catalytic core required for switching the specificity of the nucleotide binding pocket from CTP- to ATP-recognition. Despite this important function, the loop sequence varies strongly between individual class II CCA-adding enzymes. To investigate whether this loop operates as a discrete functional entity or whether it depends on the sequence context of the enzyme, we introduced reciprocal loop replacements in several enzymes. Surprisingly, many of these replacements are incompatible with enzymatic activity and inhibit ATP-incorporation. A phylogenetic analysis revealed the existence of conserved loop families. Loop replacements within families did not interfere with enzymatic activity, indicating that the loop function depends on a sequence context specific for individual enzyme families. Accordingly, modeling experiments suggest specific interactions of loop positions with important elements of the protein, forming a lever-like structure. Hence, although being part of the enzyme’s catalytic core, the loop region follows an extraordinary evolutionary path, independent of other highly conserved catalytic core elements, but depending on specific sequence features in the context of the individual enzymes.