Occurrence and Function of Hoogsteen Base Pairs in Nucleic Acids

Nucleic acids (DNA and RNA) play essential roles in the central dogma of biology for the storage and transfer of genetic information. The unique chemical and conformational structures of nucleic acids – the double helix composed of complementary Watson-Crick base pairs, provide the structural basis to carry out their biological functions. DNA double helix can dynamically accommodate Watson-Crick and Hoogsteen base-pairing, in which the purine base is flipped by ~180° degrees to adopt syn rather than anti conformation as in Watson-Crick base pairs. There is growing evidence that Hoogsteen base pairs play important roles in DNA replication, recognition, damage or mispair accommodation and repair. Here, we constructed a database for existing Hoogsteen base pairs in DNA duplexes by a structure-based survey from the Protein Data Bank, and structural analyses based on the resulted Hoogsteen structures revealed that Hoogsteen base pairs occur in a wide variety of biological contexts and can induce DNA kinking towards the major groove. As there were documented difficulties in modeling Hoogsteen or Watson-Crick by crystallography, we collaborated with the Richardsons’ lab and identified potential Hoogsteen base pairs that were mis-modeled as Watson-Crick base pairs which suggested that Hoogsteen can be more prevalent than it was thought to be. We developed solution NMR method combined with the site-specific isotope labeling to characterize the formation of, or conformational exchange with Hoogsteen base pairs in large DNA-protein complexes under solution conditions, in the absence of the crystal packing force. We showed that there are enhanced chemical exchange, potentially between Watson-Crick and Hoogsteen, at a sharp kink site in the complex formed by DNA and the Integration Host Factor protein. In stark contrast to B-form DNA, we found that Hoogsteen base pairs are strongly disfavored in A-form RNA duplex. Chemical modifications N1-methyl adenosine and N1-methyl guanosine that block Watson-Crick base-pairing, can be absorbed as Hoogsteen base pairs in DNA, but rather potently destabilized A-form RNA and caused helix melting. The intrinsic instability of Hoogsteen base pairs in A-form RNA endows the N1-methylation as a functioning post-transcriptional modification that was known to facilitate RNA folding, translation and potentially play roles in the epitranscriptome. On the other hand, the dynamic property of DNA that can accommodate Hoogsteen base pairs could be critical to maintaining the genome stability.

Delivery and Scavenging of Nucleic Acids by Polycationic Polymers

Electrostatic interaction is a strong force that attracts positively and negatively charged molecules to each other. Such an interaction is formed between positively charged polycationic polymers and negatively charged nucleic acids. In this dissertation, the electrostatic attraction between polycationic polymers and nucleic acids is exploited for applications in oral gene delivery and nucleic acid scavenging. An enhanced nanoparticle for oral gene delivery of a human Factor IX (hFIX) plasmid is developed using the polycationic polysaccharide, chitosan (Ch), in combination with protamine sulfate (PS) to treat hemophilia B. For nucleic acid scavenging purposes, the development of an effective nucleic acid scavenging nanofiber platform is described for dampening hyper-inflammation and reducing the formation of biofilms.

Non-viral gene therapy may be an attractive alternative to chronic protein replacement therapy. Orally administered non-viral gene vectors have been investigated for more than one decade with little progress made beyond the initial studies. Oral administration has many benefits over intravenous injection including patient compliance and overall cost; however, effective oral gene delivery systems remain elusive. To date, only chitosan carriers have demonstrated successful oral gene delivery due to chitosan’s stability via the oral route. In this study, we increase the transfection efficiency of the chitosan gene carrier by adding protamine sulfate to the nanoparticle formulation. The addition of protamine sulfate to the chitosan nanoparticles results in up to 42x higher in vitro transfection efficiency than chitosan nanoparticles without protamine sulfate. Therapeutic levels of hFIX protein are detected after oral delivery of Ch/PS/phFIX nanoparticles in 5/12 mice in vivo, ranging from 3 -132 ng/mL, as compared to levels below 4 ng/mL in 1/12 mice given Ch/phFIX nanoparticles. These results indicate the protamine sulfate enhances the transfection efficiency of chitosan and should be considered as an effective ternary component for applications in oral gene delivery.

Dying cells release nucleic acids (NA) and NA-complexes that activate the inflammatory pathways of immune cells. Sustained activation of these pathways contributes to chronic inflammation related to autoimmune diseases including systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease. Studies have shown that certain soluble, cationic polymers can scavenge extracellular nucleic acids and inhibit RNA-and DNA-mediated activation of Toll-like receptors (TLRs) and inflammation. In this study, the cationic polymers are incorporated onto insoluble nanofibers, enabling local scavenging of negatively charged pro-inflammatory species such as damage-associated molecular pattern (DAMP) molecules in the extracellular space, reducing cytotoxicity related to unwanted internalization of soluble cationic polymers. In vitro data show that electrospun nanofibers grafted with cationic polymers, termed nucleic acid scavenging nanofibers (NASFs), can scavenge nucleic acid-based agonists of TLR 3 and TLR 9 directly from serum and prevent the production of NF-ĸB, an immune system activating transcription factor while also demonstrating low cytotoxicity. NASFs formed from poly (styrene-alt-maleic anhydride) conjugated with 1.8 kDa branched polyethylenimine (bPEI) resulted in randomly aligned fibers with diameters of 486±9 nm. NASFs effectively eliminate the immune stimulating response of NA based agonists CpG (TLR 9) and poly (I:C) (TLR 3) while not affecting the activation caused by the non-nucleic acid TLR agonist pam3CSK4. Results in a more biologically relevant context of doxorubicin-induced cell death in RAW cells demonstrates that NASFs block ~25-40% of NF-ĸβ response in Ramos-Blue cells treated with RAW extracellular debris, ie DAMPs, following doxorubicin treatment. Together, these data demonstrate that the formation of cationic NASFs by a simple, replicable, modular technique is effective and that such NASFs are capable of modulating localized inflammatory responses.

