Random mutagenesis of Thermus aquaticus DNA polymerase I: concordance of immutable sites in vivo with the crystal structure.

Expression of Thermus aquaticus (Taq) DNA polymerase I (pol I) in Escherichia, coli complements the growth defect caused by a temperature-sensitive mutation in the host pol I. We replaced the nucleotide sequence encoding amino acids 659-671 of the O-helix of Taq DNA pol I, corresponding to the substrate binding site, with an oligonucleotide containing random nucleotides. Functional Taq pol I mutants were selected based on colony formation at the nonpermissive temperature. By using a library with 9% random substitutions at each of 39 positions, we identified 61 active Taq pol I mutants, each of which contained from one to four amino acid substitutions. Some amino acids, such as alanine-661 and threonine-664, were tolerant of several or even many diverse replacements. In contrast, no replacements or only conservative replacements were identified at arginine-659, lysine-663, and tyrosine-671. By using a library with totally random nucleotides at five different codons (arginine-659, arginine-660, lysine-663, phenylalanine-667, and glycine-668), we confirmed that arginine-659 and lysine-663 were immutable, and observed that only tyrosine substituted for phenylalanine-667. The two immutable residues and the two residues that tolerate only highly conservative replacements lie on the side of O-helix facing the incoming deoxynucleoside triphosphate, as determined by x-ray analysis. Thus, we offer a new approach to assess concordance of the active conformation of an enzyme, as interpreted from the crystal structure, with the active conformation inferred from in vivo function.

Prediction of the stability of DNA triplexes.

We present rules that allow one to predict the stability of DNA pyrimidine.purine.pyrimidine (Y.R.Y) triple helices on the basis of the sequence. The rules were derived from van't Hoff analysis of 23 oligonucleotide triplexes tested at a variety of pH values. To predict the enthalpy of triplex formation (delta H degrees), a simple nearest-neighbor model was found to be sufficient. However, to accurately predict the free energy of the triplex (delta G degrees), a combination model consisting of five parameters was needed. These parameters were (i) the delta G degrees for helix initiation, (ii) the delta G degrees for adding a T-A.T triple, (iii) the delta G degrees for adding a C(+)-G.C triple, (iv) the penalty for adjacent C bases, and (v) the pH dependence of the C(+)-G.C triple's stability. The fitted parameters are highly consistent with thermodynamic data from the basis set, generally predicting both delta H degrees and delta G degrees to within the experimental error. Examination of the parameters points out several interesting features. The combination model predicts that C(+) -G.C. triples are much more stabilizing than T-A.T triples below pH 7.0 and that the stability of the former increases approximately equal to 1 kcal/mol per pH unit as the pH is decreased. Surprisingly though, the most stable sequence is predicted to be a CT repeat, as adjacent C bases partially cancel the stability of one another. The parameters successfully predict tm values from other laboratories, with some interesting exceptions.

Sequence specificity of triplex DNA formation: Analysis by a combinatorial approach, restriction endonuclease protection selection and amplification.

We have devised a combinatorial method, restriction endonuclease protection selection and amplification (REPSA), to identify consensus ligand binding sequences in DNA. In this technique, cleavage by a type IIS restriction endonuclease (an enzyme that cleaves DNA at a site distal from its recognition sequence) is prevented by a bound ligand while unbound DNA is cleaved. Since the selection step of REPSA is performed in solution under mild conditions, this approach is amenable to the investigation of ligand-DNA complexes that are either insufficiently stable or not readily separable by other methods. Here we report the use of REPSA to identify the consensus duplex DNA sequence recognized by a G/T-rich oligodeoxyribonucleotide under conditions favoring purine-motif triple-helix formation. Analysis of 47 sequences indicated that recognition between 13 bases on the oligonucleotide 3' end and the duplex DNA was sufficient for triplex formation and indicated the possible existence of a new base triplet, G.AT. This information should help identify appropriate target sequences for purine-motif triplex formation and demonstrates the power of REPSA for investigating ligand-DNA interactions.

In vivo protein-DNA interactions at human DNA replication origin.

