A Method to Find Sequentially Separated Motifs in Biological Sequences (SSMBS)

Sequence motifs occurring in a particular order in proteins or DNA have been proved to be of biological interest. In this paper, a new method to locate the occurrences of up to five user-defined motifs in a specified order in large proteins and in nucleotide sequence databases is proposed. It has been designed using the concept of quantifiers in regular expressions and linked lists for data storage. The application of this method includes the extraction of relevant consensus regions from biological sequences. This might be useful in clustering of protein families as well as to study the correlation between positions of motifs and their functional sites in DNA sequences.

The 5-Methyl Group in Thymine Dynamically Influences the Structure of A-Tracts in DNA at the Local and Global Level

DNA sequences containing a stretch of several A:T basepairs without a 5'-TA-3' step are known as A-tracts and have been the subject of extensive investigation because of their unique structural features such as a narrow minor groove and their crucial role in several biological processes. One of the aspects under investigation has been the influence of the 5-methyl group of thymine on the properties of A-tracts. Detailed molecular dynamics simulation studies of the sequences d(CGCAAAUUUGCG) and d(CGCAAATTTGCG) indicate that the presence of the 5-methyl group in thymine increases the frequency of a narrow minor groove conformation, which could facilitate its specific recognition by proteins, and reduce its susceptibility to cleavage by DNase I. The bias toward a wider minor groove in the absence of the thymine 5-methyl group is a static structural feature. Our results also indicate that the presence of the thymine 5-methyl group is necessary for calibrating the backbone conformation and the basepair and dinucleotide step geometry of the core A-tract as well as the flanking CA/TG and the neighboring GC/GC steps, as observed in free and protein-bound DNA. As a consequence, it also fine-tunes the curvature of the longer DNA fragment in which the A-tract is embedded.

Interaction between the Z-type DNA duplex and 1,3-propanediamine:crystal structure of d(CACGTG)(2) at 1.2 angstrom resolution

The crystal structure of a hexamer duplex d(CACGTG)(2) has been determined and refined to an R-factor of 18.3% using X-ray data up to 1.2 angstrom resolution. The sequence crystallizes as a left-handed Z-form double helix with Watson-Crick base pairing. There is one hexamer duplex, a spermine molecule, 71 water molecules, and an unexpected diamine (Z-5, 1,3-propanediamine, C3H10N2)) in the asymmetric unit. This is the high-resolution non-disordered structure of a Z-DNA hexamer containing two AT base pairs in the interior of a duplex with no modifications such as bromination or methylation on cytosine bases. This structure does not possess multivalent cations such as cobalt hexaammine that are known to stabilize Z-DNA. The overall duplex structure and its crystal interactions are similar to those of the pure-spermine form of the d(CGCGCG)(2) structure. The spine of hydration in the minor groove is intact except in the vicinity of the T5A8 base pair. The binding of the Z-5 molecule in the minor grove of the d(CACGTG)(2) duplex appears to have a profound effect in conferring stability to a Z-DNA conformation via electrostatic complementarity and hydrogen bonding interactions. The successive base stacking geometry in d(CACGTG)(2) is similar to the corresponding steps in d(CG)(3). These results suggest that specific polyamines such as Z-5 could serve as powerful inducers of Z-type conformation in unmodified DNA sequences with AT base pairs. This structure provides a molecular basis for stabilizing AT base pairs incorporated into an alternating d(CG) sequence.

A novel approach to design of cis-acting DNA structural element for regulation of gene expression in vivo

Taking advantage of the degeneracy of the genetic code we have developed a novel approach to introduce, within a gene, DNA sequences capable of adopting unusual structures and to investigate the role of such sequences in regulation of gene expression in vivo. We used a computer program that generates alternative codon sequences for the same amino-acid sequence to convert a stretch of nucleotides into an inverted-repeat sequence with the potential to adopt cruciform structure. This approach was used to replace a 51-base-pair EcoRI-HindIII segment in the N-terminal region of the beta-galactosidase gene in plasmid pUC19 with a 51-bp synthetic oligonucleotide sequence with the potential to adopt a cruciform structure with 18 bp in the stem region. In selecting the 51-bp sequence, care was taken to include those codons that are preferred in E. coli. E. coli DH5-alpha cells harbouring the plasmid containing the redesigned sequence showed drastic reduction in expression of the beta-galactosidase gene compared to cells harbouring the plasmid with the native sequence. This approach demonstrates the possibility of introducing DNA secondary-structure elements to alter regulation of gene expression in vivo.

