DNA and RNA sequence determination based on phosphorothioate chemistry

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp, and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

High-throughput sequencing of EcoWI restriction fragments maps the genome-wide landscape of phosphorothioate modification at base resolution

Phosphorothioation (PT), in which a non-bridging oxygen is replaced by a sulfur, is one of the rare modifications discovered in bacteria and archaea that occurs on the sugar-phosphate backbone as opposed to the nucleobase moiety of DNA. While PT modification is widespread in the prokaryotic kingdom, how PT modifications are distributed in the genomes and their exact roles in the cell remain to be defined. In this study, we developed a simple and convenient technique called EcoWI-seq based on a modification-dependent restriction endonuclease to identify genomic positions of PT modifications. EcoWI-seq shows similar performance than other PT modification detection techniques and additionally, is easily scalable while requiring little starting material. As a proof of principle, we applied EcoWI-seq to map the PT modifications at base resolution in the genomes of both the Salmonella enterica cerro 87 and E. coli expressing the dnd+ gene cluster. Specifically, we address whether the partial establishment of modified PT positions is a stochastic or deterministic process. EcoWI-seq reveals a systematic usage of the same subset of target sites in clones for which the PT modification has been independently established.

Large number of bacteria have modified their DNA mainly as part of a strategy to resist virus infection. Most of the modifications are chemical variations on the canonical bases A, T, C or G with phosphorothioate (PT) being a rare exception of a modification that happens on the backbone of the DNA. Interestingly, this PT modification was first chemically synthesized for specific biotechnological processes before scientists discovered that bacteria and archaea naturally perform this modification using their enzymes. The exact roles of phosphorothioation in bacteria and archaea is still under investigation. To enable further investigation of PT modifications, we designed EcoWI-seq, a method to identify the exact positions of these modifications in bacterial genomes. Notably, we applied the EcoWI-seq to several strains of E. coli for which PT modification has been induced by cloning into these strains, the necessary genes for making such modification. We found that these strains, despite being independently made, followed a precise pattern of PT modification with always the same sites being modified. This result indicates a deterministic process for the establishment of PT modification.

Phosphorothioate-Based Site-Specific Labeling of Large RNAs for Structural and Dynamic Studies

Pulsed electron-electron double resonance (PELDOR) spectroscopy, X-ray scattering interferometry (XSI), and single-molecule Förster resonance energy transfer (smFRET) are molecular rulers that provide inter- or intramolecular pair-wise distance distributions in the nanometer range, thus being ideally suitable for structural and dynamic studies of biomolecules including RNAs. The prerequisite for such applications requires site-specific labeling of biomolecules with spin labels, gold nanoparticles, and fluorescent tags, respectively. Recently, site-specific labeling of large RNAs has been achieved by a combination of transcription of an expanded genetic alphabet containing A-T/G-C base pairs and NaM-TPT3 unnatural base pair (UBP) with post-transcriptional modifications at UBP bases by click chemistry or amine-NHS ester reactions. However, due to the bulky sizes of functional groups or labeling probes used, such strategies might cause structural perturbation and decrease the accuracy of distance measurements. Here, we synthesize an α-thiophosphorylated variant of rTPT3TP (rTPT3αS), which allows for post-transcriptional site-specific labeling of large RNAs at the internal α-phosphate backbone via maleimide-modified probes. Subsequent PELDOR, XSI, and smFRET measurements result in narrower distance distributions than labeling at the TPT3 base. The presented strategy provides a new route to empower the molecular rulers for structural and dynamic studies of large RNA and its complex.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

Origin of iodine preferential attack at sulfur in phosphorothioate and subsequent P-O or P-S bond dissociation

Iodine-induced cleavage at phosphorothioate DNA (PT-DNA) is characterized by extremely high sensitivity (∼1 phosphorothioate link per 106 nucleotides), which has been used for detecting and sequencing PT-DNA in bacteria. Despite its foreseeable potential for wide applications, the cleavage mechanism at the PT-modified site has not been well established, and it remains unknown as to whether or not cleavage of the bridging P-O occurs at every PT-modified site. In this work, we conducted accurate ωB97X-D calculations and high-performance liquid chromatography-mass spectrometry to investigate the process of PT-DNA cleavage at the atomic and molecular levels. We have found that iodine chemoselectively binds to the sulfur atom of the phosphorothioate link via a strong halogen-chalcogen interaction (a type of halogen bond, with binding affinity as high as 14.9 kcal/mol) and thus triggers P-O bond cleavage via phosphotriester-like hydrolysis. Additionally, aside from cleavage of the bridging P-O bond, the downstream hydrolyses lead to unwanted P-S/P-O conversions and a loss of the phosphorothioate handle. The mechanism we outline helps to explain specific selectivity at the PT-modified site but also predicts the dynamic stoichiometry of P-S and P-O bond breaking. For instance, Tris is involved in the cascade derivation of S-iodo-phosphorothioate to S-amino-phosphorothioate, suppressing the S-iodo-phosphorothioate hydrolysis to a phosphate diester. However, hydrolysis of one-third of the Tris-O-grafting phosphotriester results in unwanted P-S/P-O conversions. Our study suggests that bacterial DNA phosphorothioation may more frequently occur than previous bioinformatic estimations have predicted from iodine-induced deep sequencing data.

Directed evolution and selection of biostable L-DNA aptamers with a mirror-image DNA polymerase