Protein-DNA interactions were studied in vivo at the region containing a human DNA replication origin, located at the 3' end of the lamin B2 gene and partially overlapping the promoter of another gene, located downstream. DNase I treatment of nuclei isolated from both exponentially growing and nonproliferating HL-60 cells showed that this region has an altered, highly accessible, chromatin structure. High-resolution analysis of protein-DNA interactions in a 600-bp area encompassing the origin was carried out by the in vivo footprinting technique based on the ligation-mediated polymerase chain reaction. In growing HL-60 cells, footprints at sequences homologous to binding sites for known transcription factors (members of the basic-helix-loop-helix family, nuclear respiratory factor 1, transcription factor Sp1, and upstream binding factor) were detected in the region corresponding to the promoter of the downstream gene. Upon conversion of cells to a nonproliferative state, a reduction in the intensity of these footprints was observed that paralleled the diminished transcriptional activity of the genomic area. In addition to these protections, in close correspondence to the replication initiation site, a prominent footprint was detected that extended over 70 nucleotides on one strand only. This footprint was absent from nonproliferating HL-60 cells, indicating that this specific protein-DNA interaction might be involved in the process of origin activation.

Antiproliferative properties of the USF family of helix-loop-helix transcription factors.

USF is a family of transcription factors characterized by a highly conserved basic-helix-loop-helix-leucine zipper (bHLH-zip) DNA-binding domain. Two different USF genes, termed USF1 and USF2, are ubiquitously expressed in both humans and mice. The USF1 and USF2 proteins contain highly divergent transcriptional activation domains but share extensive homologies in the bHLH-zip region and recognize the same CACGTG DNA motifs. Although the DNA-binding and transcriptional activities of these proteins have been characterized, the biological function of USF is not well understood. Here, focus- and colony-formation assays were used to investigate the potential involvement of USF in the regulation of cellular transformation and proliferation. Both USF1 and USF2 inhibited the transformation of rat embryo fibroblasts mediated by Ras and c-Myc, a bHLH-zip transcription factor that also binds CACGTG motifs. DNA binding was required but not fully sufficient for inhibition of Myc-dependent transformation by USF, since deletion mutants containing only the DNA-binding domains of USF1 or USF2 produced partial inhibition. While the effect of USF1 was selective for Myc-dependent transformation, wild-type USF2 exerted in addition a strong inhibition of E1A-mediated transformation and a strong suppression of HeLa cell colony formation. These results suggest that members of the USF family may serve as negative regulators of cellular proliferation in two ways, one by antagonizing the transforming function of Myc, the other through a more general growth-inhibitory effect.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.

A single-stranded DNA binding protein binds the origin of replication of the duplex kinetoplast DNA.

Replication of the kinetoplast DNA (kDNA) minicircle of trypanosomatids initiates at a conserved 12-nt sequence, 5'-GGGGTTGGTGTA-3', termed the universal minicircle sequence (UMS). A sequence-specific single-stranded DNA-binding protein from Crithidia fasciculata binds the heavy strand of the 12-mer UMS. Whereas this UMS-binding protein (UMSBP) does not bind a duplex UMS dodecamer, it binds the double-stranded kDNA minicircle as well as a duplex minicircle fragment containing the origin-associated UMS. Binding of the minicircle origin region by the single-stranded DNA binding protein suggested the local unwinding of the DNA double helix at this site. Modification of thymine residues at this site by KMnO4 revealed that the UMS resides within an unwound or otherwise sharply distorted DNA at the minicircle origin region. Computer analysis predicts the sequence-directed curving of the minicircle origin region. Electrophoresis of a minicircle fragment containing the origin region in polyacrylamide gels revealed a significantly lower electrophoretic mobility than expected from its length. The fragment anomalous electrophoretic mobility is displayed only in its native conformation and is dependent on temperature and gel porosity, indicating the local curving of the DNA double helix. We suggest that binding of UMSBP at the minicircle origin of replication is possible through local unwinding of the DNA double helix at the UMS site. It is hypothesized here that this local melting is initiated through the untwisting of unstacked dinucleotide sequences at the bent origin site.

Meso1, a basic-helix-loop-helix protein involved in mammalian presomitic mesoderm development.