Probing of unusual DNA structures in topologically constrained form V DNA: use of restriction enzymes as structural probe

The ability of DNA sequences to adopt unusual structures under the superhelical torsional stress has been studied. Sequences that are forced to adopt unusual conformation in topologically constrained pBR322 form V DNA (Lk=0) were mapped using restriction enzymes as probes. Restriction enzymes such as BamHI, Pstl, Aval and HindIII could not cleave their recognition sequences. The removal of topological constraint relieved this inhibition. The influence of neighbouring sequences on the ability of a given sequence to adopt unusual DNA structure, presumably left handed Z conformation, was studied through single hit analysis. Using multiple cut restriction enzymes such as Narl and Fspl, it could be shown that under identical topological strain, the extent of structural alteration is greatly influenced by the neighbouring sequences. In the light of the variety of sequences and locations that could be mapped to adopt non-6 conformation in pBR322 form V DNA, restriction enzymes appear as potential structural probes for natural DNA sequences.

DNA Binding, Coprotease, and Strand Exchange Activities of Mycobacterial RecA Proteins: Implications for Functional Diversity among RecA Nucleoprotein Filaments

One of the fundamental questions concerning homologous recombination is how RecA or its homologues recognize several DNA sequences with high affinity and catalyze all the diverse biological activities. In this study, we show that the extent of single-stranded DNA binding and strand exchange (SE) promoted by mycobacterial RecA proteins with DNA substrates having various degrees of GC content was comparable with that observed for Escherichia coli RecA. However, the rate and extent of SE promoted by these recombinases showed a strong negative correlation with increasing amounts of sequence divergence embedded at random across the length of the donor strand. Conversely, a positive correlation was seen between SE efficiency and the degree of sequence divergence in the recipient duplex DNA. The extent of heteroduplex formation was not significantly affected when both the pairing partners contained various degrees of sequence divergence, although there was a moderate decrease in the case of mycobacterial RecA proteins with substrates containing larger amounts of sequence divergence. Whereas a high GC content had no discernible effect on E. coli RecA coprotease activity, a negative correlation was apparent between mycobacterial RecA proteins and GC content. We further show clear differences in the extent of SE promoted by E. coli and mycobacterial RecA proteins in the presence of a wide range of ATP:ADP ratios. Taken together, our findings disclose the existence of functional diversity among E. coli and mycobacterial RecA nucleoprotein filaments, and the milieu of sequence divergence (i.e., in the donor or recipient) exerts differential effects on heteroduplex formation, which has implications for the emergence of new genetic variants.

Expression of microinjected foreign DNA in silkworm,bombyx mori

As a prelude to achieving transgenesis in Bombyx mori, conditions have been established for successful microinjection of cloned foreign genes into the silk worm eggs. A sharpened metallic needle is used to pierce the thick chorion layer of the eggsheil, approaching through a droplet of DNA solution deposited on its surface. The microinjection is carried out within 2-2.5 h after oviposition and the injected eggs show 3-5% hatchability and 80-90% survival. Such larvae continuously expressed the microinjected cloned reporter gene, beta-galactosidase, placed under the control of a constitutively expressed cytoplasmic actin A3 gene promoter from B. mori. The expression is seen in different tissues, viz. the fat body, tracheae and the silk glands, till the late larval instars. The microinjected DNA sequences are retained in the adult G(o) moths.

Interaction of G-Quadruplexes with Nonintercalating Duplex-DNA Minor Groove Binding Ligands