Mirror-image aptamers made from chirally inverted nucleic acids are nuclease-resistant and exceptionally biostable, opening up opportunities for unique applications. However, the directed evolution and selection of mirror-image aptamers directly from large randomized L-DNA libraries has, to our knowledge, not been demonstrated previously. Here, we developed a 'mirror-image selection' scheme for the directed evolution and selection of biostable L-DNA aptamers with a mirror-image DNA polymerase. We performed iterative rounds of enrichment and mirror-image polymerase chain reaction (PCR) amplification of L-DNA sequences that bind native human thrombin, in conjunction with denaturing gradient gel electrophoresis (DGGE) to isolate individual aptamers and L-DNA sequencing-by-synthesis to determine their sequences. Based on the selected L-DNA aptamers, we designed biostable thrombin sensors and inhibitors, which remained functional in physiologically relevant nuclease-rich environments, even in the presence of human serum that rapidly degraded D-DNA aptamers. Mirror-image selection of biostable L-DNA aptamers directly from large randomized L-DNA libraries greatly expands the range of biomolecules that can be targeted, broadening their applications as biostable sensors, therapeutics and basic research tools.

Bioorthogonal information storage in L-DNA with a high-fidelity mirror-image Pfu DNA polymerase

Natural DNA is exquisitely evolved to store genetic information. The chirally inverted L-DNA, possessing the same informational capacity but resistant to biodegradation, may serve as a robust, bioorthogonal information repository. Here we chemically synthesize a 90-kDa high-fidelity mirror-image Pfu DNA polymerase that enables accurate assembly of a kilobase-sized mirror-image gene. We use the polymerase to encode in L-DNA an 1860 paragraph by Louis Pasteur that first proposed a mirror-image world of biology. We realize chiral steganography by embedding a chimeric D-DNA/L-DNA key molecule in a D-DNA storage library, which conveys a false or secret message depending on the chirality of reading. Furthermore, we show that a trace amount of an L-DNA barcode preserved in water from a local pond remains amplifiable and sequenceable for 1 year, whereas a D-DNA barcode under the same conditions could not be amplified after 1 day. These next-generation mirror-image molecular tools may transform the development of advanced mirror-image biology systems and pave the way for the realization of the mirror-image central dogma and exploration of their applications.

Quantification and mapping of DNA modifications

Apart from the four canonical nucleobases, DNA molecules carry a number of natural modifications. Substantial evidence shows that DNA modifications can regulate diverse biological processes. Dynamic and reversible modifications of DNA are critical for cell differentiation and development. Dysregulation of DNA modifications is closely related to many human diseases. The research of DNA modifications is a rapidly expanding area and has been significantly stimulated by the innovations of analytical methods. With the recent advances in methods and techniques, a series of new DNA modifications have been discovered in the genomes of prokaryotes and eukaryotes. Deciphering the biological roles of DNA modifications depends on the sensitive detection, accurate quantification, and genome-wide mapping of modifications in genomic DNA. This review provides an overview of the recent advances in analytical methods and techniques for both the quantification and genome-wide mapping of natural DNA modifications. We discuss the principles, advantages, and limitations of these developed methods. It is anticipated that new methods and techniques will resolve the current challenges in this burgeoning research field and expedite the elucidation of the functions of DNA modifications.

Site-Specific Scissors Based on Myeloperoxidase for Phosphorothioate DNA

The tool box of site-specific cleavage for nucleic acid has been an increasingly attractive subject. Especially, the recent emergence of the orthogonally activatable DNA device is closely related to the site-specific scission. However, most of these cleavage strategies are based on exogenous assistance, such as laser irradiation. Endogenous strategies are highly desirable for the orthogonally regulatable DNA machine to explore the crucial intracellular biological process and cell signal network. Here, we found that the accurate site-specific cleavage reaction of phosphorothioate (PT) modified DNA by using myeloperoxidase (MPO). A scissors-like mechanism by which MPO breaks PT modification through chloride oxidation has been revealed. Furthermore, we have successfully applied the scissors to activate PT-modified hairpin-DNA machines to produce horseradish peroxidase (HRP)-mimicking DNAzyme or initiate hybridization chain reaction (HCR) amplification. Since MPO plays an important role in the pathway related to oxidative stress in cells, through the HCR amplification activated by this tool box, the oxidative stress in living cells has been robustly imaged. This work proposes an accurate and endogenous site-specific cleavage tool for the research of biostimuli and the construction of DNA molecular devices.

Phosphodiester modifications in mRNA poly(A) tail prevent deadenylation without compromising protein expression

Chemical modifications enable preparation of mRNAs with augmented stability and translational activity. In this study, we explored how chemical modifications of 5',3'-phosphodiester bonds in the mRNA body and poly(A) tail influence the biological properties of eukaryotic mRNA. To obtain modified and unmodified in vitro transcribed mRNAs, we used ATP and ATP analogs modified at the α-phosphate (containing either O-to-S or O-to-BH3 substitutions) and three different RNA polymerases-SP6, T7, and poly(A) polymerase. To verify the efficiency of incorporation of ATP analogs in the presence of ATP, we developed a liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for quantitative assessment of modification frequency based on exhaustive degradation of the transcripts to 5'-mononucleotides. The method also estimated the average poly(A) tail lengths, thereby providing a versatile tool for establishing a structure-biological property relationship for mRNA. We found that mRNAs containing phosphorothioate groups within the poly(A) tail were substantially less susceptible to degradation by 3'-deadenylase than unmodified mRNA and were efficiently expressed in cultured cells, which makes them useful research tools and potential candidates for future development of mRNA-based therapeutics.