To identify genes involved in the regulation of early mammalian development, we have developed a dominant-negative mutant basic-helix-loop-helix (bHLH) protein probe for interaction cloning and have isolated a member of the bHLH family of transcription factors, Meso1. Meso1-E2A heterodimers are capable of binding to oligonucleotide probes that contain a bHLH DNA recognition motif. In mouse embryos, Meso1 is expressed prior to MyoD1 family members. Meso1 expression is first detected at the neural plate stage of development in the paraxial mesoderm of the head and in presomitic mesodermal cells prior to their condensation into somites. Our findings suggest that Meso1 may be a key regulatory gene involved in the early events of vertebrate mesoderm differentiation.

Synthesis of a thymidyl pentamer of deoxyribonucleic guanidine and binding studies with DNA homopolynucleotides.

Replacement of the phosphodiester linkages of the polyanion DNA with guanidine linkers provides the polycation deoxynucleic guanidine (DNG). The synthesis of pentameric thymidyl DNA is provided. This polycationic DNG species binds with unprecedented affinity and with base-pair specificity to negatively charged poly(dA) to provide both double and triple helices. The dramatic stability of these hybrid structures is shown by their denaturation temperatures (Tm). For example, the double helix of the pentameric thymidyl DNG and poly(dA) does not dissociate in boiling water (ionic strength = 0.12), whereas the Tm for pentameric thymidyl DNA associated with poly(dA) is approximately 13 degrees C (ionic strength = 0.12). The effect of ionic strength on Tm for DNG complexes with DNA shows an opposite correlation compared with double-stranded DNA and is much more dramatic than for double-stranded DNA.

Transplantation of a 17-amino acid alpha-helical DNA-binding domain into an antibody molecule confers sequence-dependent DNA recognition.

Recombinant antibodies capable of sequence-specific interactions with nucleic acids represent a class of DNA- and RNA-binding proteins with potential for broad application in basic research and medicine. We describe the rational design of a DNA-binding antibody, Fab-Ebox, by replacing a variable segment of the immunoglobulin heavy chain with a 17-amino acid domain derived from TFEB, a class B basic helix-loop-helix protein. DNA-binding activity was studied by electrophoretic mobility-shift assays in which Fab-Ebox was shown to form a specific complex with DNA containing the TFEB recognition motif (CACGTG). Similarities were found in the abilities of TFEB and Fab-Ebox to discriminate between oligodeoxyribonucleotides containing altered recognition sequences. Comparable interference of binding by methylation of cytosine residues indicated that Fab-Ebox and TFEB both contact DNA through interactions along the major groove of double-stranded DNA. The results of this study indicate that DNA-binding antibodies of high specificity can be developed by using the modular nature of both immunoglobulins and transcription factors.

Biochemical characterization and DNA binding properties of the MocR family member GabR, a PLP-dependent transcription factor