The enzyme telomerase synthesizes the G-rich DNA strands of the telomere and its activity is often associated with cancer. The telomerase may be therefore responsible for the ability of a cancer cell-to escape apoptosis. The G-rich DNA sequences often adopt tetra-stranded structure, known as the G-quadruplex DNA (G4-DNA). The stabilization of the telomeric DNA into the G4-DNA structures by small molecules has been the focus of many researchers for the design and development of new anticancer agents. The compounds which stabilize the G-quadruplex in the telomere inhibit the telomerase activity. Besides telomeres, the G4-DNA forming sequences are present in the genomic regions of biological significance including the transcriptional regulatory and promoter regions of several oncogenes. Inducing a G-quadruplex structure within the G-rich promoter sequences is a potential way of achieving selective gene regulation. Several G-quadruplex stabilizing ligands are known. Minor groove binding ligands (MGBLs) interact with the double-helical DNA through the minor grooves sequence-specifically and interfere with several DNA associated processes. These MGBLs when suitably modified switch their preference sometimes from the duplex DNA to G4-DNA and stabilize the G4-DNA as well. Herein, we focus on the recent advances in understanding the G-quadruplex structures, particularly made by the human telomeric ends, and review the results of various investigations of the interaction of designed organic ligands with the G-quadruplex DNA while highlighting the importance of MGBL-G-quadruplex interactions.

Single-Molecule DNA Analysis Reveals That Yeast Hop1 Protein Promotes DNA Folding and Synapsis: Implications for Condensation of Meiotic Chromosomes

During meiosis, long-range interaction between homologous chromosomes is thought to be crucial for homology recognition, exchange of DNA strands, and production of normal haploid gametes. However, little is known about the identity of the proteins involved and the actual molecular mechanism(s) by which chromosomes recognize and recombine with their appropriate homologous partners. Single-molecule analyses have the potential to provide insights into our understanding of this fascinating and long-standing question. Using atomic force microscopy and magnetic tweezers techniques, we discovered that Hop1 protein, a key structural component of Saccharomyces cerevisiae synaptonemal complex, exhibits the ability to bridge noncontiguous DNA segments into intramolecular stem-loop structures in which the DNA segments appear to be fully synapsed within the filamentous protein stems. Additional evidence suggests that Hop1 folds DNA into rigid protein DNA filaments and higher-order nucleoprotein structures. Importantly, Hop1 promotes robust intra- and intermolecular synapsis between double-stranded DNA molecules, suggesting that juxtaposition of DNA sequences may assist in strand exchange between homologues by recombination-associated proteins. Finally, the evidence from ensemble experiments is consistent with the notion that Hop1 causes rigidification of DNA molecules. These results provide the first direct evidence for long-range protein-mediated DNA DNA synapsis, independent of crossover recombination, which is presumed to occur during meiotic recombination.

Microhomology-mediated end joining is the principal mediator of double-strand break repair during mitochondrial DNA lesions

Mitochondrial DNA (mtDNA) deletions are associated with various mitochondrial disorders. The deletions identified in humans are flanked by short, directly repeated mitochondrial DNA sequences; however, the mechanism of such DNA rearrangements has yet to be elucidated. In contrast to nuclear DNA (nDNA), mtDNA is more exposed to oxidative damage, which may result in double-strand breaks (DSBs). Although DSB repair in nDNA is well studied, repair mechanisms in mitochondria are not characterized. In the present study, we investigate the mechanisms of DSB repair in mitochondria using in vitro and ex vivo assays. Whereas classical NHEJ (C-NHEJ) is undetectable, microhomology-mediated alternative NHEJ efficiently repairs DSBs in mitochondria. Of interest, robust microhomology-mediated end joining (MMEJ) was observed with DNA substrates bearing 5-, 8-, 10-, 13-, 16-, 19-, and 22-nt microhomology. Furthermore, MMEJ efficiency was enhanced with an increase in the length of homology. Western blotting, immunoprecipitation, and protein inhibition assays suggest the involvement of CtIP, FEN1, MRE11, and PARP1 in mitochondrial MMEJ. Knock-down studies, in conjunction with other experiments, demonstrated that DNA ligase III, but not ligase IV or ligase I, is primarily responsible for the final sealing of DSBs during mitochondrial MMEJ. These observations highlight the central role of MMEJ in maintenance of mammalian mitochondrial genome integrity and is likely relevant for deletions observed in many human mitochondrial disorders.

