Recognition of DNA Termini by the C-Terminal Region of the Ku80 and the DNA-Dependent Protein Kinase Catalytic Subunit

DNA double strand breaks (DSBs) can be generated by endogenous cellular processes or exogenous agents in mammalian cells. These breaks are highly variable with respect to DNA sequence and structure and all are recognized in some context by the DNA-dependent protein kinase (DNA-PK). DNA-PK is a critical component necessary for the recognition and repair of DSBs via non-homologous end joining (NHEJ). Previously studies have shown that DNA-PK responds differentially to variations in DSB structure, but how DNA-PK senses differences in DNA substrate sequence and structure is unknown. Here we explore the enzymatic mechanisms by which DNA-PK is activated by various DNA substrates and provide evidence that the DNA-PK is differentially activated by DNA structural variations as a function of the C-terminal region of Ku80. Discrimination based on terminal DNA sequence variations, on the other hand, is independent of the Ku80 C-terminal interactions and likely results exclusively from DNA-dependent protein kinase catalytic subunit interactions with the DNA. We also show that sequence differences in DNA termini can drastically influence DNA repair through altered DNA-PK activation. These results indicate that even subtle differences in DNA substrates influence DNA-PK activation and ultimately the efficiency of DSB repair.     

Polymorphic crystal structures of an all‐AT DNA dodecamer

In this work, we explore the influence of different solvents and ions on the crystallization behavior of an all‐AT dodecamer d(AATAAATTTATT)₂ In all cases, the oligonucleotides are found as continuous columns of stacked duplexes. The spatial organization of such columns is variable; consequently we have obtained seven different crystal forms. The duplexes can be made to crystallize in either parallel or crossed columns. Such versatility in the formation of a variety of crystal forms is characteristic for this sequence. It had not been previously reported for any other sequence. In all cases, the oligonucleotide duplexes have been found to crystallize in the B form. The crystallization conditions determine the organization of the crystal, although no clear local interactions have been detected. Mg²⁺ and Ni²⁺ can be used in order to obtain compact crossed structures. DNA–DNA interactions in the crystals of our all‐AT duplexes present crossovers which are different from those previously reported for mixed sequence oligonucleotides. Our results demonstrate that changes in the ionic atmosphere and the crystallization solvent have a strong influence on the DNA–DNA interactions. Similar ionic changes will certainly influence the biological activity of DNA. Modulation of the crystal structure by ions should also be explored in DNA crystal engineering. Liquid crystals with a peculiar macroscopic shape have also been observed. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 123–133, 2015.     

Profiling cellular protein complexes by proximity ligation with dual tag microarray readout

Patterns of protein interactions provide important insights in basic biology, and their analysis plays an increasing role in drug development and diagnostics of disease. We have established a scalable technique to compare two biological samples for the levels of all pairwise interactions among a set of targeted protein molecules. The technique is a combination of the proximity ligation assay with readout via dual tag microarrays. In the proximity ligation assay protein identities are encoded as DNA sequences by attaching DNA oligonucleotides to antibodies directed against the proteins of interest. Upon binding by pairs of antibodies to proteins present in the same molecular complexes, ligation reactions give rise to reporter DNA molecules that contain the combined sequence information from the two DNA strands. The ligation reactions also serve to incorporate a sample barcode in the reporter molecules to allow for direct comparison between pairs of samples. The samples are evaluated using a dual tag microarray where information is decoded, revealing which pairs of tags that have become joined. As a proof-of-concept we demonstrate that this approach can be used to detect a set of five proteins and their pairwise interactions both in cellular lysates and in fixed tissue culture cells. This paper provides a general strategy to analyze the extent of any pairwise interactions in large sets of molecules by decoding reporter DNA strands that identify the interacting molecules.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

DNA structural variations in the E. coli tyrT promoter