Development of Methods Derived from Iodine-Induced Specific Cleavage for Identification and Quantitation of DNA Phosphorothioate Modifications

DNA phosphorothioate (PT) modification is a novel modification that occurs on the DNA backbone, which refers to a non-bridging phosphate oxygen replaced by sulfur. This exclusive DNA modification widely distributes in bacteria but has not been found in eukaryotes to date. PT modification renders DNA nuclease tolerance and serves as a constitute element of bacterial restriction-modification (R-M) defensive system and more biological functions are awaiting exploration. Identification and quantification of the bacterial PT modifications are thus critical to better understanding their biological functions. This work describes three detailed methods derived from iodine-induced specific cleavage-an iodine-induced cleavage assay (ICA), a deep sequencing of iodine-induced cleavage at PT site (ICDS) and an iodine-induced cleavage PT sequencing (PT-IC-Seq)-for the investigation of PT modifications. Using these approaches, we have identified the presence of PT modifications and quantized the frequency of PT modifications in bacteria. These characterizations contributed to the high-resolution genomic mapping of PT modifications, in which the distribution of PT modification sites on the genome was marked accurately and the frequency of the specific modified sites was reliably obtained. Here, we provide time-saving and less labor-consuming methods for both of qualitative and quantitative analysis of genomic PT modifications. The application of these methodologies will offer great potential for better understanding the biology of the PT modifications and open the door to future further systematical study.

Methodologies for detecting environmentally induced DNA damage and repair

Environmental DNA damaging agents continuously challenge the integrity of the genome by introducing a variety of DNA lesions. The DNA damage caused by environmental factors will lead to mutagenesis and subsequent carcinogenesis if they are not removed efficiently by repair pathways. Methods for detection of DNA damage and repair can be applied to identify, visualize, and quantify the DNA damage formation and repair events, and they enable us to illustrate the molecular mechanisms of DNA damage formation, DNA repair pathways, mutagenesis, and carcinogenesis. Ever since the discovery of the double helical structure of DNA in 1953, a great number of methods have been developed to detect various types of DNA damage and repair. Rapid advances in sequencing technologies have facilitated the emergence of a variety of novel methods for detecting environmentally induced DNA damage and repair at the genome-wide scale during the last decade. In this review, we provide a historical overview of the development of various damage detection methods. We also highlight the current methodologies to detect DNA damage and repair, especially some next generation sequencing-based methods.

PTO-QuickStep: A Fast and Efficient Method for Cloning Random Mutagenesis Libraries

QuickStep is a cloning method that allows seamless point integration of a DNA sequence at any position within a target plasmid using only Q5 High-Fidelity DNA Polymerase and DpnI endonuclease. This efficient and cost-effective method consists of two steps: two parallel asymmetric PCRs, followed by a megaprimer-based whole-plasmid amplification. To further simplify the workflow, enhance the efficiency, and increase the uptake of QuickStep, we replaced the asymmetric PCRs with a conventional PCR that uses phosphorothioate (PTO) oligos to generate megaprimers with 3' overhangs. The ease and speed of PTO-QuickStep were demonstrated through (1) right-first-time cloning of a 1.8 kb gene fragment into a pET vector and (2) creating a random mutagenesis library for directed evolution. Unlike most ligation-free random mutagenesis library creation methods (e.g., megaprimer PCR of whole plasmid [MEGAWHOP]), PTO-QuickStep does not require the gene of interest to be precloned into an expression vector to prepare a random mutagenesis library. Therefore, PTO-QuickStep is a simple, reliable, and robust technique, adding to the ever-expanding molecular toolbox of synthetic biology and expediting protein engineering via directed evolution.

Novel Iodine-induced Cleavage Real-time PCR Assay for Accurate Quantification of Phosphorothioate Modified Sites in Bacterial DNA

DNA Phosphorothioate (PT), replacing a non-bridging phosphate oxygen atom with a sulfur atom, is one kind of common DNA modification in bacteria. Whole genome scale description of the location and frequency of PT modification is the key to understand its biological function. Herein we developed a novel method, named with iodine-induced cleavage quantitative real-time PCR (IC-qPCR), to evaluate the frequency of PT modification at a given site in bacterial DNA. The efficiency, dynamic range, sensitivity, reproducibility and accuracy of IC-qPCR were well tested and verified employing an E. coli B7A strain as example. The amplification efficiency of IC-qPCR assay ranged from 91% to 99% with a high correlation coefficient ≥0.99. The limit of quantification was determined as low as 10 copies per reaction for the 607710 and 1818096 sites, and 5 copies for the 302695 and 4120753 sites. Based on the developed IC-qPCR method, the modification frequency of four PTs in E. coli B7A was determined with high accuracy, and the results showed that the PT modification was partial and that the modification frequency varied among investigated PT sites. All these results showed that IC-qPCR was suitable for evaluating the PT modification, which would be helpful to further understand the biological function of PT modification.

DNA phosphorothioate modification-a new multi-functional epigenetic system in bacteria