Binding site size limitations of imidazole-pyrrole polyamides for recognition in the minor groove of DNA

The discovery that the three ring polyamide Im-Py-Py-Dp containing imidazole (Im) and pyrrole (Py) carboxamides binds the DNA sequence 5'-(A,T)G(A,T)C(A,T)-3' as an antiparallel dimer offers a new model for the design of ligands for specific recognition of sequences in the minor groove containing both G,C and A,T base pairs. In Chapter 2, experiments are described in which the sequential addition of five N- methylpyrrolecarboxamides to the imidazole-pyrrole polyamide Im-Py-Py-Dp affords a series of six homologous polyamides, Im-(Py)_2-7-Dp, that differ in the size of their binding site, apparent first order binding affinity, and sequence specificity. These results demonstrate that DNA sequences up to nine base pairs in length can be specifically recognized by imidazole-pyrrole polyamides containing three to seven rings by 2:1 polyamide-DNA complex formation in the minor groove. Recognition of a nine base pair site defines the new lower limit of the binding site size that can be recognized by polyamides containing exclusively imidazole and pyrrolecarboxamides. The results of this study should provide useful guidelines for the design of new polyamides that bind longer DNA sites with enhanced affinity and specificity.

In Chapter 3 the design and synthesis of the hairpin polyamide Im-Py-Im-Py-γ-Im- Py-Im-Py-Dp is described. Quantitative DNase I footprint titration experiments reveal that Im-Py-Im-Py-γ-Im-Py-Im-Py-Dp binds six base pair 5'-(A,T)GCGC(A,T)-3' sequences with 30-fold higher affinity than the unlinked polyamide Im-Py-Im-Py-Dp. The hairpin polyamide does not discriminate between A•T and T•A at the first and sixth positions of the binding site as three sites 5'-TGCGCT-3', 5'-TGCGCA-3', and 5 'AGCGCT- 3' are bound with similar affinity. However, Im-Py-Im-Py-γ-Im-Py-Im-PyDp is specific for and discriminates between G•C and C•G base pairs in the 5'-GCGC-3' core as evidenced by lower affinities for the mismatched sites 5'-AACGCA-3', 5'- TGCGTT-3', 5'-TGCGGT-3', and 5'-ACCGCT-3'.

In Chapter 4, experiments are described in which a kinetically stable hexa-aza Schiff base La³⁺ complex is covalently attached to a Tat(49-72) peptide which has been shown to bind the HIV-1 TAR RNA sequence. Although these metallo-peptides cleave TAR site-specifically in the hexanucleotide loop to afford products consistent with hydrolysis, a series of control experiments suggests that the observed cleavage is not caused by a sequence-specifically bound Tat(49-72)-La(L)³⁺ peptide.

Sterochemically modified polyamides for recognition in the minor groove of DNA

The design of synthetic molecules that recognize specific sequences of DNA is an ongoing challenge in molecular medicine. Cell-permeable small molecules targeting predetermined DNA sequences offer a potential approach for offsetting the abnormal effects of misregulated gene-expression. Over the past twenty years, Professor Peter B. Dervan has developed a set of pairing rules for the rational design of minor groove binding polyamides containing pyrrole (Py), imidazole (Im), and hydroxypyrrole (Hp). Polyamides have illustrated the capability to permeate cells and inhibit transcription of specific genes in vivo. This provides impetus to identify structural elements that expand the repetoire of polyamide motifs with recognition properties comparable to naturally occurring DNA binding proteins. Through the introduction of chiral amino acids, we have developed chiral polyamides with stereochemically regulated binding characteristics. In addition, chiral substituents have facilitated the development of new polyamide motifs that broaden binding site sizes targetable by this class of ligands.

Structure-, stereochemistry-, and metal-regulated DNA binding/cleaving molecules