GabR è un fattore di trascrizione chimerico appartenente alla famiglia dei MocR/GabR, costituito da un dominio N-terminale elica-giro-elica di legame al DNA e un dominio effettore e/o di oligomerizzazione al C-terminale. I due domini sono connessi da un linker flessibile di 29 aminoacidi. Il dominio C-terminale è strutturalmente omologo agli enzimi aminotransferasici fold-type I, i quali, utilizzando il piridossal-5’-fosfato (PLP) come cofattore, sono direttamente coinvolti nel metabolismo degli aminoacidi. L’interazione contemporanea di PLP e acido γ-aminobutirrico (GABA) a GabR fa sì che questa promuova la trascrizione di due geni, gabT e gabD, implicati nel metabolismo del GABA. GabR cristallizza come un omodimero con una configurazione testa-coda. Il legame con la regione promotrice gabTD avviene attraverso il riconoscimento specifico di due sequenze dirette e ripetute (ATACCA), separate da uno spacer di 34 bp. In questo studio sono state indagate le proprietà biochimiche, strutturali e di legame al DNA della proteina GabR di Bacillus subtilis. L’analisi spettroscopica dimostra che GabR interagisce con il PLP formando l’aldimina interna, mentre in presenza di GABA si ottiene l’aldimina esterna. L’interazione fra il promotore gabTD e le forme holo e apo di GabR è stata monitorata mediante Microscopia a Forza atomica (AFM). In queste due condizioni di legame è stata stimata una Kd di circa 40 ηM. La presenza di GABA invece, determinava un incremento di circa due volte della Kd, variazioni strutturali nei complessi GabR-DNA e una riduzione del compattamento del DNA alla proteina, indipendentemente dalla sequenza del promotore in esame. Al fine di valutare il ruolo delle caratteristiche topologiche del promotore, sono state inserite cinque e dieci bp all’interno della regione spacer che separa le due sequenze ripetute dirette riconosciute da GabR. I significativi cambiamenti topologici riscontrati nel frammento aggiunto di cinque bp si riflettono anche sulla forte riduzione dell’affinità di legame verso la proteina. Al contrario, l’inserzione di 10 bp provoca solamente l’allontanamento delle sequenze ripetute dirette. L’assenza quindi di cambiamenti significativi nella topologia di questo promotore fa sì che l’affinità di legame per GabR rimanga pressoché inalterata rispetto al promotore non mutato. L’analisi del potenziale elettrostatico superficiale di GabR mostra la presenza di una fascia carica positivamente che si estende lungo un’intera faccia della proteina. Per verificare l’importanza di questa caratteristica di GabR nel meccanismo di interazione al DNA, sono stati preparati ed indagati i mutanti R129Q e K362-366Q, in cui la carica positiva superficiale risultava indebolita. L’affinità di legame dei mutanti di GabR per il DNA era inferiore rispetto alla proteina non mutata, in particolar modo nel mutante K362-366Q. Le evidenze acquisite suggeriscono che la curvatura intrinseca del promotore ed il corretto orientamento delle sequenze sulla doppia elica, più della distanza che le separa, siano critici per sostenere l’interazione con GabR. Oltre a questo, la superficie positiva di GabR è richiesta per accomodare la curvatura del DNA sul corpo della proteina. Alla luce di questo, l’interazione GabR-gabTD è un esempio di come il riconoscimento specifico di sequenze, la topologia del DNA e le caratteristiche strutturali della proteina siano contemporaneamente necessarie per sostenere un’interazione proteina-DNA stabile.

Antitumour imidazotetrazines: probes for the major groove of DNA

The antitumour imidazotetrazinones are believed to act as prodrugs for the triazene series of alkylating agents, showing a marked pteference for the alkylation of the middle guanine residue in a run of three or more contiguous guanines. However, the. exact nature of the interactions of imidazotetrazinones within the micro~environment of DNA are; as yet unknown. In order to examine such interactions a three pronged approach involving molecular modelling, synthetic chemistry and biological analysis has been undertaken during the course of this project. . Molecular modelling studies have shown that for the 8-carboxamido substituted imidazotetrazinones antitumour activity is dependent upon the. presence of a free NH group which can be involved in the formation of both intramolecular and intermolecular hydrogen bonds, and the presence of a non-bulky substituent with a small negative potential . volume. Modelling studies involving the docking of .mitozolomide into the major groove of DNA in the region of a triguanine sequence has shown that a number of hydrogen bonding interactions are feasible. A series of 8-substituted carboxamide derivatives of mitozolomide have been synthesised via the 8-acid chloride and 8-carboxylic acid derivatives including a number of peptide analogues. The peptide derivatives were based upon the key structural features of the helix-turn-helix motif of DNA-binding proteins with a view to developing agents that are capable of binding to DNA with greater selectivity. An examination of the importance of intramolecular hydrogen bonding in influencing the antitumour activity:of :the imidazotetrazinones has led to the synthesis of the novel pyrimido[4',5' :4,3]pyrazolo[5,1-d]-1,2,3,5-tetrazine ring system. In general, in vitro cytotoxicity assays showed that the new derivatives were less active against the TLX5 lymphoma cell line. than the parent compound mitozolomide despite an increased potential for hydrogen bonding interactions. Due to the high reactivity of the: tetrazinone ring system it is difficult to study the interactions between the imidazotetrazinones and DNA. Consequently a number of structural analogues that are stable under physiological conditions have been. prepared based upon the 1,2,3 triazin-4(3H)-one ring system fused with both benzene and pyrazole rings. Although the 3-methylbenzotriazinones failed to antagonise the cytotoxic activity of temozolomide encouraging results with a 3-methylpyrazolotriazinone may suggest the existence of an imidazotetrazinone receptor site within DNA. The potential of guanine rich sequences to promote the alkylating selectivity of imidazotetrazinones by acting as a catalyst for ring cleavage and thereby generation of the alkylating agent was examined. Experiments involving the monitoring: of the rate of breakdown of mitozolomide incubated in the presence of synthetic oIigonucleotides did not reveal any catalytic effect resulting from the DNA. However, it was noted that the breakdown of mitozolomide was dependent upon the type of buffer used in the incubations and this may indeed mask any catalysis by the oligonucleotides.