Synthetic phosphorothioate (PT) internucleotide linkages, in which a nonbridging oxygen is replaced by a sulphur atom, share similar physical and chemical properties with phosphodiesters but confer enhanced nuclease tolerance on DNA/RNA, making PTs a valuable biochemical and pharmacological tool. Interestingly, PT modification was recently found to occur naturally in bacteria in a sequence-selective and RP configuration-specific manner. This oxygen-sulphur swap is catalysed by the gene products of dndABCDE, which constitute a defence barrier with DndFGH in some bacterial strains that can distinguish and attack non-PT-modified foreign DNA, resembling DNA methylation-based restriction-modification (R-M) systems. Despite their similar defensive mechanisms, PT- and methylation-based R-M systems have evolved to target different consensus contexts in the host cell because when they share the same recognition sequences, the protective function of each can be impeded. The redox and nucleophilic properties of PT sulphur render PT modification a versatile player in the maintenance of cellular redox homeostasis, epigenetic regulation and environmental fitness. The widespread presence of dnd systems is considered a consequence of extensive horizontal gene transfer, whereas the lability of PT during oxidative stress and the susceptibility of PT to PT-dependent endonucleases provide possible explanations for the ubiquitous but sporadic distribution of PT modification in the bacterial world.

DNA-encircled lipid bilayers

Lipid bilayers and lipid-associated proteins play crucial roles in biology. As in vivo studies and manipulation are inherently difficult, membrane-mimetic systems are useful for the investigation of lipidic phases, lipid-protein interactions, membrane protein function and membrane structure in vitro. In this work, we describe a route to leverage the programmability of DNA nanotechnology and create DNA-encircled bilayers (DEBs). DEBs are made of multiple copies of an alkylated oligonucleotide hybridized to a single-stranded minicircle, in which up to two alkyl chains per helical turn point to the inside of the toroidal DNA ring. When phospholipids are added, a bilayer is observed to self-assemble within the ring such that the alkyl chains of the oligonucleotides stabilize the hydrophobic rim of the bilayer to prevent formation of vesicles and support thermotropic lipid phase transitions. The DEBs are completely free of protein and can be synthesized from commercially available components using routine equipment. The diameter of DEBs can be varied in a predictable manner. The well-established toolbox from structural DNA nanotechnology, will ultimately enable the rational design of DEBs so that their size, shape or functionalization can be adapted to the specific needs of biophysical investigations of lipidic phases and the properties of membrane proteins embedded into DEB nanoparticle bilayers.

A Novel Saturation Mutagenesis Approach: Single Step Characterization of Regulatory Protein Binding Sites in RNA Using Phosphorothioates

Gene regulation plays an important role in development. Numerous DNA- and RNA-binding proteins bind their target sequences with high specificity to control gene expression. These regulatory proteins control gene expression either at the level of DNA (transcription) or at the level of RNA (pre-mRNA splicing, polyadenylation, mRNA transport, decay, and translation). Identification of regulatory sequences helps understand not only how a gene is switched on or off, but also which downstream genes are regulated by a particular regulatory protein. Here, we describe a one-step approach that allows saturation mutagenesis of a protein binding site in RNA. It involves doping DNA template with non-wild-type nucleotides within the binding site, synthesis of separate RNAs with each phosphorothioate nucleotide, and isolation of the bound fraction following incubation with protein. Interference from non-wild-type nucleotides results in their preferential exclusion from the protein-bound fraction. This is monitored by gel electrophoresis following selective chemical cleavage with iodine of phosphodiester bonds containing phosphorothioates (phosphorothioate mutagenesis or PTM). This single-step saturation mutagenesis approach is applicable to the characterization of any protein binding site in RNA.

Oxidation of phosphorothioate DNA modifications leads to lethal genomic instability

Genomic modification with sulfur as phosphorothioate (PT) is widespread among prokaryotes, including human pathogens. Apart from its physiological functions, the redox and nucleophilic properties of PT sulfur suggest effects on bacterial fitness in stressful environments. Here we show that PTs are dynamic and labile DNA modifications that cause genomic instability during oxidative stress. Using coupled isotopic labeling-mass spectrometry, we observed sulfur replacement in PTs at a rate of ~2%/h in unstressed Escherichia coli and Salmonella enterica. While PT levels were unaffected by exposure to hydrogen peroxide (H₂O₂) or hypochlorous acid (HOCl), PT turnover increased to 3.8–10%/h for HOCl and was unchanged for H₂O₂, consistent with repair of HOCl-induced sulfur damage. PT-dependent HOCl sensitivity extended to cytotoxicity and DNA strand-breaks, which occurred at orders-of-magnitude lower doses of HOCl than H₂O₂. The genotoxicity of HOCl in PT-containing bacteria suggests reduced fitness in competition with HOCl-producing organisms and during human infections.

Decoupling stability and release in disulfide bonds with antibody-small molecule conjugates

Disulfide bonds provide a bioactivatable connection with applications in imaging and therapy. The circulation stability and intracellular release of disulfides are problematically coupled in that increasing stability causes a corresponding decrease in cleavage and payload release. However, an antibody offers the potential for a reversible stabilization. We examined this by attaching a small molecule directly to engineered cysteines in an antibody. At certain sites this unhindered disulfide was stable in circulation yet cellular internalization and antibody catabolism generated a disulfide catabolite that was rapidly reduced. We demonstrated that this stable connection and facile release is applicable to a variety of payloads. The ability to reversibly stabilize a labile functional group with an antibody may offer a way to improve targeted probes and therapeutics.

In vitro analysis of phosphorothioate modification of DNA reveals substrate recognition by a multiprotein complex