In the cell, the binding of proteins to specific sequences of double helical DNA is essential for controlling the processes of protein synthesis (at the level of DNA transcription) and cell proliferation (at the level of DNA replication). In the laboratory, the sequence-specific DNA binding/cleaving properties of restriction endonuclease enzymes (secreted by microorganisms to protect them from foreign DNA molecules) have helped to fuel a revolution in molecular biology. The strength and specificity of a protein:DNA interaction depend upon structural features inherent to the protein and DNA sequences, but it is now appreciated that these features (and therefore protein:DNA complexation) may be altered (regulated) by other protein:DNA complexes, or by environmental factors such as temperature or the presence of specific organic molecules or inorganic ions. It is also now appreciated that molecules much smaller than proteins (including antibiotics of molecular weight less than 2000 and oligonucleotides) can bind to double-helical DNA in sequence-specific fashion. Elucidation of structural motifs and microscopic interactions responsible for the specific molecular recognition of DNA leads to greater understanding of natural processes and provides a basis for the design of novel sequence-specific DNA binding molecules. This thesis describes the synthesis and DNA binding/cleaving characteristics of molecules designed to probe structural, stereochemical, and environmental factors that regulate sequence-specific DNA recognition.

Chapter One introduces the DNA minor groove binding antibiotics Netropsin and Distamycin A, which are di- and tri(N-methylpyrrolecarboxamide) peptides, respectively. The method of DNA affinity cleaving, which has been employed to determine DNA binding properties of designed synthetic molecules is described. The design and synthesis of a series of Netropsin dimers linked in tail-to-tail fashion (by oxalic, malonic, succinic, or fumaric acid), or in head-to-tail fashion (by glycine, β-alanine, and γ-aminobutanoic acid (Gaba)) are presented. These Bis(Netropsin)s were appended with the iron-chelating functionality EDTA in order to make use of the technique of DNA affinity cleaving. Bis(Netropsin)-EDTA compounds are analogs of penta(N-methylpyrrolecarboxamide)-EDTA (P5E), which may be considered a head-to-tail Netropsin dimer linked by Nmethylpyrrolecarboxamide. Low- and high-resolution analysis of pBR322 DNA affinity cleaving by the iron complexes of these molecules indicated that small changes in the length and nature of the linker had significant effects on DNA binding/cleaving efficiency (a measure of DNA binding affinity). DNA binding/cleaving efficiency was found to decrease with changes in the linker in the order β-alanine > succinamide > fumaramide > N-methylpyrrolecarboxamide > malonamide >glycine, γ-aminobutanamide > oxalamide. In general, the Bis(Netropsin)-EDTA:Fe compounds retained the specificity for seven contiguous A:T base pairs characteristic of P5E:Fe binding. However, Bis(Netropsin)Oxalamide- EDTA:Fe exhibited decreased specificity for A:T base pairs, and Bis(Netropsin)-Gaba-EDT A:Fe exhibited some DNA binding sites of less than seven base pairs. Bis(Netropsin)s linked with diacids have C₂-symmmetrical DNA binding subunits and exhibited little DNA binding orientation preference. Bis(Netropsin)s linked with amino acids lack C₂-symmetrical DNA binding subunits and exhibited higher orientation preferences. A model for the high DNA binding orientation preferences observed with head-to-tail DNA minor groove binding molecules is presented.