An understandable way to clinically apply the NASF is as a wound bandage. Chronic wounds are a serious clinical problem that is attributed to an extended period of inflammation as well as the presence of biofilms. An NASF bandage can potentially have two benefits in the treatment of chronic wounds by reducing the inflammation and preventing biofilm formation. NASF can prevent biofilm formation by reducing the NA present in the wound bed, therefore removing large components of what the bacteria use to develop their biofilm matrix, the extracellular polymeric substance, without which the biofilm cannot develop. The NASF described above is used to show the effect of the nucleic acid scavenging technology on in vitro and in vivo biofilm formation of P. aeruginosa, S. aureus, and S. epidermidis biofilms. The in vitro studies demonstrated that the NASFs were able to significantly reduce the biofilm formation in all three bacterial strains. In vivo studies of the NASF on mouse wounds infected with biofilm show that the NASF retain their functionality and are able to scavenge DNA, RNA, and protein from the wound bed. The NASF remove DNA that are maintaining the inflammatory state of the open wound and contributing to the extracellular polymeric substance (EPS), such as mtDNA, and also removing proteins that are required for bacteria/biofilm formation and maintenance such as chaperonin, ribosomal proteins, succinyl CoA-ligase, and polymerases. However, the NASF are not successful at decreasing the wound healing time because their repeated application and removal disrupts the wound bed and removes proteins required for wound healing such as fibronectin, vibronectin, keratin, and plasminogen. Further optimization of NASF treatment duration and potential combination treatments should be tested to reduce the unwanted side effects of increased wound healing time.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

Comparative analyses of seven algorithms for copy number variant identification from single nucleotide polymorphism arrays.

Determination of copy number variants (CNVs) inferred in genome wide single nucleotide polymorphism arrays has shown increasing utility in genetic variant disease associations. Several CNV detection methods are available, but differences in CNV call thresholds and characteristics exist. We evaluated the relative performance of seven methods: circular binary segmentation, CNVFinder, cnvPartition, gain and loss of DNA, Nexus algorithms, PennCNV and QuantiSNP. Tested data included real and simulated Illumina HumHap 550 data from the Singapore cohort study of the risk factors for Myopia (SCORM) and simulated data from Affymetrix 6.0 and platform-independent distributions. The normalized singleton ratio (NSR) is proposed as a metric for parameter optimization before enacting full analysis. We used 10 SCORM samples for optimizing parameter settings for each method and then evaluated method performance at optimal parameters using 100 SCORM samples. The statistical power, false positive rates, and receiver operating characteristic (ROC) curve residuals were evaluated by simulation studies. Optimal parameters, as determined by NSR and ROC curve residuals, were consistent across datasets. QuantiSNP outperformed other methods based on ROC curve residuals over most datasets. Nexus Rank and SNPRank have low specificity and high power. Nexus Rank calls oversized CNVs. PennCNV detects one of the fewest numbers of CNVs.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey.

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage.

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.

Rad18 confers hematopoietic progenitor cell DNA damage tolerance independently of the Fanconi Anemia pathway in vivo.

In cultured cancer cells the E3 ubiquitin ligase Rad18 activates Trans-Lesion Synthesis (TLS) and the Fanconi Anemia (FA) pathway. However, physiological roles of Rad18 in DNA damage tolerance and carcinogenesis are unknown and were investigated here. Primary hematopoietic stem and progenitor cells (HSPC) co-expressed RAD18 and FANCD2 proteins, potentially consistent with a role for Rad18 in FA pathway function during hematopoiesis. However, hematopoietic defects typically associated with fanc-deficiency (decreased HSPC numbers, reduced engraftment potential of HSPC, and Mitomycin C (MMC) -sensitive hematopoiesis), were absent in Rad18(-/-) mice. Moreover, primary Rad18(-/-) mouse embryonic fibroblasts (MEF) retained robust Fancd2 mono-ubiquitination following MMC treatment. Therefore, Rad18 is dispensable for FA pathway activation in untransformed cells and the Rad18 and FA pathways are separable in hematopoietic cells. In contrast with responses to crosslinking agents, Rad18(-/-) HSPC were sensitive to in vivo treatment with the myelosuppressive agent 7,12 Dimethylbenz[a]anthracene (DMBA). Rad18-deficient fibroblasts aberrantly accumulated DNA damage markers after DMBA treatment. Moreover, in vivo DMBA treatment led to increased incidence of B cell malignancy in Rad18(-/-) mice. These results identify novel hematopoietic functions for Rad18 and provide the first demonstration that Rad18 confers DNA damage tolerance and tumor-suppression in a physiological setting.