A wide variety of prokaryotes possess DNA modifications consisting of sequence-specific phosphorothioates (PT) inserted by members of a five-gene cluster. Recent genome mapping studies revealed two unusual features of PT modifications: short consensus sequences and partial modification of a specific genomic site in a population of bacteria. To better understand the mechanism of target selection of PT modifications that underlies these features, we characterized the substrate recognition of the PT-modifying enzymes termed DptC, D and E in a cell extract system from Salmonella. The results revealed that double-stranded oligodeoxynucleotides underwent de novo PT modification in vitro, with the same modification pattern as in vivo, i. e., G_psAAC/G_psTTC motif. Unexpectedly, in these in vitro analyses we observed no significant effect on PT modification by sequences flanking GAAC/GTTC motif, while PT also occurred in the GAAC/GTTC motif that could not be modified in vivo. Hemi-PT DNA also served as substrate of the PT-modifying enzymes, but not single-stranded DNA. The PT-modifying enzymes were then found to function as a large protein complex, with all of three subunits in tetrameric conformations. This study provided the first demonstration of in vitro DNA PT modification by PT-modifying enzymes that function as a large protein complex.

Genomic mapping of phosphorothioates reveals partial modification of short consensus sequences

Bacterial phosphorothioate (PT) DNA modifications are incorporated by Dnd proteins A-E and often function with DndF-H as a restriction-modification (R-M) system, as in Escherichia coli B7A. However, bacteria such as Vibrio cyclitrophicus FF75 lack dndF-H, which points to other PT functions. Here we report two novel, orthogonal technologies to map PTs across the genomes of B7A and FF75 with >90% agreement: single molecule, real-time sequencing and deep sequencing of iodine-induced cleavage at PT (ICDS). In B7A, we detect PT on both strands of GpsAAC/GpsTTC motifs, but with only 12% of 40,701 possible sites modified. In contrast, PT in FF75 occurs as a single-strand modification at CpsCA, again with only 14% of 160,541 sites modified. Single-molecule analysis indicates that modification could be partial at any particular genomic site even with active restriction by DndF-H, with direct interaction of modification proteins with GAAC/GTTC sites demonstrated with oligonucleotides. These results point to highly unusual target selection by PT-modification proteins and rule out known R-M mechanisms.

Combinatorial engineering of 1-deoxy-D-xylulose 5-phosphate pathway using cross-lapping in vitro assembly (CLIVA) method

The ability to assemble multiple fragments of DNA into a plasmid in a single step is invaluable to studies in metabolic engineering and synthetic biology. Using phosphorothioate chemistry for high efficiency and site specific cleavage of sequences, a novel ligase independent cloning method (cross-lapping in vitro assembly, CLIVA) was systematically and rationally optimized in E. coli. A series of 16 constructs combinatorially expressing genes encoding enzymes in the 1-deoxy-D-xylulose 5-phosphate (DXP) pathway were assembled using multiple DNA modules. A plasmid (21.6 kb) containing 16 pathway genes, was successfully assembled from 7 modules with high efficiency (2.0 x 10(3) cfu/ µg input DNA) within 2 days. Overexpressions of these constructs revealed the unanticipated inhibitory effects of certain combinations of genes on the production of amorphadiene. Interestingly, the inhibitory effects were correlated to the increase in the accumulation of intracellular methylerythritol cyclodiphosphate (MEC), an intermediate metabolite in the DXP pathway. The overexpression of the iron sulfur cluster operon was found to modestly increase the production of amorphadiene. This study demonstrated the utility of CLIVA in the assembly of multiple fragments of DNA into a plasmid which enabled the rapid exploration of biological pathways.

Thio effects and an unconventional metal ion rescue in the genomic hepatitis delta virus ribozyme

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+). Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd(2+). Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg(2+) or Cd(2+), supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ~1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd(2+) of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd(2+) inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg(2+) binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2'-hydroxyl, supporting the experimental findings.

Convergent evolution of adenosine aptamers spanning bacterial, human, and random sequences revealed by structure-based bioinformatics and genomic SELEX

Aptamers are structured macromolecules in vitro evolved to bind molecular targets, whereas in nature they form the ligand-binding domains of riboswitches. Adenosine aptamers of a single structural family were isolated several times from random pools, but they have not been identified in genomic sequences. We used two unbiased methods, structure-based bioinformatics and human genome-based in vitro selection, to identify aptamers that form the same adenosine-binding structure in a bacterium, and several vertebrates, including humans. Two of the human aptamers map to introns of RAB3C and FGD3 genes. The RAB3C aptamer binds ATP with dissociation constants about 10 times lower than physiological ATP concentration, while the minimal FGD3 aptamer binds ATP only cotranscriptionally.

Analysis of catalytic RNA structure and function by nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) is a quick and efficient method to define concurrently, yet singly, the importance of specific functional groups at particular nucleotide residues to the structure and function of an RNA. NAIM can be utilized on virtually any RNA with an assayable function. The method hinges on the ability to successfully incorporate, within an RNA transcript, various 5'-O-(1-thio)nucleoside analogs randomly via in vitro transcription. This could be achieved by using wild-type or Y639F mutant T7 RNA polymerase, thereby creating a pool of analog doped RNAs. The pool when subjected to a selection step to separate the active transcripts from the inactive ones leads to the identification of functional groups that are crucial for RNA activity. The technique can be used to study ribozyme structure and function via monitoring of cleavage or ligation reactions, define functional groups critical for RNA folding, RNA-RNA interactions, and RNA interactions with proteins, metals, or other small molecules. All major classes of catalytic RNAs have been probed by NAIM. This is a generalized approach that should provide the scientific community with the tools to better understand RNA structure-activity relationships.