Chapter Two describes the design, synthesis, and DNA binding properties of a series of chiral molecules: Bis(Netropsin)-EDTA compounds with linkers derived from (R,R)-, (S,S)-, and (RS,SR)-tartaric acids, (R,R)-, (S,S)-, and (RS,SR)-tartaric acid acetonides, (R)- and (S)-malic acids, N ,N-dimethylaminoaspartic acid, and (R)- and (S)-alanine, as well as three constitutional isomers in which an N-methylpyrrolecarboxamide (P1) subunit and a tri(N-methylpyrrolecarboxamide)-EDTA (P3-EDTA) subunit were linked by succinic acid, (R ,R)-, and (S ,S)-tartaric acids. DNA binding/cleaving efficiencies among this series of molecules and the Bis(Netropsin)s described in Chapter One were found to decrease with changes in the linker in the order β-alanine > succinamide > P1-succinamide-P3 > fumaramide > (S)-malicamide > N-methylpyrrolecarboxamide > (R)-malicamide > malonamide > N ,N-dimethylaminoaspanamide > glycine = Gaba = (S,S)-tartaramide = P1-(S,S)-tanaramide-P3 > oxalamide > (RS,SR)-tartaramide = P1- (R,R)-tanaramide-P3 > (R,R)-tartaramide (no sequence-specific DNA binding was detected for Bis(Netropsin)s linked by (R)- or (S)-alanine or by tartaric acid acetonides). The chiral molecules retained DNA binding specificity for seven contiguous A:T base pairs. From the DNA affinity cleaving data it could be determined that: 1) Addition of one or two substituents to the linker of Bis(Netropsin)-Succinamide resulted in stepwise decreases in DNA binding affinity; 2) molecules with single hydroxyl substituents bound DNA more strongly than molecules with single dimethylamino substituents; 3) hydroxyl-substituted molecules of (S) configuration bound more strongly to DNA than molecules of (R) configuration. This stereochemical regulation of DNA binding is proposed to arise from the inherent right-handed twist of (S)-enantiomeric Bis(Netropsin)s versus the inherent lefthanded twist of (R)-enantiomeric Bis(Netropsin)s, which makes the (S)-enantiomers more complementary to the right-handed twist of B form DNA.

Chapter Three describes the design and synthesis of molecules for the study of metalloregulated DNA binding phenomena. Among a series of Bis(Netropsin)-EDTA compounds linked by homologous tethers bearing four, five, or six oxygen atoms, the Bis(Netropsin) linked by a pentaether tether exhibited strongly enhanced DNA binding/cleaving in the presence of strontium or barium cations. The observed metallospecificity was consistent with the known affinities of metal cations for the cyclic hexaether 18-crown-6 in water. High-resolution DNA affinity cleaving analysis indicated that DNA binding by this molecule in the presence of strontium or barium was not only stronger but of different sequence-specificity than the (weak) binding observed in the absence of metal cations. The metalloregulated binding sites were consistent with A:T binding by the Netropsin subunits and G:C binding by a strontium or barium:pentaether complex. A model for the observed positive metalloregulation and novel sequence-specificity is presented. The effects of 44 different cations on DNA affinity cleaving by P5E:Fe were examined. A series of Bis(Netropsin)-EDTA compounds linked by tethers bearing two, three, four, or five amino groups was also synthesized. These molecules exhibited strong and specific binding to A:T rich regions of DNA. It was found that the iron complexes of these molecules bound and cleaved DNA most efficiently at pH 6.0-6.5, while P5E:Fe bound and cleaved most efficiently at pH 7.5-8.0. Incubating the Bis(Netropsin) Polyamine-EDTA:Fe molecules with K₂PdCl₄ abolished their DNA binding/cleaving activity. It is proposed that the observed negative metalloregulation arises from kinetically inert Bis(Netropsin) Polyamine:Pd(II) complexes or aggregates, which are sterically unsuitable for DNA complexation. Finally, attempts to produce a synthetic metalloregulated DNA binding protein are described. For this study, five derivatives of a synthetic 52 amino acid residue DNA binding/cleaving protein were produced. The synthetic mutant proteins carried a novel pentaether ionophoric amino acid residue at different positions within the primary sequence. The proteins did not exhibit significant DNA binding/cleaving activity, but they served to illustrate the potential for introducing novel amino acid residues within DNA binding protein sequences, and for the development of the tricyclohexyl ester of EDTA as a superior reagent for the introduction of EDT A into synthetic proteins.