X-ray studies have established that the structure of a right-handed, Watson-Crick double helix can change from place to place along its length as a function of base sequence. The base pairs transmit deformations out to the phosphate backbone, where they can then be recognized by proteins and other DNA-binding reagents. Here we have examined at single-bond resolution the interactions of three commonly used nucleases (DNAase I, DNAase II, and copper-phenanthroline) with a DNA of natural origin, the 160 bp tyrT promoter. All three of these reagents seem sensitive to DNA backbone geometry rather than base sequence per se. Their sequence-dependent patterns of cleavage provide evidence for structural polymorphism of several sorts: global variation in helix groove width, global variation in radial asymmetry, and local variation in phosphate accessibility. These findings explain how sequence zones of a certain base composition, or purine-pyrimidine asymmetry, can influence the recognition of DNA by protein molecules.     

In vitro site-directed mutagenesis with synthetic DNA oligonucleotides yields unexpected deletions and insertions at high frequency.

We have used in vitro site-directed mutagenesis with synthetic DNA oligonucleotides to introduce single nucleotide mutations in yeast mtDNA. In addition to the expected DNA alterations we also recovered with high frequency mutants with large deletions and insertions which arose through interaction with the synthetic DNA fragment. Characterization of a number of these by DNA sequence analysis has permitted reconstruction of the mutagenic events. In all cases, the DNA fragment had base paired with non-adjacent DNA sequences sometimes more than 1000 nucleotides apart from each other on the target strand. The products of such interactions cannot be avoided due to the non-stringent annealing conditions during complementary DNA strand synthesis. However, deliberate mispairing can be directed precisely, as shown by our ability to specifically delete the 1143-bp intron from the yeast mitochondrial gene coding for large ribosomal RNA with a synthetic DNA fragment consisting of the sequence of the exon borders flanking the intron.     

Examining the contribution of a dA+dT element to the conformation of Escherichia coli integration host factor-DNA complexes.

DNA binding proteins that induce structural changes in DNA are common in both prokaryotes and eukaryotes. Integration host factor (IHF) is a multi-functional DNA binding and bending protein of Escherichia coli that can mediate protein-protein and protein-DNA interactions by bending DNA. Previously we have shown that the presence of a dA+dT element 5'-proximal to an IHF consensus sequence can affect the binding of IHF to a particular site. In this study the contribution of various sequence elements to the formation of IHF-DNA complexes was examined. We show that IHF bends DNA more when it binds to a site containing a dA+dT element upstream of its core consensus element than to a site lacking a dA+dT element. We demonstrate that IHF can be specifically crosslinked to DNA with binding sites either containing or lacking this dA+dT element. These results indicate the importance of flanking DNA and a dA+dT element in the binding and bending of a site by IHF.     

DNA-induced alpha-helix capping in conserved linker sequences is a determinant of binding affinity in Cys(2)-His(2) zinc fingers

High-affinity, sequence-specific DNA binding by Cys(2)-His(2) zinc finger proteins is mediated by both specific protein-base interactions and non-specific contacts between charged side-chains and the phosphate backbone. In addition, in DNA complexes of multiple zinc fingers, protein-protein interactions between the finger units contribute to the binding affinity. We present NMR evidence for another contribution to high- affinity binding, a highly specific DNA-induced helix capping involving residues in the linker sequence between fingers. Capping at the C terminus of the alpha-helix in each zinc finger, incorporating a consensus TGEKP linker sequence that follows each finger, provides substantial binding energy to the DNA complexes of zinc fingers 1-3 of TFIIIA (zf1-3) and the four zinc fingers of the Wilms' tumor suppressor protein (wt1-4). The same alpha-helix C-capping motif is observed in the X-ray structures of four other protein-DNA complexes. The structures of each of the TGEKP linkers in these complexes can be superimposed on the linker sequences in the zf1-3 complex, revealing a remarkable similarity in both backbone and side-chain conformations. The canonical linker structures from the zinc-finger-DNA complexes have been compared to the NMR structure of the TGEKP linker connecting fingers 1 and 2 in zf1-3 in the absence of DNA. This comparison reveals that additional stabilization likely arises in the DNA complexes from hydrogen bonding between the backbone amide of E3 and the side-chain O(gamma) of T1 in the linker. We suggest that these DNA-induced C-capping interactions provide a means whereby the multiple-finger complex, which must necessarily be domain-flexible in the unbound state as it searches for the correct DNA sequence, can be "snap-locked" in place once the correct DNA sequence is encountered. These observations provide a rationale for the high conservation of the TGEKP linker sequences in Cys(2)-His(2) zinc finger proteins.     