Analysis of helicase-RNA interactions using nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) is a combinatorial approach that probes individual atoms and functional groups in an RNA molecule and identifies those that are important for a specific biochemical function. Here, we show how NAIM can be adapted to reveal functionally important atoms and groups on RNA substrates of helicases. We explain how NAIM can be used to investigate translocation and unwinding mechanisms of helicases and discuss the advantages and limitations of this powerful chemogenetic approach.

Site-specific labeling of DNA and RNA using an efficiently replicated and transcribed class of unnatural base pairs

Site-specific labeling of enzymatically synthesized DNA or RNA has many potential uses in basic and applied research, ranging from facilitating biophysical studies to the in vitro evolution of functional nucleic acids and the construction of various nanomaterials and biosensors. As part of our efforts to expand the genetic alphabet, we have developed a class of unnatural base pairs, exemplified by d5SICS-dMMO2 and d5SICS-dNaM, which are efficiently replicated and transcribed, and which may be ideal for the site-specific labeling of DNA and RNA. Here, we report the synthesis and analysis of the ribo- and deoxyribo-variants, (d)5SICS and (d)MMO2, modified with free or protected propargylamine linkers that allow for the site-specific modification of DNA or RNA during or after enzymatic synthesis. We also synthesized and evaluated the α-phosphorothioate variant of d5SICSTP, which provides a route to backbone thiolation and an additional strategy for the postamplification site-specific labeling of DNA. The deoxynucleotides were characterized via steady-state kinetics and PCR, while the ribonucleosides were characterized by the transcription of both a short, model RNA as well as full length tRNA. The data reveal that while there are interesting nucleotide and polymerase-specific sensitivities to linker attachment, both (d)MMO2 and (d)5SICS may be used to produce DNA or RNA site-specifically modified with multiple, different functional groups with sufficient efficiency and fidelity for practical applications.

OmniChange: the sequence independent method for simultaneous site-saturation of five codons

Focused mutant library generation methods have been developed to improve mainly "localizable" enzyme properties such as activity and selectivity. Current multi-site saturation methods are restricted by the gene sequence, require subsequent PCR steps and/or additional enzymatic modifications. Here we report, a multiple site saturation mutagenesis method, OmniChange, which simultaneously and efficiently saturates five independent codons. As proof of principle, five chemically cleaved DNA fragments, each carrying one NNK-degenerated codon, were generated and assembled to full gene length in a one-pot-reaction without additional PCR-amplification or use of restriction enzymes or ligases. Sequencing revealed the presence of up to 27 different codons at individual positions, corresponding to 84.4% of the theoretical diversity offered by NNK-degeneration. OmniChange is absolutely sequence independent, does not require a minimal distance between mutated codons and can be accomplished within a day.

Probing essential nucleobase functional groups in aptamers and deoxyribozymes by nucleotide analogue interference mapping of DNA

Nucleotide analogue interference mapping of DNA (dNAIM) is here introduced as a new nonenzymatic interference-based approach that enables high-throughput identification of essential nucleobase functional groups in DNA aptamers and in the catalytic core of deoxyribozymes. Nucleobase-modified ribonucleotides are statistically incorporated into DNA by solid-phase synthesis, employing the 2'-OH group as a chemical tag for analysis of interference effects. This method is exemplified on an AMP-binding DNA aptamer and was further used to identify indispensable nucleobase functional groups for DNA-catalyzed RNA-ligation by the Mg(2+)-dependent deoxyribozymes 7S11 and 9DB1. dNAIM should prove broadly useful for facile structural probing of functional DNA for which active and inactive variants can be separated based on catalytic or ligand-binding activities.

Protection-free one-pot synthesis of 2'-deoxynucleoside 5'-triphosphates and DNA polymerization

By differentiating the functional groups on nucleosides, we have designed and developed a one-pot synthesis of deoxyribonucleoside 5'-triphosphates without any protection on the nucleosides. A facile synthesis is achieved by generating an in situ phosphitylating reagent that reacts selectively with the 5'-hydroxyl groups of the unprotected nucleosides. The synthesized triphosphates are of high quality and can be effectively incorporated into DNAs by DNA polymerase. This novel approach is straightforward and cost-effective for triphosphate synthesis.

Binding site for Xenopus ribosomal protein L5 and accompanying structural changes in 5S rRNA

The structure of the eukaryotic L5-5S rRNA complex was investigated in protection and interference experiments and is compared with the corresponding structure (L18-5S rRNA) in the Haloarcula marismortui 50S subunit. In close correspondence with the archaeal structure, the contact sites for the eukaryotic ribosomal protein are located primarily in helix III and loop C and secondarily in loop A and helix V. While the former is unique to L5, the latter is also a critical contact site for transcription factor IIIA (TFIIIA), accounting for the mutually exclusive binding of these two proteins to 5S RNA. The binding of L5 causes structural changes in loops B and C that expose nucleotides that contact the Xenopus L11 ortholog in H. marismortui. This induced change in the structure of the RNA reveals the origins of the cooperative binding to 5S rRNA that has been observed for the bacterial counterparts of these proteins. The native structure of helix IV and loop D antagonizes binding of L5, indicating that this region of the RNA is dynamic and also influenced by the protein. Examination of the crystal structures of Thermus thermophilus ribosomes in the pre- and post-translocation states identified changes in loop D and in the surrounding region of 23S rRNA that support the proposal that 5S rRNA acts to transmit information between different functional domains of the large subunit.