Chapter Four describes the discovery and characterization of a new DNA binding/cleaving agent, [SalenMn(III)]OAc. This metal complex produces single- and double-strand cleavage of DNA, with specificity for A:T rich regions, in the presence of oxygen atom donors such as iodosyl benzene, hydrogen peroxide, or peracids. Maximal cleavage by [SalenMn(III)]OAc was produced at pH 6-7. A comparison of DNA singleand double-strand cleavage by [SalenMn(III)]⁺ and other small molecules (Methidiumpropyl-EDTA:Fe, Distamycin-EDTA:Fe, Neocarzinostatin, Bleomycin:Fe) is presented. It was found that DNA cleavage by [SalenMn(III)]⁺ did not require the presence of dioxygen, and that base treatment of DNA subsequent to cleavage by [SalenMn(III)]⁺ afforded greater cleavage and alterations in the cleavage patterns. Analysis of DNA products formed upon DNA cleavage by [SalenMn(III)] indicated that cleavage was due to oxidation of the sugar-phosphate backbone of DNA. Several mechanisms consistent with the observed products and reaction requirements are discussed.

Chapter Five describes progress on some additional studies. In one study, the DNA binding/cleaving specificities of Distamycin-EDTA derivatives bearing pyrrole N-isopropyl substituents were found to be the same as those of derivatives bearing pyrrole N-methyl substituents. In a second study, the design of and synthetic progress towards a series of nucleopeptide activators of transcription are presented. Five synthetic plasmids designed to test for activation of in vitro run-off transcription by DNA triple helix-forming oligonucleotides or nucleopeptides are described.

Chapter Six contains the experimental documentation of the thesis work.

DNA-mediated oxidation of transcription factor p53

Transcription factor p53 is the most commonly altered gene in human cancer. As a redox-active protein in direct contact with DNA, p53 can directly sense oxidative stress through DNA-mediated charge transport. Electron hole transport occurs with a shallow distance dependence over long distances through the π-stacked DNA bases, leading to the oxidation and dissociation of DNA-bound p53. The extent of p53 dissociation depends upon the redox potential of the response element DNA in direct contact with each p53 monomer. The DNA sequence dependence of p53 oxidative dissociation was examined by electrophoretic mobility shift assays using radiolabeled oligonucleotides containing both synthetic and human p53 response elements with an appended anthraquinone photooxidant. Greater p53 dissociation is observed from DNA sequences containing low redox potential purine regions, particularly guanine triplets, within the p53 response element. Using denaturing polyacrylamide gel electrophoresis of irradiated anthraquinone-modified DNA, the DNA damage sites, which correspond to locations of preferred electron hole localization, were determined. The resulting DNA damage preferentially localizes to guanine doublets and triplets within the response element. Oxidative DNA damage is inhibited in the presence of p53, however, only at DNA sites within the response element, and therefore in direct contact with p53. From these data, predictions about the sensitivity of human p53-binding sites to oxidative stress, as well as possible biological implications, have been made. On the basis of our data, the guanine pattern within the purine region of each p53-binding site determines the response of p53 to DNA-mediated oxidation, yielding for some sequences the oxidative dissociation of p53 from a distance and thereby providing another potential role for DNA charge transport chemistry within the cell.

To determine whether the change in p53 response element occupancy observed in vitro also correlates in cellulo, chromatin immunoprecipition (ChIP) and quantitative PCR (qPCR) were used to directly quantify p53 binding to certain response elements in HCT116N cells. The HCT116N cells containing a wild type p53 were treated with the photooxidant [Rh(phi)2bpy]³⁺, Nutlin-3 to upregulate p53, and subsequently irradiated to induce oxidative genomic stress. To covalently tether p53 interacting with DNA, the cells were fixed with disuccinimidyl glutarate and formaldehyde. The nuclei of the harvested cells were isolated, sonicated, and immunoprecipitated using magnetic beads conjugated with a monoclonal p53 antibody. The purified immounoprecipiated DNA was then quantified via qPCR and genomic sequencing. Overall, the ChIP results were significantly varied over ten experimental trials, but one trend is observed overall: greater variation of p53 occupancy is observed in response elements from which oxidative dissociation would be expected, while significantly less change in p53 occupancy occurs for response elements from which oxidative dissociation would not be anticipated.