Establishment of distinct MyoD, E2A, and twist DNA binding specificities by different basic region-DNA conformations

Basic helix-loop-helix (bHLH) proteins perform a wide variety of biological functions. Most bHLH proteins recognize the consensus DNA sequence CAN NTG (the E-box consensus sequence is underlined) but acquire further functional specificity by preferring distinct internal and flanking bases. In addition, induction of myogenesis by MyoD-related bHLH proteins depends on myogenic basic region (BR) and BR-HLH junction residues that are not essential for binding to a muscle-specific site, implying that their BRs may be involved in other critical interactions. We have investigated whether the myogenic residues influence DNA sequence recognition and how MyoD, Twist, and their E2A partner proteins prefer distinct CAN NTG sites. In MyoD, the myogenic BR residues establish specificity for particular CAN NTG sites indirectly, by influencing the conformation through which the BR helix binds DNA. An analysis of DNA binding by BR and junction mutants suggests that an appropriate BR-DNA conformation is necessary but not sufficient for myogenesis, supporting the model that additional interactions with this region are important. The sequence specificities of E2A and Twist proteins require the corresponding BR residues. In addition, mechanisms that position the BR allow E2A to prefer distinct half-sites as a heterodimer with MyoD or Twist, indicating that the E2A BR can be directed toward different targets by dimerization with different partners. Our findings indicate that E2A and its partner bHLH proteins bind to CAN NTG sites by adopting particular preferred BR-DNA conformations, from which they derive differences in sequence recognition that can be important for functional specificity.     

The use of site-directed mutagenesis, transient transfection, and radioligand binding

Receptor-ligand interactions have traditionally been evaluated using a number of biochemical techniques including radioligand binding, photoaffinity labeling, crosslinking, and chemical modification. In modern biochemistry, these approaches have largely been superseded by site-directed mutagenesis in the study of protein function, owing in part to a better understanding of the chemical properties of oligonucleotides and to the ease with which mutant clones can now be generated. The Altered Sites II in vitro Mutagenesis System from the Promega Corporation employs oligonucleotides containing two mismatches to introduce specific nucleotide substitutions in the nucleic acid sequence of a target DNA. One of these mismatches will alter the primary sequence of a given protein, whereas the second will give rise to a silent restriction site that is used to screen for mutants. Transient transfection of tsA201 cells with mutant cDNA constructs using calcium phosphate as a carrier for plasmid DNA permits expression of recombinant receptors that can be characterized using radioligand binding assays. In this article, we focus on site-directed mutagenesis, heterologous expression in eukaryotic cells, and radioligand binding as a methodology to enable the characterization of receptor-ligand interactions.     

Assignment of position-specific error probability to primary DNA sequence data.

DNA sequence predicted from polyacrylamide gel-based technologies is inaccurate because of variations in the quality of the primary data due to limitations of the technology, and to sequence-specific variations due to nucleotide interactions within the DNA molecule and with the gel. The ability to recognize the probability of error in the primary data will be useful in reconstructing the target sequence of a DNA sequencing project, and in estimating the accuracy of the final sequence. This paper describes the use of linear discriminant analysis to assign position-specific probabilities of incorrect, over- and under-prediction of nucleotides for each predicted nucleotide position in primary sequence data generated by a gel-based DNA sequencing technology. Using this method, most of the error potential in primary sequence data can be assigned to a limited number of discrete positions. The use of probability values in the sequence reconstruction process, and in estimating the accuracy of consensus sequence determination is described.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

The study of interactions between DNA and PcrA DNA helicase by using targeted molecular dynamic simulations