Twenty years hunting for sulfur in DNA

Here we tell a 20-year long story. It began with an easily overlooked DNA degradation (Dnd) phenomenon during electrophoresis and eventually led to the discovery of an unprecedented DNA sulfur modification governed by five dnd genes. This unusual DNA modification, called phosphorothioation, is the first physiological modification identified on the DNA backbone, in which the nonbridging oxygen is replaced by sulfur in a sequence selective and stereo-specific manner. Homologous dnd gene clusters have been identified in diverse and distantly related bacteria and thus have drawn immediate attention of the entire microbial scientific community. Here, we summarize the progress in chemical, genetic, enzymatic, bioinformatical and analytical aspects of this novel postreplicative DNA modification. We also discuss perspectives on the physiological functions of the DNA phosphorothioate modification in bacteria and their implications.

Effect of cross-link structure on DNA interstrand cross-link repair synthesis

DNA interstrand cross-links (ICLs) are products of chemotherapeutic agents and cellular metabolic processes that block both replication and transcription. If left unrepaired, ICLs are extremely toxic to cells, and ICL repair mechanisms contribute to the survival of certain chemotherapeutic resistance tumors. A critical step in ICL repair involves unhooking the cross-link. In the absence of a homologous donor sequence, the resulting gap can be filled in by a repair synthesis step involving bypass of the cross-link remnant. Here, we examine the effect of cross-link structure on the ability of unhooked DNA substrates to undergo repair synthesis in mammalian whole cell extracts. Using 32P incorporation assays, we found that repair synthesis occurs efficiently past the site of damage when a DNA substrate containing a single N4C-ethyl-N4C cross-link is incubated in HeLa or Chinese hamster ovary cell extracts. This lesion, which can base pair with deoxyguanosine, is readily bypassed by both Escherichia coli DNA polymerase I and T7 DNA polymerase in a primer extension assay. In contrast, bypass was not observed in the primer extension assay or in mammalian cell extracts when DNA substrates containing a N3T-ethyl-N3T or N1I-ethyl-N3T cross-link, whose linkers obstruct the hydrogen bond face of the bases, were used. A modified phosphorothioate sequencing method was used to analyze the ICL repair patches created in the mammalian cell extracts. In the case of the N4C-ethyl-N4C substrate, the repair patch spanned the site of the cross-link, and the lesion was bypassed in an error-free manner. However, although the N3T-ethyl-N3T and N1I-ethyl-N3T substrates were unhooked in the extracts, bypass was not detected. These and our previous results suggest that although the chemical structure of an ICL may not affect initial cross-link unhooking, it can play a significant role in subsequent processing of the cross-link. Understanding how the physical and chemical differences of ICLs affect repair may provide a better understanding of the cytotoxic and mutagenic potential of specific ICLs.

Unusual DNA-DNA cross-links between a photolyase deoxyribozyme, UV1C, and its bound oligonucleotide substrate

UV1C is a photolyase deoxyribozyme that repairs thymine dimers in a DNA oligonucleotide substrate. We report that treatment with iodine generates specific DNA-DNA cross-links between UV1C and a bound substrate analogue, LDPs, in which a single phosphate at the photoreactivation site has been replaced with a phosphorothioate. Although iodine has been reported to generate lysine-cysteine cross-links within a protein, the formation of DNA-DNA cross-links is both unexpected and novel. We have used different mapping procedures to identify a number of bases located in loops of the G-quadruplex fold of UV1C as the sites for cross-linking with LDPs. Mutation of one cross-linking adenine, in particular, leads to a major reduction in UVIC's catalytic activity. A map of these contact cross-linking sites enables us to refine an earlier structural-topological model for the folded UV1C.LDPs complex. The surprising facility with which these novel contact cross-links can be generated between a nucleic acid enzyme and its substrate's reaction site opens up a powerful new approach to mapping the active sites of different ribozymes and deoxyribozymes as well as enabling the dissection of key contacts within RNA-protein complexes.

Multiple metal-binding cores are required for metalloregulation by M-box riboswitch RNAs

Riboswitches are regulatory RNAs that control downstream gene expression in response to direct association with intracellular metabolites or metals. Typically, riboswitch aptamer domains bind to a single small-molecule metabolite. In contrast, an X-ray crystallographic structural model for the M-box riboswitch aptamer revealed the absence of an organic metabolite ligand but the presence of at least six tightly associated magnesiums. This observation agrees well with the proposed role of the M-box riboswitch in functioning as a sensor of intracellular magnesium, although additional nonspecific metal interactions are also undoubtedly required for these purposes. To gain greater functional insight into the metalloregulatory capabilities of M-box RNAs, we sought to determine whether all or a subset of the RNA-chelated magnesium ions were required for riboswitch function. To accomplish this task, each magnesium-binding site was simultaneously yet individually perturbed through random incorporation of phosphorothioate nucleotide analogues, and RNA molecules were investigated for their ability to fold in varying levels of magnesium. These data revealed that all of the magnesium ions observed in the structural model are important for magnesium-dependent tertiary structure formation. Additionally, these functional data revealed a new core of potential metal-binding sites that are likely to assist formation of key tertiary interactions and were previously unobserved in the structural model. It is clear from these data that M-box RNAs require specific binding of a network of metal ions for partial fulfillment of their metalloregulatory functions.