The chemical oxidation of transcription factor p53 via DNA CT was also investigated with respect to the protein at the amino acid level. Transcription factor p53 plays a critical role in the cellular response to stress stimuli, which may be modulated through the redox modulation of conserved cysteine residues within the DNA-binding domain. Residues within p53 that enable oxidative dissociation are herein investigated. Of the 8 mutants studied by electrophoretic mobility shift assay (EMSA), only the C275S mutation significantly decreased the protein affinity (KD) for the Gadd45 response element. EMSA assays of p53 oxidative dissociation promoted by photoexcitation of anthraquinone-tethered Gadd45 oligonucleotides were used to determine the influence of p53 mutations on oxidative dissociation; mutation to C275S severely attenuates oxidative dissociation while C277S substantially attenuates dissociation. Differential thiol labeling was used to determine the oxidation states of cysteine residues within p53 after DNA-mediated oxidation. Reduced cysteines were iodoacetamide labeled, while oxidized cysteines participating in disulfide bonds were ¹³C₂D₂-iodoacetamide labeled. Intensities of respective iodoacetamide-modified peptide fragments were analyzed using a QTRAP 6500 LC-MS/MS system, quantified with Skyline, and directly compared. A distinct shift in peptide labeling toward ¹³C₂D₂-iodoacetamide labeled cysteines is observed in oxidized samples as compared to the respective controls. All of the observable cysteine residues trend toward the heavy label under conditions of DNA CT, indicating the formation of multiple disulfide bonds potentially among the C124, C135, C141, C182, C275, and C277. Based on these data it is proposed that disulfide formation involving C275 is critical for inducing oxidative dissociation of p53 from DNA.

Programming chemical kinetics: engineering dynamic reaction networks with DNA strand displacement

Over the last century, the silicon revolution has enabled us to build faster, smaller and more sophisticated computers. Today, these computers control phones, cars, satellites, assembly lines, and other electromechanical devices. Just as electrical wiring controls electromechanical devices, living organisms employ "chemical wiring" to make decisions about their environment and control physical processes. Currently, the big difference between these two substrates is that while we have the abstractions, design principles, verification and fabrication techniques in place for programming with silicon, we have no comparable understanding or expertise for programming chemistry.

In this thesis we take a small step towards the goal of learning how to systematically engineer prescribed non-equilibrium dynamical behaviors in chemical systems. We use the formalism of chemical reaction networks (CRNs), combined with mass-action kinetics, as our programming language for specifying dynamical behaviors. Leveraging the tools of nucleic acid nanotechnology (introduced in Chapter 1), we employ synthetic DNA molecules as our molecular architecture and toehold-mediated DNA strand displacement as our reaction primitive.

Abstraction, modular design and systematic fabrication can work only with well-understood and quantitatively characterized tools. Therefore, we embark on a detailed study of the "device physics" of DNA strand displacement (Chapter 2). We present a unified view of strand displacement biophysics and kinetics by studying the process at multiple levels of detail, using an intuitive model of a random walk on a 1-dimensional energy landscape, a secondary structure kinetics model with single base-pair steps, and a coarse-grained molecular model that incorporates three-dimensional geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Our findings are consistent with previously measured or inferred rates for hybridization, fraying, and branch migration, and provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

In Chapters 3 and 4, we identify and overcome the crucial experimental challenges involved in using our general DNA-based technology for engineering dynamical behaviors in the test tube. In this process, we identify important design rules that inform our choice of molecular motifs and our algorithms for designing and verifying DNA sequences for our molecular implementation. We also develop flexible molecular strategies for "tuning" our reaction rates and stoichiometries in order to compensate for unavoidable non-idealities in the molecular implementation, such as imperfectly synthesized molecules and spurious "leak" pathways that compete with desired pathways.

We successfully implement three distinct autocatalytic reactions, which we then combine into a de novo chemical oscillator. Unlike biological networks, which use sophisticated evolved molecules (like proteins) to realize such behavior, our test tube realization is the first to demonstrate that Watson-Crick base pairing interactions alone suffice for oscillatory dynamics. Since our design pipeline is general and applicable to any CRN, our experimental demonstration of a de novo chemical oscillator could enable the systematic construction of CRNs with other dynamic behaviors.