Discovering properties of new DNA-binding activity of proteins

Protein-DNA interactions are an essential feature in the genetic activities of life, and the ability to predict and manipulate such interactions has applications in a wide range of fields. This Thesis presents the methods of modelling the properties of protein-DNA interactions. In particular, it investigates the methods of visualising and predicting the specificity of DNA-binding Cys2His2 zinc finger interaction. The Cys2His2 zinc finger proteins interact via their individual fingers to base pair subsites on the target DNA. Four key residue positions on the a- helix of the zinc fingers make non-covalent interactions with the DNA with sequence specificity. Mutating these key residues generates combinatorial possibilities that could potentially bind to any DNA segment of interest. Many attempts have been made to predict the binding interaction using structural and chemical information, but with only limited success. The most important contribution of the thesis is that the developed model allows for the binding properties of a given protein-DNA binding to be visualised in relation to other protein-DNA combinations without having to explicitly physically model the specific protein molecule and specific DNA sequence. To prove this, various databases were generated, including a synthetic database which includes all possible combinations of the DNA-binding Cys2His2 zinc finger interactions. NeuroScale, a topographic visualisation technique, is exploited to represent the geometric structures of the protein-DNA interactions by measuring dissimilarity between the data points. In order to verify the effect of visualisation on understanding the binding properties of the DNA-binding Cys2His2 zinc finger interaction, various prediction models are constructed by using both the high dimensional original data and the represented data in low dimensional feature space. Finally, novel data sets are studied through the selected visualisation models based on the experimental DNA-zinc finger protein database. The result of the NeuroScale projection shows that different dissimilarity representations give distinctive structural groupings, but clustering in biologically-interesting ways. This method can be used to forecast the physiochemical properties of the novel proteins which may be beneficial for therapeutic purposes involving genome targeting in general.

Statistical Analysis of Data on Linker Histones/DNA Interactions

2000 Mathematics Subject Classification: 62P10, 92C40

Interaction of Multiwalled Carbon Nanotube Produces Structural Alteration of DNA

Multiwalled carbon nanotube (MWCNT) has been found to produce structural changes in Calf Thymus-DNA (CT-DNA). The interaction or binding of the multi-walled carbon nanotubes (MWCNT) was investigated in order to discover if it brings about any significant changes of the DNA double helix using CD spectra of the CT-DNA at two concentration levels of MWCNT representing an increasing MWCNT/DNA molar ratio. In addition, spectrophotometric titrations between MWCNT and CT-DNA were carried out in order to utilize spectral changes as a means of detecting specific binding modes of either intercalation or degradation of DNA. Interactions of MWCNT induced significant changes in the CD spectra of the B-form of natural DNA. The intensities of the positive CD band at 280 nm decreased significantly. This decrease was found to be concentration-dependent. Following spectrophotometric titrations; specific subtle conformational changes were observed with a molar ratio combination of 2:1 between MWCNT and CT-DNA and these were characterized by a formation constant of the order of 103 M-1 and a negative Gibbs free energy suggesting that MWCNT avidly binds to DNA. Thermodynamic considerations revealed that electrostatic interactions between the DNA base pairs and the MWCNT are taking place accounting for the negative free energy change, positive enthalpy change with a small entropy change. The results obtained in the study of the binding interactions of MWCNT with DNA confirm that a cytogenetic effect of MWCNT with DNA is a possibility in vivo.