DNA helicases are important enzymes involved in all aspects of nucleic acid metabolism, ranging from DNA replication and repair to recombination, rescue of stalled replication and translation. DNA helicases are molecular motors. Through conformational changes caused by ATP hydrolysis and binding, they move along the template double helix, break the hydrogen bonds between the two strands and separate the template chains, so that the genetic information can be accessed. In this paper, targeted molecular dynamic simulations were performed to study the important interactions between DNA and PcrA DNA helicase, which can not be observed from the crystal structures. The key residues on PcrA DNA helicase that have strong interactions with both double stranded DNA (ds-DNA) and single stranded DNA (ss-DNA) have been identified, and it was found that such interactions mostly exist between the protein and DNA backbone, which indicates that the translocation of PcrA is independent of the DNA sequence. The simulations indicate that the ds-DNA is separated upon ATP rebinding, rather than ATP hydrolysis, which suggests that the two strokes in the mechanism have two different major roles. Firstly, in the power stroke (ATP hydrolysis), most of the translocations of the bases from one pocket to the next occur. In the relaxation stroke (ATP binding), most of the ‘work’ is being done to ‘melt’ the DNA at the separation fork. Therefore, we propose a mechanism whereby the translocation of the ss-DNA is powered by ATP hydrolysis and the separation of the ds-DNA is powered by ATP binding.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Understanding the sequence-dependence of DNA groove dimensions: implications for DNA interactions

Background

The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves.

Methodology/Principal Findings

The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR ³¹P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The ³¹P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties.

Conclusions/Significance

The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions.     

Interaction of a protooncogene product, Myb with DNAs.

The DNA-binding domain of Myb consists of three imperfect tandem repeats and the third one which is essential for sequence-specific binding was established to have a helix-turn-helix-related motif. DNA sequences recognized by Myb have been reported to contain TAACPy sequence. Here we have examined the details of Myb-binding sequence. Using DNAs with a single mutation on the various sites of two specific DNAs and some fragments of the DNA-binding domain of Myb, we have found that (i) in a specific DNA which contains only one AAC sequence, each AAC nucleotide is found to be essential for the specific binding of Myb, while any other mutations cause no serious binding loss, (ii) in a specific DNA which contains two AAC sequences separately, one AAC is not so important in the binding, and (iii) for the specific binding with DNA, at least both repeats 2 and 3 of Myb are required. These findings suggest that repeat 3 containing a helix-turn-helix-related structure recognizes the core AAC sequence and repeat 2 supports this recognition by interactions with phosphate groups of DNA.     

Binding of EcoP15I DNA methyltransferase to DNA reveals a large structural distortion within the recognition sequence

EcoP15I DNA methyltransferase, a member of the type III restriction-modification system, binds to the sequence 5'-CAGCAG-3' transferring a methyl group from S-adenosyl-l-methionine to the second adenine base. We have investigated protein-DNA interactions in the methylase-DNA complex by three methods. Determination of equilibrium dissociation constants indicated that the enzyme had higher affinity for DNA containing mismatches at the target base within the recognition sequence. Potassium permanganate footprinting studies revealed that there was a hyper-reactive permanganate cleavage site coincident with adenine that is the target base for methylation. More importantly, to detect DNA conformational alterations within the enzyme-DNA complexes, we have used a fluorescence-based assay. When EcoP15I DNA methyltransferase bound to DNA containing 2-aminopurine substitutions within the cognate sequence, an eight to tenfold fluorescent enhancement resulting from enzymatic flipping of the target adenine base was observed. Furthermore, fluorescence spectroscopy analysis showed that the changes attributable to structural distortion were specific for only the bases within the recognition sequence. More importantly, we observed that both the adenine bases in the recognition site appear to be structurally distorted to the same extent. While the target adenine base is probably flipped out of the DNA duplex, our results also suggest that fluorescent enhancements could be derived from protein-DNA interactions other than base flipping. Taken together, our results support the proposed base flipping mechanism for adenine methyltransferases.