Identification of metabolite-riboswitch interactions using nucleotide analog interference mapping and suppression

Riboswitches are RNA elements capable of modulating gene expression through interaction with cellular metabolites. One member of the riboswitch family, the glmS riboswitch, is unique among riboswitches in that it modulates gene expression by undergoing self-cleavage in the presence of its metabolite glucosamine-6-phosphate (GlcN6P). In order to investigate the interactions between the glmS RNA and GlcN6P we performed nucleotide analog interference mapping (NAIM) and suppression (NAIS). These techniques have been previously used to identify important functional groups in and tertiary contacts necessary for self-splicing and self-cleaving by catalytic RNAs, RNA-protein complexes, RNA folding, and RNA-metal ion interactions. Described here are the details of NAIM and NAIS experiments we have utilized to investigate RNA-ligand interactions between the glmS riboswitch and GlcN6P. These techniques can be employed to study a wide variety of RNA-small molecule interactions.

In vivo analysis of ribonucleoprotein complexes using nucleotide analog interference mapping

Multicomponent RNA-protein complexes are essential for eukaryotic gene expression. Some, like the spliceosome, have been studied successfully in vitro using biochemical and structural approaches, but many have not been reconstituted in cell-free systems. Nucleotide analog interference mapping (NAIM) can report detailed atomic information about requirements for ribonucleoprotein particle assembly and function in living cells, providing a method to study complexes in a cellular context at a level of detail comparable to many biochemical assays. The method relies on incorporation of phosphorothioate-tagged nucleotide analogs during in vitro transcription, followed by a selection for the active population of molecules and analysis of the selected RNA sequence composition. Xenopus oocytes provide a cellular environment for selecting active molecules based on particle assembly or function. Functional group analysis of complexes assembled in vivo provides predictive models for further investigation either in vivo or in vitro as well as benchmarks for evaluating and refining biochemical and structural models.

Nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) is a powerful chemogenetic technique that rapidly identifies chemical groups essential for RNA function. Using a series of phosphorothioate-tagged nucleotide analogs, each carrying different modifications of nucleobase or backbone functionalities, it is possible to simultaneously, yet individually, assess the contribution of particular functional groups to an RNA's activity at every position within the molecule. In contrast to traditional mutagenesis, which modifies RNA on the nucleobase level, the smallest mutable unit in a NAIM analysis is a single atom, providing a detailed description of interactions at critical nucleotides. Because the method introduces modified nucleotides by in vitro transcription, NAIM offers a straightforward and efficient approach to study any RNA that has a selectable function, and it can be applied to RNAs of nearly any length.

Probing RNA structure and function by nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) can be used to simultaneously, yet individually, identify structurally or catalytically important functional groups within an RNA molecule. Phosphorothioate-tagged nucleotides and nucleotide analogs are randomly incorporated into an RNA of interest by in vitro transcription. The phosphorothioate tag marks the site of substitution and identifies sites at which the modification affects the structure or function of the RNA molecule. This technique has been expanded to include identification of hydrogen bonding pairs (NAIS), ionizable functional groups, metal ion ligands, and the energetics of protein binding (QNAIM). The analogs, techniques, and data analysis used in NAIM are described here.

A conserved element that stabilizes the group II intron active site

The internal loop at the base of domain 3 (D3) is one of the most conserved and catalytically important elements of a group II intron. However, the location and molecular nature of its tertiary interaction partners has remained unknown. By employing a combination of site-directed photo-cross-linking and nucleotide analog interference suppression (NAIS), we show that the domain 3 internal loop (D3IL) interacts with the epsilon-epsilon' duplex, which is an active-site element located near the 5'-splice site in D1. Our data also suggest that the D3IL may interact with the bulge of D5, which is a critical active site component. The results of this and other recent studies indicate that the D3IL participates in a complex network of tertiary interactions involving epsilon-epsilon', the bulge of D5 and J23, and that it helps to optimize active site architecture by supporting interactions among these catalytic motifs. Our results are consistent with the role of D3 as a catalytic effector that enhances intron reactivity through active site stabilization.

Structural inference of native and partially folded RNA by high-throughput contact mapping

The biological behaviors of ribozymes, riboswitches, and numerous other functional RNA molecules are critically dependent on their tertiary folding and their ability to sample multiple functional states. The conformational heterogeneity and partially folded nature of most of these states has rendered their characterization by high-resolution structural approaches difficult or even intractable. Here we introduce a method to rapidly infer the tertiary helical arrangements of large RNA molecules in their native and non-native solution states. Multiplexed hydroxyl radical (.OH) cleavage analysis (MOHCA) enables the high-throughput detection of numerous pairs of contacting residues via random incorporation of radical cleavage agents followed by two-dimensional gel electrophoresis. We validated this technology by recapitulating the unfolded and native states of a well studied model RNA, the P4-P6 domain of the Tetrahymena ribozyme, at subhelical resolution. We then applied MOHCA to a recently discovered third state of the P4-P6 RNA that is stabilized by high concentrations of monovalent salt and whose partial order precludes conventional techniques for structure determination. The three-dimensional portrait of a compact, non-native RNA state reveals a well ordered subset of native tertiary contacts, in contrast to the dynamic but otherwise similar molten globule states of proteins. With its applicability to nearly any solution state, we expect MOHCA to be a powerful tool for illuminating the many functional structures of large RNA molecules and RNA/protein complexes.